h2lib-tests 13.1.2601__py3-none-any.whl → 13.1.3102__py3-none-any.whl

h2lib_tests/test_ellipsys_couplings.py

@@ -1,8 +1,9 @@
-from h2lib._h2lib import MultiH2Lib, H2LibThread
+from h2lib._h2lib import MultiH2Lib, H2LibThread, H2Lib
 import numpy as np
 from wetb.hawc2.htc_file import HTCFile
 from numpy import newaxis as na
 from h2lib_tests.test_files import tfp
+import pytest
 
 
 class Ellipsys():
@@ -23,17 +24,38 @@ class Ellipsys():
         pass
 
 
+def test_set_aerosections_windspeeds_without_init(capfd):
+
+    with H2Lib(suppress_output=0) as h2:
+        el = Ellipsys()
+        htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT_no_aerodrag.htc')
+        htc.set_name(f'wt0')
+        htc.save()
+
+        h2.init(htc_path=f'htc/wt0.htc', model_path=tfp + 'DTU_10_MW')
+        wt_pos = np.array([0, 0, 0])
+
+        el.step()
+        pos_gl_xyz = np.array(h2.get_aerosections_position(), order='F') + wt_pos[na, na, :]
+        uvw = np.asfortranarray(el.get_uvw(pos_gl_xyz))
+
+        with pytest.raises(Exception, match='H2LibThread process died before or while executing set_aerosections_windspeed'):
+            h2.set_aerosections_windspeed(uvw)
+        out = capfd.readouterr().out
+        assert "Please call init_ad or init_al be" in out
+
+
 def test_ellipsys_dummy_workflow_1wt():
 
     N = 1
-    with MultiH2Lib(N, suppress_output=True) as h2:
+    with MultiH2Lib(N, suppress_output=0) as h2:
         el = Ellipsys()
         htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT_no_aerodrag.htc')
         for i in range(N):
             htc.set_name(f'wt{i}')
             htc.save()
 
-        h2.init(htc_path=f'htc/wt0.htc', model_path=tfp + 'DTU_10_MW')
+        h2.init_AD(htc_path=f'htc/wt0.htc', model_path=tfp + 'DTU_10_MW', tiploss_shen_c2=28)
         wt_pos = np.array([0, 0, 0])
 
         while True:
@@ -51,15 +73,14 @@ def test_ellipsys_dummy_workflow_1wt():
 def test_ellipsys_dummy_workflow():
 
     N = 4
-    with MultiH2Lib(N, suppress_output=True) as h2:
+    with MultiH2Lib(N, suppress_output=1) as h2:
         el = Ellipsys()
         htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT_no_aerodrag.htc')
         for i in range(N):
             htc.set_name(f'wt{i}')
             htc.save()
 
-        h2.init(htc_path=[f'htc/wt{i}.htc' for i in range(N)],
-                model_path=tfp + 'DTU_10_MW')
+        h2.init_AD(htc_path=[f'htc/wt{i}.htc' for i in range(N)], model_path=tfp + 'DTU_10_MW', tiploss_shen_c2=28)
         wt_pos = np.array([[0, 0, 0], [0, 500, 0], [0, 1000, 0], [0, 1500, 0]])
 
         while True:

h2lib_tests/test_h2rotor.py

@@ -1,3 +1,4 @@
+import os
 from h2lib._h2lib import H2Lib, MultiH2Lib
 
 from numpy import testing as npt
@@ -9,10 +10,11 @@ from h2lib_tests.test_files import tfp
 import matplotlib.pyplot as plt
 import numpy as np
 from wetb.hawc2.at_time_file import AtTimeFile
+import h5py
 
 
 def get_h2(htc_path='htc/DTU_10MW_RWT.htc'):
-    h2 = H2Lib(suppress_output=True)
+    h2 = H2Lib(suppress_output=1)
     h2.read_input(htc_path=htc_path, model_path=tfp + 'DTU_10_MW')
     return h2
 
@@ -65,17 +67,17 @@ def test_induction(h2):
 
 def test_rotor_orientation_multi_instance():
     dtu10 = DTU10MW()
+    dtu10.simulation.visualization = 'visu', .5, 1.5
     dtu10.output.buffer = 1
     dtu10.output.data_format = 'gtsdf64'
     dtu10.set_name('tmp_5_0')
-
     dtu10.save()
 
     tilt_ref, yaw_ref = 6, 10
     dtu10.set_tilt_cone_yaw(tilt=tilt_ref, cone=0, yaw=yaw_ref)
     dtu10.set_name('tmp_6_10')
     dtu10.save()
-    with MultiH2Lib(2, suppress_output=True) as mh2:
+    with MultiH2Lib(2, suppress_output=1) as mh2:
         mh2.read_input(['htc/tmp_5_0.htc', 'htc/tmp_6_10.htc'], model_path=tfp + "DTU_10_MW")
         # h2.suppress_output = False
         s_id = mh2.add_sensor('aero power')[0]
@@ -98,6 +100,12 @@ def test_rotor_orientation_multi_instance():
         npt.assert_allclose(data[:, 0], res[1:-1, 0] % 360 - 180, rtol=0.002)  # azi
         npt.assert_array_almost_equal(data[:, 10], res[1:-1, 1])  # power
         npt.assert_array_almost_equal(data[:, 15:18], res[1:-1, 2:5])  # rotor position
