C:/cmcintos/defOrgs/source/physical/abc/Physics.h

#if defined (_MSC_VER) && (_MSC_VER >= 1000)
#pragma once
#endif
#ifndef _INC_PHYSICS
#define _INC_PHYSICS

#include "Geometric.h"
#include "Deformation.h"
#include "stdafx.h"
#include "itkSmartPointer.h"
#include <itkLightObject.h> 
#include <vector>//Not to be confused with VNL vectors
 


namespace mial
{


        template<class Type, int nDims, class MType = vnl_matrix<Type>, class VType = vnl_vector<Type> > class Physics:public itk::LightObject
{

public:

        //Smartpointers
        typedef Physics Self;
        typedef itk::SmartPointer<Self>  Pointer;
        typedef itk::SmartPointer<const Self>  ConstPointer;
        typedef itk::WeakPointer<const Self>  ConstWeakPointer;
        //itkNewMacro(Self);

    //******************************************
    //Define the types
    //******************************************

        
    typedef Deformation<Type,nDims,MType,VType> DeformationType;


    struct Error{
        std::string msg;
        int deformationNumber;
        typename DeformationType::Error * deformationError;
    };

    typedef MType                                           MatrixType;
    typedef VType                                           VectorType;

    typedef Geometric<Type, nDims,MatrixType,VectorType>  GeometryType;


    //***************************************
    //Define the methods
    //***************************************  

    ~Physics();


        virtual bool runDeformation(const std::string defName,typename DeformationType::deformationIn* const i,std::stringstream * const s = NULL)=0;
        
        virtual bool addDeformation(Deformation<Type,nDims> * def);

    virtual bool simulate() =0;


    virtual void setExternalForces(void * img)=0;

        void setGeometry(GeometryType* a){ if(a !=NULL){geom = a;}else{Error e; e.msg="Cannot add null geometric layer"; throw & e;}};

        double getTime(){return time;};

protected:

    Physics();

    int numDeformations;

        typename GeometryType::Pointer geom;

        //use stl vectors
        std::vector<typename Deformation<Type,nDims>::Pointer> deformationsList;

        double time;
};
} // end namespace mial
#include "Physics.cxx"
#endif /* _INC_PHYSICS*/

---