pnkraemer / diffeqzoo Goto Github PK

What do we do with one-dimensional problems? Implement them as scalars (0-dimensional arrays), 1-dimensional arrays, or both?

https://github.com/williamgilpin/dysts

http://www.3d-meier.de/tut19/Seite1.html

Must have:

• A reference to the original paper (or whatever one can find)
• A description of what it models, and why it is interesting

Nice to have:

• Meanings of parameter names
• Maths
• Example usage of the function

Each IVP constructor should behave similarly (with problem-specific differences), and this should be explained in the module docstring of each collection.

The notebooks are there and working, but

• use the old API (with namedtuple outputs)
• Have no text
• Are not rendered with plots in the docs

Would be great to give them some care

The following ODE problems would be great to have. (This issue will keep being edited.)

Initial value problems

• IVP: Brusselator:
  • #82
  • #83
• IVP: Discretised PDE
  • #84
  • #85
  • #86
• #87
• #88
• #89
• #90
• #91
• #92
• #93
• #94
• #95
• #96
• #97
• #98
• #99

https://docs.juliahub.com/GeometricProblems/4fHs6/0.1.1/

They also list a whole bunch of interesting problems.

For non-trivial boundary conditions, one cannot infer the dimension of the problem from the shape of the initial condition.

Mention the design choices somewhere (implementation as raw as possible even if it implies a slightly annoying API, implement in a way that makes it easy to use for benchmarks and papers, take second-order, non-autonomous, mass-matric, problems seriously, etc.)

• should it be a subpackage, with a more granular functionality? for BVPs, vector field conversion remains the same, but the function does not apply anymore.
• The jacobian conversion is incorrect. It does not convert the output appropriately.
• The functionality is not used anywhere (e.g., in the notebooks), so I have no idea how useful it is.
• Mark somewhere that it is very experimental (even more than the rest). Complete conversion is, in the end, something the user has to do.