-
-
Notifications
You must be signed in to change notification settings - Fork 58
Small Grant Request+Completion: DAE Problem Benchmarks #168
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Conversation
|
|
||
| ## DAE Problem Benchmarks (\$100 / Benchmark) | ||
|
|
||
| Completed by **Jayant Pranjal** |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This one isn't merged yet.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Merged now 🎉 so this is done
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks😊
|
What is the new request here? |
|
I was thinking of writing another dae benchmark(Wheelset problem). |
|
Okay sounds good to me. Note that because it's an Index-2 problem, you'll likely need to use ModelingToolkit to derive the Index-1 problem for the mass matrix and DAE problem form, and then also have a separate implementation of the index-2 problem form in direct code, and then work out which solvers end up being able work if the correct index-2 initialization is given. So that will have a few more complications from before, but it's not really well-characterized how the methods will do here (Rodas5P and radau should be able to handle index-2 IIRC, https://www.sciencedirect.com/science/article/abs/pii/S0168927496000232), so it'll be very useful to have a benchmark which finally stresses this behavior and gives a good categorization of index-2 handling, along with benchmarks to see if there is performance lost by ModelingToolkit's reduction to index-1. Definitely approved here under the understanding the benchmark will include this as a required element in the benchmarking. |
|
Requires two other approvals. |
Reference PR: SciML/SciMLBenchmarks.jl#1303