PyBaMM
PyBaMM copied to clipboard
Add IDA solver for initial conditions root-finding
Description
Currently the IDA solver uses the Casadi solver for initial conditions root finding. The Casadi solver fails for any index 2 DAE systems. However, despite IDA capable of solving index 2 DAE systems, if initial root finding is disabled the solution fails (likely due to poorly defined initial conditions).
Motivation
No response
Possible Implementation
No response
Additional context
No response
Hi @Dharshannan, please could you provide some example code which reproduces this issue?