forked from Kitware/VTK
-
Notifications
You must be signed in to change notification settings - Fork 0
/
vtkRungeKutta2.h
99 lines (87 loc) · 3.33 KB
/
vtkRungeKutta2.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
/*=========================================================================
Program: Visualization Toolkit
Module: vtkRungeKutta2.h
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
/**
* @class vtkRungeKutta2
* @brief Integrate an initial value problem using 2nd
* order Runge-Kutta method.
*
*
* This is a concrete sub-class of vtkInitialValueProblemSolver.
* It uses a 2nd order Runge-Kutta method to obtain the values of
* a set of functions at the next time step.
*
* @sa
* vtkInitialValueProblemSolver vtkRungeKutta4 vtkRungeKutta45 vtkFunctionSet
*/
#ifndef vtkRungeKutta2_h
#define vtkRungeKutta2_h
#include "vtkCommonMathModule.h" // For export macro
#include "vtkInitialValueProblemSolver.h"
class VTKCOMMONMATH_EXPORT vtkRungeKutta2 : public vtkInitialValueProblemSolver
{
public:
vtkTypeMacro(vtkRungeKutta2, vtkInitialValueProblemSolver);
void PrintSelf(ostream& os, vtkIndent indent) override;
/**
* Construct a vtkRungeKutta2 with no initial FunctionSet.
*/
static vtkRungeKutta2* New();
using Superclass::ComputeNextStep;
///@{
/**
* Given initial values, xprev , initial time, t and a requested time
* interval, delT calculate values of x at t+delT (xnext).
* delTActual is always equal to delT.
* Since this class can not provide an estimate for the error error
* is set to 0.
* maxStep, minStep and maxError are unused.
* This method returns an error code representing the nature of
* the failure:
* OutOfDomain = 1,
* NotInitialized = 2,
* UnexpectedValue = 3
*/
int ComputeNextStep(double* xprev, double* xnext, double t, double& delT, double maxError,
double& error, void* userData) override
{
double minStep = delT;
double maxStep = delT;
double delTActual;
return this->ComputeNextStep(
xprev, nullptr, xnext, t, delT, delTActual, minStep, maxStep, maxError, error, userData);
}
int ComputeNextStep(double* xprev, double* dxprev, double* xnext, double t, double& delT,
double maxError, double& error, void* userData) override
{
double minStep = delT;
double maxStep = delT;
double delTActual;
return this->ComputeNextStep(
xprev, dxprev, xnext, t, delT, delTActual, minStep, maxStep, maxError, error, userData);
}
int ComputeNextStep(double* xprev, double* xnext, double t, double& delT, double& delTActual,
double minStep, double maxStep, double maxError, double& error, void* userData) override
{
return this->ComputeNextStep(
xprev, nullptr, xnext, t, delT, delTActual, minStep, maxStep, maxError, error, userData);
}
int ComputeNextStep(double* xprev, double* dxprev, double* xnext, double t, double& delT,
double& delTActual, double minStep, double maxStep, double maxError, double& error,
void* userData) override;
///@}
protected:
vtkRungeKutta2();
~vtkRungeKutta2() override;
private:
vtkRungeKutta2(const vtkRungeKutta2&) = delete;
void operator=(const vtkRungeKutta2&) = delete;
};
#endif