FE_C_Q_M8_2D.h

// ***********************************************************************
// M8 element, conforming, 2D
// ***********************************************************************

// number of degrees of freedom
static int C_Q_M8_2D_NDOF = 47;

// number of dofs on the closure of the joints
static int C_Q_M8_2D_JointDOF = 9;

// which local dofs are on the joints
static int C_Q_M8_2D_J0[9] = {  0,  1,  2,  3,  4,  5,  6,  7,  8 };
static int C_Q_M8_2D_J1[9] = {  8,  9, 10, 11, 12, 13, 14, 15, 16 };
static int C_Q_M8_2D_J2[9] = { 16, 17, 18, 19, 20, 21, 22, 23, 24 };
static int C_Q_M8_2D_J3[9] = { 24, 25, 26, 27, 28, 29, 30, 31,  0 };

static int *C_Q_M8_2D_J[4] = { C_Q_M8_2D_J0, C_Q_M8_2D_J1,
                               C_Q_M8_2D_J2, C_Q_M8_2D_J3 };

// number of inner dofs
static int C_Q_M8_2D_NInner = 15;

// array containing the numbers for the inner dofs
static int C_Q_M8_2D_Inner[15] = { 32, 33, 34, 35, 36, 37, 38, 39,
                                   40, 41, 42, 43, 44, 45, 46 };

// number of outer dofs
static int C_Q_M8_2D_NOuter = 32;

// array containing the numbers for the outer dofs
static int C_Q_M8_2D_Outer[32] = { 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, 25, 26, 27, 28, 29,
                                  30, 31 };

static char C_Q_M8_2D_String[] = "C_Q_M8_2D";

TFEDesc2D *FE_C_Q_M8_2D_Obj=new TFEDesc2D(C_Q_M8_2D_String, C_Q_M8_2D_NDOF,
                                C_Q_M8_2D_JointDOF, C_Q_M8_2D_J,
                                C_Q_M8_2D_NInner, C_Q_M8_2D_Inner,
                                C_Q_M8_2D_NOuter, C_Q_M8_2D_Outer);