Application of the surface-current method to the analysis of saturated electromagnetic devices