+        assert os.path.isfile(mh2.model_path[0] + '/visualization/tmp_5_0.hdf5')
+
+        with h5py.File(mh2.model_path[0] + '/visualization/tmp_5_0.hdf5') as f:
+            assert f.attrs['time_start'][0] == 0.5
+            assert f.attrs['time_stop'][0] == 1.5
+        assert os.path.isfile(mh2.model_path[1] + '/visualization/tmp_6_10.hdf5')
 
 
 def test_rotor_avg_windspeed():
@@ -115,14 +123,14 @@ def test_rotor_avg_windspeed():
 
 
 def test_aerosections():
-    plot = False
+    plot = 0
     h2 = get_h2(htc_path='htc/DTU_10MW_RWT_no_aerodrag.htc')
     # blade 1, global coo, r>30
     pos_ids = [h2.add_sensor(f'aero position 3 1 {xyz} 30')[0] for xyz in [1, 2, 3]]
     wsp_ids = [h2.add_sensor(f'aero windspeed 3 1 {xyz} 30')[0] for xyz in [1, 2, 3]]
     frc_ids = [h2.add_sensor(f'aero secforce 1 {xyz} 30 3')[0] for xyz in [1, 2, 3]]
     mom_ids = [h2.add_sensor(f'aero secmoment 1 {xyz} 30 3')[0] for xyz in [1, 2, 3]]
-    h2.init()
+    h2.init_AL(epsilon_smearing=5)
 
     a = h2.get_aerosections_position()
     if plot:
@@ -146,7 +154,7 @@ def test_aerosections():
     npt.assert_array_almost_equal(a[0, i], [h2.get_sensor_values(id) for id in pos_ids])
     uvw = a * 0
     uvw[:, :, 0] = 6
-    uvw[0, i, 0] = 12
+    h2.set_aerosections_windspeed(uvw)
     h2.run(3)
     npt.assert_array_equal(h2.get_sensor_values(wsp_ids), [0, 6, 0])
     npt.assert_array_almost_equal(h2.get_sensor_values(frc_ids), h2.get_aerosections_forces()[0, i] / 1000)
@@ -157,10 +165,11 @@ def test_aerosections():
         plt.figure()
         plt.plot(frc_before[:, :, 1].T)
 
-    npt.assert_array_almost_equal(h2.get_rotor_avg_wsp(), [0, 6, 0], 10)
-    h2.set_aerosections_windspeed(uvw)
-    # rotor avg freestream wsp unknown after setting wsp incl induction at aero section
+    # rotor avg freestream wsp unknown after init_AD
     npt.assert_array_equal(h2.get_rotor_avg_wsp(), [np.nan, np.nan, np.nan])
+    uvw[0, i, 0] = 12
+    h2.set_aerosections_windspeed(uvw)
+
     h2.step()
     frc_after = h2.get_aerosections_forces()
     mom_after = h2.get_aerosections_moments()
@@ -172,7 +181,7 @@ def test_aerosections():
     assert frc_before[0, i, 1] * 2 < frc_after[0, i, 1]
 
     # rest is similar (within 7N/m, max Fxyz along blade is [331 , 378, 288]
-    frc_after[0, i, :] = frc_before[0, i, :]
+    frc_after[0, i, :] = frc_before[0, i, :]  # reset frc at 30m to previous value. Now everything is similar
     npt.assert_allclose(frc_before, frc_after, atol=7)
     h2.close()
 
@@ -227,7 +236,7 @@ def test_compare_aerosection_coupling():
 
     dtu10.set_name('tmp1')
     dtu10.save()
-    with H2Lib(suppress_output=True) as h2:
+    with H2Lib(suppress_output=1) as h2:
         h2.init(dtu10.filename, dtu10.modelpath)
         h2.run(T)
 
@@ -236,8 +245,8 @@ def test_compare_aerosection_coupling():
     dtu10.set_name('tmp2')
     dtu10.save()
 
-    with H2Lib(suppress_output=True) as h2:
-        h2.init(dtu10.filename, dtu10.modelpath)
+    with H2Lib(suppress_output=1) as h2:
+        h2.init_AD(htc_path=dtu10.filename, model_path=dtu10.modelpath, tiploss_method=0)
         last_pos_gl_xyz = np.array(h2.get_aerosections_position(), order='F')
         while h2.time < T:
             pos_gl_xyz = np.array(h2.get_aerosections_position(), order='F')

h2lib_tests/test_lin.py

@@ -35,6 +35,7 @@ def test_system_not_linearized_3(h2_dtu_10mw_only_blade):
 
 
 def test_linearize(h2_dtu_10mw_only_blade):
+    h2_dtu_10mw_only_blade.structure_reset()
     n_tdofs, n_rdofs = h2_dtu_10mw_only_blade.linearize()
     assert n_rdofs == 26 * 6  # = number of nodes * 6
     assert n_tdofs == 26 * 6 * 2  # times 2 because of speed.
