mountaindust / Planktos

ABM framework for dispersal modeling
https://planktos.readthedocs.io
GNU General Public License v3.0
5 stars 3 forks source link

Upgrade SDE Euler solver to Milstein #33

Open mountaindust opened 3 years ago

mountaindust commented 3 years ago

See note in the docstring, and Milstein Method on Wikipedia. For now, we will assume that \sigma'(x) is zero for most applications, making this method equivalent to the current one.

mountaindust commented 3 years ago

Needed for porous media models. It's worth noting in the documentation somewhere that D for porous media should be the effective diffusion constant divided by volumetric water content. See this paper.