A Volume Integral Formulation Based on Facet Elements for Nonlinear Magnetostatic Problems