h2lib-tests 13.1.505__py3-none-any.whl → 13.1.506__py3-none-any.whl

h2lib_tests/test_static_solver.py

@@ -0,0 +1,98 @@
+# -*- coding: utf-8 -*-
+"""
+Test the static solver.
+
+@author: ricriv
+"""
+
+# Here we are relying on the default behavior of pytest, which is to execute
+# the tests in the same order that they are specified.
+# If one day this will not be the case anymore, we can enforce the order by
+# using the solution proposed at: https://stackoverflow.com/a/77793427/3676517
+
+import pytest
+
+import numpy as np
+from numpy import testing as npt
+
+from h2lib._h2lib import H2Lib
+from h2lib_tests.test_files import tfp
+
+from .test_write_htc import (
+    write_dtu10mw_only_blade,
+    write_dtu10mw_only_blade_low_max_iter,
+)
+
+
+@pytest.fixture(scope="module")
+def h2_dtu_10mw_only_blade(write_dtu10mw_only_blade):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_blade.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="module")
+def h2_dtu10mw_only_blade_low_max_iter(write_dtu10mw_only_blade_low_max_iter):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_blade_low_max_iter.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    yield h2
+    h2.close()
+
+
+def test_solver_static_update_no_init(h2_dtu_10mw_only_blade):
+    with pytest.raises(RuntimeError, match="STATIC_SOLVER_NOT_INITIALIZED"):
+        h2_dtu_10mw_only_blade.solver_static_update()
+
+
+def test_solver_static_solve_no_init(h2_dtu_10mw_only_blade):
+    with pytest.raises(RuntimeError, match="STATIC_SOLVER_NOT_INITIALIZED"):
+        h2_dtu_10mw_only_blade.solver_static_solve()
+
+
+def test_solver_static_init(h2_dtu_10mw_only_blade):
+
+    # First execution is fine.
+    h2_dtu_10mw_only_blade.solver_static_init()
+
+    # The next should automatically deallocate the static solver and initialize it again.
+    h2_dtu_10mw_only_blade.solver_static_init()
+
+
+def test_solver_static_update(h2_dtu_10mw_only_blade):
+    # This should happen after test_solver_static_init().
+    h2_dtu_10mw_only_blade.solver_static_update()
+
+
+def test_solver_static_solve(h2_dtu_10mw_only_blade):
+    # This should happen after test_solver_static_update().
+    h2_dtu_10mw_only_blade.solver_static_solve()
+
+
+def test_solver_static_delete(h2_dtu_10mw_only_blade):
+    h2_dtu_10mw_only_blade.solver_static_delete()
+
+
+def test_static_solver_run(h2_dtu_10mw_only_blade):
+    # Add a sensor for the blade root moment, in this case only due to gravity.
+    id = h2_dtu_10mw_only_blade.add_sensor("mbdy momentvec blade1 1 1 blade1")
+
+    # Run the static solver.
+    h2_dtu_10mw_only_blade.solver_static_run()
+
+    # Do 1 step to get the output.
+    h2_dtu_10mw_only_blade.step()
+    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.07213851e04, -4.55385871e-02, -3.94623708e01])
+    )
+
+
+def test_static_solver_run_fail(h2_dtu10mw_only_blade_low_max_iter):
+    with pytest.raises(RuntimeError, match="STATIC_SOLVER_DID_NOT_CONVERGE"):
+        h2_dtu10mw_only_blade_low_max_iter.solver_static_run()

h2lib_tests/test_write_htc.py

@@ -57,3 +57,70 @@ def write_dtu10mw_only_tower_encrypted(write_dtu10mw_only_tower):
     # 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()
+
+
+@pytest.fixture(scope="module")
+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.set_name("DTU_10MW_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 as many bodies as elements.
+    htc.new_htc_structure.main_body__7.nbodies = 26
+
+    # Reduce simulation time to 1 time step.
+    htc.simulation.time_stop = 0.01
+    htc.simulation.log_deltat.delete()
+
+    # Do not use static solver, since it will be done during the test.
+    htc.simulation.solvertype = 2
+    htc.simulation.solvertype.comments = ""
+    htc.simulation.initial_condition = 1
+
+    # No output, as we will use add_sensor().
+    # htc.add_section("output")
+    # htc.output.data_format = "gtsdf"
+    # htc.output.buffer = 10000
+
+    # Save the new file.
+    htc.save()
+
+    return htc
+
+
+@pytest.fixture(scope="module")
+def write_dtu10mw_only_blade_low_max_iter(write_dtu10mw_only_blade):
+    # Start from the write_dtu10mw_only_blade and thenreduce the number of max iterations,
+    # so that the static solver will not have time to converge.
+    htc = write_dtu10mw_only_blade.copy()
+    htc.set_name("DTU_10MW_RWT_only_blade_low_max_iter")
+    htc.simulation.max_iterations = 1
+    htc.save()
+
+    return htc

{h2lib_tests-13.1.505.dist-info → h2lib_tests-13.1.506.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.1
 Name: h2lib-tests
-Version: 13.1.505
-Summary: Tests and test_files for test h2lib (13.1.5+4-g20eca1a)
+Version: 13.1.506
+Summary: Tests and test_files for test h2lib (13.1.8)
 Download-URL: 
 Author: Mads M Pedersen
 Author-email: 

{h2lib_tests-13.1.505.dist-info → h2lib_tests-13.1.506.dist-info}/RECORD

@@ -6,9 +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_static_solver.py,sha256=Ik_CgGiFVibEPq0dGzSfSYH3gv_dTm-qcz6JKmJiH_A,3228
 h2lib_tests/test_topology_h2lib.py,sha256=cClNtAIc1wR_gZa7eDO4ZySvhScqCEIHVMKZEhk29QE,4980
-h2lib_tests/test_write_htc.py,sha256=t78STsVoQ2qR5WB-jkYhaMOy0jWHHsDfuO2nfyPj4ag,2242
+h2lib_tests/test_write_htc.py,sha256=3KMHffkR8A6C1RXODxq_O522fEBfJwp3YpjDmUmiNjU,4774
 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
@@ -87,7 +87,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.505.dist-info/METADATA,sha256=UF9ZgIIC_UKSdlCODuGG3Ao3XX3plr0A2tn_-Pnebz0,360
-h2lib_tests-13.1.505.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
-h2lib_tests-13.1.505.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
-h2lib_tests-13.1.505.dist-info/RECORD,,
+h2lib_tests-13.1.506.dist-info/METADATA,sha256=wvzBrJm0RO6QivFp3_Hx4Pol3PqGCk_9WyxwgOtyuao,349
+h2lib_tests-13.1.506.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
+h2lib_tests-13.1.506.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.1.506.dist-info/RECORD,,

h2lib_tests/test_static_solver_wrapper.py

@@ -1,20 +0,0 @@
-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()