#ifndef _Org_LevelSetSchedule_txx
#define _Org_LevelSetSchedule_txx
#include "Org_LevelSetSchedule.h"

namespace itk
{

	template <class TInputImage, class TOutputImage, class TExternalForceImage, class DataType, int nDims>
	Org_LevelSetSchedule<TInputImage, TOutputImage, TExternalForceImage, DataType, nDims>
		::Org_LevelSetSchedule()
	{
		
		//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();
		
		
		typename Beh_UniformScale<DataType, nDims>::Pointer beh2 = Beh_UniformScale<DataType,nDims>::New();
		
		typename Def_UniformScaleLevelSet<DataType,nDims,typename PhysLayerType::InternalImageType>::Pointer def2 = Def_UniformScaleLevelSet<DataType,nDims,typename PhysLayerType::InternalImageType>::New();								
		
		
		this->setCognitiveLayer(cgL);
		
		this->addBehaviour(beh2);
		
		this->addDeformation(def2);

	}

}
#endif