Method GetTraceTrans

comp/l2hofespace.cpp:2116–2185

  }

  void L2HighOrderFESpace ::
  GetTraceTrans (const FESpace & tracespace, const BaseVector & in, BaseVector & out, bool avg,
                 LocalHeap & lh) const
  {
    static Timer t("GetTraceTrans"); RegionTimer reg(t);

    out = 0.0;
    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();


    // size_t ne = ma->GetNE();

    for (auto elclass_inds : table)
      {
        if (elclass_inds.Size() == 0) continue;

        ElementId ei(VOL,elclass_inds[0]);
        auto & felx = GetFE (ei, lh);
        // auto & trafo = GetMeshAccess()->GetTrafo(ei, lh);

        Matrix<> trace_op_x = GetTraceMatrix (felx);

        Matrix<> temp_x(elclass_inds.Size(), trace_op_x.Width());
        Matrix<> temp_trace(elclass_inds.Size(), trace_op_x.Height());

        ParallelForRange
          (Range(elclass_inds),
           [&] (IntRange myrange)
           {
             Array<DofId> dofs;

             // int tid = TaskManager::GetThreadId();
             {
               for (auto i : myrange)
                 {
                   tracespace.GetDofNrs(ElementId(VOL,elclass_inds[i]), dofs);
                   in.GetIndirect(dofs, temp_trace.Row(i));
                 }
             }

             {
               // temp_trace.Rows(myrange) = 0;
               // AddABt(temp_x.Rows(myrange), trace_op_x, temp_trace.Rows(myrange));
               temp_x.Rows(myrange) = temp_trace.Rows(myrange) * trace_op_x;