Dear all, I have a parameter defined as
struct parameter{ double wave_num; bool partial; int count; }; I pass this parameter and the Jacobian to gsl_odeiv2_system and every thing works fine the first time around. But if I immediately call the solver again the integration fails. I am pretty sure it is a bug in step_bsimp because the error does not appear with other solvers. For example gsl_odeiv2_step_msbdf or for that matter even if I use gsl_odeiv2_step_rk8pd (which is not really recommended for stiff problems) every thing works just fine. So I am reporting that with a parameter defined like above and on the second and third calls to gsl_odeiv2_step_bsimp it fails. gsl_odeiv2_step_msbdf works just fine, which lead me to conclude that bsimp has a bug. Can someone please say something about which is the preffered stiff solver? bsimp or msbdf? _______________________________________________ Help-gsl mailing list Help-gsl@gnu.org https://lists.gnu.org/mailman/listinfo/help-gsl