FE_N_T_RT2_2D.h

// ***********************************************************************
// P1 Raviart-Thomas vector element, nonconforming , 2D
// History:  15.11.2011 implementation (Ulrich)
// ***********************************************************************

// number of degrees of freedom
static int N_T_RT2_2D_NDOF = 15;

// number of dofs on the closure of the joints
static int N_T_RT2_2D_JointDOF = 3;

// which local dofs are on the joints
static int N_T_RT2_2D_J0[3] = { 0, 1, 2};
static int N_T_RT2_2D_J1[3] = { 3, 4, 5};
static int N_T_RT2_2D_J2[3] = { 6, 7, 8};

 
static int *N_T_RT2_2D_J[3] = { N_T_RT2_2D_J0,
                                N_T_RT2_2D_J1,
                                N_T_RT2_2D_J2
                              };

// number of inner dofs
static int N_T_RT2_2D_NInner = 6;

// array containing the numbers for the inner dofs
static int N_T_RT2_2D_Inner[6] = {9, 10, 11, 12, 13, 14};

// number of outer dofs (dofs on edges)
static int N_T_RT2_2D_NOuter = 9;

// array containing the numbers for the outer dofs
static int N_T_RT2_2D_Outer[9] = { 0, 1, 2, 3, 4, 5, 6, 7, 8};

static char N_T_RT2_2D_String[] = "N_T_RT2_2D";

TFEDesc2D *FE_N_T_RT2_2D_Obj=new TFEDesc2D(N_T_RT2_2D_String, N_T_RT2_2D_NDOF,
                                        N_T_RT2_2D_JointDOF, N_T_RT2_2D_J,
                                        N_T_RT2_2D_NInner, N_T_RT2_2D_Inner,
                                        N_T_RT2_2D_NOuter, N_T_RT2_2D_Outer);