h2lib-tests 13.1.3101__py3-none-any.whl → 13.1.3102__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='H2LibThread process died before or while executing set_aerosections_windspeed'):
+            h2.set_aerosections_windspeed(uvw)
+        out = capfd.readouterr().out
+        assert "Please call init_ad or init_al be" 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_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-13.1.3101.dist-info → h2lib_tests-13.1.3102.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.3102
+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.3102.dist-info}/RECORD

@@ -2,11 +2,11 @@ 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_ellipsys_couplings.py,sha256=cAo036WqDkwksCX0EXdUUoJ9mJHjpkUPuIlKLW7GnzY,3284
 h2lib_tests/test_h2lib.py,sha256=opgVOkZWOD9CrnTDfwRrBaD-weBsNnutm6FeOijDa88,14197
-h2lib_tests/test_h2rotor.py,sha256=ZMhTdGvRYAORULd18X2woEDFrGB29DU-K2J2VA6sE-8,11262
+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
@@ -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.3102.dist-info/METADATA,sha256=FgPhvGVItQe7EYP_eXjmrMXOnxfT3yXxTI4ALJgN99I,419
+h2lib_tests-13.1.3102.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
+h2lib_tests-13.1.3102.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.1.3102.dist-info/RECORD,,