/*!\file Pengrid.h
 * \brief: header file for pengrid object */

#ifndef _PENGRID_H_
#define _PENGRID_H_

#include "./Load.h"
#include "./Node.h"
#include "./Element.h"
#include "./Matpar.h"

class Element;
class Pengrid: public Load{

	private: 

		int		id;
		
		Hook hnode;  //hook to 1 node
		Hook hmatpar; //hook to 1 matpar

		Parameters* parameters; //pointer to solution parameters
		Inputs*  inputs;
	
		/*internals: */
		int active;
		int zigzag_counter;

	public:

		/*FUNCTION constructors, destructors {{{1*/
		Pengrid();
		Pengrid(int pengrid_id,int pengrid_node_id int pengrid_matpar_id);
		Pengrid(int pengrid_id,Hook* pengrid_hnode, Hook* pengrid_hmatpar, Parameters* pengrid_parameters, Inputs* pengrid_inputs);
		Pengrid(int i, IoModel* iomodel);
		~Pengrid();
		/*}}}*/
		/*FUNCTION object management {{{1*/
		void  Configure(void* elements,void* nodes,void* materials);
		Object* copy();
		void  DeepEcho();
		void  Demarshall(char** pmarshalled_dataset);
		void  Echo();
		int   Enum();
		int   GetId(); 
		char* GetName();
		void  Marshall(char** pmarshalled_dataset);
		int   MarshallSize();
		int   MyRank();
		/*}}}*/
		/*FUNCTION element numerical routines {{{1*/
		void  DistributeNumDofs(int* numdofspernode,int analysis_type,int sub_analysis_type);
		void  CreateKMatrix(Mat Kgg,int analysis_type,int sub_analysis_type);
		void  CreatePVector(Vec pg, int analysis_type,int sub_analysis_type);
		void  PenaltyCreateKMatrix(Mat Kgg,double kmax,int analysis_type,int sub_analysis_type);
		void  PenaltyCreatePVector(Vec pg,double kmax,int analysis_type,int sub_analysis_type);
		void  PenaltyCreateKMatrixDiagnosticStokes(Mat Kgg,double kmax,int analysis_type,int sub_analysis_type);
		void  GetDofList(int* doflist,int* pnumberofdofspernode);
		void  PenaltyCreateKMatrixThermal(Mat Kgg,double kmax,int analysis_type,int sub_analysis_type);
		void  PenaltyCreateKMatrixMelting(Mat Kgg,double kmax,int analysis_type,int sub_analysis_type);
		void  MatparConfiguration(Matpar* matpar,int matpar_offset);
		void  PenaltyCreatePVectorThermal(Vec pg, double kmax,int analysis_type,int sub_analysis_type);
		void  PenaltyCreatePVectorMelting(Vec pg, double kmax,int analysis_type,int sub_analysis_type);
		void  PenaltyConstrain(int* punstable,int analysis_type,int sub_analysis_type);
		void  PenaltyConstrainThermal(int* punstable,int analysis_type,int sub_analysis_type);
		
		/*updates:*/
		void  UpdateFromDakota(void* inputs);
		void  UpdateInputs(double* solution, int analysis_type, int sub_analysis_type);
		/*}}}*/

};

#endif  /* _PENGRID_H_ */