Details of daessc solver

Matlab has a lot of ODE solvers available and each of them is properly documented.

However, when it comes to the "daessc" solver for DAE systems, I have not been able to find the mathematical details.

I am mostly interested in what happens after the index-reduction process has taken place, and the solver needs to integrate a 1-index DAE: how does this integration happen? Does the solver use an implicit method? Is it single-step or multistep?

Any link to some proper documentation would be appreciated

[轩]

I am also looking forward to the details.

Answer 1

[Arun]

Hey Federico,

I understand that you seek information about the DAESSC solver implementation in MATLAB.

It is multistep implicit method of low order, based on Backwards Euler and NDF2 formula. The solver algorithms were “tuned” towards robust solving of DAE’s with dominant algebraic part. Mathematical details of the NDF formulas can be found in “The MATLAB ODE Suite”. SIAM Journal of Scientific Computing, 1997, vol 18, No. 1.

Here is the link to the paper with mathematical details: https://doi.org/10.1137/S1064827594276424

Hope this helps.