Abstract
AbstractWe apply Patankar Runge–Kutta methods to y′ = M(y)y and focus on the case where M(y) is a graph Laplacian as the resulting scheme will preserve positivity and total mass. The second order Patankar Heun method is tested using four test problems (stiff and non‐stiff) cast into this form. The local error is estimated and the step size is chosen adaptively. Concerning accuracy and efficiency, the results are comparable to those obtained with a traditional L‐stable, second order Rosenbrock method.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have