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

pyvale/examples/ex1_2_thermal2d.py

@@ -1,108 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-import matplotlib.pyplot as plt
-import mooseherder as mh
-import pyvale as pyv
-
-
-def main() -> None:
-    """pyvale example: Point sensors on a 2D thermal simulation
-    ----------------------------------------------------------------------------
-    - Explanation of the usage of "get_measurements()" and "calc_measurements()"
-
-    NOTES:
-    - A sensor measurement is defined as:
-        measurement = truth + systematic error + random error.
-    - Calling the "get" methods of the sensor array will retrieve the results
-      for the current experiment.
-    - Calling the "calc" methods will generate a new
-      experiment by sampling/calculating the systematic and random errors.
-    """
-    data_path = pyv.DataSet.thermal_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    field_key = list(sim_data.node_vars.keys())[0] # type: ignore
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    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)
-    sens_data = pyv.SensorData(positions=sens_pos)
-
-    tc_array = pyv.SensorArrayFactory \
-        .thermocouples_basic_errs(sim_data,
-                                  sens_data,
-                                  field_key,
-                                  spat_dims=2)
-
-
-    measurements = tc_array.get_measurements()
-
-    print("\n"+80*"-")
-    print("For a sensor: measurement = truth + sysematic error + random error")
-    print(f"measurements.shape = {measurements.shape} = "+
-          "(n_sensors,n_field_components,n_timesteps)\n")
-    print("The truth, systematic error and random error arrays have the same "+
-          "shape.")
-
-    print(80*"-")
-    print("Looking at the last 5 time steps (measurements) of sensor 0:")
-    pyv.print_measurements(tc_array,
-                              (0,1),
-                              (0,1),
-                              (measurements.shape[2]-5,measurements.shape[2]))
-    print(80*"-")
-    print("If we call the \"calc_measurements()\" method then the errors are "+
-          "(re)calculated or sampled.")
-    measurements = tc_array.calc_measurements()
-
-    pyv.print_measurements(tc_array,
-                              (0,1),
-                              (0,1),
-                              (measurements.shape[2]-5,measurements.shape[2]))
-
-
-    (_,ax) = pyv.plot_time_traces(tc_array,field_key)
-    ax.set_title("Exp 1: called calc_measurements()")
-
-    print(80*"-")
-    print("If we call the \"get_measurements()\" method then the errors are the "+
-          "same:")
-    measurements = tc_array.get_measurements()
-
-    pyv.print_measurements(tc_array,
-                              (0,1),
-                              (0,1),
-                              (measurements.shape[2]-5,measurements.shape[2]))
-
-    (_,ax) = pyv.plot_time_traces(tc_array,field_key)
-    ax.set_title("Exp 2: called get_measurements()")
-
-    print(80*"-")
-    print("If we call the \"calc_measurements()\" method again we generate/sample"+
-           "new errors:")
-    measurements = tc_array.calc_measurements()
-
-    pyv.print_measurements(tc_array,
-                              (0,1),
-                              (0,1),
-                              (measurements.shape[2]-5,measurements.shape[2]))
-
-    (_,ax) = pyv.plot_time_traces(tc_array,field_key)
-    ax.set_title("Exp 3: called calc_measurements()")
-
-    print(80*"-")
-
-    plot_on = True
-    if plot_on:
-        plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex1_3_thermal2d.py