@@ -57,16 +58,14 @@ def test_sys_eig_no_damping(h2_dtu_10mw_only_blade):
     # Test against: result at the time of writing.
     npt.assert_allclose(
         natural_frequencies,
-        np.array([0.610409, 0.930466, 1.739094, 2.761632]),
+        np.array([0.610409, 0.930443, 1.739081, 2.761946]),
         rtol=1e-6,
     )
 
 
 def test_sys_eig_no_damping_wrong_n_modes(h2_dtu_10mw_only_blade):
     with pytest.raises(ValueError, match="TOO_MANY_MODES_REQUESTED"):
-        h2_dtu_10mw_only_blade.get_system_eigenvalues_and_eigenvectors(
-            1000, 156
-        )
+        h2_dtu_10mw_only_blade.get_system_eigenvalues_and_eigenvectors(1000, 156)
 
 
 def test_sys_eig_no_damping_wrong_ny(h2_dtu_10mw_only_blade):
@@ -77,10 +76,8 @@ def test_sys_eig_no_damping_wrong_ny(h2_dtu_10mw_only_blade):
 def test_sys_eig_no_damping_eigv(h2_dtu_10mw_only_blade):
     n_modes = 4
     n_rdofs = 156
-    eig_val, eig_vec = (
-        h2_dtu_10mw_only_blade.get_system_eigenvalues_and_eigenvectors(
-            n_modes, n_rdofs, include_damping=False
-        )
+    eig_val, eig_vec = h2_dtu_10mw_only_blade.get_system_eigenvalues_and_eigenvectors(
+        n_modes, n_rdofs, include_damping=False
     )
     assert eig_val.size == n_modes
     assert eig_vec.shape == (n_modes, n_rdofs)
@@ -88,7 +85,7 @@ def test_sys_eig_no_damping_eigv(h2_dtu_10mw_only_blade):
     assert eig_vec.dtype == np.float64
     npt.assert_allclose(
         eig_val,
-        np.array([3.83531, 5.846293, 10.927047, 17.351849]),
+        np.array([3.835311, 5.846144, 10.92697, 17.353821]),
     )
 
 
@@ -98,20 +95,18 @@ def test_sys_eig_with_damping(h2_dtu_10mw_only_blade):
     )
     # Test against: result at the time of writing.
     npt.assert_allclose(
-        freq, np.array([0.610409, 0.930468, 1.739098, 2.761702]), rtol=1e-6
+        freq, np.array([0.610409, 0.930444, 1.739086, 2.762015]), rtol=1e-6
     )
     npt.assert_allclose(
-        damp, np.array([0.004826, 0.004758, 0.013395, 0.014194]), atol=1e-6
+        damp, np.array([0.004826, 0.004758, 0.013395, 0.014198]), atol=1e-6
     )
 
 
 def test_sys_eig_with_damping_eigv(h2_dtu_10mw_only_blade):
     n_modes = 4
     n_rdofs = 156
