Hello,
On 12/22/2011 12:07 AM, Farkas, Illes wrote:
5s ) later there is an update when gsl_odeiv2_driver_apply returns the
FAILURE value: -1. After logging the current time directly from the
"function" and "jacobian" (used by the gsl_odeiv2_system, which is
driver by the driver), I found that this particular update fails,
because the time step is halved again and again until it reaches the
limit of numerical precision.
Sounds like you're driving into a singularity.
Try different steppers (msbdf, bsimp). If they're all stopping about
there, the problem might be with your system. I guess you either have a
bug in rhs function or jacobian code, or your system really has a
singluarity in it. Or maybe you just have inappropriate tolerances in
relation to your system?
BR,
Tuomo
--
tuomo.keskit...@iki.fi
http://iki.fi/tuomo.keskitalo
_______________________________________________
Help-gsl mailing list
Help-gsl@gnu.org
https://lists.gnu.org/mailman/listinfo/help-gsl
