VTK
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Static Protected Attributes | List of all members
vtkRungeKutta45 Class Reference

Integrate an initial value problem using 5th order Runge-Kutta method with adaptive stepsize control. More...

#include <vtkRungeKutta45.h>

Inheritance diagram for vtkRungeKutta45:
[legend]
Collaboration diagram for vtkRungeKutta45:
[legend]

Public Types

typedef
vtkInitialValueProblemSolver 
Superclass
 
- Public Types inherited from vtkInitialValueProblemSolver
enum  ErrorCodes { OUT_OF_DOMAIN = 1, NOT_INITIALIZED = 2, UNEXPECTED_VALUE = 3 }
 
typedef vtkObject Superclass
 

Public Member Functions

virtual int IsA (const char *type)
 
vtkRungeKutta45NewInstance () const
 
void PrintSelf (ostream &os, vtkIndent indent)
 
virtual int ComputeNextStep (double *xprev, double *xnext, double t, double &delT, double maxError, double &error)
 
virtual int ComputeNextStep (double *xprev, double *dxprev, double *xnext, double t, double &delT, double maxError, double &error)
 
virtual int ComputeNextStep (double *xprev, double *xnext, double t, double &delT, double &delTActual, double minStep, double maxStep, double maxError, double &error)
 
virtual int ComputeNextStep (double *xprev, double *dxprev, double *xnext, double t, double &delT, double &delTActual, double minStep, double maxStep, double maxError, double &error)
 
- Public Member Functions inherited from vtkInitialValueProblemSolver
vtkInitialValueProblemSolverNewInstance () const
 
virtual int IsAdaptive ()
 
virtual void SetFunctionSet (vtkFunctionSet *functionset)
 
virtual vtkFunctionSetGetFunctionSet ()
 

Static Public Member Functions

static int IsTypeOf (const char *type)
 
static vtkRungeKutta45SafeDownCast (vtkObjectBase *o)
 
static vtkRungeKutta45New ()
 
- Static Public Member Functions inherited from vtkInitialValueProblemSolver
static int IsTypeOf (const char *type)
 
static
vtkInitialValueProblemSolver
SafeDownCast (vtkObjectBase *o)
 

Protected Member Functions

virtual vtkObjectBase * NewInstanceInternal () const
 
 vtkRungeKutta45 ()
 
 ~vtkRungeKutta45 ()
 
virtual void Initialize ()
 
int ComputeAStep (double *xprev, double *dxprev, double *xnext, double t, double &delT, double &error)
 
- Protected Member Functions inherited from vtkInitialValueProblemSolver
 vtkInitialValueProblemSolver ()
 
 ~vtkInitialValueProblemSolver ()
 

Protected Attributes

double * NextDerivs [6]
 
- Protected Attributes inherited from vtkInitialValueProblemSolver
vtkFunctionSetFunctionSet
 
double * Vals
 
double * Derivs
 
int Initialized
 
int Adaptive
 

Static Protected Attributes

static double A [5]
 
static double B [5][5]
 
static double C [6]
 
static double DC [6]
 

Detailed Description

Integrate an initial value problem using 5th order Runge-Kutta method with adaptive stepsize control.

This is a concrete sub-class of vtkInitialValueProblemSolver. It uses a 5th order Runge-Kutta method with stepsize control to obtain the values of a set of functions at the next time step. The stepsize is adjusted by calculating an estimated error using an embedded 4th order Runge-Kutta formula: Press, W. H. et al., 1992, Numerical Recipes in Fortran, Second Edition, Cambridge University Press Cash, J.R. and Karp, A.H. 1990, ACM Transactions on Mathematical Software, vol 16, pp 201-222

See Also
vtkInitialValueProblemSolver vtkRungeKutta4 vtkRungeKutta2 vtkFunctionSet
Tests:
vtkRungeKutta45 (Tests)

Definition at line 43 of file vtkRungeKutta45.h.

Member Typedef Documentation

Definition at line 46 of file vtkRungeKutta45.h.

Constructor & Destructor Documentation

vtkRungeKutta45::vtkRungeKutta45 ( )
protected
vtkRungeKutta45::~vtkRungeKutta45 ( )
protected

Member Function Documentation

static int vtkRungeKutta45::IsTypeOf ( const char *  type)
static
virtual int vtkRungeKutta45::IsA ( const char *  type)
virtual

Reimplemented from vtkInitialValueProblemSolver.

static vtkRungeKutta45* vtkRungeKutta45::SafeDownCast ( vtkObjectBase *  o)
static
virtual vtkObjectBase* vtkRungeKutta45::NewInstanceInternal ( ) const
protectedvirtual

Reimplemented from vtkInitialValueProblemSolver.

vtkRungeKutta45* vtkRungeKutta45::NewInstance ( ) const
void vtkRungeKutta45::PrintSelf ( ostream &  os,
vtkIndent  indent 
)
virtual

Reimplemented from vtkInitialValueProblemSolver.

static vtkRungeKutta45* vtkRungeKutta45::New ( )
static

