C:/cmcintos/defOrgs/source/organism/Org_EulerSchedule.h

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

#include "itkOrganism.h"

#include "itkImageToImageFilter.h"
#include "itkImage.h"
#include "itkNumericTraits.h"

#include "itkImageRegionIterator.h"

//#include "itkImage.h"
#include <iostream>
#include "Phys_Euler.h"
#include "Ctrl_ScheduleDriven.h"
#include "Beh_TranslateAll.h"
#include "Def_Translation.h"
#include "Beh_UniformScale.h"
#include "Def_UniformScale.h"
#include "Sense_Gradient.h"
#include "Geom_MeshSpatialObject.h"

using namespace mial;
namespace itk
{

        template <class TInputImage, class TOutputImage,class TExternalForceImage, class DataType, int nDims>
        class ITK_EXPORT Org_EulerSchedule : public ItkOrganism<TInputImage,TOutputImage,TExternalForceImage,DataType,nDims>
        {
        public:
                itkStaticConstMacro(InputImageDimension, unsigned int,TInputImage::ImageDimension);

                itkStaticConstMacro(OutputImageDimension, unsigned int,TOutputImage::ImageDimension);

                typedef Org_EulerSchedule Self;
                typedef SmartPointer<Self> Pointer;
                typedef SmartPointer<const Self>  ConstPointer;
                itkNewMacro(Self);

                itkTypeMacro(Org_EulerSchedule, ItkOrganism);

                typedef Geom_MeshSpatialObject<float,3> GeometricType;
                typedef Phys_Euler<float,TExternalForceImage,3> PhysLayerType;
                typedef Sense_Gradient<float,TInputImage,TExternalForceImage,3> gradientSensorType;

                virtual bool setSchedule(std::string scheduleFileName){return cgL->setSchedule(scheduleFileName);};

                virtual bool setTopology(std::string fName){
                        return geomLayer->readTopologyFromFile(fName);};

                virtual void setUp()
                {
                        input->sigma = 1.0;
                        input->imageIn = this->GetInput();
                        gradientSensor->run(input);
                        gradientSensor->run(input);
                        typename gradientSensorType::sensorOut::Pointer output = (typename gradientSensorType::sensorOut *) ( gradientSensor->getOutput()).GetPointer();
                        physLayer->setExternalForces((void *) &(output->imageOut));
                }


        protected:
                Ctrl_ScheduleDriven<float, 3>::Pointer cgL;
                typename GeometricType::Pointer geomLayer;
                typename PhysLayerType::Pointer physLayer;
                typename gradientSensorType::sensorIn::Pointer input;
                typename gradientSensorType::Pointer gradientSensor;

                Org_EulerSchedule();
                virtual ~Org_EulerSchedule() {}
                //      void PrintSelf(std::ostream& os, Indent indent) const;

        };
}//end mial namespace
#include "Org_EulerSchedule.cxx"
#endif /* _INC_Org_EulerSchedule*/

---