-    eig_val, eig_vec = (
-        h2_dtu_10mw_only_blade.get_system_eigenvalues_and_eigenvectors(
-            n_modes, n_rdofs, include_damping=True
-        )
+    eig_val, eig_vec = h2_dtu_10mw_only_blade.get_system_eigenvalues_and_eigenvectors(
+        n_modes, n_rdofs, include_damping=True
     )
     assert eig_val.size == n_modes
     assert eig_vec.shape == (n_modes, 2 * n_rdofs)
@@ -121,10 +116,10 @@ def test_sys_eig_with_damping_eigv(h2_dtu_10mw_only_blade):
         eig_val,
         np.array(
             [
-                -0.01851 - 3.835266j,
-                -0.027817 - 5.846237j,
-                -0.146367 - 10.926094j,
-                -0.246296 - 17.350536j,
+                -0.01851 - 3.835268j,
+                -0.027814 - 5.846088j,
+                -0.146364 - 10.926018j,
+                -0.246401 - 17.352505j,
             ]
         ),
         atol=1e-6,
@@ -152,6 +147,7 @@ def test_get_system_matrices(h2_dtu_10mw_only_blade):
 
 
 def test_sys_eig_encrypted(h2_dtu_10mw_only_tower_encrypted):
+    h2_dtu_10mw_only_tower_encrypted.structure_reset()
     n_tdofs, n_rdofs = h2_dtu_10mw_only_tower_encrypted.linearize()
     n_modes = 4
     freq, damp = h2_dtu_10mw_only_tower_encrypted.do_system_eigenanalysis(
@@ -185,7 +181,7 @@ def test_sys_eig_encrypted(h2_dtu_10mw_only_tower_encrypted):
     )
 
 
-def test_test_get_system_matrices_encrypted(h2_dtu_10mw_only_tower_encrypted):
+def test_get_system_matrices_encrypted(h2_dtu_10mw_only_tower_encrypted):
     n_tdofs, n_rdofs = h2_dtu_10mw_only_tower_encrypted.linearize()
     with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
         h2_dtu_10mw_only_tower_encrypted.get_system_matrices(n_tdofs, n_rdofs)

h2lib_tests/test_mpi.py

@@ -1,24 +1,49 @@
-import h2lib
-from h2lib._h2lib import H2LibProcess, MultiH2Lib, set_LD_LIBRARY_PATH
 import os
+import sys
 import time
+import traceback
 
+from h2lib_tests.test_files import tfp
 from numpy import testing as npt
 import pytest
+from wetb.hawc2.htc_file import HTCFile
 
+import h2lib
+from h2lib._h2lib import H2LibProcess, MultiH2Lib, set_LD_LIBRARY_PATH
+from h2lib_tests.test_files.my_test_cls import MyTest
+from multiclass_interface import mpi_interface
 from multiclass_interface.mpi_interface import MPIClassInterface
 from multiclass_interface.multiprocess_interface import ProcessClass, MultiProcessClassInterface
-
-from h2lib_tests.test_files import tfp
 import numpy as np
-from h2lib_tests.test_files.my_test_cls import MyTest
-from multiclass_interface import mpi_interface
+import shutil
+
 
 set_LD_LIBRARY_PATH
 
+mpi_available = shutil.which('mpirun') is not None
+
+
+def mpirun(f):
+    if 'string' in traceback.format_stack()[0]:
+        return f
+    else:
+        n = f.__name__
+
+        def wrap():
+            if mpi_available:
+                exe = sys.executable.replace("\\", "/")
+                cmd = f'''mpirun -n 4 {exe} -c "from h2lib_tests.test_mpi import {n}; {n}()"'''
+                # print(cmd)
+                assert os.system(cmd) == 0, n
+            else:
+                pytest.xfail('mpirun not available')
+        wrap.__name__ = n
+        return wrap
 
-def mpitest_mpi_MyTest():
 
+@mpirun
+def test_mpi_MyTest():
+    mpi_interface.activate_mpi()
     N = 4
     try:
         with MPIClassInterface(MyTest, [(i,) for i in range(N)]) as m:
@@ -38,28 +63,31 @@ def mpitest_mpi_MyTest():
             with pytest.raises(Exception, match='Cannot make subset of SubsetMPIClassInterface'):
                 m[:3][1]
 
-            print("done, test_mpi_MyTest")
     except ChildProcessError:
         pass
 
 
-def mpitest_mpi_ProcessClass():
+@mpirun
+def test_mpi_ProcessClass():
+    mpi_interface.activate_mpi()
 
     with ProcessClass(MyTest) as cls:
         myTest = cls(1)
         assert myTest.get_id() == 1
-    print("done, test_mpi_ProcessClass")
 
 
-def mpitest_mpi_H2LibProcess():
+@mpirun
+def test_mpi_H2LibProcess():
+    mpi_interface.activate_mpi()
     with H2LibProcess(suppress_output=False) as h2:
         assert h2lib.__version__.replace("+", "-").startswith(h2.get_version().strip()
                                                               ), (h2.get_version().strip(), h2lib.__version__)
-    print("done, test_mpi_H2LibProcess")
 
 
-def mpitest_MultiH2Lib():
-    with MultiH2Lib(3) as mh2:
+@mpirun
+def test_MultiH2Lib():
+    mpi_interface.activate_mpi()
+    with MultiH2Lib(3, suppress_output=1) as mh2:
         assert all([h2lib.__version__.replace("+", "-").startswith(v.strip()) for v in mh2.get_version()])
         assert len(mh2.get_version()) == 3
         if mpi_interface.size > 1:
@@ -93,93 +121,72 @@ def mpitest_MultiH2Lib():
         with pytest.raises(Exception, match='Cannot close SubsetMPIClassInterface. Please close all instances at once'):
             mh2[:2].close()
         with pytest.raises(Exception, match='Cannot make subset of SubsetMPIClassInterface'):
-            print(mh2[:2][1].getState())
-        print('done, test_MPIH2Lib')
+            mh2[:2][1].getState()
 
 
-def ellipsys_mpi_dummy_workflow():
-    from multiclass_interface import mpi_interface
+def get_htc_lst(N=4):
+    htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT_no_aerodrag.htc')
+    for i in range(N):
+        htc.set_name(f'wt{i}')
+        htc.save()
+    return [f'htc/wt{i}.htc' for i in range(N)]
+
+
+@mpirun
+def test_ellipsys_mpi_dummy_workflow():
+    mpi_interface.activate_mpi(collective_mpi=False)
     from h2lib_tests.test_ellipsys_couplings import Ellipsys
-    from wetb.hawc2.htc_file import HTCFile
     from numpy import newaxis as na
-    mpi_interface.LOOP_UNTIL_CLOSE = False
-
     rank = mpi_interface.rank
 
     N = 4
 
-    with MultiH2Lib(N, suppress_output=True) as h2:
+    with MultiH2Lib(N, suppress_output=1) as h2:
         el = Ellipsys()
-        htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT_no_aerodrag.htc')
-        for i in range(N):
-            htc.set_name(f'wt{i}')
-            htc.save()
-        h2.read_input(htc_path=[f'htc/wt{i}.htc' for i in range(N)],
+        h2.read_input(htc_path=get_htc_lst(),
                       model_path=tfp + 'DTU_10_MW')
         wt_pos = np.array([[0, 0, 0], [0, 500, 0], [0, 1000, 0], [0, 1500, 0]])
-        h2.init()
-
+        h2.init_AD()
         t = 0
         while True:
             t = el.step()
             h2.run(t)
-            pos_gl_xyz = np.array(h2.get_aerosections_position(), order='F') + wt_pos[:, na, na, :]
-
-            if rank == 0:
-                assert pos_gl_xyz.shape == (4, 3, 50, 3), pos_gl_xyz.shape  # wt, blades, aero_sections, xyz
-                uvw = np.asfortranarray(el.get_uvw(pos_gl_xyz))
-            else:
-                uvw = None
+            pos_gl_xyz = np.array(h2.get_aerosections_position(), order='F') + wt_pos[[rank], na, na, :]
+            assert pos_gl_xyz.shape == (1, 3, 50, 3), pos_gl_xyz.shape  # 1(current wt), blades, aero_sections, xyz
+            uvw = np.asfortranarray(el.get_uvw(pos_gl_xyz))
             h2.set_aerosections_windspeed(uvw)
             frc_gl_xyz = h2.get_aerosections_forces()
-            if rank == 0:
-                el.set_fxyz(pos_gl_xyz, frc_gl_xyz)
+            assert np.shape(frc_gl_xyz) == (1, 3, 50, 3)
+            el.set_fxyz(pos_gl_xyz, frc_gl_xyz)
             if t == 1:
                 break
 
-    mpi_interface.LOOP_UNTIL_CLOSE = True
 
-
-def mpitest_all():
+@mpirun
+def test_ellipsys_mpi_dummy_workflow_collective():
     mpi_interface.activate_mpi()
-    mpi_interface.exit_mpi_on_close = False
+    from h2lib_tests.test_ellipsys_couplings import Ellipsys
+    from numpy import newaxis as na
     rank = mpi_interface.rank
-    from h2lib_tests.test_ellipsys_couplings import test_ellipsys_dummy_workflow
-    from h2lib_tests.test_h2lib import test_parallel_context_manager, test_parallel
-    from h2lib_tests.test_h2rotor import test_rotor_orientation_multi_instance
-    from h2lib_tests.test_multiprocessinterface import test_attribute, test_missing_attribute, test_execption, test_setattr, test_setattr_method
-    mpi_interface.TERMINATE_ON_CLOSE = False
-
-    def run(f, *args):
-        try:
-            if rank == 0:
-                print(f"Rank {rank} start {f.__name__}", flush=True)
-            f(*args)
-        except ChildProcessError:
-            pass
-        finally:
-            if rank == 0:
-                print(f"Rank {rank} done {f.__name__}", flush=True)
-    for f in [mpitest_mpi_MyTest,
-              mpitest_mpi_ProcessClass, mpitest_mpi_H2LibProcess,
-              mpitest_MultiH2Lib,
-
-              test_parallel,
-              test_parallel_context_manager,
-              test_rotor_orientation_multi_instance,
-              ellipsys_mpi_dummy_workflow,
-
-              ]:
-        run(f)
-
-    try:
-        with MPIClassInterface(MyTest, [(1,), (2,)]) as mpici:
-            for f in [test_attribute, test_missing_attribute, test_execption, test_setattr, test_setattr_method]:
-                break
-                run(f, mpici)
+    N = 4
 
-    except ChildProcessError:
-        pass
-    mpi_interface.TERMINATE_ON_CLOSE = True
-    test_ellipsys_dummy_workflow()
-    print(f"Rank {rank} Done test_all")
+    with MultiH2Lib(N, suppress_output=1) as h2:
+        with h2.release_mpi_workers():
+            el = Ellipsys()
+            h2.read_input(htc_path=get_htc_lst(),
+                          model_path=tfp + 'DTU_10_MW')
+            wt_pos = np.array([[0, 0, 0], [0, 500, 0], [0, 1000, 0], [0, 1500, 0]])
+            h2.init_AD()
+            t = 0
+            while True:
+                t = el.step()
+                h2.run(t)
+                pos_gl_xyz = np.array(h2.get_aerosections_position(), order='F') + wt_pos[[rank], na, na, :]
+                assert pos_gl_xyz.shape == (1, 3, 50, 3), pos_gl_xyz.shape  # 1(current wt), blades, aero_sections, xyz
+                uvw = np.asfortranarray(el.get_uvw(pos_gl_xyz))
+                h2.set_aerosections_windspeed(uvw)
+                frc_gl_xyz = h2.get_aerosections_forces()
+                assert np.shape(frc_gl_xyz) == (1, 3, 50, 3)
+                el.set_fxyz(pos_gl_xyz, frc_gl_xyz)
+                if t == 1:
+                    break

h2lib_tests/test_static_solver.py

@@ -70,8 +70,7 @@ def test_static_solver_run_1(h2_dtu_10mw_only_blade):
     val = h2_dtu_10mw_only_blade.get_sensor_values(id)
     # Test against: initial_condition 2; followed by time simulaiton.
     npt.assert_allclose(
-        val, np.array([-1.072139e+04, -3.856623e-02, -3.949042e+01]),
-        rtol=1e-6
+        val, np.array([-1.072136e04, -3.978756e-02, -4.066154e01]), rtol=1e-6
     )
 
 
@@ -89,7 +88,7 @@ def test_static_solver_run_2(h2_dtu_10mw_only_blade_rotate_base):
     h2_dtu_10mw_only_blade_rotate_base.step()
     val = h2_dtu_10mw_only_blade_rotate_base.get_sensor_values(id)
     # Test against: result at the time of writing.
-    npt.assert_allclose(val, np.array([10879.120647, 383.441758, 989.673212]))
+    npt.assert_allclose(val, np.array([10879.057846, 383.58397, 991.216685]))
 
 
 def test_static_solver_run_3(h2_dtu_10mw_only_blade_rotate_relative):
@@ -100,15 +99,13 @@ def test_static_solver_run_3(h2_dtu_10mw_only_blade_rotate_relative):
     )
 
     # Run the static solver.
-    h2_dtu_10mw_only_blade_rotate_relative.solver_static_run(
-        reset_structure=True
-    )
+    h2_dtu_10mw_only_blade_rotate_relative.solver_static_run(reset_structure=True)
 
     # Do 1 step to get the output.
     h2_dtu_10mw_only_blade_rotate_relative.step()
     val = h2_dtu_10mw_only_blade_rotate_relative.get_sensor_values(id)
     # Test against: result at the time of writing.
-    npt.assert_allclose(val, np.array([10879.074041, 383.440111, 989.67438]))
+    npt.assert_allclose(val, np.array([10879.011239, 383.582323, 991.217846]))
 
 
 def test_static_solver_run_4(h2_dtu_10mw_only_blade_rotate_bearing3):
@@ -119,17 +116,13 @@ def test_static_solver_run_4(h2_dtu_10mw_only_blade_rotate_bearing3):
     )
 
     # Run the static solver.
