#include <gsl_astep.h>
Default template arguments
param_t
- (no default)func_t
- ode_functvec_t
- ovector_basealloc_vec_t
- ovectoralloc_t
- ovector_alloc Definition at line 251 of file gsl_astep.h.
Public Member Functions | |
virtual int | astep_derivs (double &x, double &h, double xmax, size_t n, vec_t &y, vec_t &dydx, vec_t &yerr, param_t &pa, func_t &derivs) |
Make an adaptive integration step of the system derivs with derivatives. | |
virtual int | astep_full (double x, double &h, double xmax, double &x_out, size_t n, vec_t &y, vec_t &dydx, vec_t &yout2, vec_t &yerr, vec_t &dydx_out, param_t &pa, func_t &derivs) |
Make an adaptive integration step of the system derivs . | |
virtual int | astep (double &x, double &h, double xmax, size_t n, vec_t &y, vec_t &dydx_out, vec_t &yerr, param_t &pa, func_t &derivs) |
Make an adaptive integration step of the system derivs . | |
Data Fields | |
gsl_ode_control< vec_t > | con |
Control specification. | |
Protected Member Functions | |
int | evolve_apply (double t0, double &t, double &h, double t1, size_t nvar, vec_t &y, vec_t &dydx, vec_t &yout2, vec_t &yerr, vec_t &dydx_out2, param_t &pa, func_t &derivs) |
Apply the evolution for the next adaptive step. | |
Protected Attributes | |
alloc_t | ao |
Memory allocator for objects of type alloc_vec_t . | |
alloc_vec_t | yout |
Temporary storage for yout. | |
alloc_vec_t | dydx_int |
Internal storage for dydx. | |
double | last_step |
The size of the last step. | |
unsigned long int | count |
The number of steps. | |
unsigned long int | failed_steps |
The number of failed steps. | |
size_t | msize |
The size of the allocated vectors. |
virtual int astep | ( | double & | x, | |
double & | h, | |||
double | xmax, | |||
size_t | n, | |||
vec_t & | y, | |||
vec_t & | dydx_out, | |||
vec_t & | yerr, | |||
param_t & | pa, | |||
func_t & | derivs | |||
) | [inline, virtual] |
This attempts to take a step of size h
from the point x
of an n-dimensional
system derivs
starting with y
. On exit, x
and y
contain the new values at the end of the step, h
contains the size of the step, dydx_out
contains the derivative at the end of the step, and yerr
contains the estimated error at the end of the step.
If the system derivs
or the base stepper return a non-zero value, the adaptive step is aborted and y
is unmodified. The error handler is never called.
Implements adapt_step.
Definition at line 471 of file gsl_astep.h.
virtual int astep_derivs | ( | double & | x, | |
double & | h, | |||
double | xmax, | |||
size_t | n, | |||
vec_t & | y, | |||
vec_t & | dydx, | |||
vec_t & | yerr, | |||
param_t & | pa, | |||
func_t & | derivs | |||
) | [inline, virtual] |
This attempts to take a step of size h
from the point x
of an n-dimensional
system derivs
starting with y
and given the initial derivatives dydx
. On exit, x
, y
and dydx
contain the new values at the end of the step, h
contains the size of the step, dydx
contains the derivative at the end of the step, and yerr
contains the estimated error at the end of the step.
If the base stepper returns a non-zero value, the step is aborted and y
is unmodified. The error handler is never called.
Implements adapt_step.
Definition at line 381 of file gsl_astep.h.
virtual int astep_full | ( | double | x, | |
double & | h, | |||
double | xmax, | |||
double & | x_out, | |||
size_t | n, | |||
vec_t & | y, | |||
vec_t & | dydx, | |||
vec_t & | yout2, | |||
vec_t & | yerr, | |||
vec_t & | dydx_out, | |||
param_t & | pa, | |||
func_t & | derivs | |||
) | [inline, virtual] |
This function performs an adaptive integration step with the n-dimensional
system derivs
and parameter pa
. It Begins at x
with initial stepsize h
, ensuring that the step goes no farther than xmax
. At the end of the step, the size of the step taken is h
and the new value of x
is in x_out
. Initially, the function values and derivatives should be specified in y
and dydx
. The function values, derivatives, and the error at the end of the step are given in yout
, yerr
, and dydx_out
. Unlike in odestep
objects, the objects y
, yout
, dydx
, and dydx_out
must all be distinct.
If the base stepper returns a non-zero value, the step is aborted. The error handler is never called.
This adaptive stepper function is faster than astep() or astep_derivs() because it does not require any copying of vectors.
Implements adapt_step.
Definition at line 435 of file gsl_astep.h.
Documentation generated with Doxygen and provided under the GNU Free Documentation License. See License Information for details.
Project hosting provided by
,
O2scl Sourceforge Project Page