“…Partial regularity for particular weak solutions is obtained in [21,23]. Numerical simulations are described in [4,6,28,29,30,31,39]. The main difficulties are the quasilinear character of the equations, the preservation of the saturation constraint by the numerical schemes, and the non local character of the demagnetizing field in 3d.…”