element_formulation.h

vs.lib is a math library in C++ with a set of linear algebra and integrable / differentiable objects

#ifndef __ELEMENT_FORMULATION
#define __ELEMENT_FORMULATION

//==============================================================================
// Part C. Element Formulation
//==============================================================================

class Element_Type_Register { public: Element_Type_Register() {} };
class Element_Formulation {
	Global_Discretization& the_global_discretization;
   Element_Formulation* rep_ptr;
   int reference_count;
protected:
	int element_no, material_type_no, nen, nsd, ndf;
	C0 xl,               // nodal coordinates
      gl,               // nodal fixed variables
      ul,               // nodal free variables
   	stiff, force,     // calculated in user supplied Element Formulation Constructor
      the_reaction,     // - stiff * gl
      the_lhs, the_rhs, // to be mapped to Matrix_Representation
      the_element_nodal_value;
   virtual C0& __lhs();
   virtual C0& __reaction();
   virtual C0& __rhs();
   virtual C0& __element_nodal_value();
public:
	static Element_Formulation *type_list;
	Element_Formulation *next;
	Element_Formulation(Element_Type_Register);
   Element_Formulation();
	Element_Formulation(int en, Global_Discretization& gd);
   virtual ~Element_Formulation();
   Element_Formulation(const Element_Formulation&);
   Element_Formulation& operator=(const Element_Formulation&);
   Element_Formulation& create(int en, Global_Discretization& gd);
   virtual Element_Formulation* make(int en, Global_Discretization& gd) { return 0; }
   int& count() { return reference_count; }
   Element_Formulation* &rep() { return rep_ptr; }
	// to be called by Matrix_Representation assembly memember function
	C0& lhs() { return rep_ptr->__lhs(); }
	C0& reaction() { return rep_ptr->__reaction(); }
	C0& rhs() { return rep_ptr->__rhs(); }
   C0& element_nodal_value() { return rep_ptr->__element_nodal_value(); }
};
#endif