h2lib-tests 13.2.701__py3-none-any.whl → 13.2.801__py3-none-any.whl

h2lib_tests/conftest.py

@@ -4,12 +4,15 @@ from wetb.hawc2.htc_file import HTCFile
 from wetb.hawc2.mainbody import MainBody
 from h2lib_tests.test_files import tfp
 from h2lib._h2lib import H2Lib
+from hawc2models import IEA22MW
 
 
+# %% Write HAWC2 models of the DTU 10 MW.
+
 @pytest.fixture(scope="session")
 def write_dtu10mw_only_tower():
     # Start from DTU_10MW_RWT and delete everything except the tower.
-    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc = HTCFile(f"{tfp}DTU_10_MW/htc/DTU_10MW_RWT.htc")
     htc.set_name("DTU_10MW_RWT_only_tower")
     for key1 in htc["new_htc_structure"].keys():
         if key1.startswith("main_body"):
@@ -69,7 +72,7 @@ def write_dtu10mw_only_tower_encrypted(write_dtu10mw_only_tower):
 @pytest.fixture(scope="session")
 def write_dtu10mw_only_blade():
     # Start from DTU_10MW_RWT and delete everything except the blade.
-    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc = HTCFile(f"{tfp}DTU_10_MW/htc/DTU_10MW_RWT.htc")
     htc.set_name("DTU_10MW_RWT_only_blade")
     for key1 in htc["new_htc_structure"].keys():
         if key1.startswith("main_body"):
@@ -177,7 +180,7 @@ def write_dtu10mw_only_blade_rotate_relative():
     # Start from DTU_10MW_RWT and delete everything except the blade and hub.
     # The blade now rotates because of the relative rotation.
     # Because of the fix1 constraint, the blade will not rotate after time 0.
