#include <gsl_rk8pd.h>
This a fast version of gsl_rk8pd, which is a stepper for a fixed number of ODEs. It ignores the error values returned by the derivs
argument. The argument n
to step() should always be equal to the template parameter N
, and the vector parameters to step must have space allocated for at least N
elements. No error checking is performed to ensure that this is the case.
Definition at line 399 of file gsl_rk8pd.h.
Public Member Functions | |
virtual int | step (double x, double h, size_t n, vec_t &y, vec_t &dydx, vec_t &yout, vec_t &yerr, vec_t &dydx_out, param_t &pa, func_t &derivs) |
Perform an integration step. | |
Protected Attributes | |
alloc_t | ao |
Memory allocator for objects of type alloc_vec_t . | |
Storage for the intermediate steps | |
alloc_vec_t | k2 |
alloc_vec_t | k3 |
alloc_vec_t | k4 |
alloc_vec_t | k5 |
alloc_vec_t | k6 |
alloc_vec_t | k7 |
alloc_vec_t | ytmp |
alloc_vec_t | k8 |
alloc_vec_t | k9 |
alloc_vec_t | k10 |
alloc_vec_t | k11 |
alloc_vec_t | k12 |
alloc_vec_t | k13 |
Storage for the coefficients | |
double | Abar [13] |
double | A [12] |
double | ah [10] |
double | b21 |
double | b3 [2] |
double | b4 [3] |
double | b5 [4] |
double | b6 [5] |
double | b7 [6] |
double | b8 [7] |
double | b9 [8] |
double | b10 [9] |
double | b11 [10] |
double | b12 [11] |
double | b13 [12] |
virtual int step | ( | double | x, | |
double | h, | |||
size_t | n, | |||
vec_t & | y, | |||
vec_t & | dydx, | |||
vec_t & | yout, | |||
vec_t & | yerr, | |||
vec_t & | dydx_out, | |||
param_t & | pa, | |||
func_t & | derivs | |||
) | [inline, virtual] |
Perform an integration step.
Given initial value of the n-dimensional function in y
and the derivative in dydx
(which must generally be computed beforehand) at the point x
, take a step of size h
giving the result in yout
, the uncertainty in yerr
, and the new derivative in dydx_out
using function derivs
to calculate derivatives. The parameters yout
and y
and the parameters dydx_out
and dydx
may refer to the same object.
n
should be equal to the template parameter N
. Implements odestep.
Definition at line 601 of file gsl_rk8pd.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