Message-ID: <3a68c50c-3880-b2bd-99b1-e1a3b8707090@mnhn.fr>
Date: 2019-03-26T16:54:01Z
From: Alain Danet
Subject: [R-sig-dyn-mod] simecol - communication between solver and equation slot
In-Reply-To: <cc491068-8a2f-3293-8ce1-30eccbbc3151@mnhn.fr>
Dear all,
I simulate ecological models (as ODE systems) using the simecol
<https://cran.r-project.org/web/packages/simecol/index.html>
package framework which makes very easy the use and the sharing of
ecological models through the object class SimObj (See here
<https://cran.r-project.org/web/packages/simecol/vignettes/a-simecol-introduction.pdf>).
I want to implement a steady-state which stops the simulation once
the derivatives become very low.
According to this vignette
<https://cran.r-project.org/web/packages/rootSolve/vignettes/rootSolve.pdf>
and this example
<https://unamatematicaseltigre.blogspot.com/2016/08/how-to-simulate-bouncing-ball-in-r-more.html>,
you can implement it easily.
You have just to supply a custom solver which check the values of
the derivatives.
The problem is that the custom solver does not look to be able
reach the |equations| slot of the SimObj.
I would like to keep this nice functionality of the equation slots
to switch easily between different types of functional response.
I have developed a reproducible example here:
https://stackoverflow.com/questions/49819582/r-simecol-and-steady-state-rootfun-cannot-access-to-the-simobj-equations-slot
Thank you for your help.
Best regards,
Alain.
--
Alain Danet - PhD
UMR CESCO:http://cesco.mnhn.fr/en Mus?um National d'Histoire Naturelle - Paris
https://alaindanet.github.io/
/!\ New emailaddress:alain.danet at mnhn.fr /!\