-    h2_dtu_10mw_only_blade_rotate_bearing3.solver_static_run(
-        reset_structure=True
-    )
+    h2_dtu_10mw_only_blade_rotate_bearing3.solver_static_run(reset_structure=True)
 
     # Do 1 step to get the output.
     h2_dtu_10mw_only_blade_rotate_bearing3.step()
     val = h2_dtu_10mw_only_blade_rotate_bearing3.get_sensor_values(id)
     # Test against: result at the time of writing.
-    npt.assert_allclose(
-        val, np.array([-3094.986918, 115414.095937, 325.296806])
-    )
+    npt.assert_allclose(val, np.array([-3097.095312, 115414.360165, 366.321024]))
 
 
 def test_static_solver_run_no_reset(h2_dtu_10mw_only_blade):
@@ -191,9 +184,7 @@ def test_structure_reset():
         # The deflection must have changed.
         h2.step()
         val_actual = h2.get_sensor_values(id)
-        npt.assert_raises(
-            AssertionError, npt.assert_allclose, val_actual, val_desired
-        )
+        npt.assert_raises(AssertionError, npt.assert_allclose, val_actual, val_desired)
 
         # Reset the structure and check that we match the reference.
         h2.structure_reset()

h2lib_tests/test_topology_h2lib.py

