h2lib 13.0.408__cp38-cp38-win_amd64.whl → 13.0.410__cp38-cp38-win_amd64.whl

h2lib/tests/test_files/DTU_10_MW/log/dtu_10mw_rwt_ver09.log

@@ -1,210 +0,0 @@
- Logfile: ./log/dtu_10mw_rwt_ver09.log is open for log outputs
- ***********************************************************************
- *  Build information for HAWC2MB
- *  Aeroelastic tool HAWC2MB
- *  Intel, version         2021 ,     20201112
- *  WINDOWS 64-bit
- *  DEBUG version
- ***********************************************************************
- *  GIT-TAG        = 13.0.3-1-g750149b-dirty
- *  GIT-BRANCH     = master
- *  BUILD_TYPE     = 
- *  BUILDER        = mmpe
- *  COMPUTER_NAME  = DTU-5CG8462XTX
- *  BUILD_DATE     = 03-08-2023
- ***********************************************************************
-________________________________________________________________________
-  Log file output
-                                                       Time : 13:20:28
-                                                       Date : 21:08.2023
-________________________________________________________________________
- Newmark commands read with succes
- Simulation commands read with succes
- Reading data of main body : tower
- Succes opening ./data/DTU_10MW_RWT_Tower_st.dat
- timoschenko input commands read with succes
- topologi_c2def_inputs read with succes
- Topologi main body tower commands read with succes
- Reading data of main body : towertop
- Succes opening ./data/DTU_10MW_RWT_Towertop_st.dat
- timoschenko input commands read with succes
- topologi_c2def_inputs read with succes
- Topologi main body towertop commands read with succes
- Reading data of main body : shaft
- Succes opening ./data/DTU_10MW_RWT_Shaft_st.dat
- timoschenko input commands read with succes
- topologi_c2def_inputs read with succes
- Topologi main body shaft commands read with succes
- Reading data of main body : hub1
- Succes opening ./data/DTU_10MW_RWT_Hub_st.dat
- timoschenko input commands read with succes
- topologi_c2def_inputs read with succes
- Topologi main body hub1 commands read with succes
- Reading data of main body : hub2
- Topologi main body hub2 commands read with succes
- Reading data of main body : hub3
- Topologi main body hub3 commands read with succes
- Reading data of main body : blade1
- Succes opening ./data/DTU_10MW_RWT_Blade_st.dat
- timoschenko input commands read with succes
- topologi_c2def_inputs read with succes
- Topologi main body blade1 commands read with succes
- Reading data of main body : blade2
- Topologi main body blade2 commands read with succes
- Reading data of main body : blade3
- Topologi main body blade3 commands read with succes
- Base orientation input commands read with succes
- relative orientation input commands read with succes
- relative orientation input commands read with succes
- relative orientation input commands read with succes
- relative orientation input commands read with succes
- relative orientation input commands read with succes
- relative orientation input commands read with succes
- relative orientation input commands read with succes
- relative orientation input commands read with succes
- Orientation input commands read with succes
- Fix0 constraint input commands read with succes
- Fix1 constraint input commands read with succes
- Bearing1 constraint input commands read with succes
- Fix1 constraint input commands read with succes
- Fix1 constraint input commands read with succes
- Fix1 constraint input commands read with succes
- bearing2 constraint input commands read with succes
- bearing2 constraint input commands read with succes
- bearing2 constraint input commands read with succes
- constraint input commands read with succes
- Topologi commands read with succes
- Tower shadow (potential2 flow) commands read with succes
- Wind commands read with succes
- aerodrag element commands read with succes
- aerodrag element commands read with succes
- Aerodrag commands read with succes
- Aerodynamic commands read with succes
- output commands read with succes
- Output commands read
- Dll type2 input commands read with succes
- output commands read with succes
- Output commands read
- Actions commands read
- Dll type2 input commands read with succes
- output commands read with succes
- Output commands read
- Actions commands read
- Dll type2 input commands read with succes
- output commands read with succes
- Output commands read
- Actions commands read
- Dll type2 input commands read with succes
- output commands read with succes
- Output commands read
- Dll type2 input commands read with succes
- DLL commands read with succes
- output commands read with succes
- Output commands read
- Initialization of structure
- Initializing of aero rotor...
- Initialization of rotor aerodynamics
- Succes opening ./data/DTU_10MW_RWT_ae.dat
- Succes opening ./data/DTU_10MW_RWT_pc.dat
- Initialization of rotor induction
- Initialization of wind
- Initialization of external type2 DLL
- External DLL ./control/dtu_we_controller.dll is attempted to open
-DLL loaded with success ./control/dtu_we_controller_64.dll
-Using ./control/dtu_we_controller_64.dll, version: 0.5-2-g885d71c
- DLL subroutine init init_regulation_advanced is called
- In initialization call of ./control/dtu_we_controller.dll Output is 
-  0.000000000000000E+000
-*** WARNING *** Symbol not found in dll: message
- *** INFO *** The DLL subroutine message could not be loaded - bypassed!
- Initialization of external type2 DLL
- External DLL ./control/generator_servo.dll is attempted to open
-DLL loaded with success ./control/generator_servo_64.dll
-Using ./control/generator_servo_64.dll, version: 0.2
- DLL subroutine init init_generator_servo is called
- In initialization call of ./control/generator_servo.dll Output is 
-  0.000000000000000E+000
-*** WARNING *** Symbol not found in dll: message
- *** INFO *** The DLL subroutine message could not be loaded - bypassed!
- Initialization of external type2 DLL
- External DLL ./control/mech_brake.dll is attempted to open
-DLL loaded with success ./control/mech_brake_64.dll
-Using ./control/mech_brake_64.dll, version: 0.2
- DLL subroutine init init_mech_brake is called
- In initialization call of ./control/mech_brake.dll Output is 
-  0.000000000000000E+000
-*** WARNING *** Symbol not found in dll: message
- *** INFO *** The DLL subroutine message could not be loaded - bypassed!
- Initialization of external type2 DLL
- External DLL ./control/servo_with_limits.dll is attempted to open
-DLL loaded with success ./control/servo_with_limits_64.dll
-Using ./control/servo_with_limits_64.dll, version: 1.2
- DLL subroutine init init_servo_with_limits is called
- In initialization call of ./control/servo_with_limits.dll Output is 
-  0.000000000000000E+000
-*** WARNING *** Symbol not found in dll: message
- *** INFO *** The DLL subroutine message could not be loaded - bypassed!
- Initialization of external type2 DLL
- External DLL ./control/towerclearance_mblade.dll is attempted to open
-DLL loaded with success ./control/towerclearance_mblade_64.dll
-Using ./control/towerclearance_mblade_64.dll, version: 1.0.1-1-g8e11be3
- DLL subroutine init initialize is called
- In initialization call of ./control/towerclearance_mblade.dll Output is 
-  0.000000000000000E+000
-*** WARNING *** Symbol not found in dll: message
- *** INFO *** The DLL subroutine message could not be loaded - bypassed!
- Creating link between structure and aerodynamics
- Creating link between structure and aerodrag
- Initialization of Aerodrag
- Starting simulation
-Constraint residual is .6752E-13 < .10E-06. Max residual in constraint 50
-Dynamic stall method:  2 used for entire rotor
-    ... initialization of unified dynamic stall model.
-DS model deactivated: blade  1 at 0mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  1 at 88mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  1 at 355mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  1 at 797mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  1 at 1414mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  1 at 2203mm. Max f_sep=  4.24, dclda=  0.00, thickness= 99.77. See manual to change limit values
-DS model deactivated: blade  1 at 3160mm. Max f_sep=  4.24, dclda=  0.01, thickness= 98.69. See manual to change limit values
-DS model deactivated: blade  1 at 4282mm. Max f_sep=  4.24, dclda=  0.01, thickness= 97.41. See manual to change limit values
-DS model deactivated: blade  1 at 5565mm. Max f_sep=  4.24, dclda=  0.02, thickness= 94.83. See manual to change limit values
-DS model deactivated: blade  1 at 7002mm. Max f_sep=  4.24, dclda=  0.04, thickness= 90.40. See manual to change limit values
-DS model deactivated: blade  1 at 8588mm. Max f_sep=  4.24, dclda=  0.07, thickness= 84.54. See manual to change limit values
-DS model deactivated: blade  1 at 10317mm. Max f_sep=  4.24, dclda=  0.10, thickness= 77.35. See manual to change limit values
-DS model deactivated: blade  1 at 12181mm. Max f_sep=  4.24, dclda=  0.14, thickness= 69.22. See manual to change limit values
-DS model deactivated: blade  1 at 14172mm. Max f_sep=  4.24, dclda=  0.18, thickness= 60.80. See manual to change limit values
-DS model deactivated: blade  1 at 16283mm. Max f_sep=  1.12, dclda=  2.62, thickness= 53.05. See manual to change limit values
-DS model deactivated: blade  2 at 0mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  2 at 88mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  2 at 355mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  2 at 797mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  2 at 1414mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  2 at 2203mm. Max f_sep=  4.24, dclda=  0.00, thickness= 99.77. See manual to change limit values
-DS model deactivated: blade  2 at 3160mm. Max f_sep=  4.24, dclda=  0.01, thickness= 98.69. See manual to change limit values
-DS model deactivated: blade  2 at 4282mm. Max f_sep=  4.24, dclda=  0.01, thickness= 97.41. See manual to change limit values
-DS model deactivated: blade  2 at 5565mm. Max f_sep=  4.24, dclda=  0.02, thickness= 94.83. See manual to change limit values
-DS model deactivated: blade  2 at 7002mm. Max f_sep=  4.24, dclda=  0.04, thickness= 90.40. See manual to change limit values
-DS model deactivated: blade  2 at 8588mm. Max f_sep=  4.24, dclda=  0.07, thickness= 84.54. See manual to change limit values
-DS model deactivated: blade  2 at 10317mm. Max f_sep=  4.24, dclda=  0.10, thickness= 77.35. See manual to change limit values
-DS model deactivated: blade  2 at 12181mm. Max f_sep=  4.24, dclda=  0.14, thickness= 69.22. See manual to change limit values
-DS model deactivated: blade  2 at 14172mm. Max f_sep=  4.24, dclda=  0.18, thickness= 60.80. See manual to change limit values
-DS model deactivated: blade  2 at 16283mm. Max f_sep=  1.12, dclda=  2.62, thickness= 53.05. See manual to change limit values
-DS model deactivated: blade  3 at 0mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  3 at 88mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  3 at 355mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  3 at 797mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  3 at 1414mm. Max f_sep=  0.00, dclda=  0.00, thickness=100.00. See manual to change limit values
-DS model deactivated: blade  3 at 2203mm. Max f_sep=  4.24, dclda=  0.00, thickness= 99.77. See manual to change limit values
-DS model deactivated: blade  3 at 3160mm. Max f_sep=  4.24, dclda=  0.01, thickness= 98.69. See manual to change limit values
-DS model deactivated: blade  3 at 4282mm. Max f_sep=  4.24, dclda=  0.01, thickness= 97.41. See manual to change limit values
-DS model deactivated: blade  3 at 5565mm. Max f_sep=  4.24, dclda=  0.02, thickness= 94.83. See manual to change limit values
-DS model deactivated: blade  3 at 7002mm. Max f_sep=  4.24, dclda=  0.04, thickness= 90.40. See manual to change limit values
-DS model deactivated: blade  3 at 8588mm. Max f_sep=  4.24, dclda=  0.07, thickness= 84.54. See manual to change limit values
-DS model deactivated: blade  3 at 10317mm. Max f_sep=  4.24, dclda=  0.10, thickness= 77.35. See manual to change limit values
-DS model deactivated: blade  3 at 12181mm. Max f_sep=  4.24, dclda=  0.14, thickness= 69.22. See manual to change limit values
-DS model deactivated: blade  3 at 14172mm. Max f_sep=  4.24, dclda=  0.18, thickness= 60.80. See manual to change limit values
-DS model deactivated: blade  3 at 16283mm. Max f_sep=  1.12, dclda=  2.62, thickness= 53.05. See manual to change limit values
- Global time =   1.000000000000000E-002  Iter =            1
- Creating file folder res
- Creating file folder res

