gsl_odeiv_control *
gsl_odeiv_control_scaled_new(double eps_abs, double eps_rel,
double a_y, double a_dydt,
const double scale_abs[],
size_t dim)
{
gsl_odeiv_control * c =
gsl_odeiv_control_alloc (gsl_odeiv_control_scaled);
int status = gsl_odeiv_control_init (c, eps_abs, eps_rel, a_y, a_dydt);
if (status != GSL_SUCCESS)
{
gsl_odeiv_control_free (c);
GSL_ERROR_NULL ("error trying to initialize control", status);
}
{
sc_control_state_t * s = (sc_control_state_t *) c->state;
s->scale_abs = (double *)malloc(dim * sizeof(double));
if (s->scale_abs == 0)
{
free (s);
GSL_ERROR_NULL ("failed to allocate space for scale_abs",
GSL_ENOMEM);
}
memcpy(s->scale_abs, scale_abs, dim * sizeof(double));
}
return c;
}
gsl_odeiv_control *
gsl_odeiv_control_standard_new(double eps_abs, double eps_rel,
double a_y, double a_dydt)
{
gsl_odeiv_control * c =
gsl_odeiv_control_alloc (gsl_odeiv_control_standard);
int status = gsl_odeiv_control_init (c, eps_abs, eps_rel, a_y, a_dydt);
if (status != GSL_SUCCESS)
{
gsl_odeiv_control_free (c);
GSL_ERROR_NULL ("error trying to initialize control", status);
}
return c;
}