@@ -1,110 +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: Point sensors on a 2D thermal simulation
-    ----------------------------------------------------------------------------
-    - Full construction of a point sensor array from scratch
-    - Explanation of the different types of error models
-    - There are flags throughout the example allowing the user to toggle on/off
-      parts of the sensor array construction
-    """
-    data_path = pyv.DataSet.thermal_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
-
-    use_auto_descriptor = "blank"
-    if use_auto_descriptor == "factory":
-        descriptor = pyv.SensorDescriptorFactory.temperature_descriptor()
-    elif use_auto_descriptor == "manual":
-        descriptor = pyv.SensorDescriptor()
-        descriptor.name = "Temperature"
-        descriptor.symbol = "T"
-        descriptor.units = r"^{\circ}C"
-        descriptor.tag = "TC"
-    else:
-        descriptor = pyv.SensorDescriptor()
-
-    field_key = "temperature" # ("disp_x","disp_y")
-    t_field = pyv.FieldScalar(sim_data,
-                                 field_key=field_key,
-                                 spat_dims=2)
-
-    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 = False
-    if use_sim_time:
-        sample_times = None
-    else:
-        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
-
-    sensor_data = pyv.SensorData(positions=sens_pos,
-                                   sample_times=sample_times)
-
-    tc_array = pyv.SensorArrayPoint(sensor_data,
-                                       t_field,
-                                       descriptor)
-
-    errors_on = {"indep_sys": True,
-                 "rand": True,
-                 "dep_sys": True}
-
-    error_chain = []
-    if errors_on["indep_sys"]:
-        error_chain.append(pyv.ErrSysOffset(offset=-5.0))
-        error_chain.append(pyv.ErrSysUniform(low=-5.0,
-                                                high=5.0))
-
-    if errors_on["rand"]:
-        error_chain.append(pyv.ErrRandNormPercent(std_percent=1.0))
-        error_chain.append(pyv.ErrRandUnifPercent(low_percent=-1.0,
-                                            high_percent=1.0))
-
-    if errors_on["dep_sys"]:
-        error_chain.append(pyv.ErrSysDigitisation(bits_per_unit=2**8/100))
-        error_chain.append(pyv.ErrSysSaturation(meas_min=0.0,meas_max=300.0))
-
-    if len(error_chain) > 0:
-        error_integrator = pyv.ErrIntegrator(error_chain,
-                                                  sensor_data,
-                                                  tc_array.get_measurement_shape())
-        tc_array.set_error_integrator(error_integrator)
-
-
-    measurements = tc_array.get_measurements()
-
-    print("\n"+80*"-")
-    print("For a sensor: measurement = truth + sysematic error + random error")
-    print(f"measurements.shape = {measurements.shape} = "+
-          "(n_sensors,n_field_components,n_timesteps)\n")
-    print("The truth, systematic error and random error arrays have the same "+
-          "shape.")
-
-    print(80*"-")
-    print("Looking at the last 5 time steps (measurements) of sensor 0:")
-    pyv.print_measurements(tc_array,
-                              (0,1),
-                              (0,1),
-                              (measurements.shape[2]-5,measurements.shape[2]))
-    print(80*"-")
-
-    pyv.plot_time_traces(tc_array,field_key)
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex1_5_thermal2d.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
-import mooseherder as mh
-import pyvale as pyv
-
-
-def main() -> None:
-    """pyvale example: point sensors on a 2D thermal simulation
-    ----------------------------------------------------------------------------
-    """
-    data_path = pyv.DataSet.thermal_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    field_key = list(sim_data.node_vars.keys())[0] # type: ignore
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    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)
-
-    sample_times = np.linspace(0.0,np.max(sim_data.time),50) # | None
-
-    sensor_data = pyv.SensorData(positions=sens_pos,
-                                    sample_times=sample_times)
-
-    tc_array = pyv.SensorArrayFactory \
-        .thermocouples_no_errs(sim_data,
-                               sensor_data,
-                               field_key,
-                               spat_dims=2)
-
-    #===========================================================================
-    # Examples of full error library
-
-    #---------------------------------------------------------------------------
-    # Standard independent systematic errors
-    err_chain = []
-    err_chain.append(pyv.ErrSysOffset(offset=-1.0))
-    err_chain.append(pyv.ErrSysOffsetPercent(offset_percent=-1.0))
-
-    err_chain.append(pyv.ErrSysUniform(low=-2.0,
-                                        high=2.0))
-    err_chain.append(pyv.ErrSysUniformPercent(low_percent=-2.0,
-                                                    high_percent=2.0))
-
-    err_chain.append(pyv.ErrSysNormal(std=1.0))
-    err_chain.append(pyv.ErrSysNormPercent(std_percent=2.0))
-
-    sys_gen = pyv.GeneratorTriangular(left=-1.0,
-                                         mode=0.0,
-                                         right=1.0)
-    err_chain.append(pyv.ErrSysGenerator(sys_gen))
-
-    #---------------------------------------------------------------------------
-    err_chain.append(pyv.ErrRandNormal(std = 2.0))
-    err_chain.append(pyv.ErrRandNormPercent(std_percent=2.0))
-
-    err_chain.append(pyv.ErrRandUniform(low=-2.0,high=2.0))
-    err_chain.append(pyv.ErrRandUnifPercent(low_percent=-2.0,
-                                               high_percent=2.0))
-
-    rand_gen = pyv.GeneratorTriangular(left=-5.0,
-                                          mode=0.0,
-                                          right=5.0)
-    err_chain.append(pyv.ErrRandGenerator(rand_gen))
-
-    #---------------------------------------------------------------------------
-    err_chain.append(pyv.ErrSysDigitisation(bits_per_unit=2**8/100))
-    err_chain.append(pyv.ErrSysSaturation(meas_min=0.0,meas_max=300.0))
-
-    err_int = pyv.ErrIntegrator(err_chain,
-                                     sensor_data,
-                                     tc_array.get_measurement_shape())
-    tc_array.set_error_integrator(err_int)
-
-
-    #===========================================================================
-
-    measurements = tc_array.calc_measurements()
-    print(80*'-')
-    sens_num = 4
-    print('The last 5 time steps (measurements) of sensor {sens_num}:')
-    pyv.print_measurements(tc_array,
-                              (sens_num-1,sens_num),
-                              (0,1),
-                              (measurements.shape[2]-5,measurements.shape[2]))
-    print(80*'-')
-
-    pyv.plot_time_traces(tc_array,field_key)
-    plt.show()
-
-
-if __name__ == '__main__':
-    main()

pyvale/examples/ex2_1_thermal3d .py

@@ -1,84 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-from pathlib import Path
-import matplotlib.pyplot as plt
-import mooseherder as mh
-import pyvale as pyv
-
-
-def main() -> None:
-    """pyvale example: thermocouples on a 3D divertor monoblock heatsink
-    ----------------------------------------------------------------------------
-    """
-    #data_path = pyv.DataSet.thermal_3d_path()
-    data_path = Path.cwd()/"src"/"pyvale"/"simcases"/"case16_out.e"
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    field_name = 'temperature'
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    n_sens = (1,4,1)
-    x_lims = (12.5,12.5)
-    y_lims = (0.0,33.0)
-    z_lims = (0.0,12.0)
-    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
-
-    sens_data = pyv.SensorData(positions=sens_pos)
-
-    tc_array = pyv.SensorArrayFactory \
-        .thermocouples_basic_errs(sim_data,
-                                  sens_data,
-                                  field_name,
-                                  spat_dims=3)
-
-    measurements = tc_array.get_measurements()
-    print(f'\nMeasurements for sensor at top of block:\n{measurements[-1,0,:]}\n')
-
-    pv_plot = pyv.plot_point_sensors_on_sim(tc_array,field_name)
-
-    # Set this to 'interactive' to get an interactive 3D plot of the simulation
-    # and labelled sensor locations, set to 'save_fig' to create a vector
-    # graphic using a specified camera position.
-    pv_plot_mode = 'interactive'
-
-    if pv_plot_mode == 'interactive':
-        pv_plot.camera_position = [(59.354, 43.428, 69.946),
-                                   (-2.858, 13.189, 4.523),
-                                   (-0.215, 0.948, -0.233)]
-        pv_plot.show()
-
-        print(80*"=")
-        print('Camera positions = ')
-        print(pv_plot.camera_position)
-        print(80*"="+"\n")
-
-    if pv_plot_mode == 'save_fig':
-        # Determined manually by moving camera and then dumping camera position
-        # to console after window close - see 'interactive above'
-        pv_plot.camera_position = [(59.354, 43.428, 69.946),
-                                   (-2.858, 13.189, 4.523),
-                                   (-0.215, 0.948, -0.233)]
-        save_render = Path('src/examples/monoblock_thermal_sim_view.svg')
-        pv_plot.save_graphic(save_render) # only for .svg .eps .ps .pdf .tex
-        pv_plot.screenshot(save_render.with_suffix('.png'))
-
-    # Set this to 'interactive' to get a matplotlib.pyplot with the sensor
-    # traces plotted over time. Set to 'save_fig' to save an image of the plot
-    # to file.
-    trace_plot_mode = 'interactive'
-
-    (fig,_) = pyv.plot_time_traces(tc_array,field_name)
-
-    if trace_plot_mode == 'interactive':
-        plt.show()
-    if trace_plot_mode == 'save_fig':
-        save_traces = Path('src/examples/monoblock_thermal_traces.png')
-        fig.savefig(save_traces, dpi=300, format='png', bbox_inches='tight')
-
-
-if __name__ == '__main__':
-    main()

pyvale/examples/ex2_2_thermal3d.py

@@ -1,51 +0,0 @@
-'''
-================================================================================
-Example: 3d thermocouples on a monoblock
-
-pyvale: the python validation engine
-License: MIT
-Copyright (C) 2025 The Computer Aided Validation Team
-================================================================================
-'''
-import mooseherder as mh
-import pyvale as pyv
-
-
-def main() -> None:
-    """pyvale example: thermocouples on a 3D divertor monoblock heatsink
-    ----------------------------------------------------------------------------
-    """
-    data_path = pyv.DataSet.thermal_3d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    field_name = 'temperature'
-
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    n_sens = (1,4,1)
-    x_lims = (12.5,12.5)
-    y_lims = (0,31.0)
-    z_lims = (0.0,12.5)
-    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
-
-    sens_data = pyv.SensorData(positions=sens_pos)
-
-    tc_array = pyv.SensorArrayFactory() \
-        .thermocouples_basic_errs(sim_data,
-                                  sens_data,
-                                  field_name,
-                                  spat_dims=3)
-
-
-    measurements = tc_array.get_measurements()
-    print(f'\nMeasurements for sensor at top of block:\n{measurements[-1,0,:]}\n')
-
-    pv_plot = pyv.plot_point_sensors_on_sim(tc_array,field_name)
-    pv_plot.camera_position = [(59.354, 43.428, 69.946),
-                                (-2.858, 13.189, 4.523),
-                                (-0.215, 0.948, -0.233)]
-    pv_plot.show()
-
-
-if __name__ == '__main__':
-    main()

pyvale/examples/ex2_3_thermal3d.py

@@ -1,106 +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: thermocouples on a 3D divertor monoblock heatsink
-    ----------------------------------------------------------------------------
-    """
-    data_path = pyv.DataSet.thermal_3d_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
-
-    use_auto_descriptor = "manual"
-    if use_auto_descriptor == "factory":
-        descriptor = pyv.SensorDescriptorFactory.temperature_descriptor()
-    elif use_auto_descriptor == "manual":
-        descriptor = pyv.SensorDescriptor()
-        descriptor.name = "Temperature"
-        descriptor.symbol = "T"
-        descriptor.units = r"^{\circ}C"
-        descriptor.tag = "TC"
-    else:
-        descriptor = pyv.SensorDescriptor()
-
-
-    field_key = "temperature"
-    t_field = pyv.FieldScalar(sim_data,
-                                 field_key=field_key,
-                                 spat_dims=3)
-
-    n_sens = (1,4,1)
-    x_lims = (12.5,12.5)
-    y_lims = (0.0,33.0)
-    z_lims = (0.0,12.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),80)
-
-    sens_data = pyv.SensorData(positions=sens_pos,
-                                  sample_times=sample_times)
-
-    tc_array = pyv.SensorArrayPoint(sens_data,
-                                       t_field,
-                                       descriptor)
-
-    errors_on = {"indep_sys": True,
-                 "rand": True,
-                 "dep_sys": True}
-
-    error_chain = []
-    if errors_on["indep_sys"]:
-        error_chain.append(pyv.ErrSysOffset(offset=-5.0))
-        error_chain.append(pyv.ErrSysUniform(low=-10.0,
-                                                   high=10.0))
-
-    if errors_on["rand"]:
-        error_chain.append(pyv.ErrRandNormPercent(std_percent=5.0))
-        error_chain.append(pyv.ErrRandUnifPercent(low_percent=-5.0,
-                                                   high_percent=5.0))
-
-    if errors_on["dep_sys"]:
-        error_chain.append(pyv.ErrSysDigitisation(bits_per_unit=1/20))
-        error_chain.append(pyv.ErrSysSaturation(meas_min=0.0,meas_max=800.0))
-
-    if len(error_chain) > 0:
-        error_integrator = pyv.ErrIntegrator(
-            error_chain,
-            sens_data,
-            tc_array.get_measurement_shape(),
-        )
-        tc_array.set_error_integrator(error_integrator)
-
-    measurements = tc_array.get_measurements()
-    print(f"\nMeasurements for sensor at top of block:\n{measurements[-1,0,:]}\n")
-
-
-    (fig,_) = pyv.plot_time_traces(tc_array,field_key)
-    plt.show()
-
-    save_fig = True
-    if save_fig:
-        image_path = Path.cwd() / "example_output"
-        image_file = image_path / "monoblock_thermal_traces_syserrs.png"
-        if not image_path.is_dir():
-            image_path.mkdir()
-
-        fig.savefig(image_file, dpi=300, format="png", bbox_inches="tight")
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex3_1_displacement2d.py

