BoxConstantForceField.inl
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025 #ifndef SOFA_COMPONENT_BOXCONSTANTFORCEFIELD_INL
00026 #define SOFA_COMPONENT_BOXCONSTANTFORCEFIELD_INL
00027
00028
00029 #include <sofa/component/forcefield/BoxConstantForceField.h>
00030
00031 namespace sofa
00032 {
00033
00034 namespace component
00035 {
00036
00037 namespace forcefield
00038 {
00039
00040 using namespace sofa::helper;
00041
00042
00043 template<class DataTypes>
00044 BoxConstantForceField<DataTypes>::BoxConstantForceField()
00045 : _box( initData( &_box, Vec6(0,0,0,1,1,1), "box", "DOFs in the box defined by xmin,ymin,zmin, xmax,ymax,zmax undergo a constant force") )
00046 , _force( initData( &_force, Deriv(), "force", "the constant force to apply"))
00047 {
00048 }
00049
00050 template<class DataTypes>
00051 void BoxConstantForceField<DataTypes>::init()
00052 {
00053 ConstantForceField<DataTypes>::init();
00054 vector <unsigned> indices;
00055 const Vec6& b=_box.getValue();
00056 this->mstate->getIndicesInSpace( indices, b[0],b[3],b[1],b[4],b[2],b[5] );
00057 for(int i = 0; i < (int)indices.size(); i++)
00058 {
00059 this->points.beginEdit()->push_back(indices[i]);
00060 this->forces.beginEdit()->push_back(_force.getValue());
00061 }
00062 }
00063 }
00064 }
00065 }
00066
00067
00068 #endif
00069