A new vector element in the volume integral equation method for nonlinear magnetostatics