h2lib/tests/test_files/__init__.py

@@ -1,2 +0,0 @@
-import os
-tfp = os.path.dirname (__file__) + "/"

h2lib/tests/test_files/minimal/htc/minimal.htc

@@ -1,32 +0,0 @@
-begin simulation;
-  time_stop    1000.0;
-  solvertype   1 ;    (newmark)
-  on_no_convergence continue ;
-  convergence_limits 1E3 1.0 1E-7 ;
-  log_deltat 1;
-  begin newmark;
-    deltat    0.1;  
-  end newmark;
-end simulation;
-;
-;---------------------------------------------------------------------------------------------------------------------------------------------------------------- 
-begin wind ;
-  density                 1.225 ;
-  wsp                     5.0   ;
-  tint                    0.1   ;
-  horizontal_input        0     ;
-  windfield_rotations     0.0 0 0.0 ;    yaw, tilt, rotation
-  center_pos0             0 2 0 ; hub heigth
-  shear_format            0 0.2 ; 
-  turb_format             0     ;  0=none, 1=mann,2=flex
-  tower_shadow_method     0     ;  0=none, 1=potential flow, 2=jet
-  scale_time_start       0.0 ; 
-end wind;
-;
-begin output;
-  filename res/minimal;
-  data_format gtsdf64;
-  buffer 1;
-  general time;
-end output;
-exit;

