h2lib-tests 13.1.2301__py3-none-any.whl → 13.1.2302__py3-none-any.whl

h2lib_tests/conftest.py

@@ -105,6 +105,9 @@ def write_dtu10mw_only_blade():
     htc.simulation.solvertype.comments = ""
     htc.simulation.initial_condition = 1
 
+    # 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.
@@ -206,6 +209,9 @@ def write_dtu10mw_only_blade_rotate_relative():
     htc.simulation.solvertype.comments = ""
     htc.simulation.initial_condition = 1
 
+    # Set low convergence limits.
+    htc.simulation.convergence_limits = [1e2, 1e-5, 1e-07]
+
     # Save the new file.
     htc.save()
 

h2lib_tests/test_h2rotor.py

@@ -157,7 +157,10 @@ 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
+    npt.assert_array_equal(h2.get_rotor_avg_wsp(), [np.nan, np.nan, np.nan])
     h2.step()
     frc_after = h2.get_aerosections_forces()
     mom_after = h2.get_aerosections_moments()

h2lib_tests/test_static_solver.py

@@ -63,58 +63,65 @@ def test_static_solver_run_1(h2_dtu_10mw_only_blade):
     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()
+    h2_dtu_10mw_only_blade.solver_static_run(reset_structure=True)
 
     # 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])
+        val, np.array([-1.072139e+04, -3.856623e-02, -3.949042e+01]),
+        rtol=1e-6
     )
 
 
 def test_static_solver_run_2(h2_dtu_10mw_only_blade_rotate_base):
     # Apply centrifugal loading with the base command.
     # Add a sensor for the blade root force, in this case only due to centrifugal force.
-    id = h2_dtu_10mw_only_blade_rotate_base.add_sensor("mbdy forcevec blade1 1 1 blade1")
+    id = h2_dtu_10mw_only_blade_rotate_base.add_sensor(
+        "mbdy forcevec blade1 1 1 blade1"
+    )
 
     # Run the static solver.
-    h2_dtu_10mw_only_blade_rotate_base.solver_static_run()
+    h2_dtu_10mw_only_blade_rotate_base.solver_static_run(reset_structure=True)
 
     # Do 1 step to get the output.
     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.129998, 383.43168, 989.68145])
-    )
+    npt.assert_allclose(val, np.array([10879.120647, 383.441758, 989.673212]))
 
 
 def test_static_solver_run_3(h2_dtu_10mw_only_blade_rotate_relative):
     # Apply centrifugal loading with the relative command.
     # Add a sensor for the blade root force, in this case only due to centrifugal force.
-    id = h2_dtu_10mw_only_blade_rotate_relative.add_sensor("mbdy forcevec blade1 1 1 blade1")
+    id = h2_dtu_10mw_only_blade_rotate_relative.add_sensor(
+        "mbdy forcevec blade1 1 1 blade1"
+    )
 
     # Run the static solver.
-    h2_dtu_10mw_only_blade_rotate_relative.solver_static_run()
+    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.083576, 383.430039, 989.67935])
-    )
+    npt.assert_allclose(val, np.array([10879.074041, 383.440111, 989.67438]))
 
 
 def test_static_solver_run_4(h2_dtu_10mw_only_blade_rotate_bearing3):
     # Apply centrifugal loading with the bearing3 command.
     # Add a sensor for the blade root moment, in this case only due to centrifugal force.
-    id = h2_dtu_10mw_only_blade_rotate_bearing3.add_sensor("mbdy momentvec blade1 1 1 blade1")
+    id = h2_dtu_10mw_only_blade_rotate_bearing3.add_sensor(
+        "mbdy momentvec blade1 1 1 blade1"
+    )
 
     # Run the static solver.
-    h2_dtu_10mw_only_blade_rotate_bearing3.solver_static_run()
+    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()
@@ -123,3 +130,73 @@ def test_static_solver_run_4(h2_dtu_10mw_only_blade_rotate_bearing3):
     npt.assert_allclose(
         val, np.array([-3094.986918, 115414.095937, 325.296806])
     )
