Text this: Solution of three-dimensional, anisotropic, nonlinear problems of magnetostatics using two scalar potentials, finite and infinite multipolar elements and automatic mesh generation