#ifndef _DEF_TRANSLATION_TXX
#define _DEF_TRANSLATION_TXX

#include "Def_Translation.h"

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

	
	template<class DataType,int nDims,class MType, class VType>
	bool Def_Translation<DataType, nDims,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_TranslateAll";
			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_TranslateAll";
			throw e;
		}
		VType magVec(nDims);
		magVec.fill(0);
		//vnl_Matrix tmp = org.nodesF;
		int numNodes = orgUp->nodesF->rows();
		for( unsigned int i=0; i<nDims; i++ )
		{
			magVec(i) = in->amplitude[i]; //Read in the translation magnitude for each dimension
			//magVec(i) = 1000;
		}
		//Apply the force to the model
		for( unsigned int i=0; i<numNodes; i++ )
		{
			//orgUp->nodesFDef->set_row( i, magVec );
			orgUp->nodesV->set_row(i,magVec);
		}

		return false;
	}

} // end namespace mial

#endif