FE_N_H_Q4_3D.h

// ***********************************************************************
// Q4 element, nonconforming, 3D
// ***********************************************************************

// number of degrees of freedom
static int N_H_Q4_3D_NDOF = 70;

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

// which local dofs are on the joints
static int N_H_Q4_3D_J0[10] = { 0,  6, 12, 18, 24, 30, 36, 42, 48, 54 };
static int N_H_Q4_3D_J1[10] = { 1,  7, 13, 19, 25, 31, 37, 43, 49, 55 };
static int N_H_Q4_3D_J2[10] = { 2,  8, 14, 20, 26, 32, 38, 44, 50, 56 };
static int N_H_Q4_3D_J3[10] = { 3,  9, 15, 21, 27, 33, 39, 45, 51, 57 };
static int N_H_Q4_3D_J4[10] = { 4, 10, 16, 22, 28, 34, 40, 46, 52, 58 };
static int N_H_Q4_3D_J5[10] = { 5, 11, 17, 23, 29, 35, 41, 47, 53, 59 };

static int *N_H_Q4_3D_J[6] = { N_H_Q4_3D_J0, N_H_Q4_3D_J1,
                               N_H_Q4_3D_J2, N_H_Q4_3D_J3,
                               N_H_Q4_3D_J4, N_H_Q4_3D_J5};

// number of inner dofs
static int N_H_Q4_3D_NInner = 10;

// array containing the numbers for the inner dofs
static int N_H_Q4_3D_Inner[10] = { 60, 61, 62, 63, 64, 65, 66, 67, 68, 69 };

static char N_H_Q4_3D_String[] = "N_H_Q4_3D";

TFEDesc3D *FE_N_H_Q4_3D_Obj=new TFEDesc3D(N_H_Q4_3D_String, N_H_Q4_3D_NDOF, 
                                N_H_Q4_3D_JointDOF,
                                N_H_Q4_3D_J, N_H_Q4_3D_NInner, N_H_Q4_3D_Inner);