#if !defined(VGEOMETRY)
#define VGEOMETRY

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000

#include <vnl/vnl_matrix.h>

#include "Geometric.h"

//vType can be any data type that has the following operations defined:
//vector[element] = number. 

namespace mial{

template <class dType, int nDims,class MType = vnl_matrix<dType>, class VType = vnl_vector<dType> >
class Geom_vGeometry : public Geometric<dType, nDims,MType,VType>  
{
public:
	typedef MType MatrixType;
	typedef VType VectorType;
	private:
	
		MatrixType nodes;
	
	public:
	
		Geom_vGeometry();
	//	virtual ~vGeometry();
	
		virtual MatrixType getVectorNodePositions(){return nodes;};
	
		//virtual bool setVectNodePositions(vType a);
	
		virtual bool setNodePosition(int index,dType* pos);
	
		virtual void getNodePosition(int index, dType *);
		
	
		virtual void writeNodesToFile(std::string fileName);

	
		virtual void readNodesFromFile(std::string fileName);
		
	
		virtual unsigned int numNodes(){ return (unsigned int)nodes.rows(); }
};
}//end mial
#include "Geom_vGeometry.cxx"
#endif