pyvale 2025.4.1__py3-none-any.whl → 2025.5.2__py3-none-any.whl

pyvale/examples/ex3_3_displacement2d.py

@@ -1,101 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-from pathlib import Path
-import numpy as np
-import matplotlib.pyplot as plt
-from scipy.spatial.transform import Rotation
-import mooseherder as mh
-import pyvale as pyv
-
-def main() -> None:
-    """pyvale example: displacement sensors on a 2D plate with a hole
-    ----------------------------------------------------------------------------
-    - Demonstrates sensor rotation for vector fields
-    """
-    data_path = pyv.DataSet.mechanical_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    descriptor = pyv.SensorDescriptorFactory.displacement_descriptor()
-
-    spat_dims = 2
-    field_key = 'disp'
-    components = ('disp_x','disp_y')
-    disp_field = pyv.FieldVector(sim_data,field_key,components,spat_dims)
-
-    n_sens = (2,2,1)
-    x_lims = (0.0,100.0)
-    y_lims = (0.0,150.0)
-    z_lims = (0.0,0.0)
-    sensor_positions = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
-
-    use_sim_time = False
-    if use_sim_time:
-        sample_times = None
-    else:
-        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
-
-    sens_data_norot = pyv.SensorData(positions=sensor_positions,
-                                        sample_times=sample_times)
-
-    disp_sens_norot = pyv.SensorArrayPoint(sens_data_norot,
-                                              disp_field,
-                                              descriptor)
-
-    meas_norot = disp_sens_norot.get_measurements()
-
-    sens_angles = sensor_positions.shape[0] * \
-        (Rotation.from_euler("zyx", [45, 0, 0], degrees=True),)
-
-    sens_data_rot = pyv.SensorData(positions=sensor_positions,
-                                      sample_times=sample_times,
-                                      angles=sens_angles)
-
-    disp_sens_rot = pyv.SensorArrayPoint(sens_data_rot,
-                                            disp_field,
-                                            descriptor)
-
-
-    angle_offset = np.zeros_like(sensor_positions)
-    angle_offset[:,0] = 1.0 # only rotate about z in 2D
-    angle_error_data = pyv.ErrFieldData(ang_offset_zyx=angle_offset)
-
-    sys_err_rot = pyv.ErrSysField(disp_field,angle_error_data)
-
-    sys_err_int = pyv.ErrIntegrator([sys_err_rot],
-                                         sens_data_rot,
-                                         disp_sens_rot.get_measurement_shape())
-    disp_sens_rot.set_error_integrator(sys_err_int)
-
-    meas_rot = disp_sens_rot.get_measurements()
-
-
-    print(80*'-')
-    sens_num = 4
-    print('The last 5 time steps (measurements) of sensor {sens_num}:')
-    pyv.print_measurements(disp_sens_rot,
-                              (sens_num-1,sens_num),
-                              (0,1),
-                              (meas_rot.shape[2]-5,meas_rot.shape[2]))
-    print(80*'-')
-
-    plot_field = 'disp_x'
-    if plot_field == 'disp_x':
-        pv_plot = pyv.plot_point_sensors_on_sim(disp_sens_rot,'disp_x')
-        pv_plot.show(cpos="xy")
-    elif plot_field == 'disp_y':
-        pv_plot = pyv.plot_point_sensors_on_sim(disp_sens_rot,'disp_y')
-        pv_plot.show(cpos="xy")
-
-    pyv.plot_time_traces(disp_sens_norot,plot_field)
-    pyv.plot_time_traces(disp_sens_rot,plot_field)
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex3_4_displacement2d.py

