“…The derivative and function evaluations in dynamical and static micromagnetic computations are very expensive, mostly due to the nonlocal stray field component. However, the use of a stray field algorithm that scales optimally or quasi-optimially with problem size such as the algebraic multigrid method [46,42,13], the fast multipole method [50,4,33] and hierarchical matrices [36], fast Fourier transform based methods (FFT) [50,16,49] or non-uniform FFT algorithms [29,15,20] is not sufficient to obtain a micromagnetic solver that scales linearly with problem size. Owing to the exchange interactions the micromagnetic equations can be considered as stiff [9] and the number of time steps in a dynamical solver or the number of iterations in a static solver increase with increasing problem size.…”