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

pyvale/core/visualtraceplotter.py

@@ -1,256 +0,0 @@
-"""
-================================================================================
-pyvale: the python validation engine
-License: MIT
-Copyright (C) 2025 The Computer Aided Validation Team
-================================================================================
-"""
-from typing import Any
-
-import numpy as np
-import matplotlib.pyplot as plt
-from pyvale.core.sensorarraypoint import SensorArrayPoint
-from pyvale.core.visualopts import (PlotOptsGeneral,
-                                    TraceOptsSensor,
-                                    TraceOptsExperiment)
-from pyvale.core.experimentsimulator import ExperimentSimulator
-
-#TODO: Docstrings
-
-def plot_time_traces(sensor_array: SensorArrayPoint,
-                     component: str,
-                     trace_opts: TraceOptsSensor | None = None,
-                     plot_opts: PlotOptsGeneral | None = None
-                     ) -> tuple[Any,Any]:
-
-    #---------------------------------------------------------------------------
-    field = sensor_array.field
-    comp_ind = sensor_array.field.get_component_index(component)
-    samp_time = sensor_array.get_sample_times()
-    measurements = sensor_array.get_measurements()
-    n_sensors = sensor_array.sensor_data.positions.shape[0]
-    descriptor = sensor_array.descriptor
-    sensors_perturbed = sensor_array.get_sensor_data_perturbed()
-
-    #---------------------------------------------------------------------------
-    if plot_opts is None:
-        plot_opts = PlotOptsGeneral()
-
-    if trace_opts is None:
-        trace_opts = TraceOptsSensor()
-
-    if trace_opts.sensors_to_plot is None:
-        sensors_to_plot = range(n_sensors)
-    else:
-        sensors_to_plot = trace_opts.sensors_to_plot
-
-    #---------------------------------------------------------------------------
-    # Figure canvas setup
-    fig, ax = plt.subplots(figsize=plot_opts.single_fig_size_landscape,
-                           layout='constrained')
-    fig.set_dpi(plot_opts.resolution)
-
-    #---------------------------------------------------------------------------
-    # Plot simulation and truth lines
-    if trace_opts.sim_line is not None:
-        sim_time = field.get_time_steps()
-        sim_vals = field.sample_field(sensor_array.sensor_data.positions,
-                                      None,
-                                      sensor_array.sensor_data.angles)
-
-        for ii,ss in enumerate(sensors_to_plot):
-            ax.plot(sim_time,
-                    sim_vals[ss,comp_ind,:],
-                    trace_opts.sim_line,
-                    lw=plot_opts.lw,
-                    ms=plot_opts.ms,
-                    color=plot_opts.colors[ii % plot_opts.n_colors])
-
-    if trace_opts.truth_line is not None:
-        truth = sensor_array.get_truth()
-        for ii,ss in enumerate(sensors_to_plot):
-            ax.plot(samp_time,
-                    truth[ss,comp_ind,:],
-                    trace_opts.truth_line,
-                    lw=plot_opts.lw,
-                    ms=plot_opts.ms,
-                    color=plot_opts.colors[ii % plot_opts.n_colors])
-
-    sensor_tags = descriptor.create_sensor_tags(n_sensors)
-    for ii,ss in enumerate(sensors_to_plot):
-        sensor_time = samp_time
-        if sensors_perturbed is not None:
-            if sensors_perturbed.sample_times is not None:
-                sensor_time = sensors_perturbed.sample_times
-
-        ax.plot(sensor_time,
-                measurements[ss,comp_ind,:],
-                trace_opts.meas_line,
-                label=sensor_tags[ss],
-                lw=plot_opts.lw,
-                ms=plot_opts.ms,
-                color=plot_opts.colors[ii % plot_opts.n_colors])
-
-    #---------------------------------------------------------------------------
-    # Axis / legend labels and options
-    ax.set_xlabel(trace_opts.time_label,
-                fontsize=plot_opts.font_ax_size, fontname=plot_opts.font_name)
-    ax.set_ylabel(descriptor.create_label(comp_ind),
-                fontsize=plot_opts.font_ax_size, fontname=plot_opts.font_name)
-
-    if trace_opts.time_min_max is None:
-        min_time = np.min((np.min(samp_time),np.min(sensor_time)))
-        max_time = np.max((np.max(samp_time),np.max(sensor_time)))
-        ax.set_xlim((min_time,max_time)) # type: ignore
-    else:
-        ax.set_xlim(trace_opts.time_min_max)
-
-    if trace_opts.legend:
-        ax.legend(prop={"size":plot_opts.font_leg_size},loc='best')
-
-    plt.grid(True)
-    plt.draw()
-
-    return (fig,ax)
-
-
-def plot_exp_traces(exp_sim: ExperimentSimulator,
-                    component: str,
-                    sens_array_num: int,
-                    sim_num: int,
-                    trace_opts: TraceOptsExperiment | None = None,
-                    plot_opts: PlotOptsGeneral | None = None) -> tuple[Any,Any]:
-    #---------------------------------------------------------------------------
-    n_sensors = exp_sim.sensor_arrays[sens_array_num].sensor_data.positions.shape[0]
-    descriptor = exp_sim.sensor_arrays[sens_array_num].descriptor
-    comp_ind = exp_sim.sensor_arrays[sens_array_num].field.get_component_index(component)
-    samp_time = exp_sim.sensor_arrays[sens_array_num].get_sample_times()
-    num_sens = exp_sim.sensor_arrays[sens_array_num].get_measurement_shape()[0]
-
-    #---------------------------------------------------------------------------
-    if trace_opts is None:
-        trace_opts = TraceOptsExperiment()
-
-    if plot_opts is None:
-        plot_opts = PlotOptsGeneral()
-
-    #---------------------------------------------------------------------------
-    exp_data = exp_sim.get_data()
-    exp_stats = exp_sim.get_stats()
-
-    if trace_opts.sensors_to_plot is None:
-        sensors_to_plot = range(num_sens)
-    else:
-        sensors_to_plot = trace_opts.sensors_to_plot
-
-    #---------------------------------------------------------------------------
-    # Figure canvas setup
-    fig, ax = plt.subplots(figsize=plot_opts.single_fig_size_landscape,
-                           layout='constrained')
-    fig.set_dpi(plot_opts.resolution)
-
-    #---------------------------------------------------------------------------
-    # Plot all simulated experimental points
-    if trace_opts.plot_all_exp_points:
-        for ii,ss in enumerate(sensors_to_plot):
-            for ee in range(exp_sim.num_exp_per_sim):
-                ax.plot(samp_time,
-                        exp_data[sens_array_num][sim_num,ee,ss,comp_ind,:],
-                        "+",
-                        lw=plot_opts.lw,
-                        ms=plot_opts.ms,
-                        color=plot_opts.colors[ii % plot_opts.n_colors])
-
-    #---------------------------------------------------------------------------
-    sensor_tags = descriptor.create_sensor_tags(n_sensors)
-    for ii,ss in enumerate(sensors_to_plot):
-        if trace_opts.centre == "median":
-            ax.plot(samp_time,
-                    exp_stats[sens_array_num].median[sim_num,ss,comp_ind,:],
-                    trace_opts.exp_mean_line,
-                    label=sensor_tags[ss],
-                    lw=plot_opts.lw,
-                    ms=plot_opts.ms,
-                    color=plot_opts.colors[ii % plot_opts.n_colors])
-        else:
-            ax.plot(samp_time,
-                    exp_stats[sens_array_num].mean[sim_num,ss,comp_ind,:],
-                    trace_opts.exp_mean_line,
-                    label=sensor_tags[ss],
-                    lw=plot_opts.lw,
-                    ms=plot_opts.ms,
-                    color=plot_opts.colors[ii % plot_opts.n_colors])
-
-        if trace_opts is not None:
-            upper = np.zeros_like(exp_stats[sens_array_num].min)
-            lower = np.zeros_like(exp_stats[sens_array_num].min)
-
-            if trace_opts.fill_between == 'max':
-                upper = exp_stats[sens_array_num].min
-                lower = exp_stats[sens_array_num].max
-
-            elif trace_opts.fill_between == 'quartile':
-                upper = exp_stats[sens_array_num].q25
-                lower = exp_stats[sens_array_num].q75
-
-            elif trace_opts.fill_between == '2std':
-                upper = exp_stats[sens_array_num].mean + \
-                        2*exp_stats[sens_array_num].std
-                lower = exp_stats[sens_array_num].mean - \
-                        2*exp_stats[sens_array_num].std
-
-            elif trace_opts.fill_between == '3std':
-                upper = exp_stats[sens_array_num].mean + \
-                        3*exp_stats[sens_array_num].std
-                lower = exp_stats[sens_array_num].mean - \
-                        3*exp_stats[sens_array_num].std
-
-            ax.fill_between(samp_time,
-                upper[sim_num,ss,comp_ind,:],
-                lower[sim_num,ss,comp_ind,:],
-                color=plot_opts.colors[ii % plot_opts.n_colors],
-                alpha=0.2)
-
-    #---------------------------------------------------------------------------
-    # Plot simulation and truth line
-    if trace_opts.sim_line is not None:
-        sim_time = exp_sim.sensor_arrays[sens_array_num].field.get_time_steps()
-        sim_vals = exp_sim.sensor_arrays[sens_array_num].field.sample_field(
-                    exp_sim.sensor_arrays[sens_array_num].positions)
-
-        for ss in sensors_to_plot:
-            ax.plot(sim_time,
-                    sim_vals[ss,comp_ind,:],
-                    trace_opts.sim_line,
-                    lw=plot_opts.lw,
-                    ms=plot_opts.ms)
-
-    if trace_opts.truth_line is not None:
-        truth = exp_sim.sensor_arrays[sens_array_num].get_truth()
-        for ss in sensors_to_plot:
-            ax.plot(samp_time,
-                    truth[ss,comp_ind,:],
-                    trace_opts.truth_line,
-                    lw=plot_opts.lw,
-                    ms=plot_opts.ms,
-                    color=plot_opts.colors[ii % plot_opts.n_colors])
-
-    #---------------------------------------------------------------------------
-    # Axis / legend labels and options
-    ax.set_xlabel(trace_opts.time_label,
-                fontsize=plot_opts.font_ax_size, fontname=plot_opts.font_name)
-    ax.set_ylabel(descriptor.create_label(comp_ind),
-                fontsize=plot_opts.font_ax_size, fontname=plot_opts.font_name)
-
-    if trace_opts.time_min_max is None:
-        ax.set_xlim((np.min(samp_time),np.max(samp_time))) # type: ignore
-    else:
-        ax.set_xlim(trace_opts.time_min_max)
-
-    if trace_opts.legend:
-        ax.legend(prop={"size":plot_opts.font_leg_size},loc='best')
-
-    plt.grid(True)
-    plt.draw()
-
-    return (fig,ax)

