Method IntegrationRuleTP

fem/intrule.cpp:1579–1597

  template <>
  IntegrationRuleTP<1> :: IntegrationRuleTP (const ElementTransformation & eltrans,
                                             IVec<1> vorder, bool compute_duffy, bool compute_points) 
  {
    int order = vorder[0];
    irx = &SelectIntegrationRule (ET_SEGM, order);
    int nip = irx->GetNIP();

    SetSize(nip);
    dxdxi_duffy.SetSize(nip);

    for (int i = 0; i < irx->GetNIP(); i++)
      {
        (*this)[i] = IntegrationPoint ((*irx)[i](0), 0, 0, 
                                       (*irx)[i].Weight());
      }
        
    for (int i = 0; i < nip; i++)
      dxdxi_duffy[i] = 1;
  }


  template <>