MagnetoPyElastica, an extension of PyElastica, is an open-source project for simulating magnetic Cosserat rods interacting with external magnetic environments.
I am wondering if MagnetoPyElastica only consider magnetic torque $\tau = m \times B$, where $m$ represents magnetization and $B$ is the externel magnetic field, or if it also accounts for the force induced by non-uniform magnetic field $F = \nabla (m\cdot B)$? I noticed that in magnetic_forces.py, the MagneticForces class only applies torque to the rod.