h2lib/tests/test_files/minimal/htc/minimal_mann_turb.htc

@@ -1,41 +0,0 @@
-begin simulation;
-  time_stop    140.0;
-  solvertype   1 ;    (newmark)
-  on_no_convergence continue ;
-  convergence_limits 1E3 1.0 1E-7 ;
-  begin newmark;
-    deltat    0.1;  
-  end newmark;
-end simulation;
-;
-;---------------------------------------------------------------------------------------------------------------------------------------------------------------- 
-begin wind ;
-  density                 1.225 ;
-  wsp                     5.0   ;
-  tint                    0.1   ;
-  horizontal_input        0     ;
-  windfield_rotations     0.0 0.0 0.0 ;    yaw, tilt, rotation
-  center_pos0             46.5 2 -30 ; hub heigth
-  shear_format            0 0.2 ; 
-  turb_format             1     ;  0=none, 1=mann,2=flex
-  tower_shadow_method     0     ;  0=none, 1=potential flow, 2=jet
-  scale_time_start       0.0 ; 
-  begin mann;
-      create_turb_parameters 33.6 .1 3.9 1 0;
-      filename_u turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_u;
-      filename_v turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_v;
-      filename_w turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_w;
-      box_dim_u  512 2;
-      box_dim_v  32 3;
-      box_dim_w  16 4;
-      reverse_mann_box 0;
-      dont_scale 1;
-  end mann;
-end wind;
-;
-begin output;
-  filename res/minimal_mann_turb;
-  data_format gtsdf64;
-  wind free_wind 1 0 0 0;
-end output;
-exit;

