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 
44 class VTKCOMMONMATH_EXPORT vtkRungeKutta45 : public vtkInitialValueProblemSolver
45 {
46 public:
48  void PrintSelf(ostream& os, vtkIndent indent) override;
49 
53  static vtkRungeKutta45* New();
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 
108 protected:
109  vtkRungeKutta45();
110  ~vtkRungeKutta45() override;
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 
125 private:
126  vtkRungeKutta45(const vtkRungeKutta45&) = delete;
127  void operator=(const vtkRungeKutta45&) = delete;
128 };
129 
130 #endif
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...
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 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...
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
a simple class to control print indentation
Definition: vtkIndent.h:39
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...
Integrate an initial value problem using 5th order Runge-Kutta method with adaptive stepsize control...
static vtkObject * New()
Create an object with Debug turned off, modified time initialized to zero, and reference counting on...
Integrate a set of ordinary differential equations (initial value problem) in time.