+
+
+def test_static_solver_run_no_reset(h2_dtu_10mw_only_blade):
+    # Use gravity to deflect the clamped blade.
+
+    # Run the static solver.
+    h2_dtu_10mw_only_blade.solver_static_run(reset_structure=True)
+
+    # Run it again, without resetting the structure.
+    # The static solver must exit immediately and output the same residuals.
+    _, resq_1, resg_1, resd_1 = h2_dtu_10mw_only_blade.check_convergence()
+    h2_dtu_10mw_only_blade.solver_static_run(reset_structure=False)
+    _, resq_2, resg_2, resd_2 = h2_dtu_10mw_only_blade.check_convergence()
+
+    npt.assert_allclose(
+        np.array([resq_2, resg_2, resd_2]), np.array([resq_2, resg_2, resd_2])
+    )
+
+
+def test_static_solver_run_with_reset(h2_dtu_10mw_only_blade):
+    # Use gravity to deflect the clamped blade.
+
+    # Run the static solver.
+    h2_dtu_10mw_only_blade.solver_static_run(reset_structure=True)
+
+    # Run it again, but first reset the structure.
+    # The static solver must follow the same convergence history.
+    # Note that here we are only checking the last value of the residuals.
+    _, resq_1, resg_1, resd_1 = h2_dtu_10mw_only_blade.check_convergence()
+    h2_dtu_10mw_only_blade.solver_static_run(reset_structure=True)
+    _, resq_2, resg_2, resd_2 = h2_dtu_10mw_only_blade.check_convergence()
+
+    npt.assert_allclose(
+        np.array([resq_2, resg_2, resd_2]), np.array([resq_2, resg_2, resd_2])
+    )
+
+
+def test_structure_reset():
+    # Apply centrifugal loading with the base command.
+    # We cannot use h2_dtu_10mw_only_blade_rotate_base because we need to add a sensor,
+    # and this must be done before init.
+
+    with H2Lib(suppress_output=True) as h2:
+        model_path = f"{tfp}DTU_10_MW/"
+        htc_path = "htc/DTU_10MW_RWT_only_blade_rotate_base.htc"
+        id = h2.add_sensor(
+            "mbdy statevec_new blade1 c2def global absolute 85.0 1.0 0.0 0.0"
+        )
+        h2.init(htc_path=htc_path, model_path=model_path)
+        h2.stop_on_error(False)
+
+        # Do 1 step to get the output. This is our reference.
+        h2.step()
+        val_desired = h2.get_sensor_values(id)
+
+        # Run static solver.
+        h2.solver_static_run(reset_structure=False)
+
+        # 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
+        )
+
+        # Reset the structure and check that we match the reference.
+        h2.structure_reset()
+        h2.step()
+        val_actual = h2.get_sensor_values(id)
+        npt.assert_allclose(val_actual, val_desired)

h2lib_tests/test_topology_h2lib.py

@@ -172,16 +172,27 @@ def test_set_orientation_base_without_reset_orientation(
     )
 
 
-def test_set_orientation_base_speed(
-    h2_dtu_10mw_only_blade
-):
+def test_set_orientation_base_speed(h2_dtu_10mw_only_blade):
+    # Add a sensor for the blade root force, in this case only due to centrifugal force.
+    id = h2_dtu_10mw_only_blade.add_sensor("mbdy forcevec blade1 1 1 blade1")
+
     # Set speed.
     h2_dtu_10mw_only_blade.set_orientation_base(
         main_body="blade1",
         reset_orientation=False,
         mbdy_ini_rotvec_d1=np.array([0.0, 1.0, 0.0, 1.0]),
     )