h2lib/tests/test_h2lib.py

@@ -1,217 +0,0 @@
-from wetb.gtsdf import gtsdf
-import matplotlib.pyplot as plt
-import numpy as np
-from numpy import newaxis as na
-import pytest
-from h2lib._h2lib import H2Lib, MultiH2Lib
-from numpy import testing as npt
-from h2lib.tests.test_files import tfp
-from wetb.hawc2.htc_file import HTCFile
-import time
-from tqdm import tqdm
-
-
-def test_echo_version():
-    with H2Lib() as h2:
-        h2.suppress_output = False
-        h2.echo_version()
-        h2.suppress_output = True
-
-
-def test_sensors():
-    with H2Lib() as h2:
-        model_path = tfp + "minimal/"
-        htc = HTCFile(model_path + 'htc/minimal.htc')
-        htc.output.add_sensor("general", 'variable', [1, 5, 'mysensor', 'myunit', 'mydesc'])
-        htc.set_name('tmp')
-        htc.save()
-        h2.read_input('htc/tmp.htc', model_path)
-        i1 = h2.add_sensor('wind free_wind 1 0 0 -119;')  # adds three sensors
-        assert i1 == 3
-        i2 = h2.add_sensor('general step 0.5 1 2;')
-        assert i2 == 4
-        i3 = h2.add_sensor('general variable 1 5 myname myunit mydesc;')
-        assert i3 == 5
-        h2.init()
-        t = h2.step()
-        assert h2.get_sensor_info(1) == ['WSP gl. coo.,Vx', 'm/s',
-                                         'Free wind speed Vx, gl. coo, of gl. pos    0.00,   0.00,-119.00']
-        assert h2.get_sensor_info(2) == ['WSP gl. coo.,Vy', 'm/s',
-                                         'Free wind speed Vy, gl. coo, of gl. pos    0.00,   0.00,-119.00']
-        assert h2.get_sensor_info(3) == ['WSP gl. coo.,Vz', 'm/s',
-                                         'Free wind speed Vz, gl. coo, of gl. pos    0.00,   0.00,-119.00']
-        assert h2.get_sensor_info(4) == ['Step', '-', 'Step value']
-        assert h2.get_sensor_info(5) == ['myname', 'myunit', 'mydesc']
-
-        r = []
-        for _ in range(10):
-            r.append([t, h2.get_sensor_value(4), h2.get_sensor_value(5)])
-            if t == 0.6:
-                h2.set_variable_sensor_value(1, 2)
-            t = h2.step()
-        r = np.array(r)
-        npt.assert_array_equal(r[:, 1], np.where(r[:, 0] <= .5, 1, 2))
-        npt.assert_array_equal(r[:, 2], np.where(r[:, 0] >= .7, 2, 5.))
-        time, data, _ = gtsdf.load(h2.model_path + 'res/tmp.hdf5')
-        npt.assert_array_equal(data[:, 0], np.where(time >= .7, 2, 5.))
-
-
-def test_set_windspeed_uvw():
-    with H2Lib() as h2:
-        pos = [0, 0, -119]
-        model_path = tfp + "minimal/"
-        h2.read_input('htc/minimal.htc', model_path)
-        Nxyz = [8, 4, 2]
-        h2.init_windfield(Nxyz, dxyz=[1., 1, 1], box_offset_yz=[0., -119])
-        h2.init()
-        u, v, w = 1, 2, 3
-        uvw = np.zeros([3] + Nxyz) + np.array([u, v, w])[:, na, na, na]
-        h2.set_windfield(uvw, 0)
-        vx, vy, vz = h2.get_wind_speed(pos)
-        assert vx == v
-        assert vy == u
-        assert vz == -w
-        npt.assert_array_equal([u, v, w], h2.get_uvw(pos))
-
-
-def test_set_windspeed_interpolation():
-    with H2Lib() as h2:
-        model_path = tfp + "minimal/"
-        h2.read_input('htc/minimal.htc', model_path)
-        Nxyz = np.array([8, 4, 2])
-        h2.init_windfield(Nxyz, dxyz=[2., 3, 4], box_offset_yz=[-4.5, 117])
-        h2.init()
-        uvw = np.array(np.meshgrid(np.arange(Nxyz[0]) + 100, np.arange(Nxyz[1]) + 200, np.arange(Nxyz[2]) + 300,
-                                   indexing='ij'),
-                       dtype=np.float64)
-        h2.set_windfield(uvw, 0)
-
-        # test u
-        y = np.linspace(0, 14, 10)
-        np.testing.assert_array_almost_equal([h2.get_wind_speed([0, y_, 0])[1] for y_ in y], np.linspace(100, 107, 10))
-
-        h2.set_windfield(uvw, -5.)
-        # box moved -6m offset
-        np.testing.assert_array_almost_equal([h2.get_wind_speed([0, y_, 0])[1]
-                                             for y_ in y - 5], np.linspace(100, 107, 10))
-
-        h2.run(3)
-        # box moved 3s * 5m/s - 5m offset = 10
-        np.testing.assert_array_almost_equal([h2.get_wind_speed([0, y_, 0])[1]
-                                             for y_ in y + 10], np.linspace(100, 107, 10))
-
-        h2.set_windfield(uvw, -5.)
-        # box moved 6m offset
-        np.testing.assert_array_almost_equal([h2.get_wind_speed([0, y_, 0])[1]
-                                             for y_ in y - 5], np.linspace(100, 107, 10))
-
-        # test v
-        x = np.linspace(-4.5, 4.5, 10)
-        np.testing.assert_array_almost_equal([h2.get_wind_speed([x_, 0, 0])[0] for x_ in x], np.linspace(200, 203, 10))
-
-        # test w
-        z = np.linspace(-117, -121, 10)
-        np.testing.assert_array_almost_equal([h2.get_wind_speed([0, 0, z_])[2] for z_ in z], -np.linspace(300, 301, 10))
-
-
-def test_mann_turb():
-    with H2Lib() as h2:
-        U = 5
-        from hipersim import MannTurbulenceField
-        mtf = MannTurbulenceField.from_hawc2(filenames=[tfp + f'minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_{uvw}' for uvw in 'uvw'],
-                                             alphaepsilon=.1, L=33.6, Gamma=3.9,
-                                             Nxyz=(512, 32, 16), dxyz=(2, 3, 4), seed=1, HighFreqComp=0)
-
-        time, data, _ = gtsdf.load(tfp + 'minimal/res/minimal_mann_turb.hdf5')
-        time, data = time[1000:], data[1000:]
-        vx, vy, vz = data.T
-        h2_uvw = [vy, vx, -vz]
-        model_path = tfp + "minimal/"
-        h2.read_input('htc/minimal.htc', model_path)
-        h2.init_windfield(Nxyz=mtf.Nxyz, dxyz=mtf.dxyz, box_offset_yz=[0, 0])
-        h2.init()
-        h2.set_windfield(uvw=mtf.uvw.astype(np.float64), box_offset_x=-(mtf.Nx - 1) * mtf.dx)
-        h2l_vxyz = []
-        for t in tqdm(time, disable=0):
-            h2.run(t)
-            h2l_vxyz.append(h2.get_wind_speed([0, 0, 0]))
-        vx, vy, vz = np.array(h2l_vxyz).T
-        h2l_uvw = [vy, vx, -vz]
-
-        x = (mtf.Nx * mtf.dx) - U * time - 2
-        mtf_uvw = mtf(x, x * 0, x * 0).T
-
-        if 0:
-            plt.plot(time, h2_uvw[0], label='h2')
-            plt.plot(time, mtf_uvw[0], label='mtf')
-            plt.plot(time, h2l_uvw[0], label='h2l')
-            plt.legend()
-            plt.show()
-
-        npt.assert_array_almost_equal(mtf_uvw, h2_uvw)
-        npt.assert_array_almost_equal(mtf_uvw, h2l_uvw)
-
-
-def test_mann_turb_small_buffer():
-    with H2Lib() as h2:
-        from hipersim import MannTurbulenceField
-        U = 5
-        mtf = MannTurbulenceField.from_hawc2(filenames=[tfp + f'minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_{uvw}' for uvw in 'uvw'],
-                                             alphaepsilon=.1, L=33.6, Gamma=3.9,
-                                             Nxyz=(512, 32, 16), dxyz=(2, 3, 4), seed=1, HighFreqComp=0)
-
-        time, data, _ = gtsdf.load(tfp + 'minimal/res/minimal_mann_turb.hdf5')
-        time, data = time[:1400], data[:1400]
-        vx, vy, vz = data.T
-        h2_uvw = [vy, vx, -vz]
-        model_path = tfp + "minimal/"
-        h2.read_input('htc/minimal.htc', model_path)
-        h2.init_windfield(Nxyz=(128, 32, 16), dxyz=mtf.dxyz, box_offset_yz=[0, 0])
-        h2.init()
-        h2l_vxyz = []
-        h2.set_windfield(uvw=mtf.uvw[:, -128:].astype(np.float64), box_offset_x=-(mtf.Nx - 1) * mtf.dx)
-        for t in tqdm(time):
-            t = np.round(t, 6)
-            h2.run(t)
-            if t % 40 == 0:
-                i = int(512 - 128 - t * U / 2)
-                h2.set_windfield(uvw=mtf.uvw[:, i:i + 128].astype(np.float64), box_offset_x=-(mtf.Nx - 1) * mtf.dx)
-            h2l_vxyz.append(h2.get_wind_speed([0, 0, 0]))
-        vx, vy, vz = np.array(h2l_vxyz).T
-        h2l_uvw = [vy, vx, -vz]
-
-        x = ((mtf.Nx - 1) * mtf.dx) - U * time
-        mtf_uvw = mtf(x, x * 0, x * 0).T
-
-        if 0:
-            plt.plot(time, h2_uvw[0], '-', label='h2')
-            plt.plot(time, mtf_uvw[0], '--', label='mtf')
-            plt.plot(time, h2l_uvw[0], '-', label='h2l')
-            plt.legend()
-            plt.show()
-
-        npt.assert_array_almost_equal(mtf_uvw, h2_uvw)
-        npt.assert_array_almost_equal(mtf_uvw, h2l_uvw)
-
-
-def test_parallel():
-    N = 3
-    mh2 = MultiH2Lib(N)
-    assert mh2.getState() == [0, 0, 0]
-    mh2.setState([4, 5, 6])
-    assert mh2.getState() == [4, 5, 6]
-    mh2.close()
-
-
-def test_parallel_context_manager():
-    N = 3
-    with MultiH2Lib(N) as mh2:
-        mh2.init_windfield(Nxyz=(128, 32, 16), dxyz=[(1, 2, 3), (4, 5, 6), (7, 8, 9)], box_offset_yz=[0, 0])
-        assert mh2.getState() == [0, 0, 0]
-        mh2.setState([4, 5, 6])
-        assert mh2.getState() == [4, 5, 6]
-        mh2[0].setState(10)
-        assert mh2.getState() == [10, 5, 6]
-        t = time.time()
-        mh2.work(2)
-        assert time.time() - t < 2.5

