Add new stiff ODE benchmark problems #1282
Merged
+589
−514
Conversation
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Implement benchmark files for Ring Modulator, E5, and Beam problems from the classic IVP test set: - RingModulator.jmd: 15-dimensional stiff system modeling electronic ring modulator circuit - E5.jmd: 4-dimensional chemical pyrolysis model with widely separated time scales - Beam.jmd: 80-dimensional elastic beam dynamics with structural stiffness These benchmarks complement the existing stiff ODE test suite and provide additional challenging test cases for evaluating solver performance across different system sizes and stiffness characteristics. Addresses issue #358 for adding PDFs from stiff ODE benchmark sets. 🤖 Generated with [Claude Code](https://claude.ai/code) Co-Authored-By: Claude <[email protected]>
Complete the missing benchmark problems from issue #358: - MedicalAkzoNobel.jmd: 6-dimensional singularly perturbed system modeling radioimmunotherapy, representing the penetration of radio-labeled antibodies into tumor tissue with small parameter ε causing stiffness - EMEP.jmd: 30-dimensional European atmospheric chemistry model focusing on ozone formation, including NOx-Ox chemistry, hydroxyl radical processes, organic compound reactions, and photochemistry with diurnal cycles These benchmarks complement the existing Ring Modulator, E5, and Beam problems to provide a complete implementation of the classic IVP test set problems mentioned in issue #358. The EMEP benchmark is distinct from the existing Pollution.jmd, representing a larger and more complex atmospheric chemistry system. 🤖 Generated with [Claude Code](https://claude.ai/code) Co-Authored-By: Claude <[email protected]>
Member
Author
|
One week of Claude and it was 3/5 after about a few thousand dollars of API calls for what would have taken about 3 hours to translate by hand, on equations where an original source was given from Fortran. I think humans still have a place in this world. |
Member
Author
|
Okay it got about 50% right on E5. Let's now fix the others. |
Member
Author
|
Okay, it actually got those other 2 wrong. It also got e5 kind of wrong? I had to fix things up by hand to get something running, but it did get the right equations (4 equations). So AI ended up 1 for 5 after trying to get it to correct itself. We need at least 100x larger models to get close to our day to day work 😅 |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Implement benchmark files for Ring Modulator, E5, and Beam problems from the classic IVP test set:
These benchmarks complement the existing stiff ODE test suite and provide additional challenging test cases for evaluating solver performance across different system sizes and stiffness characteristics.
Addresses issue #358 for adding PDFs from stiff ODE benchmark sets.
Complete the missing benchmark problems from issue #358:
MedicalAkzoNobel.jmd: 6-dimensional singularly perturbed system modeling radioimmunotherapy, representing the penetration of radio-labeled antibodies into tumor tissue with small parameter ε causing stiffness
EMEP.jmd: 30-dimensional European atmospheric chemistry model focusing on ozone formation, including NOx-Ox chemistry, hydroxyl radical processes, organic compound reactions, and photochemistry with diurnal cycles
These benchmarks complement the existing Ring Modulator, E5, and Beam problems to provide a complete implementation of the classic IVP test set problems mentioned in issue #358. The EMEP benchmark is distinct from the existing Pollution.jmd, representing a larger and more complex atmospheric chemistry system.