physioblocks 1.0.1__py3-none-any.whl → 1.0.3__py3-none-any.whl

physioblocks/__init__.py

@@ -26,7 +26,7 @@
 
 """Physioblocks package definition"""
 
-__version__ = "1.0.1"
+__version__ = "1.0.3"
 __copyright__ = "INRIA"
 __license__ = "LGPL-3.0-only"
 __authors__ = [

physioblocks/configuration/constants.py

@@ -38,11 +38,8 @@ PARAMETERS_ID = "parameters"
 # Definition of the flux-dof types couples
 FLUX_DOF_DEFINITION_ID = "flux_dof_definitions"
 
-
-# Magnitudes
-
-# The variable magnitude item label in the configuration
-MAGNITUDES = "magnitudes"
+# Options label for specific simulation types
+SIMULATION_OPTIONS = "simulation_options"
 
 # The variable magnitude item label in the configuration
 VARIABLES_MAGNITUDES = "variables_magnitudes"

physioblocks/configuration/simulation/simulations.py

@@ -35,10 +35,10 @@ from typing import Any
 from physioblocks.configuration.base import Configuration, ConfigurationError
 from physioblocks.configuration.constants import (
     INIT_VARIABLES_ID,
-    MAGNITUDES,
     NET_ID,
     OUTPUTS_FUNCTIONS_ID,
     PARAMETERS_ID,
+    SIMULATION_OPTIONS,
     SOLVER_ID,
     TIME_MANAGER_ID,
     VARIABLES_MAGNITUDES,
@@ -86,16 +86,15 @@ def load_simulation_config(
         if SOLVER_ID in config:
             solver = load(config[SOLVER_ID])
 
-        # magnitudes
-        magnitudes = None
-        if VARIABLES_MAGNITUDES in config:
-            magnitudes = load(config[VARIABLES_MAGNITUDES])
+        simulation_options = None
+        if SIMULATION_OPTIONS in config:
+            simulation_options = load(config[SIMULATION_OPTIONS])
 
         sim_factory = SimulationFactory(
             configuration_type,
             solver,
             net,
-            simulation_options={MAGNITUDES: magnitudes},
+            simulation_options=simulation_options,
         )
 
         configuration_object = sim_factory.create_simulation()
@@ -147,7 +146,7 @@ def save_simulation_config(
                 type(variable_init_values).__name__,
             )
         )
-    sim_config[VARIABLES_MAGNITUDES] = save(simulation.magnitudes)
+    sim_config[VARIABLES_MAGNITUDES] = save(simulation.state.magnitudes)
 
     # Parameters
     # Get quantities
@@ -226,6 +225,11 @@ def _configure_simulation(
         configuration_references=simulation.quantities,
     )
 
+    # magnitudes
+    if VARIABLES_MAGNITUDES in config:
+        magnitudes = load(config[VARIABLES_MAGNITUDES])
+        simulation.state.set_variables_magnitudes(magnitudes)
+
     references.update(simulation.quantities)
     references.update(simulation.models)
 

physioblocks/library/aliases/simulations/circulation_alone_forward_simulation.jsonc

@@ -49,7 +49,5 @@
             "inverses": ["aorta_distal.resistance"]
         }
     },
-    "solver": {
-        "iteration_max": 2 // This problem is linear
-    }
+    "solver": {"type": "linear_solver"}
 }

physioblocks/library/model_components/dynamics.py

@@ -902,7 +902,9 @@ class SphericalDynamicsModelComponent(ModelComponent):
         eq_system = EqSystem(1)
 
         state = State()
-        state.add_variable(SPHERICAL_DYNAMICS_STATIC_DISP_LOCAL_ID, 1)
+        state.add_variable(
+            SPHERICAL_DYNAMICS_STATIC_DISP_LOCAL_ID, 1, _STATIC_PROBLEM_DISP_MAG
+        )
         state[SPHERICAL_DYNAMICS_STATIC_DISP_LOCAL_ID].initialize(self.disp.current)
 
         static_block = _SphericalDynamicsStaticModelComponent(
@@ -942,11 +944,7 @@ class SphericalDynamicsModelComponent(ModelComponent):
             p_test <= static_block.pressure.current
             and p_test > _STATIC_PROBLEM_MIN_PRESSURE_STEP
         ):
