Iterative solvers for nonlinear magnetostatics exhibit global convergence through generalized gradient descent methods.