DAE Solver Issues in Octave

ComputationalScienceLaboratory / ODE-Test-Problems

A MATLAB suite of initial value problems

https://computationalsciencelaboratory.github.io/ODE-Test-Problems/

MIT License

9 stars 3 forks source link

DAE Solver Issues in Octave #42

Closed Steven-Roberts closed 1 year ago

Steven-Roberts commented 1 year ago

For the Ascher linear DAE, the default Octave solver ode23s is giving incorrect answers. Presumably time dependent mass matrix is not well supported. ode15s works but breaks for other problems due to Vectorized property. Need to reevaluate default stiff/DAE solver for this problem and others.

AndreyAPopov commented 1 year ago

Could we solve this by using our own Lobatto IIIC?