Method AddDualTrans

fem/tangentialfacetfe.cpp:1192–1209

  template <ELEMENT_TYPE ET>
  void TangentialFacetVolumeFE<ET> :: AddDualTrans (const SIMD_BaseMappedIntegrationRule & bmir, BareSliceMatrix<SIMD<double>> values, BareSliceVector<double> coefs) const
  {
    Switch<4-DIM>
      (bmir.DimSpace()-DIM,[this,&bmir,coefs,values](auto CODIM)
       {
         constexpr int DIMSPACE = DIM+CODIM.value;
         auto & mir = static_cast<const SIMD_MappedIntegrationRule<DIM,DIMSPACE>&> (bmir);
         for (size_t i = 0; i < mir.Size(); i++)
           {
             Vec<DIMSPACE,SIMD<double>> value = values.Col(i);
             this -> CalcDualShape2 (mir[i], mir[i].IP().FacetNr(),
                                     SBLambda([value, coefs] (size_t j, auto val)
                                              {
                                                coefs(j) += HSum(InnerProduct(value,val));
                                              }));
           }
       });
  }
  
  template class TangentialFacetFacetFE<ET_SEGM>;
  template class TangentialFacetFacetFE<ET_TRIG>;