Micromagnetics simulations require accurate approximation of the magnetiza- tion dynamics described by the Landau-Lifshitz-Gilbert equation, which is non- linear, nonlocal, and has a non-convex constraint, posing interesting ...