h2lib-tests 13.1.303__py3-none-any.whl → 13.1.501__py3-none-any.whl

h2lib_tests/test_static_solver_wrapper.py

@@ -0,0 +1,20 @@
+import pytest
+
+from h2lib._h2lib import H2Lib
+from h2lib_tests.test_files import tfp
+
+from .test_write_htc import write_dtu10mw_only_tower
+
+
+def test_init_static_solver(write_dtu10mw_only_tower):
+    with H2Lib() as h2:
+        # Load a basic model.
+        model_path = tfp + "DTU_10_MW/"
+        h2.init("htc/DTU_10MW_RWT_only_tower.htc", model_path)
+
+        # First execution is fine.
+        h2.init_static_solver()
+
+        # The next should give an error.
+        with pytest.raises(RuntimeError, match="STATIC_SOLVER_ALREADY_INITIALIZED"):
+            h2.init_static_solver()

h2lib_tests/test_topology_h2lib.py

@@ -0,0 +1,143 @@
+from h2lib._h2lib import H2Lib
+
+from numpy import testing as npt
+from h2lib_tests.test_files import tfp
+import numpy as np
+import pytest
+from .test_write_htc import (
+    write_dtu10mw_only_tower,
+    write_dtu10mw_only_tower_rotated,
+    write_dtu10mw_only_tower_encrypted
+)
+
+
+@pytest.fixture(scope="module")
+def h2_dtu_10mw_only_tower(write_dtu10mw_only_tower):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="module")
+def h2_dtu_10mw_only_tower_rotated(write_dtu10mw_only_tower_rotated):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower_rotated.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="module")
+def h2_dtu_10mw_only_tower_encrypted(write_dtu10mw_only_tower_encrypted):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower_encrypted.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    yield h2
+    h2.close()
+
+
+def test_number_of_bodies_and_constraints(
+    h2_dtu_10mw_only_tower,
+):
+    nbdy, ncst = h2_dtu_10mw_only_tower.get_number_of_bodies_and_constraints()
+    assert nbdy == 3
+    assert ncst == 9
+
+
+def test_number_of_bodies_and_constraints_encrypted(
+    h2_dtu_10mw_only_tower_encrypted,
+):
+    with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
+        h2_dtu_10mw_only_tower_encrypted.get_number_of_bodies_and_constraints()
+
+
+def test_get_number_of_elements(h2_dtu_10mw_only_tower):
+    nelem = h2_dtu_10mw_only_tower.get_number_of_elements()
+    npt.assert_array_equal(nelem, np.array([3, 3, 4]))
+
+
+def test_get_number_of_elements_encrypted(
+    h2_dtu_10mw_only_tower_encrypted,
+):
+    # This test is not really needed, since the check for confidential structure
+    # is already done by test_number_of_bodies_and_constraints_encrypted().
+    with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
+        h2_dtu_10mw_only_tower_encrypted.get_number_of_elements()
+
+
+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)
+    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]))
+    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]]),
+    )
+
+    # Test last element.
+    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_equal(
+        tes,
+        np.array([[-1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, -1.0]]),
+    )
+
+
+def test_get_timoshenko_location_body_does_not_exist(
+    h2_dtu_10mw_only_tower,
+):
+    with pytest.raises(IndexError, match="BODY_DOES_NOT_EXIST"):
+        h2_dtu_10mw_only_tower.get_timoshenko_location(ibdy=1000, ielem=0)
+
+
+def test_get_timoshenko_location_element_does_not_exist(
+    h2_dtu_10mw_only_tower,
+):
+    with pytest.raises(IndexError, match="ELEMENT_DOES_NOT_EXIST"):
+        h2_dtu_10mw_only_tower.get_timoshenko_location(ibdy=0, ielem=1000)
+
+
+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)
+
+
+def test_get_body_rotation_tensor_1(h2_dtu_10mw_only_tower):
+    amat = h2_dtu_10mw_only_tower.get_body_rotation_tensor(ibdy=0)
+    npt.assert_array_equal(amat, np.eye(3))
+
+
+def test_get_body_rotation_tensor_2(
+    h2_dtu_10mw_only_tower_rotated, write_dtu10mw_only_tower_rotated
+):
+    amat = h2_dtu_10mw_only_tower_rotated.get_body_rotation_tensor(ibdy=0)
+    _, alpha = write_dtu10mw_only_tower_rotated
+    alpha_rad = np.deg2rad(alpha)
+    sa = np.sin(alpha_rad)
+    ca = np.cos(alpha_rad)
+    npt.assert_array_almost_equal_nulp(
+        amat, np.array([[1.0, 0.0, 0.0], [0.0, ca, -sa], [0.0, sa, ca]])
+    )
+
+
+def test_get_body_rotation_tensor_body_does_not_exist(h2_dtu_10mw_only_tower):
+    with pytest.raises(IndexError, match="BODY_DOES_NOT_EXIST"):
+        h2_dtu_10mw_only_tower.get_body_rotation_tensor(ibdy=1000)
+
+
+def test_get_body_rotation_tensor_encrypted(h2_dtu_10mw_only_tower_encrypted):
+    with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
+        h2_dtu_10mw_only_tower_encrypted.get_body_rotation_tensor(ibdy=0)