@@ -1,102 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-import numpy as np
-import matplotlib.pyplot as plt
-from scipy.spatial.transform import Rotation
-import mooseherder as mh
-import pyvale as pyv
-
-def main() -> None:
-    """pyvale example: displacement sensors on a 2D plate with a hole
-    ----------------------------------------------------------------------------
-    """
-    data_path = pyv.DataSet.mechanical_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    descriptor = pyv.SensorDescriptorFactory.displacement_descriptor()
-
-    spat_dims = 2
-    field_key = "disp"
-    components = ("disp_x","disp_y")
-    disp_field = pyv.FieldVector(sim_data,field_key,components,spat_dims)
-
-    #---------------------------------------------------------------------------
-    n_sens = (2,2,1)
-    x_lims = (0.0,100.0)
-    y_lims = (0.0,150.0)
-    z_lims = (0.0,0.0)
-    sensor_positions = pyv.create_sensor_pos_array(n_sens,
-                                                      x_lims,
-                                                      y_lims,
-                                                      z_lims)
-
-    use_sim_time = False
-    if use_sim_time:
-        sample_times = None
-    else:
-        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
-
-    sensor_angles = sensor_positions.shape[0] * \
-        (Rotation.from_euler("zyx", [0, 0, 0], degrees=True),)
-
-    sensor_data = pyv.SensorData(positions=sensor_positions,
-                                  sample_times=sample_times,
-                                  angles=sensor_angles,
-                                  spatial_averager=pyv.EIntSpatialType.QUAD4PT,
-                                  spatial_dims=np.array([5.0,5.0,0.0]))
-
-    #---------------------------------------------------------------------------
-    disp_sensors = pyv.SensorArrayPoint(sensor_data,
-                                           disp_field,
-                                           descriptor)
-
-    pos_offset = -10.0*np.ones_like(sensor_positions)
-    pos_offset[:,2] = 0 # in 2d we only have offset in x and y so zero z
-
-    angle_offset = np.zeros_like(sensor_positions)
-    angle_offset[:,0] = 5.0 # only rotate about z in 2D
-
-    time_offset = 1.0*np.ones_like(disp_sensors.get_sample_times())
-
-    field_error_data = pyv.ErrFieldData(pos_offset_xyz=pos_offset,
-                                           ang_offset_zyx=angle_offset,
-                                           time_offset=time_offset)
-
-    error_chain = []
-    error_chain.append(pyv.ErrSysField(disp_field,field_error_data))
-    error_integrator = pyv.ErrIntegrator(error_chain,
-                                            sensor_data,
-                                            disp_sensors.get_measurement_shape())
-
-    disp_sensors.set_error_integrator(error_integrator)
-
-    measurements = disp_sensors.calc_measurements()
-
-    #---------------------------------------------------------------------------
-    print(80*"-")
-    sens_num = 4
-    print("The last 5 time steps (measurements) of sensor {sens_num}:")
-    pyv.print_measurements(disp_sensors,
-                              (sens_num-1,sens_num),
-                              (0,1),
-                              (measurements.shape[2]-5,measurements.shape[2]))
-    print(80*"-")
-
-    #---------------------------------------------------------------------------
-    plot_field = "disp_x"
-
-    pv_plot = pyv.plot_point_sensors_on_sim(disp_sensors,plot_field)
-    pv_plot.show(cpos="xy")
-
-    pyv.plot_time_traces(disp_sensors,plot_field)
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex4_1_strain2d.py

@@ -1,54 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-from pathlib import Path
-import numpy as np
-import matplotlib.pyplot as plt
-import mooseherder as mh
-import pyvale as pyv
-
-def main() -> None:
-    """pyvale example: strain sensors on a 2D plate with a hole
-    ----------------------------------------------------------------------------
-    """
-    data_path = pyv.DataSet.mechanical_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    n_sens = (2,3,1)
-    x_lims = (0.0,100.0)
-    y_lims = (0.0,150.0)
-    z_lims = (0.0,0.0)
-    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
-
-    use_sim_time = False
-    if use_sim_time:
-        sample_times = None
-    else:
-        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
-
-    sens_data = pyv.SensorData(positions=sens_pos,
-                                  sample_times=sample_times)
-
-    straingauge_array = pyv.SensorArrayFactory \
-                            .strain_gauges_basic_errs(sim_data,
-                                                      sens_data,
-                                                      "strain",
-                                                      spat_dims=2)
-
-    plot_field = 'strain_yy'
-    pv_plot = pyv.plot_point_sensors_on_sim(straingauge_array,plot_field)
-    pv_plot.show(cpos="xy")
-
-    pyv.plot_time_traces(straingauge_array,'strain_xx')
-    pyv.plot_time_traces(straingauge_array,'strain_yy')
-    pyv.plot_time_traces(straingauge_array,'strain_xy')
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex4_2_strain2d.py