-            sol = solver.solve(
-                state,
-                eq_system,
-                {SPHERICAL_DYNAMICS_STATIC_DISP_LOCAL_ID: _STATIC_PROBLEM_DISP_MAG},
-            )
+            sol = solver.solve(state, eq_system)
 
             if sol.converged is True:
                 state[SPHERICAL_DYNAMICS_STATIC_DISP_LOCAL_ID].initialize(sol.x[0])

physioblocks/simulation/runtime.py

@@ -96,9 +96,6 @@ class AbstractSimulation(ABC):
 
     :param eq_system: the equation system to solve at each time step
     :type eq_system: EqSystem
-
-    :param magnitudes: magnitude of the state variables
-    :type magnitudes: dict[str, float]
     """
 
     def __init__(
@@ -111,7 +108,6 @@ class AbstractSimulation(ABC):
         models: dict[str, ModelComponent],
         solver: AbstractSolver,
         eq_system: EqSystem,
-        magnitudes: dict[str, float] | None = None,
     ):
         self.factory = factory
         self.state = state
@@ -121,9 +117,6 @@ class AbstractSimulation(ABC):
         self.time_manager = time_manager
         self.solver = solver
         self.eq_system = eq_system
-        if magnitudes is None:
-            magnitudes = {}
-        self.magnitudes = self._check_magnitudes(magnitudes, state)
         self._timed_updates: dict[str, AbstractFunction] = {}
         self._output_functions_updates: dict[str, AbstractFunction] = {}
 
@@ -272,33 +265,6 @@ class AbstractSimulation(ABC):
         self.time_manager.time.initialize(self.time_manager.start)
         self.state.set_state_vector(self._initial_state)
 
-    def _check_magnitudes(
-        self, magnitudes: dict[str, float], state: State
-    ) -> dict[str, float]:
-        checked_magnitudes = {}
-
-        for variable_id in state:
-            if variable_id not in magnitudes:
-                message = str.format(
-                    "No magnitude initialized for variable {0}. Magnitude set to 1.0",
-                    variable_id,
-                )
-                _logger.warning(message)
-                checked_magnitudes[variable_id] = 1.0
-
-            elif magnitudes[variable_id] == 0.0:
-                message = str.format(
-                    "Magnitude for variable {0} is initialized to 0.0. "
-                    "Replacing with 1.0",
-                    variable_id,
-                )
-                _logger.warning(message)
-                checked_magnitudes[variable_id] = 1.0
-            else:
-                checked_magnitudes[variable_id] = magnitudes[variable_id]
-
-        return checked_magnitudes
-
     @abstractmethod
     def run(self) -> Results:
         """
@@ -406,7 +372,7 @@ class ForwardSimulation(AbstractSimulation):
                 ):
                     self.state.reset_state_vector()
 
-                    sol = self.solver.solve(self.state, self.eq_system, self.magnitudes)
+                    sol = self.solver.solve(self.state, self.eq_system)
 
                     if sol.converged is False:
                         inter_time = 0.5 * self.time_manager.current_step_size

physioblocks/simulation/solvers.py

@@ -81,31 +81,8 @@ class AbstractSolver(ABC):
         self.tolerance = tolerance
         self.iteration_max = iteration_max
 
-    def _get_state_magnitude(
-        self, state: State, magnitudes: dict[str, float] | None = None
-    ) -> NDArray[np.float64]:
-        if magnitudes is None:
-            return np.ones(
-                state.size,
-            )
-
-        mag_dict = {}
-        for var_mag_key, var_mag_value in magnitudes.items():
-            var_index = state.get_variable_index(var_mag_key)
-            mag_dict[var_index] = var_mag_value
-        sorted_mag = sorted(mag_dict.items())
-        state_mag_list = [x[1] for x in sorted_mag]
-        return np.array(
-            state_mag_list,
-        )
-
     @abstractmethod
