C:/cmcintos/defOrgs/source/geometrical/Geom_vGeometry.cxx

#ifndef _Geom_vGeometry_txx
#define _Geom_vGeometry_txx
#include "Geom_vGeometry.h"

namespace mial
{

template <class dType, int nDims,class MType,class VType>
Geom_vGeometry<dType, nDims,MType,VType>::Geom_vGeometry()
{
        //TODO fill in
}


template <class dType, int nDims,class MType, class VType>
void Geom_vGeometry<dType, nDims,MType,VType>::getNodePosition(int index, dType* container)
{
        //double a[nDims];
        for(int i=0; i<nDims; i++)
        {
                container[i] = nodes(index,i);
        }
        //return a;
}



template <class dType, int nDims,class MType, class VType >
bool Geom_vGeometry<dType, nDims,MType,VType>::setNodePosition(int index, dType* pos)
{
//      nodes(index) = 0;
        for(int i=0; i<nDims; i++)
        {
                nodes(index,i) = pos[i];
        }
        return 1;
}



template <class dType, int nDims,class MType, class VType>
void Geom_vGeometry<dType, nDims,MType,VType>::writeNodesToFile(std::string fileName)
{
        std::ofstream out(fileName.c_str());
        out << nodes.rows() << std::endl;
        for(int i =0; i<nodes.rows(); i++)
        {
                out << nodes.get_row(i) << std::endl;
        }
}

template <class dType, int nDims,class MType,class VType >
void Geom_vGeometry<dType, nDims,MType,VType>::readNodesFromFile(std::string fileName)
{
        std::ifstream in(fileName.c_str());
        int count = 0;
        int size =0;
        in >> size;
        nodes.set_size(size,nDims);
        while(count < size)
        {
                for(int i=0; i<nDims; i++)
                {
                        in >> nodes(count,i);
                }
                count++;
        }
}
}//end mial
#endif

---