@@ -1,76 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-import numpy as np
-import matplotlib.pyplot as plt
-import mooseherder as mh
-import pyvale as pyv
-
-def main() -> None:
-    """pyvale example: strain sensors on a 2D plate with a hole
-    ----------------------------------------------------------------------------
-    """
-    data_path = pyv.DataSet.mechanical_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    descriptor = pyv.SensorDescriptor()
-    descriptor.name = 'Strain'
-    descriptor.symbol = r'\varepsilon'
-    descriptor.units = r'-'
-    descriptor.tag = 'SG'
-    descriptor.components = ('xx','yy','xy')
-
-    spat_dims = 2
-    field_key = 'strain'
-    norm_components = ('strain_xx','strain_yy')
-    dev_components = ('strain_xy',)
-    strain_field = pyv.FieldTensor(sim_data,
-                                    field_key,
-                                    norm_components,
-                                    dev_components,
-                                    spat_dims)
-
-    n_sens = (2,3,1)
-    x_lims = (0.0,100.0)
-    y_lims = (0.0,150.0)
-    z_lims = (0.0,0.0)
-    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
-
-    use_sim_time = False
-    if use_sim_time:
-        sample_times = None
-    else:
-        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
-
-    sens_data = pyv.SensorData(positions=sens_pos,
-                                  sample_times=sample_times)
-
-    straingauge_array = pyv.SensorArrayPoint(sens_data,
-                                                strain_field,
-                                                descriptor)
-
-    error_chain = []
-    error_chain.append(pyv.ErrSysUniform(low=-0.1e-3,high=0.1e-3))
-    error_chain.append(pyv.ErrRandNormal(std=0.1e-3))
-    error_int = pyv.ErrIntegrator(error_chain,
-                                       sens_data,
-                                       straingauge_array.get_measurement_shape())
-    straingauge_array.set_error_integrator(error_int)
-
-    plot_field = 'strain_yy'
-    pv_plot = pyv.plot_point_sensors_on_sim(straingauge_array,plot_field)
-    pv_plot.show(cpos="xy")
-
-    pyv.plot_time_traces(straingauge_array,'strain_xx')
-    pyv.plot_time_traces(straingauge_array,'strain_yy')
-    pyv.plot_time_traces(straingauge_array,'strain_xy')
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex4_3_strain2d.py