-    def solve(
-        self,
-        state: State,
-        system: EqSystem,
-        magnitudes: dict[str, float] | None = None,
-    ) -> Solution:
+    def solve(self, state: State, system: EqSystem) -> Solution:
         """
         Child classes have to override this method
 
@@ -114,6 +91,31 @@ class AbstractSolver(ABC):
         """
 
 
+LINEAR_SOLVER = "linear_solver"
+
+
+@register_type(LINEAR_SOLVER)
+class LinearSolver(AbstractSolver):
+    """
+    Implementation of the :class:`~.AbstractSolver` class for linear problems.
+    """
+
+    def solve(self, state: State, system: EqSystem) -> Solution:
+        """
+        Solve the equation system directly
+
+        :return: the solution
+        :rtype: Solution
+        """
+        res = system.compute_residual()
+        grad = system.compute_gradient()
+
+        xn = np.linalg.solve(grad, res)
+        x = state.state_vector - xn
+
+        return Solution(x, True)
+
+
 # Type id for the Newton Solver
 NEWTON_SOLVER_TYPE_ID = "newton_solver"
 
@@ -166,12 +168,7 @@ class NewtonSolver(AbstractSolver):
         grad_rescaled = gradient * grad_mag_inv
         return res_rescaled, grad_rescaled
 
-    def solve(
-        self,
-        state: State,
-        system: EqSystem,
-        magnitudes: dict[str, float] | None = None,
-    ) -> Solution:
+    def solve(self, state: State, system: EqSystem) -> Solution:
         """
         Solve the equation system using the Newton method.
 
@@ -183,15 +180,16 @@ class NewtonSolver(AbstractSolver):
             try:
                 i = 0
                 # initialize residual and magnitude
-                state_mag = self._get_state_magnitude(state, magnitudes)
                 res = np.ones(state.state_vector.shape)
 
                 # step 0 outside ou the loop to compute the residual and gradient
                 # magnitude
                 res, grad = self._compute_residual_and_gradient(system)
-                res_mag_inv, grad_mag_inv = self._compute_res_grad_mag(grad, state_mag)
+                res_mag_inv, grad_mag_inv = self._compute_res_grad_mag(
+                    grad, state.state_magnitudes
+                )
                 res, grad = self._rescale_res_grad(res, res_mag_inv, grad, grad_mag_inv)
-                x = self._compute_new_state(state, res, grad, state_mag)
+                x = self._compute_new_state(state, res, grad, state.state_magnitudes)
                 state.update_state_vector(x)
 
                 # Begin loop at iteration 1 (0 already done)
@@ -204,7 +202,9 @@ class NewtonSolver(AbstractSolver):
                     res, grad = self._rescale_res_grad(
                         res, res_mag_inv, grad, grad_mag_inv
                     )
-                    x = self._compute_new_state(state, res, grad, state_mag)
+                    x = self._compute_new_state(
+                        state, res, grad, state.state_magnitudes
+                    )
                     state.update_state_vector(x)
                     i += 1
 

physioblocks/simulation/state.py

@@ -28,7 +28,7 @@
 Define the **State** that holds simulation variables.
 """
 
-from collections.abc import Callable, Generator, Mapping
+from collections.abc import Callable, Generator, Mapping, Sized
 from pprint import pformat
 from typing import Any
 
@@ -53,8 +53,12 @@ class State:
     _variables: dict[str, Quantity[Any]]
     """The variables ids and quantities values"""
 
+    _magnitudes: dict[str, Any]
+    """The variables magnitudes"""
+
     def __init__(self) -> None:
         self._variables = {}
+        self._magnitudes = {}
 
     @property
     def size(self) -> int:
@@ -76,6 +80,16 @@ class State:
         """
         return self._variables.copy()
 
+    @property
+    def magnitudes(self) -> dict[str, Any]:
+        """
+        Get a mapping of variables names and magnitudes.
+
+        :return: the variables names and magnitudes.
+        :rtype: dict[str, Quantity]
+        """
+        return self._magnitudes.copy()
+
     @property
     def state_vector(self) -> NDArray[np.float64]:
         """
