Circulatory System Models

CirculatorySystemModels v0.3.1

Diff since v0.3.0

Breaking changes since version 0.1.X:

1. The Double-Hill chamber DHChamber is now in terms of dV/dt rather than dp/dt which increases computation speed by an order of magnitude. The old form can be recovered using the inP argument. The parameter Ev has been removed. Please use the ventricular volume as an initial condition instead of the pressure. This can be calculated as $V = p/E_{min}$
2. Compliance and Elastance are now in terms of dV/dt. Use the new argument inP=true to recover the old form. Initial conditions need to be changed accordingly.
3. Compliance and Elastance have options for external pressure as parameter has_ep and as a connector pin to connect other pressures to has_variable_ep - the latter can take an equation like C.ep.p ~ LV.p * K.
4. Chamber has been deleted and replaced with VariableCompliance which can take a function for E(t) and has external pressure options like the other compliances and elastances.
5. ShiChamber already was formulated in terms of dV/dt. It now has the argument inP=true to allow the dp/dt form. Initial conditions need to be changed accordingly.

*NOTE: the use of inP is strongly discouraged, since it needs the first derivative of the Elastance function, is less stable, and an order of magnitude slower. If you want to define initial condtions in pressure, it is advised to calculate the corresponding initial volume.*