00001 /****************************************************************************** 00002 * SOFA, Simulation Open-Framework Architecture, version 1.0 beta 3 * 00003 * (c) 2006-2008 MGH, INRIA, USTL, UJF, CNRS * 00004 * * 00005 * This library is free software; you can redistribute it and/or modify it * 00006 * under the terms of the GNU Lesser General Public License as published by * 00007 * the Free Software Foundation; either version 2.1 of the License, or (at * 00008 * your option) any later version. * 00009 * * 00010 * This library is distributed in the hope that it will be useful, but WITHOUT * 00011 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or * 00012 * FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License * 00013 * for more details. * 00014 * * 00015 * You should have received a copy of the GNU Lesser General Public License * 00016 * along with this library; if not, write to the Free Software Foundation, * 00017 * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. * 00018 ******************************************************************************* 00019 * SOFA :: Framework * 00020 * * 00021 * Authors: M. Adam, J. Allard, B. Andre, P-J. Bensoussan, S. Cotin, C. Duriez,* 00022 * H. Delingette, F. Falipou, F. Faure, S. Fonteneau, L. Heigeas, C. Mendoza, * 00023 * M. Nesme, P. Neumann, J-P. de la Plata Alcade, F. Poyer and F. Roy * 00024 * * 00025 * Contact information: contact@sofa-framework.org * 00026 ******************************************************************************/ 00027 #ifndef SOFA_CORE_COMPONENTMODEL_BEHAVIOR_ODESOLVER_H 00028 #define SOFA_CORE_COMPONENTMODEL_BEHAVIOR_ODESOLVER_H 00029 00030 #include <sofa/core/objectmodel/BaseObject.h> 00031 #include <sofa/core/componentmodel/behavior/BaseMechanicalState.h> 00032 #include <sofa/core/componentmodel/behavior/MultiVector.h> 00033 #include <sofa/core/componentmodel/behavior/MultiMatrix.h> 00034 #include <sofa/defaulttype/BaseMatrix.h> 00035 #include <sofa/defaulttype/BaseVector.h> 00036 00037 namespace sofa 00038 { 00039 00040 namespace core 00041 { 00042 00043 namespace componentmodel 00044 { 00045 00046 namespace behavior 00047 { 00048 00066 class OdeSolver : public virtual objectmodel::BaseObject 00067 { 00068 public: 00069 typedef BaseMechanicalState::VecId VecId; 00070 00071 OdeSolver(); 00072 00073 virtual ~OdeSolver(); 00074 00079 virtual void solve (double dt) = 0; 00080 00084 virtual void propagatePositionAndVelocity(double t, BaseMechanicalState::VecId x, BaseMechanicalState::VecId v) = 0; 00085 00101 virtual double getIntegrationFactor(int inputDerivative, int outputDerivative) const = 0; 00102 00106 virtual double getSolutionIntegrationFactor(int outputDerivative) const = 0; 00107 00108 00112 virtual double getVelocityIntegrationFactor() const 00113 { 00114 return getSolutionIntegrationFactor(1); 00115 } 00116 00120 virtual double getPositionIntegrationFactor() const 00121 { 00122 return getSolutionIntegrationFactor(0); 00123 } 00124 00125 }; 00126 00127 } // namespace behavior 00128 00129 } // namespace componentmodel 00130 00131 } // namespace core 00132 00133 } // namespace sofa 00134 00135 #endif
1.5.7.1