@@ -91,8 +105,23 @@ class State:
         else:
             return np.array([])
 
-    def __array__(self) -> NDArray[Any]:
-        return self.state_vector
+    @property
+    def state_magnitudes(self) -> NDArray[np.float64]:
+        """
+        Get the state variables magnitudes.
+
+        :return: the state variables magnitudes
+        :rtype: NDArray[np.float64]
+        """
+        if len(self._variables) > 0:
+            return np.concatenate(
+                [self._magnitudes[var_id] for var_id in self._variables], axis=None
+            )
+        else:
+            return np.array([])
+
+    def __array__(self, *args: Any, **kwargs: Any) -> NDArray[Any]:
+        return np.array(self.state_vector, *args, **kwargs)
 
     def __getitem__(self, var_id: str) -> Quantity[Any]:
         """
@@ -311,7 +340,7 @@ class State:
                 quantity_part = x[var_index : var_index + quantity.size]
                 func(quantity, quantity_part)
 
-    def add_variable(self, var_id: str, var_value: Any) -> None:
+    def add_variable(self, var_id: str, var_value: Any, magnitude: Any = None) -> None:
         """
         Add a variable to the state.
 
@@ -328,6 +357,21 @@ class State:
         quantity = var_value if isinstance(var_value, Quantity) else Quantity(var_value)
         self._variables[var_id] = quantity
 
+        if magnitude is not None:
+            if quantity.size == 1 and isinstance(magnitude, float):
+                self._magnitudes[var_id] = magnitude
+            elif isinstance(magnitude, float):
+                self._magnitudes[var_id] = [magnitude] * quantity.size
+            elif isinstance(magnitude, Sized) and len(magnitude) == quantity.size:
+                self._magnitudes[var_id] = magnitude
+            else:
+                raise ValueError("Variable and magnitude sizes mismatch")
+        else:
+            if quantity.size == 1:
+                self._magnitudes[var_id] = 1.0
+            else:
+                self._magnitudes[var_id] = [1.0] * quantity.size
+
     def remove_variable(self, var_id: str) -> None:
         """
         Remove a variable from the state
@@ -338,3 +382,22 @@ class State:
         # remove the variable
         if var_id in self._variables:
             self._variables.pop(var_id)
+
+        if var_id in self._magnitudes:
+            self._magnitudes.pop(var_id)
+
+    def set_variables_magnitudes(self, magnitudes: dict[str, Any]) -> None:
+        """
+        Update all the state variables magnitudes.
+
+        :param magnitudes: the state variables magnitudes
+        :type magnitudes: dict[str, Any]
+
+        :raise KeyError: Raise a KeyError if one or more of the magnitude parameters key
+          is not in the state.
+        """
+
+        if set(magnitudes.keys()).issubset(self._variables.keys()) is False:
+            difference = set(magnitudes.keys()).difference(self._variables.keys())
+            raise KeyError(str.format("State has no variables named {0}.", difference))
+        self._magnitudes.update(magnitudes)

physioblocks/utils/gradient_test_utils.py

@@ -77,22 +77,16 @@ def gradient_test_from_file(config_file_path: str) -> bool:
     :param config_file_path: the file path to the simulation configuration file.
     :type config_file_path: str
 
-    :return: True if the gradient test is successfull, false otherwise.
+    :return: True if the gradient test is successful, false otherwise.
     """
     configuration = read_json(config_file_path)
     configuration = unwrap_aliases(configuration)
     sim: AbstractSimulation = load(configuration)
     sim._initialize()  # noqa SLF001