pyvale/examples/analyticdatagen/__init__.py

@@ -1,7 +0,0 @@
-"""
-================================================================================
-pyvale: the python validation engine
-License: MIT
-Copyright (C) 2025 The Computer Aided Validation Team
-================================================================================
-"""

pyvale/examples/ex1_1_thermal2d.py

@@ -1,89 +0,0 @@
-"""
-================================================================================
-Example: thermocouples on a 2d plate
-
-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: point sensors on a 2D thermal simulation
-    ----------------------------------------------------------------------------
-    - Quick start
-    - Basic sensor array construction using the sensor array factory
-    - Basic visualisation of sensor locations and sensor traces with the pyvale
-      wrapper for pyvista and matplotlib.
-    """
-
-    data_path = pyv.DataSet.thermal_2d_path()
-    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
-    field_key = "temperature"
-    # Scale to mm to make 3D visualisation scaling easier
-    sim_data.coords = sim_data.coords*1000.0 # type: ignore
-
-    n_sens = (3,2,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(f"\nMeasurements for last sensor:\n{measurements[-1,0,:]}\n")
-
-    pv_plot = pyv.plot_point_sensors_on_sim(tc_array,field_key)
-    # 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 = [(-7.547, 59.753, 134.52),
-                                   (41.916, 25.303, 9.297),
-                                   (0.0810, 0.969, -0.234)]
-        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 = [(-7.547, 59.753, 134.52),
-                                   (41.916, 25.303, 9.297),
-                                   (0.0810, 0.969, -0.234)]
-        save_render = Path("src/examples/plate_thermal_2d_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_key)
-
-    if trace_plot_mode == "interactive":
-        plt.show()
-    if trace_plot_mode == "save_fig":
-        save_traces = Path("src/examples/plate_thermal_2d_traces.png")
-        fig.savefig(save_traces, dpi=300, format="png", bbox_inches="tight")
-
-
-if __name__ == "__main__":
-    main()

pyvale/examples/ex1_2_thermal2d.py

@@ -1,111 +0,0 @@
-"""
-================================================================================
-Example: thermocouples on a 2d plate
-
-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,113 +0,0 @@
-"""
-================================================================================
-Example: thermocouples on a 2d plate
-
-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,105 +0,0 @@
-'''
-================================================================================
-Example: thermocouples on a 2d plate
-
-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()