Add magnetic field real-space derivatives to `MagneticField` objects

[dpanici]

TODO for later: implement a way to do this for all fields, either with AD like this or through analytic derivatives on a case-by-case basis

Originally posted by @dpanici in #1820

[dpanici]

Maybe would add with or after new magnetic field API #1360