h2lib/tests/test_h2rotor.py

@@ -1,162 +0,0 @@
-import pytest
-
-from h2lib._h2lib import H2Lib, MultiH2Lib
-from h2lib.tests.test_files import tfp
-import numpy as np
-import matplotlib.pyplot as plt
-from h2lib.tests.dtu10mw import DTU10MW
-from wetb.gtsdf import gtsdf
-from numpy import testing as npt
-from h2lib.utils import ProcessClass
-
-
-def get_h2(htc_path='htc/DTU_10MW_RWT.htc'):
-    h2 = ProcessClass(H2Lib)()
-    h2.read_input(htc_path=htc_path, model_path=tfp + 'DTU_10_MW')
-
-    return h2
-
-
-@pytest.fixture(scope='module')
-def h2():
-    h2 = get_h2()
-    h2.init()
-    h2.step()
-    yield h2
-    h2.close()
-
-
-def test_get_nxxx(h2):
-    assert h2.get_rotor_dims() == [[50, 50, 50]]
-
-
-def test_get_bem_grid_dim(h2):
-    assert h2.get_bem_grid_dim() == [16, 50]
-
-
-def test_get_bem_grid(h2):
-    azi, rad = h2.get_bem_grid()
-
-    npt.assert_array_almost_equal(np.roll(np.linspace(-np.pi, np.pi, 17)[1:], 9), azi)
-    npt.assert_array_almost_equal([3.09929662, 11.4350243, 33.76886193, 60.8641887, 82.01217305], rad[::10])
-
-
-def test_induction(h2):
-    azi, rad = h2.get_bem_grid()
-    induc_grid = h2.get_induction_polargrid()
-    induc_axisymmetric = h2.get_induction_axisymmetric()
-    if 0:
-        Azi, Rad = np.meshgrid(azi, rad)
-        ax = plt.gcf().add_subplot(111, polar=True)
-        ax.set_theta_zero_location('S')
-        ax.set_theta_direction(-1)
-        cntf = ax.contourf(Azi, Rad, induc_grid.T, 50)
-        plt.colorbar(cntf)
-        plt.figure()
-        plt.plot(rad, induc_axisymmetric)
-        plt.show()
-
-    npt.assert_array_almost_equal(np.mean(induc_grid, 0), induc_axisymmetric)
-
-
-def test_rotor_orientation_multi_instance():
-    dtu10 = DTU10MW()
-
-    dtu10.set_name('tmp_5_0')
-    dtu10.save()
-
-    tilt_ref, yaw_ref = 6, 10
-    dtu10.set_tilt_cone_yaw(tilt=tilt_ref, cone=0, yaw=yaw_ref)
-    dtu10.set_name('tmp_6_10')
-    dtu10.save()
-    with MultiH2Lib(2) as mh2:
-        mh2.read_input(['htc/tmp_5_0.htc', 'htc/tmp_6_10.htc'], model_path=tfp + "DTU_10_MW")
-        #h2.suppress_output = False
-        s_id = mh2.add_sensor('aero power')[0]
-        mh2.init()
-        yaw, tilt, _ = zip(*mh2.get_rotor_orientation(deg=True))
-        np.testing.assert_almost_equal(yaw, [0, 10])
-        np.testing.assert_almost_equal(tilt, [5, 6])
-        h2 = mh2[0]
-        res = []
-        for t in np.arange(0, 2.5, .01):
-            h2.run(t)
-            res.append([h2.get_rotor_orientation(deg=True)[2], h2.get_sensor_value(s_id)] +
-                       h2.get_rotor_position().tolist())
-
-        data = gtsdf.load(h2.model_path + '/res/tmp_5_0.hdf5')[1]
-        res = np.array(res)
-        npt.assert_allclose(data[:, 0], res[1:-1, 0] % 360 - 180, rtol=0.002)  # azi
-        npt.assert_array_almost_equal(data[:, 10], res[1:-1, 1])  # power
-        npt.assert_array_almost_equal(data[:, 15:18], res[1:-1, 2:5])  # rotor position
-
-
-def test_rotor_avg_windspeed():
-    h2 = get_h2()
-    h2.init_windfield(Nxyz=(2, 2, 2), dxyz=(200, 200, 200), box_offset_yz=(-100, -19))
-    h2.init()
-    u = np.zeros((2, 2, 2))
-    u[:, 1, :] = 10  # 0 in one side, 10 in other, 5 in avg
-    h2.set_windfield([u, u * 0, u * 0], -100)
-    h2.step()
-
-    npt.assert_almost_equal(h2.get_rotor_avg_wsp(1), 5)
-
-
-def test_aerosections():
-    plot = False
-    h2 = get_h2()
-    h2.suppress_output = False
-    # blade 1, global coo, r>30
-    pos_ids = [h2.add_sensor(f'aero position 3 1 {xyz} 30') for xyz in [1, 2, 3]]
-    wsp_ids = [h2.add_sensor(f'aero windspeed 3 1 {xyz} 30') for xyz in [1, 2, 3]]
-    frc_ids = [h2.add_sensor(f'aero secforce 1 {xyz} 30 3') for xyz in [1, 2, 3]]
-    mom_ids = [h2.add_sensor(f'aero secmoment 1 {xyz} 30 3') for xyz in [1, 2, 3]]
-    h2.init()
-
-    a = h2.get_aerosections_position()
-    if plot:
-        ax = plt.figure().add_subplot(projection='3d')
-        for b in a:
-            ax.plot(*b.T)
-        plt.show()
-
-    assert a.shape == (3, 50, 3)
-
-    r = np.sqrt(np.sum((a[0, :] - a[0, 0])**2, 1))
-    i = np.searchsorted(r, 30)
-
-    h2.step()
-    name, unit, desc = h2.get_sensor_info(pos_ids)
-    assert str(np.round(r[i], 2)) in desc
-    assert str(np.round(r[i], 1)) in name
-    assert unit == 'm'
-
-    a = h2.get_aerosections_position()
-    npt.assert_array_almost_equal(a[0, i], [h2.get_sensor_value(id) for id in pos_ids])
-    uvw = a * 0
-    uvw[:, :, 0] = 6
-    uvw[0, i, 0] = 12
-    h2.run(3)
-    npt.assert_array_equal([h2.get_sensor_value(id) for id in wsp_ids], [0, 6, 0])
-    npt.assert_array_almost_equal([h2.get_sensor_value(id) for id in frc_ids],
-                                  h2.get_aerosections_load()[0][0, i] / 1000)
-
-    frc_before = h2.get_aerosections_load()[0]
-    if plot:
-        plt.figure()
-        plt.plot(frc_before[:, :, 1].T)
-
-    h2.set_aerosections_windspeed(uvw)
-    h2.step()
-    frc_after = h2.get_aerosections_load()[0]
-    if plot:
-        plt.plot(frc_after[:, :, 1].T, '--')
-        plt.show()
-
-    # Fy in section with u=12 instead of u=6m/s more than doubled
-    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, :]
-    npt.assert_allclose(frc_before, frc_after, atol=7)