Construct a vtkRungeKutta45 with no initial FunctionSet.

virtual int vtkRungeKutta45::ComputeNextStep ( double *  xprev,
double *  xnext,
double  t,
double &  delT,
double  maxError,
double &  error 
)
inlinevirtual

Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of x at t+delTActual (xnext). Possibly delTActual != delT. This may occur because this solver supports adaptive stepsize control. It tries to change to stepsize such that the (estimated) error of the integration is less than maxError. The solver will not set the stepsize smaller than minStep or larger than maxStep (note that maxStep and minStep should both be positive, whereas delT can be negative). Also note that delT is an in/out argument. vtkRungeKutta45 will modify delT to reflect the best (estimated) size for the next integration step. An estimated value for the error is returned (by reference) in error. This is the norm of the error vector if there are more than one function to be integrated. This method returns an error code representing the nature of the failure: OutOfDomain = 1, NotInitialized = 2, UnexpectedValue = 3

Reimplemented from vtkInitialValueProblemSolver.

Definition at line 68 of file vtkRungeKutta45.h.

virtual int vtkRungeKutta45::ComputeNextStep ( double *  xprev,
double *  dxprev,
double *  xnext,
double  t,
double &  delT,
double  maxError,
double &  error 
)
inlinevirtual

Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of x at t+delTActual (xnext). Possibly delTActual != delT. This may occur because this solver supports adaptive stepsize control. It tries to change to stepsize such that the (estimated) error of the integration is less than maxError. The solver will not set the stepsize smaller than minStep or larger than maxStep (note that maxStep and minStep should both be positive, whereas delT can be negative). Also note that delT is an in/out argument. vtkRungeKutta45 will modify delT to reflect the best (estimated) size for the next integration step. An estimated value for the error is returned (by reference) in error. This is the norm of the error vector if there are more than one function to be integrated. This method returns an error code representing the nature of the failure: OutOfDomain = 1, NotInitialized = 2, UnexpectedValue = 3

Reimplemented from vtkInitialValueProblemSolver.

Definition at line 77 of file vtkRungeKutta45.h.

virtual int vtkRungeKutta45::ComputeNextStep ( double *  xprev,
double *  xnext,
double  t,
double &  delT,
double &  delTActual,
double  minStep,
double  maxStep,
double  maxError,
double &  error 
)
inlinevirtual

Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of x at t+delTActual (xnext). Possibly delTActual != delT. This may occur because this solver supports adaptive stepsize control. It tries to change to stepsize such that the (estimated) error of the integration is less than maxError. The solver will not set the stepsize smaller than minStep or larger than maxStep (note that maxStep and minStep should both be positive, whereas delT can be negative). Also note that delT is an in/out argument. vtkRungeKutta45 will modify delT to reflect the best (estimated) size for the next integration step. An estimated value for the error is returned (by reference) in error. This is the norm of the error vector if there are more than one function to be integrated. This method returns an error code representing the nature of the failure: OutOfDomain = 1, NotInitialized = 2, UnexpectedValue = 3

Reimplemented from vtkInitialValueProblemSolver.

Definition at line 87 of file vtkRungeKutta45.h.

virtual int vtkRungeKutta45::ComputeNextStep ( double *  xprev,
double *  dxprev,
double *  xnext,
double  t,
double &  delT,
double &  delTActual,
double  minStep,
double  maxStep,
double  maxError,
double &  error 
)
virtual

Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of x at t+delTActual (xnext). Possibly delTActual != delT. This may occur because this solver supports adaptive stepsize control. It tries to change to stepsize such that the (estimated) error of the integration is less than maxError. The solver will not set the stepsize smaller than minStep or larger than maxStep (note that maxStep and minStep should both be positive, whereas delT can be negative). Also note that delT is an in/out argument. vtkRungeKutta45 will modify delT to reflect the best (estimated) size for the next integration step. An estimated value for the error is returned (by reference) in error. This is the norm of the error vector if there are more than one function to be integrated. This method returns an error code representing the nature of the failure: OutOfDomain = 1, NotInitialized = 2, UnexpectedValue = 3

Implements vtkInitialValueProblemSolver.

virtual void vtkRungeKutta45::Initialize ( )
protectedvirtual

Reimplemented from vtkInitialValueProblemSolver.

int vtkRungeKutta45::ComputeAStep ( double *  xprev,
double *  dxprev,
double *  xnext,
double  t,
double &  delT,
double &  error 
)
protected

Member Data Documentation

double vtkRungeKutta45::A[5]
staticprotected

Definition at line 108 of file vtkRungeKutta45.h.

double vtkRungeKutta45::B[5][5]
staticprotected

Definition at line 109 of file vtkRungeKutta45.h.

double vtkRungeKutta45::C[6]
staticprotected

Definition at line 110 of file vtkRungeKutta45.h.

double vtkRungeKutta45::DC[6]
staticprotected

Definition at line 111 of file vtkRungeKutta45.h.

double* vtkRungeKutta45::NextDerivs[6]
protected

Definition at line 113 of file vtkRungeKutta45.h.


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