[R-sig-dyn-mod] Possible bug specifying root functions for lsodes

[McGrath, Justin M]

Thanks Thomas, that table is helpful. I wasn't aware that you could pass options to the solvers to switch between stiff and nonstiff methods.

Justin
________________________________________
From: r-sig-dynamic-models-bounces at r-project.org [r-sig-dynamic-models-bounces at r-project.org] on behalf of Thomas Petzoldt [Thomas.Petzoldt at TU-Dresden.de]
Sent: Tuesday, March 05, 2013 5:26 PM
To: Special Interest Group for Dynamic Simulation Models in R
Subject: Re: [R-sig-dyn-mod] Possible bug specifying root functions     for     lsodes

Justin,

most of the solvers, except rk and euler, are able to handle stiff
systems, see the following table:

http://desolve.r-forge.r-project.org/slides/tutorial.pdf#24

Sometimes, tuning a few options may be necessary, e.g. specific settings
of atol and rtol for individual state variables, or choosing specific
method flags for lsode.

The use of lsodes with root finding is somewhat exotic and I guess that
nobody tried this before. Anyway, we will work on it, thanks for your
report.

Thomas

_______________________________________________
R-sig-dynamic-models mailing list
R-sig-dynamic-models at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-dynamic-models