Im trying to implement the implicit Euler method in high-performance software for micromagnetic simulations, where I'm restricted to using the driving function of the ODE (Landau-Lifshitz equation) and the previous solution points. This obviously not a problem for an explicit method, since we...