[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] Checking error in ODE system

From: Pau Cervera Badia
Subject: Re: [Help-gsl] Checking error in ODE system
Date: Wed, 16 Aug 2006 18:11:51 +0200
User-agent: Thunderbird (X11/20060728)


I think the error is located at the yerr field of the gsl_odeiv_evolve pointer type that has the name evolve in your piece of code. I think you can access the value via,

double error_at_i = evolve->yerr[i];

You can see the definition of the evolution object at gsl_odeiv.h. You can see other uses of yerr in the source for the stepper methods where the errors ara estimated; see rk2.c in rk2_apply, for instance.

If someone knows a better way, I would like to know it. Thanks!

I hope that will help you Sebastien,

Sebastien Maret wrote:

I am trying to solve a system of stiff ODEs using the GSL. The part of
my code that solve the equations looks like:

     for (i = 0; i < time_steps; i++)
         while (t < tim[i])
             gsl_odeiv_evolve_apply (evolve, control, step, &syst,
                                     &t, tim[i], &h, y);

I would like to check what is the error for y[] reached at the end of each
gsl_odeiv_evolve_apply() call.

How can I do that? An example would be appreciated.


Help-gsl mailing list

Pau Cervera i Badia (e-mail address@hidden)
  Departament de Física Fonamental               Martí i Franqués, 1
  Universitat de Barcelona                 Planta 3, despatx 346 bis
                                                     08028 Barcelona
  tel: +34 934 921 155                                         Spain

  "Be careful of reading health books, you might die of a misprint."
                                                       -- Mark Twain

reply via email to

[Prev in Thread] Current Thread [Next in Thread]