-    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc = HTCFile(f"{tfp}DTU_10_MW/htc/DTU_10MW_RWT.htc")
     htc.set_name("DTU_10MW_RWT_only_blade_rotate_relative")
     for key1 in (
         "main_body",  # tower
@@ -265,7 +268,7 @@ def write_dtu10mw_only_blade_rotate_relative():
 def write_dtu10mw_only_blade_rotate_bearing3():
     # Start from DTU_10MW_RWT and delete everything except the blade and hub.
     # The blade now rotates because of the bearing3 between the blade and hub.
-    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc = HTCFile(f"{tfp}DTU_10_MW/htc/DTU_10MW_RWT.htc")
     htc.set_name("DTU_10MW_RWT_only_blade_rotate_bearing3")
     for key1 in (
         "main_body",  # tower
@@ -376,11 +379,12 @@ def write_dtu10mw_only_blade_high_gravity(
     htc.new_htc_structure.main_body__7.gravity = 100.0
     mb = MainBody(htc, "blade1")
     blade_c2def = mb.c2def
+    blade_st = mb.stFile
     # Somehow the wind and output blocks are back.
     htc["wind"].delete()
     htc["output"].delete()
     htc.save()
-    return htc, blade_c2def
+    return htc, blade_c2def, blade_st
 
 
 @pytest.fixture(scope="session")
@@ -388,17 +392,18 @@ def write_dtu10mw_only_blade_high_gravity_1_body(
     write_dtu10mw_only_blade_high_gravity,
 ):
     # Start from the write_dtu10mw_only_blade_high_gravity and then set 1 body.
-    htc_ori, _ = write_dtu10mw_only_blade_high_gravity
+    htc_ori, _, _ = write_dtu10mw_only_blade_high_gravity
     htc = htc_ori.copy()
     htc.set_name("DTU_10MW_RWT_only_blade_high_gravity_1_body")
     htc.new_htc_structure.main_body__7.nbodies = 1
     mb = MainBody(htc, "blade1")
     blade_c2def = mb.c2def
+    blade_st = mb.stFile
     # Somehow the wind and output blocks are back.
     htc["wind"].delete()
     htc["output"].delete()
     htc.save()
-    return htc, blade_c2def
+    return htc, blade_c2def, blade_st
 
 
 @pytest.fixture(scope="session")
@@ -429,6 +434,32 @@ def write_dtu10mw_only_blade_high_gravity_deformed(
     return htc, blade_c2def
 
 
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_blade_high_gravity_changed_st(
+    write_dtu10mw_only_blade,
+):
+    # Start from the write_dtu10mw_only_blade and then increase the gravity loading.
+    htc = write_dtu10mw_only_blade.copy()
+    htc.set_name("DTU_10MW_RWT_only_blade_high_gravity_changed_st")
+    htc.new_htc_structure.main_body__7.gravity = 100.0
+    # Change ST
+    mb = MainBody(htc, "blade1")
+    blade_st = mb.stFile
+    rng = np.random.default_rng(seed=582)
+    factor = rng.uniform(low=0.7, high=1.3, size=19)
+    blade_st.main_data_sets[1][1] *= factor
+    # Current directory is hawc2lib/tests/h2lib_tests
+    blade_st.save(f"{tfp}DTU_10_MW/data/changed_st.dat", precision="%28.16e")
+    htc.new_htc_structure.main_body__7.timoschenko_input.filename = "./data/changed_st.dat"
+    # Somehow the wind and output blocks are back.
+    htc["wind"].delete()
+    htc["output"].delete()
+    htc.save()
+    return htc, blade_st
+
+
+# %% Create H2Lib objects of the DTU 10 MW.
+
 @pytest.fixture(scope="session")
 def h2_dtu_10mw_only_tower(write_dtu10mw_only_tower):
     h2 = H2Lib(suppress_output=True)
@@ -592,3 +623,177 @@ def h2_dtu10mw_only_blade_high_gravity_1_body(write_dtu10mw_only_blade_high_grav
     h2.stop_on_error(False)
     yield h2
     h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_dtu10mw_only_blade_high_gravity_changed_st(write_dtu10mw_only_blade_high_gravity_changed_st):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_blade_high_gravity_changed_st.htc"
+    h2.add_sensor("mbdy statevec_new blade1 c2def global absolute 90.0 1.0 0.0 0.0")  # 1-6
+    h2.add_sensor("mbdy forcevec blade1 1 1 blade1")  # 7, 8, 9
+    h2.add_sensor("mbdy momentvec blade1 1 1 blade1")  # 10, 11, 12
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+# %% Write HAWC2 models of the IEA 22 MW.
+
+@pytest.fixture(scope="session")
+def write_iea22mw_rwt():
+    # Download HAWC2 model from GitHub repo.
+    htc = IEA22MW(folder=f"{tfp}IEA-22-280-RWT",
+                  version="5f9bc425bcbda0723a1245bb4e0c1fc5d5439ed3")  # Last commit on master from Aug 26, 2025
+    htc.make_onshore()
+    htc.save()
+    return htc
+
+
+@pytest.fixture(scope="session")
+def write_iea22mw_only_blade(write_iea22mw_rwt):
+    # Start from IEA-22-280-RWT and delete everything except the blade.
+    # htc = write_iea22mw_rwt.copy()
+    htc = HTCFile(f"{tfp}IEA-22-280-RWT/htc/iea_22mw_rwt.htc")
+    htc.set_name("IEA_22MW_RWT_only_blade")
+    for key1 in htc["new_htc_structure"].keys():
+        if key1.startswith("main_body"):
+            if "blade1" not in htc["new_htc_structure"][key1]["name"].values:
+                htc["new_htc_structure"][key1].delete()
+        if key1 == "orientation":
+            htc["new_htc_structure"][key1].delete()
+        if key1 == "constraint":
+            htc["new_htc_structure"][key1].delete()
+    htc["wind"].delete()
+    htc["aerodrag"].delete()
+    htc["aero"].delete()
+    htc["dll"].delete()
+    htc["output"].delete()
+
+    # Set the blade horizontal, to maximize gravity loading.
+    htc.new_htc_structure.add_section("orientation")
+    htc.new_htc_structure.orientation.add_section("base")
+    htc.new_htc_structure.orientation.base.mbdy = "blade1"
+    htc.new_htc_structure.orientation.base.inipos = [0.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.base["mbdy_eulerang"] = [90.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.base.mbdy_eulerang.comments = (
+        "Blade span is horizontal."
+    )
+
+    # Clamp the blade.
+    htc.new_htc_structure.add_section("constraint")
+    htc.new_htc_structure.constraint.add_section("fix0")
+    htc.new_htc_structure.constraint.fix0.mbdy = "blade1"
+
+    # Set low convergence limits.
+    htc.simulation.convergence_limits = [1e2, 1e-5, 1e-07]
+
+    # No output, as we will use add_sensor().
+
+    # Save the new file.
+    htc.save()
+
+    return htc
+
+
+@pytest.fixture(scope="session")
+def write_iea22mw_only_blade_high_gravity(write_iea22mw_only_blade):
+    # Start from the write_iea22mw_only_blade and then increase the gravity loading.
+    htc = write_iea22mw_only_blade.copy()
+    htc.set_name("IEA_22MW_RWT_only_blade_high_gravity")
+    htc.new_htc_structure.main_body__8.gravity = 100.0
+    mb = MainBody(htc, "blade1")
+    blade_c2def = mb.c2def
+    blade_st = mb.stFile
+    # Somehow the wind and output blocks are back.
+    htc["wind"].delete()
+    htc["output"].delete()
+    htc.save()
+    return htc, blade_c2def, blade_st
+
+
+@pytest.fixture(scope="session")
+def write_iea22mw_only_blade_high_gravity_1_body(
+    write_iea22mw_only_blade_high_gravity,
+):
+    # Start from the write_dtu10mw_only_blade_high_gravity and then set 1 body.
+    htc_ori, _, _ = write_iea22mw_only_blade_high_gravity
+    htc = htc_ori.copy()
+    htc.set_name("IEA_22MW_RWT_only_blade_high_gravity_1_body")
+    htc.new_htc_structure.main_body__8.nbodies = 1
+    mb = MainBody(htc, "blade1")
+    blade_c2def = mb.c2def
+    blade_st = mb.stFile
+    # Somehow the wind and output blocks are back.
+    htc["wind"].delete()
+    htc["output"].delete()
+    htc.save()
+    return htc, blade_c2def, blade_st
+
+
+@pytest.fixture(scope="session")
+def write_iea22mw_only_blade_high_gravity_changed_st(
+    write_iea22mw_only_blade,
+):
+    # Start from the write_iea22mw_only_blade and then increase the gravity loading.
+    htc = write_iea22mw_only_blade.copy()
+    htc.set_name("IEA_22MW_RWT_only_blade_high_gravity_changed_st")
+    htc.new_htc_structure.main_body__8.gravity = 100.0
+    # Change ST
+    mb = MainBody(htc, "blade1")
+    blade_st = mb.stFile
+    rng = np.random.default_rng(seed=582)
+    # We change all columns except for the curved length.
+    factor = rng.uniform(low=0.7, high=1.3, size=29)
+    blade_st.main_data_sets[1][1][:, 1:] *= factor
+    # Current directory is hawc2lib/tests/h2lib_tests
+    blade_st.save(f"{tfp}IEA-22-280-RWT/data/changed_st_fpm.dat", precision="%28.16e")
+    htc.new_htc_structure.main_body__8.timoschenko_input.filename = "./data/changed_st_fpm.dat"
+    # Somehow the wind and output blocks are back.
+    htc["wind"].delete()
+    htc["output"].delete()
+    htc.save()
+    return htc, blade_st
+
+
+# %% Create H2Lib objects of the IEA 22 MW.
+
+@pytest.fixture(scope="session")
+def h2_iea22mw_only_blade_high_gravity(write_iea22mw_only_blade_high_gravity):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}IEA-22-280-RWT/"
+    htc_path = "htc/IEA_22MW_RWT_only_blade_high_gravity.htc"
+    h2.add_sensor("mbdy statevec_new blade1 c2def global absolute 90.0 1.0 0.0 0.0")  # 1-6
+    h2.add_sensor("mbdy forcevec blade1 1 1 blade1")  # 7, 8, 9
+    h2.add_sensor("mbdy momentvec blade1 1 1 blade1")  # 10, 11, 12
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_iea22mw_only_blade_high_gravity_1_body(write_iea22mw_only_blade_high_gravity_1_body):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}IEA-22-280-RWT/"
+    htc_path = "htc/IEA_22MW_RWT_only_blade_high_gravity_1_body.htc"
+    h2.add_sensor("mbdy statevec_new blade1 c2def global absolute 90.0 1.0 0.0 0.0")
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_iea22mw_only_blade_high_gravity_changed_st(write_iea22mw_only_blade_high_gravity_changed_st):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}IEA-22-280-RWT/"
+    htc_path = "htc/IEA_22MW_RWT_only_blade_high_gravity_changed_st.htc"
+    h2.add_sensor("mbdy statevec_new blade1 c2def global absolute 90.0 1.0 0.0 0.0")  # 1-6
+    h2.add_sensor("mbdy forcevec blade1 1 1 blade1")  # 7, 8, 9
+    h2.add_sensor("mbdy momentvec blade1 1 1 blade1")  # 10, 11, 12
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()

h2lib_tests/test_distributed_sections.py

@@ -6,6 +6,7 @@ from h2lib.distributed_sections import DistributedSections, LinkType
 import numpy as np
 from numpy import testing as npt
 import matplotlib.pyplot as plt
+from wetb.hawc2.htc_file import HTCFile
 
 
 def test_distributed_sections():
@@ -109,3 +110,88 @@ def test_distributed_sections():
                                                 [-0.0, 0.7, -46.2],
                                                 [-0.0, 2.4, -92.5],
                                                 [0.0, 3.5, -115.6]], 1)
+
+
+def test_distributed_sections_static_solver():
+    with H2Lib() as h2:
+        model_path = f"{tfp}DTU_10_MW/"
+        htc_path = f"{tfp}DTU_10_MW/htc/DTU_10MW_RWT_only_tower.htc"
+        h2.init(htc_path=htc_path, model_path=model_path)
+        ds = h2.add_distributed_sections(mainbody_name='tower', section_relative_position=[0, .5, 1])
+        h2.initialize_distributed_sections()
+
+        ax = plt.figure().add_subplot(projection='3d')
+
+        def draw(label, ref):
+            pos = h2.get_mainbody_nodes_state(mainbody_nr=1, state='pos')
+            # print(label, np.round(pos[-1], 4).tolist())
+            npt.assert_array_almost_equal(pos[-1], ref, 4)
+            ax.plot(*pos.T, marker='.', label=label)
+            ax.set_zlim([0, -120])
+
+        pos = h2.get_mainbody_nodes_state(mainbody_nr=1, state='pos')
+
+        draw('initial', [0.0, 0.0, -115.63])
+
+        frc = np.zeros((3, 3))
+        frc[:, 1] = 100000
+        h2.set_distributed_section_force_and_moment(ds, sec_frc=frc, sec_mom=frc * 0)
+        h2.solver_static_run(reset_structure=True)
+        draw('set frc + static solver', [0.0, 1.8295, -115.6123])
+        c2_def = np.concatenate([pos, pos[:, :1]], 1)
+        c2_def[:, 0] = np.r_[np.arange(0, 60, 10), np.arange(50, 0, -10)]
+        h2.set_c2_def('tower', c2_def)
+        draw('set_c2_def', [10.0, 0.1184, -115.6296])
+        h2.solver_static_run(reset_structure=True)
+        draw('static solver', [10.2777, 2.8085, -115.5066])
+        h2.set_distributed_section_force_and_moment(ds, sec_frc=-frc, sec_mom=frc * 0)
+        h2.solver_static_run(reset_structure=True)
+        draw('set -frc + static solver', [10.2859, -2.8064, -115.5011])
+        if 0:
+            plt.legend()
+            plt.show()
+        else:
+            plt.close('all')
+
+
+def test_set_distributed_section_force_and_moment_coo():
+    with H2Lib() as h2:
+        model_path = f"{tfp}DTU_10_MW/"
+        htc = HTCFile(model_path + "htc/DTU_10MW_RWT_only_tower.htc")
+        htc.new_htc_structure.orientation.base.body_eulerang = 0, 0, 90
+        htc.save(model_path + "htc/DTU_10MW_RWT_only_tower_rot90.htc")
+        print(htc)
+        h2.init(htc_path=htc.filename, model_path=model_path)
+        ds = h2.add_distributed_sections(mainbody_name='tower', section_relative_position=[0, .5, 1])
+        h2.initialize_distributed_sections()
+
+        ax = plt.figure().add_subplot(projection='3d')
+
+        def draw(label, ref):
+            h2.solver_static_run()
+            pos = h2.get_mainbody_nodes_state(mainbody_nr=1, state='pos')
+            print(label, np.round(pos[-1], 4).tolist())
+            ax.plot(*pos.T, marker='.', label=label)
+            npt.assert_array_almost_equal(pos[-1], ref, 4)
+
+        draw('init', [0.0, 0.0, -115.6281])
+        frc = np.zeros((3, 3))
+        frc[:, 1] = 2e6
+        h2.set_distributed_section_force_and_moment(ds, sec_frc=frc, sec_mom=frc * 0)
+        draw('frc_y_global', [-0.0, 34.7544, -109.7119])
+
+        h2.set_distributed_section_force_and_moment(ds, sec_frc=frc, sec_mom=frc * 0, mainbody_coo_nr=1)
+        draw('frc_y_tower', [-34.7545, 0.0, -109.7119])
+
+        if 0:
+            ax.axis('equal')
+            ax.set_zlim([0, -120])
+            ax.set_xlabel('x')
+            ax.set_ylabel('y')
+            ax.set_zlabel('z')
+            ax.plot([0, 0], [0, 10], label='global y')
+            ax.plot([0, -10], [0, 0], label='tower y')
+            ax.legend()
+            plt.show()
+        else:
+            plt.close('all')

h2lib_tests/test_h2lib.py

@@ -335,13 +335,6 @@ def test_fail():
         assert h2.get_version() == version  # hawc2 still alive
         stop_code = 0
 
-        # cover get_stop_code, get_stop_message and reset_stop_code_and_message, but
-        # reset_stop_code_and_message is already called, so stop_code and stop_message are 0 and ""
-        assert h2.get_stop_code(stop_code)[0][0] == 0
-        stop_msg = " " * 1024
-        assert h2.get_stop_message(stop_msg)[0][0].strip() == ''
-        h2.reset_stop_code_and_message()
-
     with MultiH2Lib(2) as h2:
         # , match=re.escape('H2LibThread process died before or while executing fail(...)')):
         with pytest.raises(Exception, match='MyError'):

h2lib_tests/test_h2rotor.py

@@ -55,11 +55,14 @@ def test_get_mainbody_position_orientation(h2: H2LibThread):
         ax = plt.figure().add_subplot(projection='3d')
         for name, mainbody_nr in mb_dict.items():
             pos, ori = h2.get_mainbody_position_orientation(mainbody_nr, mainbody_coo_nr)
-            ax.plot(*pos, 'o', label=name)
+            c = ax.plot(*pos, 'o', label=name)[0].get_color()
+            nodes_pos = h2.get_mainbody_nodes_state(mainbody_nr, state='pos')
+            ax.plot(*nodes_pos.T, '.-', color=c)
             for (exyz), c in zip(ori, 'rgb'):
                 ax.plot(*np.array([pos, pos + exyz * 10]).T, color=c)
         plt.axis('equal')
         plt.legend()
+        plt.gca().set_zlim([0, -200])
         plt.show()
 
     pos, ori = h2.get_mainbody_position_orientation(mb_dict['blade2'], 0)
@@ -73,6 +76,19 @@ def test_get_mainbody_position_orientation(h2: H2LibThread):
     npt.assert_array_almost_equal(ori, [[0.49826698, 0.03781904, 0.86619844],
                                         [0.07556557, 0.99335236, -0.0868386],
                                         [-0.86372441, 0.10872359, 0.49209687]], 6)
+    npt.assert_array_almost_equal(
+        h2.get_mainbody_nodes_state(mainbody_nr=mb_dict['blade1'], state='pos', mainbody_coo_nr=0)[-1],
+        [-0.113689, -6.504988, -208.222845], 6)
+    npt.assert_array_almost_equal(
+        h2.get_mainbody_nodes_state(mainbody_nr=mb_dict['blade1'], state='pos', mainbody_coo_nr=mb_dict['blade1'])[-1],
+        [-0.069654, -3.326239, 86.364205], 6)
+
+    npt.assert_array_almost_equal(
+        h2.get_mainbody_nodes_state(mainbody_nr=mb_dict['blade1'], state='vel', mainbody_coo_nr=0)[-1],
+        [-17.787036, 1.8e-05, 0.061966], 6)
+    npt.assert_array_almost_equal(
+        h2.get_mainbody_nodes_state(mainbody_nr=mb_dict['blade1'], state='acc', mainbody_coo_nr=0)[-1],
+        [-0.032976, 0.08213, 11.252719], 6)
 
 
 def test_induction(h2):

h2lib_tests/test_topology_h2lib.py

@@ -238,6 +238,16 @@ def test_set_orientation_base_speed(h2_dtu_10mw_only_blade):
     )
 
 
+def test_set_c2_def_too_few_sections(h2_dtu_10mw_only_blade):
+    with pytest.raises(ValueError, match="TOO_FEW_SECTIONS_IN_C2DEF"):
+        h2_dtu_10mw_only_blade.set_c2_def("blade1", np.zeros((1, 4)))
+
+
+def test_set_c2_def_wrong_number_of_columns(h2_dtu_10mw_only_blade):
+    with pytest.raises(ValueError, match="WRONG_NUMBER_OF_COLUMNS"):
+        h2_dtu_10mw_only_blade.set_c2_def("blade1", np.zeros((2, 6)))
+
+
 def test_set_c2_def_main_body_not_found(h2_dtu_10mw_only_blade):
     with pytest.raises(ValueError, match="MAIN_BODY_NOT_FOUND"):
         h2_dtu_10mw_only_blade.set_c2_def(
@@ -246,11 +256,6 @@ def test_set_c2_def_main_body_not_found(h2_dtu_10mw_only_blade):
         )
 
 
-def test_set_c2_def_too_few_sections(h2_dtu_10mw_only_blade):
-    with pytest.raises(ValueError, match="TOO_FEW_SECTIONS_IN_C2DEF"):
-        h2_dtu_10mw_only_blade.set_c2_def("blade1", np.zeros((1, 4)))
-
-
 def test_set_c2_def_different_nsec(h2_dtu_10mw_only_blade):
     with pytest.raises(ValueError, match="DIFFERENT_NSEC"):
         h2_dtu_10mw_only_blade.set_c2_def(
@@ -288,7 +293,7 @@ def test_set_c2_def_blade_static_back_to_original(
     #  4. Revert it to the original one and check that the static solution matches the one from step 2.
 
     # Get the clamped DTU 10 MW blade subjected to high gravity loading.
-    _, blade_c2def = write_dtu10mw_only_blade_high_gravity
+    _, blade_c2def, _ = write_dtu10mw_only_blade_high_gravity
     h2 = h2_dtu10mw_only_blade_high_gravity
 
     # Sensors 1 to 6 are the displacement and rotation of the blade tip.
@@ -340,7 +345,7 @@ def test_set_c2_def_blade_static_deformed(
     npt.assert_allclose(blade_tip_actual, blade_tip_desired)
 
     # Restore c2_def.
-    _, c2def_original = write_dtu10mw_only_blade_high_gravity
+    _, c2def_original, _ = write_dtu10mw_only_blade_high_gravity
     h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_original)
 
 
@@ -351,7 +356,7 @@ def test_set_c2_def_blade_eig_back_to_original(
 ):
     # Get the clamped DTU 10 MW blade.
     # We use the fixture with high gravity because it also returns c2_def.
-    _, blade_c2def = write_dtu10mw_only_blade_high_gravity
+    _, blade_c2def, _ = write_dtu10mw_only_blade_high_gravity
     h2 = h2_dtu10mw_only_blade_high_gravity
 
     # Solve the eigenvalue problem.
@@ -383,13 +388,13 @@ def test_set_c2_def_blade_eig_back_to_original(
     npt.assert_allclose(damping_ratios_actual, damping_ratios_desired)
 
 
-def test_set_c2_def_blade_mass_back_to_original(
+def test_set_c2_def_blade_inertia_back_to_original(
     write_dtu10mw_only_blade_high_gravity_1_body,
     h2_dtu10mw_only_blade_high_gravity_1_body,
 ):
-    # Get the clamped DTU 10 MW blade.
+    # Get the clamped DTU 10 MW blade with only 1 body.
     # We use the fixture with high gravity because it also returns c2_def.
-    _, blade_c2def = write_dtu10mw_only_blade_high_gravity_1_body
+    _, blade_c2def, _ = write_dtu10mw_only_blade_high_gravity_1_body
     h2 = h2_dtu10mw_only_blade_high_gravity_1_body
 
     # Get the inertia properties.
@@ -415,7 +420,7 @@ def test_set_c2_def_blade_mass_back_to_original(
         npt.assert_allclose(inertia_actual[i], inertia_desired[i])
 
 
-def test_set_c2_def_blade_mass_deformed(
+def test_set_c2_def_blade_inertia_deformed(
     write_dtu10mw_only_blade_high_gravity,
     write_dtu10mw_only_blade_high_gravity_deformed,
     h2_dtu10mw_only_blade_high_gravity,
@@ -439,7 +444,7 @@ def test_set_c2_def_blade_mass_deformed(
             npt.assert_allclose(inertia_actual[i], inertia_desired[i], rtol=1e-6)
 
     # Restore c2_def.
-    _, c2def_original = write_dtu10mw_only_blade_high_gravity
+    _, c2def_original, _ = write_dtu10mw_only_blade_high_gravity
     h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_original)
 
 
@@ -469,10 +474,364 @@ def test_set_c2_def_blade_element_deformed(
             npt.assert_allclose(mat_actual[2], mat_desired[2], rtol=1e-10)  # damping
 
     # Restore c2_def.
-    _, c2def_original = write_dtu10mw_only_blade_high_gravity
+    _, c2def_original, _ = write_dtu10mw_only_blade_high_gravity
     h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_original)
 
 
+def test_set_st_wrong_number_of_columns(h2_dtu_10mw_only_blade):
+    with pytest.raises(ValueError, match="WRONG_NUMBER_OF_COLUMNS"):
+        h2_dtu_10mw_only_blade.set_st("blade1", np.empty((2, 6)))
+
+
+def test_set_st_main_body_not_found(h2_dtu_10mw_only_blade):
+    with pytest.raises(ValueError, match="MAIN_BODY_NOT_FOUND"):
+        h2_dtu_10mw_only_blade.set_st("tower", np.empty((2, 19)))
+
+
+def test_set_st_z_not_continuously_increasing(h2_dtu_10mw_only_blade):
+    with pytest.raises(ValueError, match="ST_Z_NOT_CONTINUOUSLY_INCREASING"):
+        h2_dtu_10mw_only_blade.set_st("blade1", np.zeros((3, 19)))
+
+
+def test_set_st_uniform_node_distribution(
+    h2_dtu_10mw_only_blade_uniform_node_distribution,
+):
+    with pytest.raises(NotImplementedError):
+        st = np.zeros((2, 19))
+        st[1, 0] = 1.0
+        h2_dtu_10mw_only_blade_uniform_node_distribution.set_st("blade1", st)
+
+
+def test_set_st_classic_timoshenko_blade_static_back_to_original(
+    write_dtu10mw_only_blade_high_gravity,
+    h2_dtu10mw_only_blade_high_gravity,
+):
+    # This test will:
+    #  1. Take the DTU 10 MW blade subjected to high gravity loading.
+    #  2. Compute the static solution.
+    #  3. Change ST a few times and re-compute the static solution.
+    #  4. Revert it to the original one and check that the static solution matches the one from step 2.
+
+    # Get the clamped DTU 10 MW blade subjected to high gravity loading.
+    _, _, blade_st = write_dtu10mw_only_blade_high_gravity
+    h2 = h2_dtu10mw_only_blade_high_gravity
+
+    # Sensors 1 to 6 are the displacement and rotation of the blade tip.
+    #   h2.get_sensor_info(1)
+
+    # Run the static solver and get the blade tip position and rotation.
+    h2.solver_static_run(reset_structure=True)
+    blade_tip_desired = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
+
+    # Change blade ST and run the static solver.
+    rng = np.random.default_rng(seed=582)
+    for _ in range(10):
+        factor = rng.uniform(low=0.7, high=1.3, size=19)
+        st_new = factor * blade_st.main_data_sets[1][1]
+        # We do not change z, to prevent ST_Z_NOT_CONTINUOUSLY_INCREASING.
+        st_new[:, 0] = blade_st.main_data_sets[1][1][:, 0]
+        h2.set_st("blade1", st_new)
+        h2.solver_static_run(reset_structure=True)
+
+        # Get new blade displacement.
+        blade_tip_actual = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
+
+        # Must differ from blade_tip_desired.
+        with npt.assert_raises(AssertionError):
+            npt.assert_allclose(blade_tip_actual, blade_tip_desired)
+
+    # Restore blade c2_def.
+    h2.set_st("blade1", blade_st.main_data_sets[1][1])
+    h2.solver_static_run(reset_structure=True)
+    blade_tip_actual = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
+    npt.assert_allclose(blade_tip_actual, blade_tip_desired)
+
+
+def test_set_st_classic_timoshenko_blade_static_changed_st(
+    write_dtu10mw_only_blade_high_gravity,
+    write_dtu10mw_only_blade_high_gravity_changed_st,
+    h2_dtu10mw_only_blade_high_gravity,
+    h2_dtu10mw_only_blade_high_gravity_changed_st,
+):
+    # Solve the static problem with the changed ST loaded directly by HAWC2.
+    h2_dtu10mw_only_blade_high_gravity_changed_st.solver_static_run(reset_structure=True)
+    blade_tip_desired = h2_dtu10mw_only_blade_high_gravity_changed_st.get_sensor_values((1, 2, 3, 4, 5, 6))
+
+    # Set ST in the original blade, thus making it equivalent to the changed one.
+    # Then, solve again the static problem.
+    _, st_changed = write_dtu10mw_only_blade_high_gravity_changed_st
+    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+    h2_dtu10mw_only_blade_high_gravity.solver_static_run(reset_structure=True)
+    blade_tip_actual = h2_dtu10mw_only_blade_high_gravity_changed_st.get_sensor_values((1, 2, 3, 4, 5, 6))
+
+    npt.assert_allclose(blade_tip_actual, blade_tip_desired)
+
+    # Restore c2_def.
+    _, _, st_riginal = write_dtu10mw_only_blade_high_gravity
+    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_riginal.main_data_sets[1][1])
+
+
+def test_set_st_classic_timoshenko_inertia_back_to_original(
+    write_dtu10mw_only_blade_high_gravity_1_body,
+    h2_dtu10mw_only_blade_high_gravity_1_body,
+):
+    # Get the clamped DTU 10 MW blade with only 1 body.
+    # We use the fixture with high gravity because it also returns st.
+    _, _, st = write_dtu10mw_only_blade_high_gravity_1_body
+    h2 = h2_dtu10mw_only_blade_high_gravity_1_body
+
+    # Get the inertia properties.
+    h2.structure_reset()
+    inertia_desired = h2.body_output_mass(0)
+
+    # Change blade density and compute inertia properties.
+    mass_original = st.m()
+    rng = np.random.default_rng(seed=582)
+    for _ in range(10):
+        # Uniformly scale to change the mass.
+        density_new = rng.uniform(low=0.5, high=2.0) * mass_original
+        # Make the blade tip heavier to change center of gravity.
+        density_new *= np.linspace(0.8, 2.5, density_new.size)
+        st.set_value(mset=1, set=1, m=density_new)
+        h2.set_st("blade1", st.main_data_sets[1][1])
+        inertia_actual = h2.body_output_mass(0)
+
+        # Must differ from the desired ones.
+        for i in range(4):  # Loop over tuple of arrays.
+            with npt.assert_raises(AssertionError):
+                npt.assert_allclose(inertia_actual[i], inertia_desired[i])
+
+    # Restore ST.
+    st.set_value(mset=1, set=1, m=mass_original)
+    h2.set_st("blade1", st.main_data_sets[1][1])
+    inertia_actual = h2.body_output_mass(0)
+    for i in range(4):  # Loop over tuple of arrays.
+        npt.assert_allclose(inertia_actual[i], inertia_desired[i])
+
+
+def test_set_st_classic_timoshenko_blade_inertia_changed_st(
+    write_dtu10mw_only_blade_high_gravity,
+    write_dtu10mw_only_blade_high_gravity_changed_st,
+    h2_dtu10mw_only_blade_high_gravity,
+    h2_dtu10mw_only_blade_high_gravity_changed_st,
+):
+    # Revert blades to the undeflected configuration.
+    h2_dtu10mw_only_blade_high_gravity.structure_reset()
+    h2_dtu10mw_only_blade_high_gravity_changed_st.structure_reset()
+
+    # Set ST in the original blade, thus making it equivalent to the changed one.
+    _, st_changed = write_dtu10mw_only_blade_high_gravity_changed_st
+    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+
+    # Compare inertia properties for all bodies.
+    nbdy, _ = h2_dtu10mw_only_blade_high_gravity_changed_st.get_number_of_bodies_and_constraints()
+    for i in range(nbdy):
+        inertia_desired = h2_dtu10mw_only_blade_high_gravity_changed_st.body_output_mass(i)
+        inertia_actual = h2_dtu10mw_only_blade_high_gravity.body_output_mass(i)
+
+        for i in range(4):  # Loop over tuple of arrays.
+            npt.assert_allclose(inertia_actual[i], inertia_desired[i])
+
+    # Restore ST.
+    _, _, st_original = write_dtu10mw_only_blade_high_gravity
+    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_original.main_data_sets[1][1])
+
+
+def test_set_st_classic_timoshenko_blade_element_changed_st(
+    write_dtu10mw_only_blade_high_gravity,
+    write_dtu10mw_only_blade_high_gravity_changed_st,
+    h2_dtu10mw_only_blade_high_gravity,
+    h2_dtu10mw_only_blade_high_gravity_changed_st,
+):
+    # Revert blades to the undeflected configuration.
+    h2_dtu10mw_only_blade_high_gravity.structure_reset()
+    h2_dtu10mw_only_blade_high_gravity_changed_st.structure_reset()
+
+    # Set ST in the original blade, thus making it equivalent to the changed one.
+    _, st_changed = write_dtu10mw_only_blade_high_gravity_changed_st
+    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+
+    # Compare element matrices for all bodies.
+    nelem = h2_dtu10mw_only_blade_high_gravity_changed_st.get_number_of_elements()
+    for ibdy in range(nelem.size):
+        for ielem in range(nelem[ibdy]):
+            mat_desired = h2_dtu10mw_only_blade_high_gravity_changed_st.body_output_element(ibdy, ielem)
+            mat_actual = h2_dtu10mw_only_blade_high_gravity.body_output_element(ibdy, ielem)
+
+            npt.assert_allclose(mat_actual[0], mat_desired[0], rtol=1e-14)  # mass
+            npt.assert_allclose(mat_actual[1], mat_desired[1], rtol=1e-14)  # stiffness
+            npt.assert_allclose(mat_actual[2], mat_desired[2], rtol=1e-14)  # damping
+
+    # Restore ST.
+    _, _, st_original = write_dtu10mw_only_blade_high_gravity
+    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_original.main_data_sets[1][1])
+
+
+def test_set_st_fpm_blade_static_back_to_original(
+    write_iea22mw_only_blade_high_gravity,
+    h2_iea22mw_only_blade_high_gravity,
+):
+    # This test will:
+    #  1. Take the IEA 22 MW blade subjected to high gravity loading.
+    #  2. Compute the static solution.
+    #  3. Change ST a few times and re-compute the static solution.
+    #  4. Revert it to the original one and check that the static solution matches the one from step 2.
+
+    # Get the clamped IEA 22 MW blade subjected to high gravity loading.
+    _, _, blade_st = write_iea22mw_only_blade_high_gravity
+    h2 = h2_iea22mw_only_blade_high_gravity
+
+    # Sensors 1 to 6 are the displacement and rotation of the blade tip.
+    #   h2.get_sensor_info(1)
+
+    # Run the static solver and get the blade tip position and rotation.
+    h2.solver_static_run(reset_structure=True)
+    blade_tip_desired = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
+
+    # Change blade ST and run the static solver.
+    rng = np.random.default_rng(seed=582)
+    for _ in range(10):
+        factor = rng.uniform(low=0.7, high=1.3, size=30)
+        st_new = factor * blade_st.main_data_sets[1][1]
+        # We do not change z, to prevent ST_Z_NOT_CONTINUOUSLY_INCREASING.
+        st_new[:, 0] = blade_st.main_data_sets[1][1][:, 0]
+        h2.set_st("blade1", st_new)
+        h2.solver_static_run(reset_structure=True)
+
+        # Get new blade displacement.
+        blade_tip_actual = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
+
+        # Must differ from blade_tip_desired.
+        with npt.assert_raises(AssertionError):
+            npt.assert_allclose(blade_tip_actual, blade_tip_desired)
+
+    # Restore blade c2_def.
+    h2.set_st("blade1", blade_st.main_data_sets[1][1])
+    h2.solver_static_run(reset_structure=True)
+    blade_tip_actual = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
+    npt.assert_allclose(blade_tip_actual, blade_tip_desired)
+
+
+def test_set_st_fpm_blade_static_changed_st(
+    write_iea22mw_only_blade_high_gravity,
+    write_iea22mw_only_blade_high_gravity_changed_st,
+    h2_iea22mw_only_blade_high_gravity,
+    h2_iea22mw_only_blade_high_gravity_changed_st,
+):
+    # Solve the static problem with the changed ST loaded directly by HAWC2.
+    h2_iea22mw_only_blade_high_gravity_changed_st.solver_static_run(reset_structure=True)
+    blade_tip_desired = h2_iea22mw_only_blade_high_gravity_changed_st.get_sensor_values((1, 2, 3, 4, 5, 6))
+
+    # Set ST in the original blade, thus making it equivalent to the changed one.
+    # Then, solve again the static problem.
+    _, st_changed = write_iea22mw_only_blade_high_gravity_changed_st
+    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+    h2_iea22mw_only_blade_high_gravity.solver_static_run(reset_structure=True)
+    blade_tip_actual = h2_iea22mw_only_blade_high_gravity_changed_st.get_sensor_values((1, 2, 3, 4, 5, 6))
+
+    npt.assert_allclose(blade_tip_actual, blade_tip_desired)
+
+    # Restore c2_def.
+    _, _, st_riginal = write_iea22mw_only_blade_high_gravity
+    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_riginal.main_data_sets[1][1])
+
+
+def test_set_st_fpm_inertia_back_to_original(
+    write_iea22mw_only_blade_high_gravity_1_body,
+    h2_iea22mw_only_blade_high_gravity_1_body,
+):
+    # Get the clamped DTU 10 MW blade with only 1 body.
+    # We use the fixture with high gravity because it also returns st.
+    _, _, st = write_iea22mw_only_blade_high_gravity_1_body
+    h2 = h2_iea22mw_only_blade_high_gravity_1_body
+
+    # Get the inertia properties.
+    h2.structure_reset()
+    inertia_desired = h2.body_output_mass(0)
+
+    # Change blade density and compute inertia properties.
+    mass_original = st.m()
+    rng = np.random.default_rng(seed=582)
+    for _ in range(10):
+        # Uniformly scale to change the mass.
+        density_new = rng.uniform(low=0.5, high=2.0) * mass_original
+        # Make the blade tip heavier to change center of gravity.
+        density_new *= np.linspace(0.8, 2.5, density_new.size)
+        st.set_value(mset=1, set=1, m=density_new)
+        h2.set_st("blade1", st.main_data_sets[1][1])
+        inertia_actual = h2.body_output_mass(0)
+
+        # Must differ from the desired ones.
+        for i in range(4):  # Loop over tuple of arrays.
+            with npt.assert_raises(AssertionError):
+                npt.assert_allclose(inertia_actual[i], inertia_desired[i])
+
+    # Restore ST.
+    st.set_value(mset=1, set=1, m=mass_original)
+    h2.set_st("blade1", st.main_data_sets[1][1])
+    inertia_actual = h2.body_output_mass(0)
+    for i in range(4):  # Loop over tuple of arrays.
+        npt.assert_allclose(inertia_actual[i], inertia_desired[i])
+
+
+def test_set_st_fpm_blade_inertia_changed_st(
+    write_iea22mw_only_blade_high_gravity,
+    write_iea22mw_only_blade_high_gravity_changed_st,
+    h2_iea22mw_only_blade_high_gravity,
+    h2_iea22mw_only_blade_high_gravity_changed_st,
+):
+    # Revert blades to the undeflected configuration.
+    h2_iea22mw_only_blade_high_gravity.structure_reset()
+    h2_iea22mw_only_blade_high_gravity_changed_st.structure_reset()
+
+    # Set ST in the original blade, thus making it equivalent to the changed one.
+    _, st_changed = write_iea22mw_only_blade_high_gravity_changed_st
+    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+
+    # Compare inertia properties for all bodies.
+    nbdy, _ = h2_iea22mw_only_blade_high_gravity_changed_st.get_number_of_bodies_and_constraints()
+    for i in range(nbdy):
+        inertia_desired = h2_iea22mw_only_blade_high_gravity_changed_st.body_output_mass(i)
+        inertia_actual = h2_iea22mw_only_blade_high_gravity.body_output_mass(i)
+
+        for i in range(4):  # Loop over tuple of arrays.
+            npt.assert_allclose(inertia_actual[i], inertia_desired[i])
+
+    # Restore ST.
+    _, _, st_original = write_iea22mw_only_blade_high_gravity
+    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_original.main_data_sets[1][1])
+
+
+def test_set_st_fpm_blade_element_changed_st(
+    write_iea22mw_only_blade_high_gravity,
+    write_iea22mw_only_blade_high_gravity_changed_st,
+    h2_iea22mw_only_blade_high_gravity,
+    h2_iea22mw_only_blade_high_gravity_changed_st,
+):
+    # Revert blades to the undeflected configuration.
+    h2_iea22mw_only_blade_high_gravity.structure_reset()
+    h2_iea22mw_only_blade_high_gravity_changed_st.structure_reset()
+
+    # Set ST in the original blade, thus making it equivalent to the changed one.
+    _, st_changed = write_iea22mw_only_blade_high_gravity_changed_st
+    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+
+    # Compare element matrices for all bodies.
+    nelem = h2_iea22mw_only_blade_high_gravity_changed_st.get_number_of_elements()
+    for ibdy in range(nelem.size):
+        for ielem in range(nelem[ibdy]):
+            mat_desired = h2_iea22mw_only_blade_high_gravity_changed_st.body_output_element(ibdy, ielem)
+            mat_actual = h2_iea22mw_only_blade_high_gravity.body_output_element(ibdy, ielem)
+
+            npt.assert_allclose(mat_actual[0], mat_desired[0], rtol=1e-14)  # mass
+            npt.assert_allclose(mat_actual[1], mat_desired[1], rtol=1e-14)  # stiffness
+            npt.assert_allclose(mat_actual[2], mat_desired[2], rtol=1e-14)  # damping
+
+    # Restore ST.
+    _, _, st_original = write_iea22mw_only_blade_high_gravity
+    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_original.main_data_sets[1][1])
+
+
 def test_set_orientation_relative_main_body_not_found(
     h2_dtu_10mw_only_blade_rotate_relative,
 ):
@@ -606,4 +965,4 @@ def test_set_orientation_relative_static(
 
 if __name__ == "__main__":
     # pytest.main([__file__])
-    pytest.main([__file__, "-k test_set_c2_def_blade_element_deformed"])
+    pytest.main([__file__, "-k test_set_st_fpm_blade_element_changed_st"])

{h2lib_tests-13.2.701.dist-info → h2lib_tests-13.2.801.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.4
 Name: h2lib-tests
-Version: 13.2.701
-Summary: Tests and test_files for test h2lib (13.2.7+4-gb779ee7)
+Version: 13.2.801
+Summary: Tests and test_files for test h2lib (13.2.8)
 Download-URL: 
 Author: Mads M Pedersen
 Author-email: 
@@ -13,6 +13,7 @@ Requires-Dist: pytest-cov
 Requires-Dist: hipersim
 Requires-Dist: wetb
 Requires-Dist: h2lib
+Requires-Dist: Hawc2Models==0.0.4
 Dynamic: author
 Dynamic: requires-dist
 Dynamic: summary

{h2lib_tests-13.2.701.dist-info → h2lib_tests-13.2.801.dist-info}/RECORD

@@ -1,16 +1,16 @@
 h2lib_tests/__init__.py,sha256=VjSqfGg8BzdmSjfSFhJh4hZbYZ_cME7xp9EWFKHQphA,61
-h2lib_tests/conftest.py,sha256=o9bkMEaB40rr_KSGovy1nuvq2Y09rKfQWxTcTGZ0MfY,21607
+h2lib_tests/conftest.py,sha256=r4p868zCfJ5MMd2hs6XvGaiqu3jYlSR5G0OKAk5WX7Y,29629
 h2lib_tests/dtu10mw.py,sha256=a7SXfyDwDQPastYKb5CgghOQcYfgO1eTwGrd-H3Enok,4374
 h2lib_tests/test_calc.py,sha256=VNLfr2J9R2Jy9xTbdZ9dfbQ4dCwr7H7nbZRI3yP69fQ,2152
-h2lib_tests/test_distributed_sections.py,sha256=JtguIddd2nmRyNsy2jqjPFB8E5cGpH2_IVErP1wGPPg,5432
+h2lib_tests/test_distributed_sections.py,sha256=-YqB2GgRWmZMophwubRjRh8y5X6ptAm4xqor9sgDQ30,9042
 h2lib_tests/test_ellipsys_couplings.py,sha256=mfSTk1IamaKETU6Be8p8W9vA-yostWJl17m5sFlEK3Q,3345
-h2lib_tests/test_h2lib.py,sha256=LEV1n0Fwa1llQwIM06Mcqr07Y6WOn-j8mNxjzlCEOl0,14219
-h2lib_tests/test_h2rotor.py,sha256=yAcG9RDD3uO4g4w8LJnJ3Md_bC-_mBn1--4KHwvlTRc,12963
+h2lib_tests/test_h2lib.py,sha256=t2jwmXzibej9DQoAkldwWX1MmEmpcV0kp2hLDRV5dbQ,13836
+h2lib_tests/test_h2rotor.py,sha256=cF3KMdzTyGoZGyihhHCJqdH_n_t2vjS1BPHDaUCk9pU,13918
 h2lib_tests/test_lin.py,sha256=yydbsMX44ym_MSqR1cbQkqil8qC4icSW_JccPSElGac,6292
 h2lib_tests/test_mpi.py,sha256=CTT160yc6uZFRr_QNFWxyOwJ-y0qHiIp8jPYs1MQpyg,7111
 h2lib_tests/test_multiprocessinterface.py,sha256=h2o4havtK6IuMXsplNjGUa3VxOnbpEYGxdrrAKQilj0,1470
 h2lib_tests/test_static_solver.py,sha256=IuUkKrEoz_EKLCmGQ94CLFgjKPvO20iVx4773JTKO-8,6503
-h2lib_tests/test_topology_h2lib.py,sha256=7ZfVjEUYBJacPcFCDi_Qb0LqW6W-0no_zTBUKbNcKQw,23858
+h2lib_tests/test_topology_h2lib.py,sha256=g6H4E5OR4b2NlCDTZkOkUds2NyA-9XSG9Vu9Cy7TJw8,40348
 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
@@ -91,7 +91,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.2.701.dist-info/METADATA,sha256=-mXXeU5VEAY9tlTcDTGdIrr3YjPISLCZsZpFATvR-cg,417
-h2lib_tests-13.2.701.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-h2lib_tests-13.2.701.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
-h2lib_tests-13.2.701.dist-info/RECORD,,
+h2lib_tests-13.2.801.dist-info/METADATA,sha256=c2VITTY0PxlaMfuXZrhOsoTeB8MVqNudNxuApGiQAPI,441
+h2lib_tests-13.2.801.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+h2lib_tests-13.2.801.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.2.801.dist-info/RECORD,,