VTK
vtkRungeKutta4.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkRungeKutta4.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 =========================================================================*/
29 #ifndef vtkRungeKutta4_h
30 #define vtkRungeKutta4_h
31 
32 #include "vtkCommonMathModule.h" // For export macro
34 
35 class VTKCOMMONMATH_EXPORT vtkRungeKutta4 : public vtkInitialValueProblemSolver
36 {
37 public:
39  void PrintSelf(ostream& os, vtkIndent indent) override;
40 
44  static vtkRungeKutta4 *New();
45 
46 
48 
61  int ComputeNextStep(double* xprev, double* xnext,
62  double t, double& delT,
63  double maxError, double& error) override
64  {
65  double minStep = delT;
66  double maxStep = delT;
67  double delTActual;
68  return this->ComputeNextStep(xprev, nullptr, xnext, t, delT, delTActual,
69  minStep, maxStep, maxError, error);
70  }
71  int ComputeNextStep(double* xprev, double* dxprev, double* xnext,
72  double t, double& delT,
73  double maxError, double& error) override
74  {
75  double minStep = delT;
76  double maxStep = delT;
77  double delTActual;
78  return this->ComputeNextStep(xprev, dxprev, xnext, t, delT, delTActual,
79  minStep, maxStep, maxError, error);
80  }
81  int ComputeNextStep(double* xprev, double* xnext,
82  double t, double& delT, double& delTActual,
83  double minStep, double maxStep,
84  double maxError, double& error) override
85  {
86  return this->ComputeNextStep(xprev, nullptr, xnext, t, delT, delTActual,
87  minStep, maxStep, maxError, error);
88  }
89  int ComputeNextStep(double* xprev, double* dxprev, double* xnext,
90  double t, double& delT, double& delTActual,
91  double minStep, double maxStep,
92  double maxError, double& error) override;
94 
95 protected:
97  ~vtkRungeKutta4() override;
98 
99  void Initialize() override;
100 
101  double* NextDerivs[3];
102 private:
103  vtkRungeKutta4(const vtkRungeKutta4&) = delete;
104  void operator=(const vtkRungeKutta4&) = delete;
105 };
106 
107 #endif
108 
109 
110 
111 
112 
113 
114 
115 
int ComputeNextStep(double *xprev, double *xnext, double t, double &delT, double &delTActual, double minStep, double maxStep, double maxError, double &error) override
Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of...
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 *dxprev, double *xnext, double t, double &delT, double maxError, double &error) 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:33
Integrate an initial value problem using 4th order Runge-Kutta method.
int ComputeNextStep(double *xprev, double *xnext, double t, double &delT, double maxError, double &error) override
Given initial values, xprev , initial time, t and a requested time interval, delT calculate values of...
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.