@@ -1,97 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-import numpy as np
-import matplotlib.pyplot as plt
-from scipy.spatial.transform import Rotation as R
-import mooseherder as mh
-import pyvale as pyv
-
-def main() -> None:
-    """pyvale example: strain sensors on a 2D plate with a hole
-    ----------------------------------------------------------------------------
-    - Demonstrates rotation of tensor fields
-    """
-    data_path = pyv.DataSet.mechanical_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    descriptor = pyv.SensorDescriptorFactory.strain_descriptor()
-
-    spat_dims = 2
-    field_key = 'strain'
-    norm_components = ('strain_xx','strain_yy')
-    dev_components = ('strain_xy',)
-    strain_field = pyv.FieldTensor(sim_data,
-                                    field_key,
-                                    norm_components,
-                                    dev_components,
-                                    spat_dims)
-
-    n_sens = (2,3,1)
-    x_lims = (0.0,100.0)
-    y_lims = (0.0,150.0)
-    z_lims = (0.0,0.0)
-    sensor_positions = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
-
-    use_sim_time = False
-    if use_sim_time:
-        sample_times = None
-    else:
-        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
-
-    sens_data_norot = pyv.SensorData(positions=sensor_positions,
-                                        sample_times=sample_times)
-
-    sg_array_norot = pyv.SensorArrayPoint(sens_data_norot,
-                                             strain_field,
-                                             descriptor)
-
-    meas_norot = sg_array_norot.get_measurements()
-
-    sens_angles = sensor_positions.shape[0] * \
-        (R.from_euler("zyx", [45, 0, 0], degrees=True),)
-
-    sens_data_rot = pyv.SensorData(positions=sensor_positions,
-                                      sample_times=sample_times,
-                                      angles=sens_angles)
-
-    sg_array_rot = pyv.SensorArrayPoint(sens_data_rot,
-                                           strain_field,
-                                           descriptor)
-
-    angle_offset = np.zeros_like(sensor_positions)
-    angle_offset[:,0] = 1.0 # only rotate about z in 2D
-    angle_error_data = pyv.ErrFieldData(ang_offset_zyx=angle_offset)
-
-    sys_err_rot = pyv.ErrSysField(strain_field,angle_error_data)
-    err_int = pyv.ErrIntegrator([sys_err_rot],
-                                     sens_data_rot,
-                                     sg_array_rot.get_measurement_shape())
-    sg_array_rot.set_error_integrator(err_int)
-
-
-    meas_rot = sg_array_rot.get_measurements()
-
-
-    print(80*'-')
-    sens_num = 4
-    print('The last 5 time steps (measurements) of sensor {sens_num}:')
-    pyv.print_measurements(sg_array_rot,
-                              (sens_num-1,sens_num),
-                              (1,2),
-                              (meas_rot.shape[2]-5,meas_rot.shape[2]))
-    print(80*'-')
-
-    plot_comp = 'strain_yy'
-    pyv.plot_time_traces(sg_array_norot,plot_comp)
-    pyv.plot_time_traces(sg_array_rot,plot_comp)
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex5_1_multiphysics2d.py

@@ -1,75 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-import numpy as np
-import matplotlib.pyplot as plt
-import mooseherder as mh
-import pyvale as pyv
-
-
-def main() -> None:
-    """pyvale example: thermo-mechanical multi-physics sensors on a 2D plate
-    ----------------------------------------------------------------------------
-    """
-    #===========================================================================
-    # Load Simulations as mooseherder.SimData objects
-    data_path = pyv.DataSet.thermomechanical_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    #===========================================================================
-    # Specify sensor locations and sample times
-    n_sens = (4,1,1)
-    x_lims = (0.0,100.0)
-    y_lims = (0.0,50.0)
-    z_lims = (0.0,0.0)
-    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
-
-    use_sim_time = True
-    if use_sim_time:
-        sample_times = None
-    else:
-        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
-
-    sens_data = pyv.SensorData(positions=sens_pos,
-                                  sample_times=sample_times)
-
-    #===========================================================================
-    # Create pyvale sensor arrays for thermal and mechanical data
-    tc_field = 'temperature'
-    tc_array = pyv.SensorArrayFactory \
-        .thermocouples_basic_errs(sim_data,
-                                  sens_data,
-                                  tc_field,
-                                  spat_dims=2)
-
-    sg_field = 'strain'
-    sg_array = pyv.SensorArrayFactory \
-        .strain_gauges_basic_errs(sim_data,
-                                  sens_data,
-                                  sg_field,
-                                  spat_dims=2)
-
-    #===========================================================================
-    # Visualise Traces
-    print(80*'-')
-    sens_num = 4
-    print('THERMAL: The last 5 time steps (measurements) of sensor {sens_num}:')
-    pyv.print_measurements(tc_array,
-                              (sens_num-1,sens_num),
-                              (0,1),
-                              (tc_array.get_measurement_shape()[2]-5,
-                               tc_array.get_measurement_shape()[2]))
-    print(80*'-')
-
-    pyv.plot_time_traces(tc_array,"temperature")
-    pyv.plot_time_traces(sg_array,"strain_xx")
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex6_1_multiphysics2d_expsim.py

