Skip to content

Commit 35ee35c

Browse files
authored
First thoughts on realted software
1 parent a588020 commit 35ee35c

2 files changed

Lines changed: 25 additions & 0 deletions

File tree

paper/paper.bib

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -190,4 +190,20 @@ @article {OeffnerTorlo2020
190190
MRNUMBER = {4064785},
191191
DOI = {10.1016/j.apnum.2020.01.025},
192192
URL = {https://doi.org/10.1016/j.apnum.2020.01.025},
193+
}
194+
195+
@article {Shampine2005,
196+
AUTHOR = {Shampine, L. F. and Thompson, S. and Kierzenka, J. A. and
197+
Byrne, G. D.},
198+
TITLE = {Non-negative solutions of {ODE}s},
199+
JOURNAL = {Applied Mathematics and Computation},
200+
VOLUME = {170},
201+
YEAR = {2005},
202+
NUMBER = {1},
203+
PAGES = {556--569},
204+
ISSN = {0096-3003,1873-5649},
205+
MRCLASS = {34C11},
206+
MRNUMBER = {2177561},
207+
DOI = {10.1016/j.amc.2004.12.011},
208+
URL = {https://doi.org/10.1016/j.amc.2004.12.011},
193209
}

paper/paper.md

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -54,6 +54,15 @@ All implemented schemes were originally introduced for conservative production-d
5454

5555
Readers interested in additional theoretical background and further properties of the implemented schemes are referred to the following papers: @kopeczmeister2019, @izgin2022stability1, @izgin2022stability2, @huang2023, @torlo2022, @izginoeffner2023
5656

57+
Existing software libraries do not have a strong focus on unconditional positivity and, to the authors' knowledge, there is no other software library which offers MPRK schemes.
58+
59+
A common strategy to obtain nonnegative solutions used in the `PositiveDomain` callback of `Differentialequtions'jl` or the commercial package `Matlab` is described by @Shampine2005. In this approach negative components of approximate solutions that have been accepted by the adaptive time stepping algorithm are set to zero.
60+
61+
Another possibility is to reduce the chosen time step size beyond accuracy considerations until a non-negative approximation is calculated. This can be achieved in `DifferentialEquations.jl` using the solver option `isoutofdomain`.
62+
63+
Some papers on MPRK schemes offer supplementary codes. However, these are mainly small scripts for the reproduction of results shown in the papers and are not intended as software libraries.
64+
65+
5766
TODO
5867

5968
@bartel2024structure

0 commit comments

Comments
 (0)