12 static double NF_N_Q_RT2_2D_Xi[] = {0};
13 static double NF_N_Q_RT2_2D_Eta[] = {0};
19 static double NF_N_Q_RT2_2D_T[] = {-0.866025403784439,0,0.866025403784439};
39 cout <<
"Raviart-Thomas elements of order 2 on quadrilaterals: "
40 <<
"Nodal functionals are not implemented properly!" << endl;
41 for(
int i=0; i<24; i++)
42 Functionals[i] = PointValues[0];
69 l = sqrt((x0-x1)*(x0-x1) + (y0-y1)*(y0-y1));
70 Functionals[0] = PointValues[0]*l;
71 Functionals[1] = PointValues[1]*l;
72 Functionals[2] = PointValues[2]*l;
76 (NF_N_Q_RT2_2D, 24, 3, 1, 3, NF_N_Q_RT2_2D_Xi, NF_N_Q_RT2_2D_Eta,
77 NF_N_Q_RT2_2D_T, NF_N_Q_RT2_2D_EvalAll, NF_N_Q_RT2_2D_EvalEdge);
store cells in an array, used by cell iterators
Definition: Collection.h:18
Definition: NodalFunctional2D.h:20
virtual TVertex * GetVertex(int Vert_i)=0
return the pointer to vertex with number i
void GetCoords(double &x, double &y, double &z) const
Definition: Vertex.h:106
information for finite element data structure
Definition: BaseCell.h:25