#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