Method GenerateMatrix

comp/facetfespace.cpp:315–339

    template <typename FEL, typename MIP, typename MAT>
    static void GenerateMatrix (const FEL & bfel, const MIP & mip,
                                MAT && mat, LocalHeap & lh)
    {
      int facetnr = mip.IP().FacetNr();
      if (facetnr >= 0)
        {
          HeapReset hr(lh);
          
          const FacetVolumeFiniteElement<D> & fel_facet = static_cast<const FacetVolumeFiniteElement<D>&> (bfel);
          auto r = fel_facet.GetFacetDofs(facetnr);
          
          FlatMatrix<> dshaperef(r.Size(), DIM_ELEMENT, lh);
          mat.AddSize(DIM_DMAT, bfel.GetNDof()) = 0.0;
          /*
          fel_facet.Facet(facetnr).CalcDShape(mip.IP(), 
                                              mat.Row(0).Range(fel_facet.GetFacetDofs(facetnr)));
          */
          fel_facet.Facet(facetnr).CalcDShape(mip.IP(), dshaperef);
          mat.Cols(r) = Trans(mip.GetJacobianInverse()) * Trans(dshaperef);
        }
      else
        {
            throw Exception("cannot evaluate facet-fe inside element");
        }
    }
  };