-    # TODO: check against h2_dtu_10mw_only_blade_rotate_base.
+    # Run static solver.
+    h2_dtu_10mw_only_blade.solver_static_run(reset_structure=True)
+
+    # 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: 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]))
+
     # Reset speed.
     h2_dtu_10mw_only_blade.set_orientation_base(
         main_body="blade1",
@@ -189,15 +200,23 @@ def test_set_orientation_base_speed(
     )
 
 
-def test_set_orientation_relative_main_body_not_found(h2_dtu_10mw_only_blade_rotate_relative):
+def test_set_orientation_relative_main_body_not_found(
+    h2_dtu_10mw_only_blade_rotate_relative,
+):
     h2_dtu_10mw_only_blade_rotate_relative.stop_on_error(False)
     with pytest.raises(ValueError, match="MAIN_BODY_NOT_FOUND"):
-        h2_dtu_10mw_only_blade_rotate_relative.set_orientation_relative("hub1", "last", "blade", 0)
+        h2_dtu_10mw_only_blade_rotate_relative.set_orientation_relative(
+            "hub1", "last", "blade", 0
+        )
 
 
-def test_set_orientation_relative_rot_not_found(h2_dtu_10mw_only_blade_rotate_relative):
+def test_set_orientation_relative_rot_not_found(
+    h2_dtu_10mw_only_blade_rotate_relative,
+):
     with pytest.raises(ValueError, match="RELATIVE_ROTATION_NOT_FOUND"):
-        h2_dtu_10mw_only_blade_rotate_relative.set_orientation_relative("hub1", "last", "blade1", "last")
+        h2_dtu_10mw_only_blade_rotate_relative.set_orientation_relative(
+            "hub1", "last", "blade1", "last"
+        )
 
 
 def test_set_orientation_relative_reset(
@@ -273,3 +292,43 @@ def test_set_orientation_relative_2(
         h2_dtu_10mw_only_blade_rotate_relative.get_body_rotation_tensor(ibdy=1)
     )
     npt.assert_array_almost_equal_nulp(amat_actual, amat_desired)
+
+
+def test_set_orientation_relative_static(
+    h2_dtu_10mw_only_blade_rotate_relative,
+):
+    # Add a sensor for the blade root force.
+    id = h2_dtu_10mw_only_blade_rotate_relative.add_sensor(
+        "mbdy forcevec blade1 1 1 blade1"
+    )
+
+    # Set arbitrary orientation and speed.
+    h2_dtu_10mw_only_blade_rotate_relative.set_orientation_relative(
+        main_body_1="hub1",
+        node_1="last",
+        main_body_2="blade1",
+        node_2=0,
+        reset_orientation=True,
+        mbdy2_eulerang_table=np.array([-80.0, 0.0, 0.0]),
+        mbdy2_ini_rotvec_d1=[0.0, 1.0, 0.0, 0.8],
+    )
+    # Run static solver.
+    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([8702.239718, 306.617339, 638.840412]))
+
+    # Reset to original value.
+    h2_dtu_10mw_only_blade_rotate_relative.set_orientation_relative(
+        main_body_1="hub1",
+        node_1="last",
+        main_body_2="blade1",
+        node_2=0,
+        mbdy2_eulerang_table=np.array([-90.0, 0.0, 0.0]),
+        reset_orientation=True,
+        mbdy2_ini_rotvec_d1=np.array([0.0, 1.0, 0.0, 1.0]),
+    )

{h2lib_tests-13.1.2301.dist-info → h2lib_tests-13.1.2302.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.2
 Name: h2lib-tests
-Version: 13.1.2301
-Summary: Tests and test_files for test h2lib (13.1.23+5-gcdb6f43)
+Version: 13.1.2302
+Summary: Tests and test_files for test h2lib (13.1.25)
 Download-URL: 
 Author: Mads M Pedersen
 Author-email: 

{h2lib_tests-13.1.2301.dist-info → h2lib_tests-13.1.2302.dist-info}/RECORD

@@ -1,15 +1,15 @@
 h2lib_tests/__init__.py,sha256=VjSqfGg8BzdmSjfSFhJh4hZbYZ_cME7xp9EWFKHQphA,61
-h2lib_tests/conftest.py,sha256=uJrYYp-R7bVpBhfdzJH0CUrshS-S1_hAzST36JJZNbo,14559
+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_h2lib.py,sha256=opgVOkZWOD9CrnTDfwRrBaD-weBsNnutm6FeOijDa88,14197
-h2lib_tests/test_h2rotor.py,sha256=SP4NTByJGJDc9yezR9QdijYlsJjMJ4BwItbhksijL1U,10597
+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_multiprocessinterface.py,sha256=h2o4havtK6IuMXsplNjGUa3VxOnbpEYGxdrrAKQilj0,1470
-h2lib_tests/test_static_solver.py,sha256=OiabTF7IwIskkvBslp32c5o4v53YHQNTTV3vOU5vVb8,4679
-h2lib_tests/test_topology_h2lib.py,sha256=AIPgO6ePVoZbjX_ECmlQu8Iog1dIBzd192hkXBvg1V0,9651
+h2lib_tests/test_static_solver.py,sha256=mGapb1vwSoe0Uqpc3-xjDegqauf2R_XMHshnvVojwmY,7509
+h2lib_tests/test_topology_h2lib.py,sha256=dCwvfgO1eHb53PCXwm-nr6iwKdsviU2tuAnU0WIrsJc,11645
 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.2301.dist-info/METADATA,sha256=OUugVXDBqUPSq4FO51fyoNuSbB_kT1ucDZA3BfpGsZU,419
-h2lib_tests-13.1.2301.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
-h2lib_tests-13.1.2301.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
-h2lib_tests-13.1.2301.dist-info/RECORD,,
+h2lib_tests-13.1.2302.dist-info/METADATA,sha256=2xAxSCnVwrOtfNcEru0fQ9fCtEhYR9ZLmvkfTvF5_zk,408
+h2lib_tests-13.1.2302.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
+h2lib_tests-13.1.2302.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.1.2302.dist-info/RECORD,,