C:/cmcintos/defOrgs/source/physical/Def_Translation.cxx

#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

---