computation proofed by tests

magpylib · OrtnerMichael · Sep 24, 2025 · Jul 21, 2025 · Jul 21, 2025 · Jul 21, 2025
commit 0d1a6173711c26139254afb4af2046511a54df3a
diff --git a/src/magpylib/__init__.py b/src/magpylib/__init__.py
@@ -10,7 +10,7 @@
 from magpylib._src.defaults.defaults_classes import default_settings as defaults
 from magpylib._src.defaults.defaults_utility import SUPPORTED_PLOTTING_BACKENDS
 from magpylib._src.display.display import show, show_context
-from magpylib._src.fields import getB, getH, getJ, getM
+from magpylib._src.fields import getB, getH, getJ, getM, getFT
 from magpylib._src.obj_classes.class_Collection import Collection
 from magpylib._src.obj_classes.class_Sensor import Sensor
 
@@ -28,6 +28,7 @@
     "getH",
     "getJ",
     "getM",
+    "getFT",
     "graphics",
     "magnet",
     "misc",

diff --git a/src/magpylib/_src/fields/__init__.py b/src/magpylib/_src/fields/__init__.py
@@ -1,6 +1,7 @@
 """_src.fields"""
 
-__all__ = ["getB", "getH", "getJ", "getM"]
+__all__ = ["getB", "getH", "getJ", "getM", "getFT"]
 
 # create interface to outside of package
 from magpylib._src.fields.field_wrap_BH import getB, getH, getJ, getM
+from magpylib._src.fields.field_FT import getFT
diff --git a/src/magpylib/_src/fields/field_FT.py b/src/magpylib/_src/fields/field_FT.py
@@ -275,12 +275,12 @@ def getFT(sources, targets, pivot="centroid", eps=1e-5, squeeze=True):
 
         force = np.einsum('ijk,ik->ij', DB, MOM) # numpy only
         #force = np.sum(DB * MOM[:, np.newaxis, :], axis=2) # array API
-        torque = np.cross(B, MOM)
+        torque = np.cross(MOM, B)
 
         # Add pivot point contribution to torque
         if pivot is not None:
             PIV = np.tile(np.repeat(pivot[mask_magnet], mesh_sizes_mag, axis=0), (n_sources, 1))
-            torque -= np.cross(POS - PIV, force)
+            torque += np.cross(POS - PIV, force)
 
         # Sum over mesh cells
         idx_starts_all = np.concatenate([idx_starts + n_mesh_mag*j for j in range(n_sources)])
@@ -330,7 +330,7 @@ def getFT(sources, targets, pivot="centroid", eps=1e-5, squeeze=True):
         # Add pivot point contribution to torque
         if pivot is not None:
             PIV = np.tile(np.repeat(pivot[mask_current], mesh_sizes_cur, axis=0), (n_sources, 1))
-            torque -= np.cross(POS - PIV, force)
+            torque += np.cross(POS - PIV, force)
 
         # Sum over mesh cells
         idx_starts_all = np.concatenate([idx_starts + n_mesh_cur*j for j in range(n_sources)])