Method AddMatrixGF

comp/bilinearform.cpp:5970–8118

  template <class SCAL>
  void S_BilinearForm<SCAL> :: AddMatrixGF (SCAL val,
                                            const BaseVector & x,
                                            BaseVector & y,
                                            bool transpose,
                                            LocalHeap & lh) const
  {
    static Timer t("BilinearForm::Apply - geomfree");
    static Timer tgetx("BilinearForm::Apply - get x");    
    static Timer tx("BilinearForm::Apply - transform x");
    static Timer ty("BilinearForm::Apply - transform y");
    static Timer taddy("BilinearForm::Apply - add y");        
    static Timer tgf("BilinearForm::Apply - geomfree gridfunction");
    static Timer tgfmult("BilinearForm::Apply - geomfree gridfunction - mult");
    static Timer tm("BilinearForm::Apply - geomfree mult");
    static Timer teval("BilinearForm::Apply - evaluate");
    RegionTimer reg(t);

    // classify:
    auto fesx = GetTrialSpace();
    auto fesy = GetTestSpace();
    if (transpose) Swap (fesx, fesy);
    auto ma = GetMeshAccess();

    /*
    Array<short> classnr(ma->GetNE());
    ma->IterateElements
      (VOL, lh, [&] (auto el, LocalHeap & llh)
       {
         classnr[el.Nr()] = 
           SwitchET<ET_TRIG,ET_TET>
           (el.GetType(),
            [el] (auto et) { return ET_trait<et.ElementType()>::GetClassNr(el.Vertices()); });
       });
    
    TableCreator<size_t> creator;
    for ( ; !creator.Done(); creator++)
      for (auto i : Range(classnr))
        creator.Add (classnr[i], i);
    Table<size_t> table = creator.MoveTable();
    */
    
    const Table<size_t> & table = ma->GetElementsOfClass();
    // auto dof_tablex = fesx->CreateDofTable(VOL);



#ifdef OLDPARALLEL
    
    for (auto elclass_inds : table)
      {
        if (elclass_inds.Size() == 0) continue;
        HeapReset hr(lh);
        
        ElementId ei(VOL,elclass_inds[0]);
        auto & felx = GetTrialSpace()->GetFE (ei, lh);
        auto & fely = GetTestSpace()->GetFE (ei, lh);
        auto & trafo = GetTrialSpace()->GetMeshAccess()->GetTrafo(ei, lh);