h2lib_tests/test_write_htc.py

@@ -0,0 +1,59 @@
+import pytest
+from wetb.hawc2.htc_file import HTCFile
+from h2lib_tests.test_files import tfp
+
+
+@pytest.fixture(scope="module")
+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.set_name("DTU_10MW_RWT_only_tower")
+    for key1 in htc["new_htc_structure"].keys():
+        if key1.startswith("main_body"):
+            if "tower" not in htc["new_htc_structure"][key1]["name"].values:
+                htc["new_htc_structure"][key1].delete()
+        if key1 == "orientation":
+            for key2 in htc["new_htc_structure"]["orientation"].keys():
+                if key2.startswith("relative"):
+                    htc["new_htc_structure"]["orientation"][key2].delete()
+        if key1 == "constraint":
+            for key2 in htc["new_htc_structure"]["constraint"].keys():
+                if key2 != "fix0":
+                    htc["new_htc_structure"]["constraint"][key2].delete()
+    htc["wind"].delete()
+    htc["aerodrag"].delete()
+    htc["aero"].delete()
+    htc["dll"].delete()
+    htc["output"].delete()
+    # Reduce simulation time.
+    htc.simulation.time_stop = 10.0
+    # Change number of bodies in the tower.
+    htc.new_htc_structure.main_body.nbodies = 3
+    # Save the new file.
+    htc.save()
+    return htc
+
+
+@pytest.fixture(scope="module")
+def write_dtu10mw_only_tower_rotated(write_dtu10mw_only_tower):
+    # Start from the DTU_10MW_RWT_only_tower and rotate the tower.
+    htc = write_dtu10mw_only_tower.copy()
+    htc.set_name("DTU_10MW_RWT_only_tower_rotated")
+    alpha = 30.0
+    htc.new_htc_structure.orientation.base.body_eulerang = [
+        alpha,
+        0.0,
+        0.0,
+    ]
+    htc.save()
+    return (htc, alpha)
+
+
+@pytest.fixture(scope="module")
+def write_dtu10mw_only_tower_encrypted(write_dtu10mw_only_tower):
+    # Start from the DTU_10MW_RWT_only_tower and then encrypt the tower.
+    htc = write_dtu10mw_only_tower.copy()
+    htc.set_name("DTU_10MW_RWT_only_tower_encrypted")
+    # Only the tower is left.
+    htc.new_htc_structure.main_body.timoschenko_input.filename = "./data/DTU_10MW_RWT_Tower_st.dat.v3.enc"
+    htc.save()

