VTK  9.2.6
vtkRungeKutta45.h
Go to the documentation of this file.
1/*=========================================================================
2
3 Program: Visualization Toolkit
4 Module: vtkRungeKutta45.h
5
6 Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
7 All rights reserved.
8 See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
9
10 This software is distributed WITHOUT ANY WARRANTY; without even
11 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
12 PURPOSE. See the above copyright notice for more information.
13
14=========================================================================*/
38#ifndef vtkRungeKutta45_h
39#define vtkRungeKutta45_h
40
41#include "vtkCommonMathModule.h" // For export macro
43
44class VTKCOMMONMATH_EXPORT vtkRungeKutta45 : public vtkInitialValueProblemSolver
45{
46public:
48 void PrintSelf(ostream& os, vtkIndent indent) override;
49
54
55 using Superclass::ComputeNextStep;
57
79 int ComputeNextStep(double* xprev, double* xnext, double t, double& delT, double maxError,
80 double& error, void* userData) override
81 {
82 double minStep = delT;
83 double maxStep = delT;
84 double delTActual;
85 return this->ComputeNextStep(
86 xprev, nullptr, xnext, t, delT, delTActual, minStep, maxStep, maxError, error, userData);
87 }
88 int ComputeNextStep(double* xprev, double* dxprev, double* xnext, double t, double& delT,
89 double maxError, double& error, void* userData) override
90 {
91 double minStep = delT;
92 double maxStep = delT;
93 double delTActual;
94 return this->ComputeNextStep(
95 xprev, dxprev, xnext, t, delT, delTActual, minStep, maxStep, maxError, error, userData);
96 }
97 int ComputeNextStep(double* xprev, double* xnext, double t, double& delT, double& delTActual,
98 double minStep, double maxStep, double maxError, double& error, void* userData) override
99 {
100 return this->ComputeNextStep(
101 xprev, nullptr, xnext, t, delT, delTActual, minStep, maxStep, maxError, error, userData);
102 }
103 int ComputeNextStep(double* xprev, double* dxprev, double* xnext, double t, double& delT,
104 double& delTActual, double minStep, double maxStep, double maxError, double& estErr,
105 void* userData) override;
107
108protected:
111
112 void Initialize() override;
113
114 // Cash-Karp parameters
115 static double A[5];
116 static double B[5][5];
117 static double C[6];
118 static double DC[6];
119
120 double* NextDerivs[6];
121
122 int ComputeAStep(double* xprev, double* dxprev, double* xnext, double t, double& delT,
123 double& delTActual, double& error, void* userData);
124
125private:
126 vtkRungeKutta45(const vtkRungeKutta45&) = delete;
127 void operator=(const vtkRungeKutta45&) = delete;
128};
129
130#endif
a simple class to control print indentation
Definition vtkIndent.h:40
Integrate a set of ordinary differential equations (initial value problem) in time.
virtual int ComputeNextStep(double *xprev, double *xnext, double t, double &delT, double maxError, double &error)
Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of...
Integrate an initial value problem using 5th order Runge-Kutta method with adaptive stepsize control.
~vtkRungeKutta45() override
int ComputeNextStep(double *xprev, double *dxprev, double *xnext, double t, double &delT, double maxError, double &error, void *userData) override
Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of...
int ComputeNextStep(double *xprev, double *xnext, double t, double &delT, double maxError, double &error, void *userData) override
Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of...
int ComputeAStep(double *xprev, double *dxprev, double *xnext, double t, double &delT, double &delTActual, double &error, void *userData)
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
static vtkRungeKutta45 * New()
Construct a vtkRungeKutta45 with no initial FunctionSet.
int ComputeNextStep(double *xprev, double *xnext, double t, double &delT, double &delTActual, double minStep, double maxStep, double maxError, double &error, void *userData) override
Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of...
void Initialize() override
int ComputeNextStep(double *xprev, double *dxprev, double *xnext, double t, double &delT, double &delTActual, double minStep, double maxStep, double maxError, double &estErr, void *userData) override
Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of...