Summary

We are in need of faster solving and simultaneously better steady state checking if the algebra solver is used.
I’ve included additional unit tests so we can be sure that each of the configurations are solving to the correct steady state.

The basic premise is to solve the ODE solution for a period of time so that it is close to the steady state solution and will be numerically stable when solving. Additionally, we’ve included the adjoint solver which produces more efficient solving -- even when using the algebra solver. This indicates that the backwards solution to the ODE is still being calculated despite being not used in the log-probability calculation.

Reason

• faster ode solving.

Drawbacks

• more parameters and finer tuning to ensure stable sampling.

Checklist:

• Updated any relevant documentation
• Add an adr doc if appropriate
• Include links to any relevant issues in the description
• Unit tests passing
• Integration tests passing