@@ -1,44 +0,0 @@
-# ================================================================================
-# pyvale: the python validation engine
-# License: MIT
-# Copyright (C) 2025 The Computer Aided Validation Team
-# ================================================================================
-
-import matplotlib.pyplot as plt
-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
-
-    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)
-
-    sens_data = pyv.SensorData(positions=sens_pos)
-
-    disp_sens_array = pyv.SensorArrayFactory \
-                            .disp_sensors_basic_errs(sim_data,
-                                                     sens_data,
-                                                     "displacement",
-                                                     spat_dims=2)
-
-    plot_field = 'disp_x'
-    pv_plot = pyv.plot_point_sensors_on_sim(disp_sens_array,plot_field)
-    pv_plot.show(cpos="xy")
-
-    pyv.plot_time_traces(disp_sens_array,'disp_x')
-    pyv.plot_time_traces(disp_sens_array,'disp_y')
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex3_2_displacement2d.py

@@ -1,76 +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: 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.SensorDescriptor()
-    descriptor.name = 'Disp.'
-    descriptor.symbol = r'u'
-    descriptor.units = r'm'
-    descriptor.tag = 'DS'
-    descriptor.components = ('x','y','z')
-
-    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,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)
-
-    disp_sens_array = pyv.SensorArrayPoint(sens_data,
-                                              disp_field,
-                                              descriptor)
-
-    error_chain = []
-    error_chain.append(pyv.ErrSysUniform(low=-0.01e-3,high=0.01e-3))
-    error_chain.append(pyv.ErrRandNormal(std=0.01e-3))
-    error_int = pyv.ErrIntegrator(error_chain,
-                                       sens_data,
-                                       disp_sens_array.get_measurement_shape())
-    disp_sens_array.set_error_integrator(error_int)
-
-
-    plot_field = 'disp_x'
-    if plot_field == 'disp_x':
-        pv_plot = pyv.plot_point_sensors_on_sim(disp_sens_array,'disp_x')
-        pv_plot.show(cpos="xy")
-    elif plot_field == 'disp_y':
-        pv_plot = pyv.plot_point_sensors_on_sim(disp_sens_array,'disp_y')
-        pv_plot.show(cpos="xy")
-
-    pyv.plot_time_traces(disp_sens_array,'disp_x')
-    pyv.plot_time_traces(disp_sens_array,'disp_y')
-    plt.show()
-
-
-if __name__ == "__main__":
-    main()