@@ -39,9 +39,7 @@ def test_get_timoshenko_location(
     h2_dtu_10mw_only_tower,
 ):
     # Test first element.
-    l, r1, r12, tes = h2_dtu_10mw_only_tower.get_timoshenko_location(
-        ibdy=0, ielem=0
-    )
+    l, r1, r12, tes = h2_dtu_10mw_only_tower.get_timoshenko_location(ibdy=0, ielem=0)
     assert l - 11.5 < 1e-14
     npt.assert_array_equal(r1, np.array([0.0, 0.0, 0]))
     npt.assert_array_almost_equal_nulp(r12, np.array([0.0, 0.0, -11.5]))
@@ -51,14 +49,10 @@ def test_get_timoshenko_location(
     )
 
     # Test last element.
-    l, r1, r12, tes = h2_dtu_10mw_only_tower.get_timoshenko_location(
-        ibdy=2, ielem=3
-    )
+    l, r1, r12, tes = h2_dtu_10mw_only_tower.get_timoshenko_location(ibdy=2, ielem=3)
     assert l - 12.13 < 1e-14
     npt.assert_array_almost_equal_nulp(r1, np.array([0.0, 0.0, -34.5]))
-    npt.assert_array_almost_equal_nulp(
-        r12, np.array([0.0, 0.0, -12.13]), nulp=3
-    )
+    npt.assert_array_almost_equal_nulp(r12, np.array([0.0, 0.0, -12.13]), nulp=3)
     npt.assert_array_equal(
         tes,
         np.array([[-1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, -1.0]]),
@@ -83,9 +77,7 @@ def test_get_timoshenko_location_encrypted(
     h2_dtu_10mw_only_tower_encrypted,
 ):
     with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
-        h2_dtu_10mw_only_tower_encrypted.get_timoshenko_location(
-            ibdy=0, ielem=0
-        )
+        h2_dtu_10mw_only_tower_encrypted.get_timoshenko_location(ibdy=0, ielem=0)
 
 
 def test_get_body_rotation_tensor_1(h2_dtu_10mw_only_tower):
@@ -121,17 +113,13 @@ def test_set_orientation_base_not_found(h2_dtu_10mw_only_tower):
         h2_dtu_10mw_only_tower.set_orientation_base(main_body="blade")
 
 
-def test_set_orientation_base_1(
-    h2_dtu_10mw_only_tower, h2_dtu_10mw_only_tower_rotated
-):
+def test_set_orientation_base_1(h2_dtu_10mw_only_tower, h2_dtu_10mw_only_tower_rotated):
     # Start from h2_dtu_10mw_only_tower and rotate the base.
     # See if it matches h2_dtu_10mw_only_tower_rotated.
     h2_dtu_10mw_only_tower.set_orientation_base(
         main_body="tower", mbdy_eulerang_table=np.array([30.0, 0.0, 0.0])
     )
-    amat_desired = h2_dtu_10mw_only_tower_rotated.get_body_rotation_tensor(
-        ibdy=0
-    )
+    amat_desired = h2_dtu_10mw_only_tower_rotated.get_body_rotation_tensor(ibdy=0)
     amat_actual = h2_dtu_10mw_only_tower.get_body_rotation_tensor(ibdy=0)
     npt.assert_array_almost_equal_nulp(amat_actual, amat_desired)
     # Reset orientation.
@@ -144,9 +132,7 @@ def test_set_orientation_base_with_reset_orientation(
     h2_dtu_10mw_only_tower_rotated.set_orientation_base(
         main_body="tower", reset_orientation=True
     )
-    amat_actual = h2_dtu_10mw_only_tower_rotated.get_body_rotation_tensor(
-        ibdy=0
-    )
+    amat_actual = h2_dtu_10mw_only_tower_rotated.get_body_rotation_tensor(ibdy=0)
     npt.assert_array_almost_equal_nulp(amat_actual, np.eye(3))
     # Reset orientation.
     h2_dtu_10mw_only_tower_rotated.set_orientation_base(
@@ -162,9 +148,7 @@ def test_set_orientation_base_without_reset_orientation(
         mbdy_eulerang_table=np.array([-30.0, 0.0, 0.0]),
         reset_orientation=False,
     )
-    amat_actual = h2_dtu_10mw_only_tower_rotated.get_body_rotation_tensor(
-        ibdy=0
-    )
+    amat_actual = h2_dtu_10mw_only_tower_rotated.get_body_rotation_tensor(ibdy=0)
     npt.assert_array_almost_equal_nulp(amat_actual, np.eye(3))
     # Reset orientation.
     h2_dtu_10mw_only_tower_rotated.set_orientation_base(
@@ -191,7 +175,7 @@ def test_set_orientation_base_speed(h2_dtu_10mw_only_blade):
 
     # Test against: result at the time of writing.
     # The result is close, but not identical, to test_static_solver_run_2.
-    npt.assert_allclose(val, np.array([10879.429172, 793.66543, 1033.395135]))
+    npt.assert_allclose(val, np.array([10879.363636, 793.815211, 1034.923463]))
 
     # Reset speed.
     h2_dtu_10mw_only_blade.set_orientation_base(
@@ -232,12 +216,12 @@ def test_set_orientation_relative_reset(
         reset_orientation=True,
     )
     # Get orientation of blade root.
-    amat_actual = (
-        h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(ibdy=1)
+    amat_actual = h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(
+        ibdy=1
     )
     # It must be the same as the hub.
-    amat_desired = (
-        h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(ibdy=0)
+    amat_desired = h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(
+        ibdy=0
     )
     npt.assert_array_almost_equal_nulp(amat_actual, amat_desired)
     # This matches a rotation around x by 180 deg.
@@ -263,8 +247,8 @@ def test_set_orientation_relative_2(
     h2_dtu_10mw_only_blade_rotate_relative,
 ):
     # Get orientation of blade root.
-    amat_desired = (
-        h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(ibdy=1)
+    amat_desired = h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(
+        ibdy=1
     )
     # Change orientation a few times.
     rng = np.random.default_rng(seed=123)
@@ -288,8 +272,8 @@ def test_set_orientation_relative_2(
         mbdy2_ini_rotvec_d1=np.array([0.0, 1.0, 0.0, 1.0]),
     )
     # Check.
-    amat_actual = (
-        h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(ibdy=1)
+    amat_actual = h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(
+        ibdy=1
     )
     npt.assert_array_almost_equal_nulp(amat_actual, amat_desired)
 
@@ -320,7 +304,7 @@ def test_set_orientation_relative_static(
     val = h2_dtu_10mw_only_blade_rotate_relative.get_sensor_values(id)
 
     # Test against: result at the time of writing.
-    npt.assert_allclose(val, np.array([8702.239718, 306.617339, 638.840412]))
+    npt.assert_allclose(val, np.array([8702.206018, 306.728782, 640.051269]))
 
     # Reset to original value.
     h2_dtu_10mw_only_blade_rotate_relative.set_orientation_relative(

{h2lib_tests-13.1.2601.dist-info → h2lib_tests-13.1.3102.dist-info}/METADATA

@@ -1,7 +1,7 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: h2lib-tests
-Version: 13.1.2601
-Summary: Tests and test_files for test h2lib (13.1.26)
+Version: 13.1.3102
+Summary: Tests and test_files for test h2lib (13.1.31+5-g08eea40)
 Download-URL: 
 Author: Mads M Pedersen
 Author-email: 

{h2lib_tests-13.1.2601.dist-info → h2lib_tests-13.1.3102.dist-info}/RECORD

@@ -2,14 +2,14 @@ h2lib_tests/__init__.py,sha256=VjSqfGg8BzdmSjfSFhJh4hZbYZ_cME7xp9EWFKHQphA,61
 h2lib_tests/conftest.py,sha256=DAxRUsOP7-3-JUTglMSuyr1qcsdrqeKEDerEtenL0lA,14753
 h2lib_tests/dtu10mw.py,sha256=a7SXfyDwDQPastYKb5CgghOQcYfgO1eTwGrd-H3Enok,4374
 h2lib_tests/test_calc.py,sha256=VNLfr2J9R2Jy9xTbdZ9dfbQ4dCwr7H7nbZRI3yP69fQ,2152
-h2lib_tests/test_ellipsys_couplings.py,sha256=GUvO43QF2YKxXby7eoBObv6Iu0rgMhtOXGcuCHO2JtM,2424
+h2lib_tests/test_ellipsys_couplings.py,sha256=cAo036WqDkwksCX0EXdUUoJ9mJHjpkUPuIlKLW7GnzY,3284
 h2lib_tests/test_h2lib.py,sha256=opgVOkZWOD9CrnTDfwRrBaD-weBsNnutm6FeOijDa88,14197
-h2lib_tests/test_h2rotor.py,sha256=b3enflMN4l4Hvy3odWLXsZ4-9VBRjVaWUSllBRzXbMM,10838
-h2lib_tests/test_lin.py,sha256=FxuDvPuu1MjWE4B7dIWC3qR-9NtjqodGPO1aurUDBwY,6362
-h2lib_tests/test_mpi.py,sha256=sDiUMnZ5j-3y7fOpUy5-s624guTejKwlP-NHRZet3iY,6847
+h2lib_tests/test_h2rotor.py,sha256=4P5bUGvwYDl8Z9hVNP_SXKa79DWhsA-ptpp0v0rT9ow,11309
+h2lib_tests/test_lin.py,sha256=KujYIy9YXTo-uxi_umRWIWZwcHu0oH8GyeVgMcMCq0o,6383
+h2lib_tests/test_mpi.py,sha256=8bMeukbWboijMgWLJIIu3Pn0QtaDadOGhhL099MfMls,6990
 h2lib_tests/test_multiprocessinterface.py,sha256=h2o4havtK6IuMXsplNjGUa3VxOnbpEYGxdrrAKQilj0,1470
-h2lib_tests/test_static_solver.py,sha256=mGapb1vwSoe0Uqpc3-xjDegqauf2R_XMHshnvVojwmY,7509
-h2lib_tests/test_topology_h2lib.py,sha256=dCwvfgO1eHb53PCXwm-nr6iwKdsviU2tuAnU0WIrsJc,11645
+h2lib_tests/test_static_solver.py,sha256=ckzcBAuVge5V5tdrac4PPYxrH5TheNTUc3PKksbB89A,7426
+h2lib_tests/test_topology_h2lib.py,sha256=VjJ1SUBHRHwuSn0EPI6KGzwAULZbPzA39rXdDGXyxkk,11510
 h2lib_tests/test_files/__init__.py,sha256=9e6ZUPb42e0wf2E1rutdcTM8hROcWFRVPXtZriU3ySw,50
 h2lib_tests/test_files/my_test_cls.py,sha256=7ZDsFkxrLfOY6q00U5Y-daxfuhATK-K5H04RP-VmQdE,850
 h2lib_tests/test_files/DTU_10_MW/control/dtu_we_controller.dll,sha256=C5T_CuAFtIuDgCXSYAoNu24yKPwj2nWOeORacJbLN9s,1134592
@@ -89,7 +89,7 @@ h2lib_tests/test_files/minimal/res/minimal_mann_turb.hdf5,sha256=Q3cs3bZyplZjBpo
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_u,sha256=byiorJmXDL6uKFbyfXthHTjJdm6ELvLR2lS202KrhRI,1048576
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_v,sha256=cxK5Rfgfm3gyJsEYi_KlmYY8DIIl_G0aizN2jt18Glc,1048576
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_w,sha256=xs61jAwhP3fIR1P5Oa8ovEt2baLoF8uCNs6pKIT8L4o,1048576
-h2lib_tests-13.1.2601.dist-info/METADATA,sha256=k-4FO0odue2--QCwk0naVXRxBtySjpp3PRP-SvoNWYI,408
-h2lib_tests-13.1.2601.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
-h2lib_tests-13.1.2601.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
-h2lib_tests-13.1.2601.dist-info/RECORD,,
+h2lib_tests-13.1.3102.dist-info/METADATA,sha256=FgPhvGVItQe7EYP_eXjmrMXOnxfT3yXxTI4ALJgN99I,419
+h2lib_tests-13.1.3102.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
+h2lib_tests-13.1.3102.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.1.3102.dist-info/RECORD,,

{h2lib_tests-13.1.2601.dist-info → h2lib_tests-13.1.3102.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.8.0)
+Generator: setuptools (78.1.0)
 Root-Is-Purelib: true
 Tag: py3-none-any