{h2lib_tests-13.1.303.dist-info → h2lib_tests-13.1.501.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.1
 Name: h2lib-tests
-Version: 13.1.303
-Summary: Tests and test_files for test h2lib (13.1.3)
+Version: 13.1.501
+Summary: Tests and test_files for test h2lib (13.1.5+2-g1058cd2)
 Download-URL: 
 Author: Mads M Pedersen
 Author-email: 

{h2lib_tests-13.1.303.dist-info → h2lib_tests-13.1.501.dist-info}/RECORD

@@ -6,6 +6,9 @@ h2lib_tests/test_h2lib.py,sha256=ZKp8pfQQpCUQ7IBXypTmvewbJwTQowSiCEt6t4R9TgY,133
 h2lib_tests/test_h2rotor.py,sha256=kaE9AFDYSPaMFSgBAbeegd7-l6Z_4_BsQaWie0k32q4,6210
 h2lib_tests/test_mpi.py,sha256=emBgXRAvvFFOsVrAziiQCUZvEF9HS5Wc-x4KqqltQP0,6835
 h2lib_tests/test_multiprocessinterface.py,sha256=h2o4havtK6IuMXsplNjGUa3VxOnbpEYGxdrrAKQilj0,1470
+h2lib_tests/test_static_solver_wrapper.py,sha256=xCWDkxtW8K2kD6l6MRgyPS1kz91IpbRA0Cc3BG7o3Es,604
+h2lib_tests/test_topology_h2lib.py,sha256=cClNtAIc1wR_gZa7eDO4ZySvhScqCEIHVMKZEhk29QE,4980
+h2lib_tests/test_write_htc.py,sha256=t78STsVoQ2qR5WB-jkYhaMOy0jWHHsDfuO2nfyPj4ag,2242
 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
@@ -28,6 +31,7 @@ h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Blade_st.dat,sha256=7stkwSf7R
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Hub_st.dat,sha256=UPIVkaThyxRsEeMIyjhXIr_n9kr1SAGQYudgj_ofJ04,916
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Shaft_st.dat,sha256=Jz_QA1hJUY6ZG4re1oFj0_7MkXu8az9jn2-WLGb0jY0,905
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Tower_st.dat,sha256=0asFvlyhL1CZltWve-PNGgjfpfdMGRCCwMf3U0K_x_c,8367
+h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Tower_st.dat.v3.enc,sha256=2MIV4QV0Y5XZ6hHsUdO95pG4b1HVlNWeuIzSM6n-88A,8422
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Towertop_st.dat,sha256=iW8RzjgQe7D3y9DEQMD74e0PWqJ4d7meYGTT6Fsl4D0,955
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_ae.dat,sha256=LkAlElRW_Fl1V4GD8hQZIgLD4IW0o3uxXL2olnEU3dI,1531
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_pc.dat,sha256=2Nc_ofcqfCB3iu8dSFRj8Ub6AhaQADdW0O8BXBWj5Uk,26888
@@ -79,7 +83,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.303.dist-info/METADATA,sha256=u93e4yUJ_aWSDJe7juA3_c-V6d-rIT_gsOgYutazCjI,349
-h2lib_tests-13.1.303.dist-info/WHEEL,sha256=R0nc6qTxuoLk7ShA2_Y-UWkN8ZdfDBG2B6Eqpz2WXbs,91
-h2lib_tests-13.1.303.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
-h2lib_tests-13.1.303.dist-info/RECORD,,
+h2lib_tests-13.1.501.dist-info/METADATA,sha256=-kMUlvOUDYWvd8sU6fjj7n0bG-MWPRHm4Sz4fb4qrcI,360
+h2lib_tests-13.1.501.dist-info/WHEEL,sha256=cVxcB9AmuTcXqmwrtPhNK88dr7IR_b6qagTj0UvIEbY,91
+h2lib_tests-13.1.501.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.1.501.dist-info/RECORD,,

{h2lib_tests-13.1.303.dist-info → h2lib_tests-13.1.501.dist-info}/WHEEL

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