Method CalcDualShape

fem/tangentialfacetfe.cpp:191–202

  template <ELEMENT_TYPE ET>
  void TangentialFacetFacetFE<ET> :: CalcDualShape (const SIMD_BaseMappedIntegrationRule & bmir, BareSliceMatrix<SIMD<double>> shapes) const
  {
    constexpr int DIMSPACE = DIM+1;
    auto & mir = static_cast<const SIMD_MappedIntegrationRule<DIM,DIMSPACE>&> (bmir);
    for (size_t i = 0; i < mir.Size(); i++)
      this->CalcDualShape2 (mir[i], 
                            SBLambda ([&] (size_t j, auto val)
                            {
                              for (size_t k = 0; k < DIMSPACE; k++)
                                shapes(j*DIMSPACE+k, i) = val(k);
                            }));
  }