Irreducible causes issues with some ODE/SDE solvers thus keep ODEs separate from algebraic equations, ex. Balloon model

[david-hofmann]

For instance it is not possible to use Euler-Maruyama solver: “This solver is not able to use mass matrices.”
Unless there is another solution, the bold signal needs to be separated from the hemodynamic response function in:

Neuroblox.jl/src/measurementmodels/fmri.jl

Line 46 in 0918db7

struct BalloonModel <: ObserverBlox