Scalable magnetic field modeling using structured kernel interpolation for Gaussian process regression