@@ -1,115 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-import numpy as np
-import matplotlib.pyplot as plt
-import mooseherder as mh
-import pyvale as pyv
-
-def main() -> None:
-    """pyvale example: thermo-mechanical multi-physics sensors on a 2D plate
-    ----------------------------------------------------------------------------
-    - Demonstrates the experiment module for running many monte-carlo style
-      experiments and statistically analysing the results.
-    """
-    # Load Simulations as mooseherder.SimData objects
-    #base_path = Path("src/pyvale/data")
-    data_paths = pyv.DataSet.thermomechanical_2d_experiment_paths()
-
-    sim_list = []
-    for pp in data_paths:
-        sim_data = mh.ExodusReader(pp).read_all_sim_data()
-        # Scale to mm to make 3D visualisation scaling easier
-        sim_data.coords = sim_data.coords*1000.0 # type: ignore
-        sim_list.append(sim_data)
-
-    #===========================================================================
-    # Create pyvale sensor arrays for thermal and mechanical data
-    sim_data = sim_list[0]
-
-    n_sens = (4,1,1)
-    x_lims = (0.0,100.0)
-    y_lims = (0.0,50.0)
-    z_lims = (0.0,0.0)
-    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
-
-    use_sim_time = True
-    if use_sim_time:
-        sample_times = None
-    else:
-        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
-
-    sens_data = pyv.SensorData(positions=sens_pos,
-                                  sample_times=sample_times)
-
-    tc_field = 'temperature'
-    tc_array = pyv.SensorArrayFactory \
-        .thermocouples_basic_errs(sim_data,
-                                  sens_data,
-                                  tc_field,
-                                  spat_dims=2,
-                                  errs_pc=1.0)
-
-    sg_field = 'strain'
-    sg_array = pyv.SensorArrayFactory \
-        .strain_gauges_basic_errs(sim_data,
-                                  sens_data,
-                                  sg_field,
-                                  spat_dims=2,
-                                  errs_pc=1.0)
-
-    sensor_arrays = [tc_array,sg_array]
-
-    #===========================================================================
-    # Create and run the simulated experiment
-    exp_sim = pyv.ExperimentSimulator(sim_list,
-                                         sensor_arrays,
-                                         num_exp_per_sim=1000)
-
-    exp_data = exp_sim.run_experiments()
-    exp_stats = exp_sim.calc_stats()
-
-    #===========================================================================
-    print(80*"=")
-    print("exp_data and exp_stats are lists where the index is the sensor array")
-    print("position in the list as field components are not consistent dims.\n")
-
-    print(80*"-")
-    print("Thermal sensor array @ exp_data[0]")
-    print(80*"-")
-    print("shape=(n_sims,n_exps,n_sensors,n_field_comps,n_time_steps)")
-    print(f"{exp_data[0].shape=}")
-    print()
-    print("Stats are calculated over all experiments (axis=1)")
-    print("shape=(n_sims,n_sensors,n_field_comps,n_time_steps)")
-    print(f"{exp_stats[0].max.shape=}")
-    print()
-    print(80*"-")
-    print("Mechanical sensor array @ exp_data[1]")
-    print(80*"-")
-    print("shape=(n_sims,n_exps,n_sensors,n_field_comps,n_time_steps)")
-    print(f"{exp_data[1].shape=}")
-    print()
-    print("shape=(n_sims,n_sensors,n_field_comps,n_time_steps)")
-    print(f"{exp_stats[1].max.shape=}")
-    print(80*"=")
-
-    #===========================================================================
-    # VISUALISE RESULTS
-    (fig,ax) = pyv.plot_exp_traces(exp_sim,
-                                      component="temperature",
-                                      sens_array_num=0,
-                                      sim_num=0)
-
-    (fig,ax) = pyv.plot_exp_traces(exp_sim,
-                                    component="strain_yy",
-                                    sens_array_num=1,
-                                    sim_num=2)
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()