#ifndef _DEF_UniformScaleLevelSet_TXX
#define _DEF_UniformScaleLevelSet_TXX

#include "Def_UniformScaleLevelSet.h"

namespace mial
{
	// constructor:
	template<class DataType,int nDims,class TDistanceImageType,class MType, class VType>
	Def_UniformScaleLevelSet<DataType,nDims,TDistanceImageType,MType,VType>
		::Def_UniformScaleLevelSet()
	{
		this->name = "Def_UniformScale";
	}

	template<class DataType,int nDims,class TDistanceImageType,class MType, class VType>
	bool Def_UniformScaleLevelSet<DataType,nDims,TDistanceImageType,MType,VType>
		::run(typename Deformation<DataType,nDims>::deformationIn* i,typename Deformation<DataType,nDims>::DefArgSet* org, std::stringstream *s)
	{
		DefArgSet* orgUp = reinterpret_cast<DefArgSet *>(org);
		typename deformationIn::Pointer in;
		if(i !=NULL && s != NULL)
		{
			Error e;
			e.msg = "Only one of struct or stream input may be provided.";
			e.name = "Def_UniformScaleLevelSet";
			throw e;
		}
		else if(i !=NULL)
			in = reinterpret_cast<deformationIn *>(i); //Typecast the input to its desired form
		else if( s!= NULL)
		{
			in = deformationIn::New();
			in->fillFromStream(*s);
		}
		else
		{
			Error e;
			e.msg = "Either struct or stream input must be provided.";
			e.name = "Def_UniformScaleLevelSet";
			throw e;
		}
		//Add the scale factor to each pixel in the distance image
		IteratorType phiIt( orgUp->distanceImg, orgUp->distanceImg->GetRequestedRegion() );
		DataType sc = in->scaleAmount;
		for ( phiIt.GoToBegin(); !phiIt.IsAtEnd(); ++phiIt)
		{
			phiIt.Set( phiIt.Get() + sc);
		}

		return false;
	}

} // end namespace mial

#endif