ComputationalScienceLaboratory / ODE-Test-Problems

A MATLAB suite of initial value problems
https://computationalsciencelaboratory.github.io/ODE-Test-Problems/
MIT License
9 stars 2 forks source link

Constrained pendulum problem #36

Open AndreyAPopov opened 2 years ago

AndreyAPopov commented 2 years ago

For use in half-explicit methods and for data assimilation.

Steven-Roberts commented 2 years ago

Long term, I would like to have a problem class for constrained mechanical systems (Solving ODEs II pg. 464), and the pendulum would simply be a preset. This would also unify the handling of reduced index forms.

AndreyAPopov commented 2 years ago

Long term, I would like to have a problem class for constrained mechanical systems (Solving ODEs II pg. 464), and the pendulum would simply be a preset. This would also unify the handling of reduced index forms.

I agree longterm. Right now we need this problem for something. We might also add a constrained shallow water (2D)

AndreyAPopov commented 2 years ago

Unless I am missing something, I think I have addressed all the issues.

Steven-Roberts commented 2 years ago

Unless I am missing something, I think I have addressed all the issues.

Added some final pedantic comments

AndreyAPopov commented 2 years ago

Alright, The last thing I want to add is the half-explicit jacobian, and Jvp, but I am not sure what to call it.

AndreyAPopov commented 2 years ago

All resolved.

AndreyAPopov commented 2 years ago

Should I remove the other DAE solvers (SDIRK and ESDIRK)?