Solution of transient hysteretic magnetic field problems with hybrid Newton-polarization methods