00001 /*========================================================================= 00002 00003 Program: Visualization Toolkit 00004 Module: $RCSfile: vtkRungeKutta4.h,v $ 00005 Language: C++ 00006 00007 00008 Copyright (c) 1993-2001 Ken Martin, Will Schroeder, Bill Lorensen 00009 All rights reserved. 00010 00011 Redistribution and use in source and binary forms, with or without 00012 modification, are permitted provided that the following conditions are met: 00013 00014 * Redistributions of source code must retain the above copyright notice, 00015 this list of conditions and the following disclaimer. 00016 00017 * Redistributions in binary form must reproduce the above copyright notice, 00018 this list of conditions and the following disclaimer in the documentation 00019 and/or other materials provided with the distribution. 00020 00021 * Neither name of Ken Martin, Will Schroeder, or Bill Lorensen nor the names 00022 of any contributors may be used to endorse or promote products derived 00023 from this software without specific prior written permission. 00024 00025 * Modified source versions must be plainly marked as such, and must not be 00026 misrepresented as being the original software. 00027 00028 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS ``AS IS'' 00029 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 00030 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 00031 ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHORS OR CONTRIBUTORS BE LIABLE FOR 00032 ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 00033 DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR 00034 SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER 00035 CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, 00036 OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE 00037 OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 00038 00039 =========================================================================*/ 00062 #ifndef __vtkRungeKutta4_h 00063 #define __vtkRungeKutta4_h 00064 00065 #include "vtkInitialValueProblemSolver.h" 00066 00067 class VTK_COMMON_EXPORT vtkRungeKutta4 : public vtkInitialValueProblemSolver 00068 { 00069 public: 00070 vtkTypeMacro(vtkRungeKutta4,vtkInitialValueProblemSolver); 00071 virtual void PrintSelf(ostream& os, vtkIndent indent); 00072 00074 static vtkRungeKutta4 *New(); 00075 00077 00081 virtual float ComputeNextStep(float* xprev, float* xnext, float t, 00082 float delT) 00083 { 00084 return this->ComputeNextStep(xprev, 0, xnext, t, delT); 00085 } 00086 virtual float ComputeNextStep(float* xprev, float* dxprev, float* xnext, 00087 float t, float delT); 00089 00091 00092 virtual vtkInitialValueProblemSolver* MakeObject() 00093 { 00094 return vtkRungeKutta4::New(); 00095 } 00097 00098 00099 protected: 00100 vtkRungeKutta4(); 00101 ~vtkRungeKutta4(); 00102 00103 virtual void Initialize(); 00104 00105 float* NextDerivs[3]; 00106 private: 00107 vtkRungeKutta4(const vtkRungeKutta4&); // Not implemented. 00108 void operator=(const vtkRungeKutta4&); // Not implemented. 00109 }; 00110 00111 #endif 00112 00113 00114 00115 00116 00117 00118 00119