/*!\file:  ElementVector.h
 * \brief container for information needed to plug element vector generated by elements 
 * into the pf global load vector. 
 * This object will hold the element vector on the g-set, the local as well as global 
 * dof lists in the f set
 */ 

#ifndef _ELEMENT_VECTOR_H_
#define _ELEMENT_VECTOR_H_

/*Headers:*/
#include "../../datastructures/datastructures.h"
#include "../../toolkits/toolkits.h"
#include "../../shared/shared.h"
class Node;
template <class doubletype> class Vector;
class Parameters;
class ElementMatrix;

class ElementVector{

	private:
		int *gglobaldoflist;
		int *fglobaldoflist;

	public:
		int         nrows;
		IssmDouble *values;

		/*ElementVector constructors, destructors*/
		ElementVector();
		ElementVector(ElementVector* pe1,ElementVector* pe2);
		ElementVector(ElementVector* pe1,ElementVector* pe2,ElementVector* pe3);
		ElementVector(Node** nodes,int numnodes,Parameters* parameters,int approximation=NoneApproximationEnum);
		~ElementVector();

		/*ElementVector specific routines*/
		void AddToGlobal(Vector<IssmDouble>* pf);
		void CheckConsistency(void);
		void Echo(void);
		void Init(ElementVector* pe);
		void InsertIntoGlobal(Vector<IssmDouble>* pf);
		void SetValue(IssmDouble scalar);
		void StaticCondensation(ElementMatrix* Ke,int numindices,int* indices);
};
#endif //#ifndef _ELEMENT_VECTOR_H_
