#ifndef _Org_EulerSchedule_txx
#define _Org_EulerSchedule_txx
#include "Org_EulerSchedule.h"

namespace itk
{

	template <class TInputImage, class TOutputImage, class TExternalForceImage, class DataType, int nDims>
	Org_EulerSchedule<TInputImage, TOutputImage, TExternalForceImage, DataType, nDims>::Org_EulerSchedule()
	{
		//Instatiate a sensor
		input = gradientSensorType::sensorIn::New();
		gradientSensor = gradientSensorType::New();
		//gradientSensor.sensorIn input;
		


		//Instantiate geomtery and physics layers
		physLayer = PhysLayerType::New();
		geomLayer = GeometricType::New();

		
		physLayer->setGeometry(geomLayer);

		this->setPhysicsLayer(physLayer);
		this->setGeometricLayer(geomLayer);
		std::cout << "Physics layer added..." << std::endl;

		//testOrg->writeNodesToFile(nodeOutputFileName);

		//TEST BEHAVIOR
		cgL = Ctrl_ScheduleDriven<DataType, nDims>::New();
		
		this->beh0 = Beh_SearchForObject<DataType,TInputImage,nDims>::New();
		typename Beh_TranslateAll<DataType, nDims>::Pointer beh1 = Beh_TranslateAll<DataType,nDims>::New();
		typename Beh_UniformScale<DataType, nDims>::Pointer beh2 = Beh_UniformScale<DataType,nDims>::New();
		typename Def_Translation<DataType, nDims>::Pointer def1 = Def_Translation<DataType,nDims>::New();								
		typename Def_UniformScale<DataType, nDims>::Pointer def2 = Def_UniformScale<DataType,nDims>::New();								
		this->setCognitiveLayer(cgL);
		this->addBehaviour(beh0);
		this->addBehaviour(beh1);
		this->addBehaviour(beh2);
		this->addDeformation(def1);
		this->addDeformation(def2);

	}

}
#endif