yousef/Applying DIRKs to DAE and developing error estimators

Mar 20th

  • Successfully implemented code that approximates solutions to the Van de Pol DAE with DIRK methods using Newton’s method.

  • Problem: unlike low order stiffly accurate DIRK where the error goes to 1e-16, the error for the high order DIRK doesn’t go below 1e-9.

Mar 29 th

  • VDP problem:

  • simple DAE against exact solution with (T = 1)

  • Robert’s chemical Run DAE

  • Pendulum