NF_D_Q_P1_2D.h

static double NF_D_Q_P1_2D_Xi[9]={ 0.774596669241483, -0.774596669241483,
                                 0.774596669241483, -0.774596669241483,
                                 0.774596669241483, -0.774596669241483,
                                 0,                  0,
                                 0 };
static double NF_D_Q_P1_2D_Eta[9]={  0.774596669241483,  0.774596669241483,
                                  -0.774596669241483, -0.774596669241483,
                                   0,                  0,
                                   0.774596669241483, -0.774596669241483,
                                   0 };

static double *NF_D_Q_P1_2D_t = NULL;

void NF_D_Q_P1_2D_EvalAll(TCollection *Coll, TBaseCell *Cell, double *PointValues,
                          double *Functionals)
{
  static double weights[9]={ 0.308641975308642, 0.308641975308642,
                             0.308641975308642, 0.308641975308642,
                             0.493827160493827, 0.493827160493827,
                             0.493827160493827, 0.493827160493827,
                             0.790123456790123 };
  Functionals[0] =  weights[0]*PointValues[0]
                   +weights[1]*PointValues[1]
                   +weights[2]*PointValues[2]
                   +weights[3]*PointValues[3]
                   +weights[4]*PointValues[4]
                   +weights[5]*PointValues[5]
                   +weights[6]*PointValues[6]
                   +weights[7]*PointValues[7]
                   +weights[8]*PointValues[8];
  Functionals[1] =  weights[0]*PointValues[0]*NF_D_Q_P1_2D_Xi[0]
                   +weights[1]*PointValues[1]*NF_D_Q_P1_2D_Xi[1]
                   +weights[2]*PointValues[2]*NF_D_Q_P1_2D_Xi[2]
                   +weights[3]*PointValues[3]*NF_D_Q_P1_2D_Xi[3]
                   +weights[4]*PointValues[4]*NF_D_Q_P1_2D_Xi[4]
                   +weights[5]*PointValues[5]*NF_D_Q_P1_2D_Xi[5]
                   +weights[6]*PointValues[6]*NF_D_Q_P1_2D_Xi[6]
                   +weights[7]*PointValues[7]*NF_D_Q_P1_2D_Xi[7]
                   +weights[8]*PointValues[8]*NF_D_Q_P1_2D_Xi[8];
  Functionals[2] =  weights[0]*PointValues[0]*NF_D_Q_P1_2D_Eta[0]
                   +weights[1]*PointValues[1]*NF_D_Q_P1_2D_Eta[1]
                   +weights[2]*PointValues[2]*NF_D_Q_P1_2D_Eta[2]
                   +weights[3]*PointValues[3]*NF_D_Q_P1_2D_Eta[3]
                   +weights[4]*PointValues[4]*NF_D_Q_P1_2D_Eta[4]
                   +weights[5]*PointValues[5]*NF_D_Q_P1_2D_Eta[5]
                   +weights[6]*PointValues[6]*NF_D_Q_P1_2D_Eta[6]
                   +weights[7]*PointValues[7]*NF_D_Q_P1_2D_Eta[7]
                   +weights[8]*PointValues[8]*NF_D_Q_P1_2D_Eta[8];

  Functionals[0] *= 0.25;
  Functionals[1] *= 0.25;
  Functionals[2] *= 0.25;
}

void NF_D_Q_P1_2D_EvalEdge(TCollection *Coll, TBaseCell *Cell, int Joint, double *PointValues,
                           double *Functionals)
{
}

/*
    TNodalFunctional2D(NodalFunctional2D id,
                       int n_allfunctionals, int n_edgefunctionals,
                       int n_pointsall, int n_pointsedge,
                       double *xi, double *eta, double *t,
                       DoubleFunctVect *evalall,
                       DoubleFunctVect *evaledge);
*/

TNodalFunctional2D *NF_D_Q_P1_2D_Obj = new TNodalFunctional2D
        (NF_D_Q_P1_2D, 3, 0, 9, 0, NF_D_Q_P1_2D_Xi, NF_D_Q_P1_2D_Eta,
         NF_D_Q_P1_2D_t, NF_D_Q_P1_2D_EvalAll, NULL);