FE_C_T_P3_3D.h

// ***********************************************************************
// P3 element, conforming, 3D
// ***********************************************************************

// number of degrees of freedom
static int C_T_P3_3D_NDOF = 20;

// number of dofs on the closure of the joints
static int C_T_P3_3D_JointDOF = 10;

// which local dofs are on the joints
static int C_T_P3_3D_J0[10] = { 0,  1,  2,  3,  4,  5,  6,  7,  8,  9 };
static int C_T_P3_3D_J1[10] = { 0, 10, 16, 19,  1, 11, 17,  2, 12,  3 };
static int C_T_P3_3D_J2[10] = { 9,  8,  6,  3, 15, 14, 12, 18, 17, 19 };
static int C_T_P3_3D_J3[10] = { 0,  4,  7,  9, 10, 13, 15, 16, 18, 19 };

static int *C_T_P3_3D_J[4] = { C_T_P3_3D_J0, C_T_P3_3D_J1,
                             C_T_P3_3D_J2, C_T_P3_3D_J3 };

// number of inner dofs
static int C_T_P3_3D_NInner = 0;

// array containing the numbers for the inner dofs (here is no inner dof)
static int *C_T_P3_3D_Inner = NULL;

static char C_T_P3_3D_String[] = "C_T_P3_3D";

TFEDesc3D *FE_C_T_P3_3D_Obj=new TFEDesc3D(C_T_P3_3D_String, C_T_P3_3D_NDOF, 
                                C_T_P3_3D_JointDOF,
                                C_T_P3_3D_J, C_T_P3_3D_NInner, C_T_P3_3D_Inner);