gsl_rkck_fast Class Template Reference

#include <gsl_rkck.h>

Inheritance diagram for gsl_rkck_fast:

odestep

Detailed Description

template<size_t N, class param_t, class func_t, class vec_t = ovector_view, class alloc_vec_t = ovector, class alloc_t = ovector_alloc>
class gsl_rkck_fast< N, param_t, func_t, vec_t, alloc_vec_t, alloc_t >

Faster Cash-Karp embedded Runge-Kutta ODE stepper (GSL).

This a faster version of gsl_rkck, 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 215 of file gsl_rkck.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 ytmp
Storage for the coefficients
double ah [5]
double b3 [2]
double b4 [3]
double b5 [4]
double b6 [5]
double ec [7]
double b21
double c1
double c3
double c4
double c6

Member Function Documentation

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.

Note:
The value of the parameter n should be equal to the template parameter N.

Reimplemented from odestep.

Definition at line 307 of file gsl_rkck.h.


The documentation for this class was generated from the following file:

Documentation generated with Doxygen and provided under the GNU Free Documentation License. See License Information for details.

Project hosting provided by SourceForge.net Logo, O2scl Sourceforge Project Page