additional tests

Author: fslanovc

tests/test_basics.py

@@ -1,6 +1,7 @@
 import numpy as np
 import magpylib as magpy
 from magpylib_force import getFT
+from scipy.spatial.transform import Rotation as R
 
 
 def test_rotation1():
@@ -54,3 +55,20 @@ def test_rotation2():
     FT = getFT(s1, [c1, c2], anchor=(0,0,0))
 
     np.testing.assert_allclose(FT[0], FT[1])
+
+
+def test_orientation():
+    """
+    test if dipole with orientation gives same result as rotated magnetic moment
+    """
+    mm, md = np.array((0.976, 4.304, 2.055)), np.array((0.878, -1.527, 2.918))
+    pm, pd = np.array((-1.248, 7.835, 9.273)), np.array((-2.331, 5.835, 0.578))
+
+    magnet = magpy.magnet.Cuboid(position=pm, dimension=(1,2,3), polarization=mm)
+    r = R.from_euler('xyz', (25, 65, 150), degrees=True)
+    dipole1 = magpy.misc.Dipole(position=pd, moment=md, orientation=r)
+    dipole2 = magpy.misc.Dipole(position=pd, moment=r.apply(md))
+
+    F = getFT(magnet, [dipole1, dipole2], anchor=(0,0,0))
+
+    np.testing.assert_allclose(F[0], F[1])

tests/test_physics.py

@@ -266,4 +266,4 @@ def test_physics_force_between_two_dipoles():
         5*r_unit*(np.dot(np.cross(r_unit, m1), np.cross(r_unit, m2)))
         )
 
-    np.testing.assert_allclose(F_num, F_ana)
+    np.testing.assert_allclose(F_num, F_ana, rtol=1e-8, atol=1e-8)

tests/test_self_consistency.py

@@ -267,3 +267,26 @@ def test_consistency_polyline_circle():
     F2,T2 = getFT(src, loop2, anchor=(0,0,0))
     assert abs(np.linalg.norm(F1-F2)/np.linalg.norm(F1+F2)) < 1e-7
     assert abs(np.linalg.norm(T1-T2)/np.linalg.norm(T1+T2)) < 1e-7
+
+
+def test_consistency_sphere_dipole():
+    """
+    force on sphere and dipole should be the same in nearly homogeneous field
+    """
+
+    src = magpy.current.Circle(diameter=10, current=123)
+    pos = (0,0,0)
+    diameter = 0.5
+    magnetization_sphere = np.array((1e6,2e6,3e6))
+    moment_dipole = magnetization_sphere*diameter**3/6*np.pi
+
+    sphere = magpy.magnet.Sphere(diameter=diameter, magnetization=magnetization_sphere, position=pos)
+    sphere.meshing = 100
+
+    dipole = magpy.misc.Dipole(position=pos, moment=moment_dipole)
+
+    FT_sphere = getFT(src, sphere, anchor=(0,0,0))
+    FT_dipole = getFT(src, dipole, anchor=(0,0,0))
+
+    np.testing.assert_allclose(FT_sphere, FT_dipole, rtol=1e-6, atol=1e-6)
+