Cogging torque fast computation method for surface permanent magnet machines based on winding function theory