h2lib-tests 13.1.3101__py3-none-any.whl → 13.1.3103__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='Please call init_AD or init_AL before set_aerosections_windspeed'):
+            h2.set_aerosections_windspeed(uvw)
+        out = capfd.readouterr().out
+        assert "Please call init_AD or init_AL before set_aerosections_windspeed" 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_h2lib.py

@@ -301,7 +301,7 @@ def test_process_died():
             h2.echo_version()
 
 
-def test_error():
+def test_error(capfd):
     model_path = tfp + "minimal/"
     htc = HTCFile(model_path + 'htc/minimal_mann_turb.htc')
     htc.wind.scale_time_start = 200
@@ -310,44 +310,39 @@ def test_error():
     htc.set_name('tmp')
     htc.save()
 
-    with H2Lib() as h2:
-        with pytest.raises(Exception, match=re.escape('H2LibThread process died before or while executing init(htc/tmp.htc, ')):
+    with H2Lib(suppress_output=0) as h2:
+        version = h2.get_version()
+        with pytest.raises(Exception, match='Turbulence scale_time_start >= simulation length'):
             h2.init('htc/tmp.htc', model_path)
+        out, err = capfd.readouterr()
+        assert "*** ERROR *** Turbulence scale_time_start >= simulation length" in out
+        assert h2.get_version() == version  # hawc2 still alive
+    with pytest.raises(ValueError):
+        h2.get_version()  # now it is closed
 
     with MultiH2Lib(2) as h2:
-        with pytest.raises(Exception, match=re.escape('H2LibThread process died before or while executing init(...)')):
+        version = h2.get_version()
+        with pytest.raises(Exception, match='Turbulence scale_time_start >= simulation length'):
             h2.init('htc/tmp.htc', model_path)
+        assert h2.get_version() == version  # hawc2 still alive
 
 
 def test_fail():
     with H2Lib() as h2:
-        with pytest.raises(Exception, match=re.escape('H2LibThread process died before or while executing fail(MyError)')):
+        version = h2.get_version()
+        with pytest.raises(Exception, match='MyError'):
             h2.fail('MyError')
+        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):
+        with pytest.raises(Exception, match='MyError'):
             h2.fail('MyError')
-
-
-def test_dont_stop(capfd):
-    # This test is the same as test_error(), but instead we disable stop.
-    model_path = tfp + "minimal/"
-    htc = HTCFile(model_path + 'htc/minimal_mann_turb.htc')
-    htc.wind.scale_time_start = 200
-    htc.wind.turb_format = 1
-    htc.wind.mann.dont_scale = 0
-    htc.set_name('tmp')
-    htc.save()
-
-    with H2Lib() as h2:
-        h2.stop_on_error(False)
-        h2.init('htc/tmp.htc', model_path)
-        out, err = capfd.readouterr()
-        assert "*** ERROR *** Turbulence scale_time_start >= simulation length" in out
-
-    with MultiH2Lib(2) as h2:
-        h2.stop_on_error(False)
-        h2.init('htc/tmp.htc', model_path)
-        out, err = capfd.readouterr()
-        assert "*** ERROR *** Turbulence scale_time_start >= simulation length" in out

h2lib_tests/test_h2rotor.py

@@ -14,7 +14,7 @@ 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
 
@@ -123,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:
@@ -154,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)
@@ -165,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()
@@ -180,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()
 
@@ -235,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)
 
@@ -244,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_mpi.py

@@ -146,7 +146,7 @@ def test_ellipsys_mpi_dummy_workflow():
         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()
@@ -176,7 +176,7 @@ def test_ellipsys_mpi_dummy_workflow_collective():
             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()

h2lib_tests/test_topology_h2lib.py

@@ -189,9 +189,7 @@ def test_set_orientation_relative_main_body_not_found(
 ):
     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(

{h2lib_tests-13.1.3101.dist-info → h2lib_tests-13.1.3103.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.4
 Name: h2lib-tests
-Version: 13.1.3101
-Summary: Tests and test_files for test h2lib (13.1.31+1-g14e98da)
+Version: 13.1.3103
+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.3101.dist-info → h2lib_tests-13.1.3103.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_h2lib.py,sha256=opgVOkZWOD9CrnTDfwRrBaD-weBsNnutm6FeOijDa88,14197
-h2lib_tests/test_h2rotor.py,sha256=ZMhTdGvRYAORULd18X2woEDFrGB29DU-K2J2VA6sE-8,11262
+h2lib_tests/test_ellipsys_couplings.py,sha256=JLos66l_zNDBY6glxwqSfDa6p0Uko5aRKAckh2uoSSM,3302
+h2lib_tests/test_h2lib.py,sha256=LEV1n0Fwa1llQwIM06Mcqr07Y6WOn-j8mNxjzlCEOl0,14219
+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=l3V6fts0NJ2Vmsm7Tm5eYbOtewUyrcQurdaLJVITALM,6984
+h2lib_tests/test_mpi.py,sha256=8bMeukbWboijMgWLJIIu3Pn0QtaDadOGhhL099MfMls,6990
 h2lib_tests/test_multiprocessinterface.py,sha256=h2o4havtK6IuMXsplNjGUa3VxOnbpEYGxdrrAKQilj0,1470
 h2lib_tests/test_static_solver.py,sha256=ckzcBAuVge5V5tdrac4PPYxrH5TheNTUc3PKksbB89A,7426
-h2lib_tests/test_topology_h2lib.py,sha256=VjJ1SUBHRHwuSn0EPI6KGzwAULZbPzA39rXdDGXyxkk,11510
+h2lib_tests/test_topology_h2lib.py,sha256=U1NOREhoJ9m7cYxb8sLsVuKSHRqjkEde6hoyl88IF4A,11486
 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.3101.dist-info/METADATA,sha256=oqZnjRxb1_2zZvCrPkhXUd9kJBvUb4AiVHf8m1UYclU,419
-h2lib_tests-13.1.3101.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
-h2lib_tests-13.1.3101.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
-h2lib_tests-13.1.3101.dist-info/RECORD,,
+h2lib_tests-13.1.3103.dist-info/METADATA,sha256=2hyzPnwoJPCSkngKuqGVx5i49g2Kl_FRGaSA5OYxuek,419
+h2lib_tests-13.1.3103.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
+h2lib_tests-13.1.3103.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.1.3103.dist-info/RECORD,,