autoreduce.utils.reduction.solve_sensitivity

solve_sensitivity(system_obj, timepoints, normalize=False, **kwargs)[source]

Solve the coupled state and local sensitivity equations.

Returns an array indexed by time point, parameter, and state. This direct path builds symbolic Jacobians and solves dS_j/dt = J(x, p) S_j + df/dp_j alongside the original ODE.