Method AddMatrixTrans

comp/bilinearform.cpp:5045–5137

  template <class SCAL>
  void S_BilinearForm<SCAL> :: AddMatrixTrans (double val,
                                               const BaseVector & x,
                                               BaseVector & y, LocalHeap & clh) const
  {
    if (geom_free_parts.Size())
      AddMatrixGF(val, x, y, true, clh);
    if (geom_free_parts.Size() == parts.Size()) return;

    if (!MixedSpaces())
      {
        for (auto vb : { VOL, BND, BBND, BBBND } )
          if (VB_parts[vb].Size())
            {
              IterateElements 
                (*fespace, vb, clh, 
                 [&] (FESpace::Element el, LocalHeap & lh)
                 {
                   // RegionTimer reg (timer_loop);                   
                   auto & fel = el.GetFE();
                   auto & trafo = el.GetTrafo();
                   auto dnums = el.GetDofs();
                   
                   FlatVector<SCAL> elvecx (dnums.Size() * fespace->GetDimension(), lh);
                   FlatVector<SCAL> elvecy (dnums.Size() * fespace->GetDimension(), lh);
                   
                   x.GetIndirect (dnums, elvecx);
                   this->fespace->TransformVec (el, elvecx, TRANSFORM_SOL);
                   
                   for (auto & bfi : VB_parts[vb])
                     {
                       if (!bfi->DefinedOn (el.GetIndex())) continue;
                       if (!bfi->DefinedOnElement (el.Nr())) continue;

                       auto & mapped_trafo = trafo.AddDeformation(bfi->GetDeformation().get(), lh);
                       
                       bfi->ApplyElementMatrixTrans (fel, mapped_trafo, elvecx, elvecy, 0, lh);

                       this->fespace->TransformVec (el, elvecy, TRANSFORM_RHS);
                       
                       elvecy *= val;
                       y.AddIndirect (dnums, elvecy, fespace->HasAtomicDofs());  // coloring
                     }
                 });
            } 
      }
    else // MixedSpaces
      {
        static Timer timer ("Apply Matrix Trans - mixed"); RegionTimer reg(timer);

        for (auto vb : { VOL, BND, BBND } )
          if (VB_parts[vb].Size())
            {        
              IterateElements 
                (*fespace, vb, clh, 
                 [&] (ElementId ei, LocalHeap & lh)
           
                 {
                   if (!fespace->DefinedOn (ei)) return;