FE_D_Q_Q4_2D.h

// ***********************************************************************
// Q4 element, discontinous, 2D, quadrilateral
// ***********************************************************************

// number of degrees of freedom
static int D_Q_Q4_2D_NDOF = 25;

// number of dofs on the closure of the joints
static int D_Q_Q4_2D_JointDOF = 0;

// which local dofs are on the joints
static int *D_Q_Q4_2D_J0 = NULL;
static int *D_Q_Q4_2D_J1 = NULL;
static int *D_Q_Q4_2D_J2 = NULL;
static int *D_Q_Q4_2D_J3 = NULL;

static int *D_Q_Q4_2D_J[4] = { D_Q_Q4_2D_J0, D_Q_Q4_2D_J1,
                               D_Q_Q4_2D_J2, D_Q_Q4_2D_J3 };

// number of inner dofs
static int D_Q_Q4_2D_NInner = 25;

// array containing the numbers for the inner dofs (here is no inner dof)
static int D_Q_Q4_2D_Inner[25] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
                                   10, 11, 12, 13, 14, 15, 16, 17,
                                   18, 19, 20, 21, 22, 23, 24 };

static char D_Q_Q4_2D_String[] = "D_Q_Q4_2D";

TFEDesc2D *FE_D_Q_Q4_2D_Obj=new TFEDesc2D(D_Q_Q4_2D_String, D_Q_Q4_2D_NDOF, 
                              D_Q_Q4_2D_JointDOF,
                              D_Q_Q4_2D_J, D_Q_Q4_2D_NInner, D_Q_Q4_2D_Inner);