3D nonlinear magnetostatic field computation with finite volume method by means of M-B iteration