C:/cmcintos/defOrgs/source/physical/Phys_Euler.h

#if defined (_MSC_VER) && (_MSC_VER >= 1000)
#pragma once
#endif
#ifndef _INC_Phys_Euler_422B64F703C8_INCLUDED
#define _INC_Phys_Euler_422B64F703C8_INCLUDED
#include "stdafx.h"
#include "Physics.h"
#include "SpringMassDeformation.h"
/*#include <vnl/MatrixType.h>
#include <vnl/VectorType.h>
#include <vnl/VectorType_fixed.h>
#include <vnl/MatrixType_fixed.h>*/
#include <itkVectorLinearInterpolateImageFunction.h> 

#include <math.h>

#define DEBUG 0
#define BOUND_CHECKING 1

namespace mial
{


        template<class DataType, class TGradientImage, int nDims,class MType = vnl_matrix<DataType>, class VType = vnl_vector<DataType> >
        class Phys_Euler: public Physics<DataType,nDims, MType,VType>
        {
        
                //***************************/
                //      Public definitions
                //***************************/
        public:
                //Smart pointers
                typedef Phys_Euler Self;
                typedef itk::SmartPointer<Self>  Pointer;
                typedef itk::SmartPointer<const Self>  ConstPointer;
                typedef itk::WeakPointer<const Self>  ConstWeakPointer;
                itkNewMacro(Self);


                typedef SpringMassDeformation<DataType,nDims,MType,VType> DeformationType; //Use spring mass deformations specificially

                typedef MType MatrixType;
                typedef VType VectorType;
                typedef TGradientImage GradientImageType;

                struct Error: public Physics<DataType,nDims,MType,VType>::Error
                {
                        // TODO fill in
                };
                //***************************/
                //      Member Variables
                //***************************/

        protected:


                typedef typename GradientImageType::Pointer GradientImageTypePointer;

                MatrixType nodes;

                MatrixType nodesV;

                MatrixType nodesF;

                MatrixType nodesFDef;

                MatrixType nodesA;

                VectorType nodesM;
                
                DataType defaultMass;

                VectorType springsRest;

                VectorType springsDamp;
                
                DataType defaultDamp;

                MatrixType springsNodes;

                VectorType springLengths;

                VectorType springsK;
                
                DataType defaultK;

                DataType defaultDrag;

                double timeStep;

                
                GradientImageTypePointer gradientPointer;

                bool imageForces;
                //***************************/
                //      Member Functions
                //***************************/
        public:


                virtual bool runDeformation(const std::string defName,typename DeformationType::deformationIn* const i,std::stringstream * const s = NULL);


                virtual  void setRestLengths(int* a, DataType* values, int n);


                virtual void setRestLengths(VectorType a, VectorType values);


                virtual  void setTimeStep(double a){timeStep = a;}      


                virtual  void setSpringLengths(int* a, DataType* values, int n);


                virtual void setSpringLengths(VectorType a, VectorType values);

                

                virtual  void setSpringsK(int* a, DataType* values, int n);


                virtual void setSpringsK(VectorType a, VectorType values);

                virtual void enableImageForces(){imageForces = true;};

                virtual void disableImageForces(){imageForces = false;};

                virtual bool simulate();

                virtual void setExternalForces(void * f){gradientPointer = *((GradientImageTypePointer *) f);};

                
        protected:
                        Phys_Euler(int numNodes =0 ,int numSprings = 0, int numPossibleDeformations = 0, int defK = 35);

        private:
                virtual void updateSpringsFromGeometric();
                
        };
} // end namespace mial 
#include "Phys_Euler.cxx"
#endif /* _INC_Phys_Euler_422B64F703C8_INCLUDED */

---