-    return gradient_test(
-        sim.eq_system,
-        sim.state,
-        sim.solver._get_state_magnitude(sim.state, sim.magnitudes),  # noqa: SLF001
-    )
+    return gradient_test(sim.eq_system, sim.state)
 
 
-def gradient_test_from_model(
-    model: ModelComponent, state: State, state_magnitude: NDArray[np.float64]
-) -> bool:
+def gradient_test_from_model(model: ModelComponent, state: State) -> bool:
     """
     Create an equation system for the given block only and perform a gradient test.
 
@@ -106,10 +100,7 @@ def gradient_test_from_model(
     :param state: the state used to determine the variables in the model.
     :type state: str
 
-    :param state_magnitude: the state variables magnitudes
-    :type state_magnitude: str
-
-    :return: True if the gradient test is successfull, false otherwise.
+    :return: True if the gradient test is successful, false otherwise.
     """
     line_index = 0
     expressions = setup.SystemExpressions()
@@ -124,14 +115,11 @@ def gradient_test_from_model(
 
     eq_system = setup.build_eq_system(expressions, state)
     model.initialize()
-    return gradient_test(eq_system, state, state_magnitude)
+    return gradient_test(eq_system, state)
 
 
 def gradient_test_from_expression(
-    expr: Expression,
-    expr_params: Any,
-    state: State,
-    state_magnitude: NDArray[np.float64],
+    expr: Expression, expr_params: Any, state: State
 ) -> bool:
     """
     Create an equation system for the given expression only and perform a gradient test.
@@ -145,20 +133,15 @@ def gradient_test_from_expression(
     :param state: the state used to determine the variables in the expression.
     :type state: str
 
-    :param state_magnitude: the state variables magnitudes
-    :type state_magnitude: str
-
-    :return: True if the gradient test is successfull, false otherwise.
+    :return: True if the gradient test is successful, false otherwise.
     :rtype: bool
     """
 
     eq_system = setup.build_eq_system([(0, expr, expr_params)], state)
-    return gradient_test(eq_system, state, state_magnitude)
+    return gradient_test(eq_system, state)
 
 
-def gradient_test(
-    eq_system: EqSystem, state: State, state_magnitude: NDArray[np.float64]
-) -> bool:
+def gradient_test(eq_system: EqSystem, state: State) -> bool:
     """
     Test the computed gradient for the equation system by comparing it to
     a gradient estimated with finite differences.
@@ -170,9 +153,6 @@ def gradient_test(
     :param state: system state
     :type state: State
 
-    :param state_magnitude: the state variables magnitudes
-    :type state_magnitude: str
-
     :return: True if the estimated and computed gradient meet tolerance,
       False otherwise.
     """
@@ -180,9 +160,9 @@ def gradient_test(
     _logger.info(str.format("State:{0}{1}", linesep, state))
     _logger.info(str.format("System:{0}{1}", linesep, eq_system))
 
-    new_state = state.state_vector + _NEW_STATE_SHIFT_FACTOR * state_magnitude
+    new_state = state.state_vector + _NEW_STATE_SHIFT_FACTOR * state.state_magnitudes
 
-    shift_1 = _SHIFT_FACTOR * state_magnitude
+    shift_1 = _SHIFT_FACTOR * state.state_magnitudes
     state.update_state_vector(new_state)
 
     res = eq_system.compute_residual()

{physioblocks-1.0.1.dist-info → physioblocks-1.0.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: physioblocks
-Version: 1.0.1
+Version: 1.0.3
 Dynamic: Author
 Dynamic: Author-email
 Dynamic: License
@@ -124,4 +124,7 @@ python -m physioblocks.launcher references/spherical_heart_sim.jsonc -v -t -s Qu
 ```
 
 
-Results will be available in the `$LAUNCHER_FOLDER_PATH$/simulations/QuickStart` series folder.
+Results will be available in the `$LAUNCHER_FOLDER_PATH$/simulations/QuickStart` series folder:
+* the csv file contains the simulation results.
+* the html allows you to visualize the results.
+* the log and json files are here for debug purposes.

{physioblocks-1.0.1.dist-info → physioblocks-1.0.3.dist-info}/RECORD

@@ -1,4 +1,4 @@
-physioblocks/__init__.py,sha256=cSi4KUq8godvq235F5W2XvSwnWGwtIdTSE3fINK6vn8,1169
+physioblocks/__init__.py,sha256=cISNrgJhzM89OfbdFh9adOQcp0TvpiqizNJtm-sBWEE,1169
 physioblocks/base/__init__.py,sha256=2h6bk0fe8pwQvJ15Buytg4SZzYOH-Ep7WqoD6pd1eOQ,986
 physioblocks/base/operators.py,sha256=uUfswQZsaMNd2bd7wsaxSGFIYLAr1hQCQ6SQfYawdBc,5689
 physioblocks/base/registers.py,sha256=Cd0hmy1XqojjPQgb_jn4mGeAMTyWWgDkl3qJHCVjCRU,3119
@@ -9,7 +9,7 @@ physioblocks/computing/quantities.py,sha256=f2wU0VDylJ8nuUtOsV0FsA1mJnqNy1AvoK6-
 physioblocks/configuration/__init__.py,sha256=XNeIev97Fc9DAy_mA3CAZ0kMif7Nkop9sCI38ZZ2wgM,1454
 physioblocks/configuration/aliases.py,sha256=o7E54IYeW0NfmCX_r3zLKURVHsQoCmTlLD5085kY0eg,6302
 physioblocks/configuration/base.py,sha256=eLxZXZdc4QkSsA7TOA0-E5ckfMZZkdgPJ6JVrFwM4nY,3627
-physioblocks/configuration/constants.py,sha256=JsXGcoxHXgaixpIIiYZDPNiKSeD7Y_4YswVm0Uaf-Hs,3056
+physioblocks/configuration/constants.py,sha256=R72050sju7X_2kRtxY9mhyZ0btYLj8RaMeOsDoZT9DY,3046
 physioblocks/configuration/functions.py,sha256=pJl5BkSwRoIydH3bSMNilu0wj7yOVJ_BVAPNS-yMl8k,23520
 physioblocks/configuration/computing/__init__.py,sha256=duNJXdaHTZ9_V9XScNXjJgPsKSc6FbHGrwNDLrcYDg8,1012
 physioblocks/configuration/computing/quantities.py,sha256=fM5L4YeeM1NQQjhbdU4e50yxZ9B2ogF0WzDIYqFDbgg,1721
@@ -17,7 +17,7 @@ physioblocks/configuration/description/__init__.py,sha256=-KyjKPheXzgKiVZ0jL-nl8
 physioblocks/configuration/description/blocks.py,sha256=FGYFR5S15jGlZnIIfvZLdMr33bXabXMi9Vj9jETLeeM,8003
 physioblocks/configuration/description/nets.py,sha256=wKAH1Qq_tQ65-bn17qV6c_n0dt__PtY9r8TnYtZsPTY,5124
 physioblocks/configuration/simulation/__init__.py,sha256=RuRdBDvwqBcyghwh1D0ud3qwdUEcw2vgHPoB8drXp6w,1126
-physioblocks/configuration/simulation/simulations.py,sha256=4W-pcSDIMBKRFYjkkHOePrcZfXiFdfCHdyI_qpFSZ_c,8795
+physioblocks/configuration/simulation/simulations.py,sha256=bOf2a7rlYcw00XsoPspL2LBSTfAgl0HC7xCeMJiQ-Zs,8969
 physioblocks/description/__init__.py,sha256=o6G4LgLuSKSXOF6Uex-nua-8MyWsiwBd_G4dQLJsA6Y,1200
 physioblocks/description/blocks.py,sha256=y0qLvdE4VBmdDmOeHFaRFsmPRg-7_JVSGs2thAAUzHc,13542
 physioblocks/description/flux.py,sha256=aKs31bO1KyfAlHhmWUP9ZS7V7SryyVe_6MZG_Sp9jCQ,4858
@@ -44,7 +44,7 @@ physioblocks/library/aliases/model_components/spherical_dynamics.json,sha256=Tcd
 physioblocks/library/aliases/model_components/velocity_law_hht.json,sha256=HKbc5_2hGuOiQKUvLGto6tBQI_Z5AbbhXp8Itwmo2ag,93
 physioblocks/library/aliases/nets/circulation_alone_net.json,sha256=M3TiXIXgaXCsdwfyA5cm3fqspViiIl0Kyb6NKMiFsQ8,874
 physioblocks/library/aliases/nets/spherical_heart_net.json,sha256=UnOSQslMJuvbT6dXXPWj-HEvqe6U4fG1CK17_sKTxM8,3057
-physioblocks/library/aliases/simulations/circulation_alone_forward_simulation.jsonc,sha256=0l4dKEvmY6-3NLB70B4zVF_OJs-bzgxAsesgr3KZ7P0,1983
+physioblocks/library/aliases/simulations/circulation_alone_forward_simulation.jsonc,sha256=F5MIEC-OI-VUcAhYU8SLesbZPAiy3Us1f4VuDd0TEiI,1948
 physioblocks/library/aliases/simulations/default_forward_simulation.jsonc,sha256=TE_82KWUYAQdMX4MpDGY-Hi_djpwYAhSB5uzYtzYqcA,136
 physioblocks/library/aliases/simulations/default_time.jsonc,sha256=LM2_EU-GvOeGXJMvXXyJeWEOv-cfy25luhcaSMj5hkw,165
 physioblocks/library/aliases/simulations/newton_method_solver.json,sha256=p0jgVCpySTR0Npo5iFZrUUys3SRcy4j_EesQw9SvoyA,79
@@ -62,7 +62,7 @@ physioblocks/library/functions/trigonometric.py,sha256=Y5UolXZTASEi4EH_GpLpUIxMZ
 physioblocks/library/functions/watchers.py,sha256=fY1SSPn-7R9cix0VMFP69H5-pmJzl4OXs1tgvjCetvM,3212
 physioblocks/library/model_components/__init__.py,sha256=JcDL2PaPTrT9mTlDUbPjUsyFfsiGh5GgroJ7J3g2uwE,998
 physioblocks/library/model_components/active_law.py,sha256=97D3fC7Dmzw8jLqDALxUJgl48EZYOx39SYjZgGJx-Eo,11634
-physioblocks/library/model_components/dynamics.py,sha256=SZaOU2_EFMRlqvJvsW9NXLTviN8NBAhy_9z_0EErPb0,35282
+physioblocks/library/model_components/dynamics.py,sha256=T-ZmfhTgtr1Tj6434GDYU-5XSztHx-dd6sf228c4cyo,35198
 physioblocks/library/model_components/rheology.py,sha256=7x5RrB_0DXMWlm0mIrTxdLlsjXVa1WPuFYYifW4fKeg,4724
 physioblocks/library/model_components/velocity_law.py,sha256=UYr9MywlLkRw0LszK0u_AamwThjobDMV9jYbaSbdDkU,5047
 physioblocks/registers/__init__.py,sha256=trxT3MHDFX8cJ-IBX2cEIm8-uZBnfs2UkV8QDclq9so,1081
@@ -72,22 +72,22 @@ physioblocks/registers/type_register.py,sha256=JUlrPfP-X_TjIOHewSTkzf7Ck3L0ivVoH
 physioblocks/simulation/__init__.py,sha256=wqKVI3OLLVWuQoUcWwKSs5ibdkTP3Da_5DOe36v1up4,1504
 physioblocks/simulation/constants.py,sha256=w9CsB4V7YtvcdW4N1p7Iu6MeV4i-kazcEPAygE7mSx0,1075
 physioblocks/simulation/functions.py,sha256=c6HFwXK8IAzH7gSkz9aMiXLqnWXn_4_9WsX4jrGN9pY,2307
-physioblocks/simulation/runtime.py,sha256=xMtUL07SSQYpskK_EnNi-qz2D1d4wlDAESWHPqMUHn8,16762
+physioblocks/simulation/runtime.py,sha256=JDRBoN5lSDUmEKohXIWPadxjO6JcFHd7JYzUlkFrzxE,15498
 physioblocks/simulation/saved_quantities.py,sha256=fKWr1wScqxHPhuCE7aDNs8lpm9oWfyLN8NZDus2knPM,4111
 physioblocks/simulation/setup.py,sha256=UmIklaHmfkxi_BIdYKOhorLoC9QT9D4oEXwC-Xfd6v0,18109
-physioblocks/simulation/solvers.py,sha256=IrLHrsLJhM97SYzeygz2-UkkNbpp2hsixCuUi0Nis98,7353
-physioblocks/simulation/state.py,sha256=jE6fTxzsnGKch1H6mErunpi7lwTb_9fvs2Hx77LIaMA,10310
+physioblocks/simulation/solvers.py,sha256=xztL-QboEO2mgtws13QE0-wHJ22vjQQzAQahYfxXVBQ,7226
+physioblocks/simulation/state.py,sha256=YOWqJSSHgTzMZEmaR4SX4qsckFZ_UQVfJ-jQ5niAUPc,12628
 physioblocks/simulation/time_manager.py,sha256=AbP_Y8mcG7-6x-3rGY5CHFQwKMqksO4HvdCJ-mU03d4,9381
 physioblocks/utils/__init__.py,sha256=EMwPHVznBryA1uOEFfY_11L9Sd9bdCU_x_huNgIaZtQ,988
 physioblocks/utils/dynamic_import_utils.py,sha256=AvGFRYPQcI2JSVpScXCgnYJDgF2uh8QI7xBewoe_9Ks,4862
 physioblocks/utils/exceptions_utils.py,sha256=2w2q1XP-UDnxCXFu4XI2k1-jLUub24oYuwdUfd6tOkE,3421
-physioblocks/utils/gradient_test_utils.py,sha256=41FVSrjzOLWmxvaiZI258QbmnQOA7M2tBm-Vd2GsEuo,11208
+physioblocks/utils/gradient_test_utils.py,sha256=ISvkFk-HPbBafwmi9Zqimo6Dpb6yJpcLyyAoTge63RA,10666
 physioblocks/utils/math_utils.py,sha256=w5MDs4cP5WiTPHe8WofFuNZZ20WDNLT0XmFGlbisJJ0,2881
 physioblocks/references/circulation_alone_sim.jsonc,sha256=qv27nzmK7jennhtCJ7wQ_f83mt9_YJSOTABjGCYMcZw,736
 physioblocks/references/spherical_heart_respiration_sim.jsonc,sha256=kKq5xf94r90O1KpkCHk1qwgy6DJ37ePo3-8buEtcu5Y,1397
 physioblocks/references/spherical_heart_sim.jsonc,sha256=kGGBqxY8dnmGKmpMv216kbneyJRGz_3VkcYn1OWavIA,1151
-physioblocks-1.0.1.dist-info/METADATA,sha256=q9nDe7s9M-P1FdTf5yuGpLSlfaY1tU96Al-T3k5IMSM,4558
-physioblocks-1.0.1.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
-physioblocks-1.0.1.dist-info/licenses/licenses/GPL-3.0-only.txt,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-physioblocks-1.0.1.dist-info/licenses/licenses/LGPL-3.0-only.txt,sha256=46mU2C5kSwOnkqkw9XQAJlhBL2JAf1_uCD8lVcXyMRg,7652
-physioblocks-1.0.1.dist-info/RECORD,,
+physioblocks-1.0.3.dist-info/METADATA,sha256=-_Pm7SVFyAzY4eHI79u3CGR93C-Yo5iWGwUJzDtM6gM,4708
+physioblocks-1.0.3.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
+physioblocks-1.0.3.dist-info/licenses/licenses/GPL-3.0-only.txt,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+physioblocks-1.0.3.dist-info/licenses/licenses/LGPL-3.0-only.txt,sha256=46mU2C5kSwOnkqkw9XQAJlhBL2JAf1_uCD8lVcXyMRg,7652
+physioblocks-1.0.3.dist-info/RECORD,,