Function IterateElements

comp/fespace.cpp:2123–2198

  }

  void IterateElements (const FESpace & fes, 
			VorB vb, 
			LocalHeap & clh, 
			const function<void(FESpace::Element,LocalHeap&)> & func)
  {
    static mutex copyex_mutex;
    const Table<int> & element_coloring = fes.ElementColoring(vb);
    
    if (auto *task_manager = GetTaskManager(); task_manager)
      {
        for (FlatArray<int> els_of_col : element_coloring)
          {
            SharedLoop2 sl(els_of_col.Range());

            task_manager -> CreateJob
              ( [&] (const TaskInfo & ti) 
                {
                  LocalHeap lh = clh.Split(ti.thread_nr, ti.nthreads);
                  ArrayMem<int,100> temp_dnums;

                  for (int mynr : sl)
                    {
                      HeapReset hr(lh);
                      FESpace::Element el(fes, 
                                          ElementId (vb, els_of_col[mynr]), 
                                          temp_dnums, lh);
                      
                      func (std::move(el), lh);
                    }

                  ProgressOutput::SumUpLocal();
                } );
          }
        return;
      }


    Exception * ex = nullptr;

    for (FlatArray<int> els_of_col : element_coloring)
      ParallelForRange( IntRange(els_of_col.Size()), [&] ( IntRange r )
      {
        LocalHeap lh = clh.Split();
        Array<int> temp_dnums;

        // lh.ClearValues();
        
        for (int i : r)
          {
            try
              {
                HeapReset hr(lh);
                FESpace::Element el(fes, ElementId (vb, els_of_col[i]), temp_dnums, lh);
                
                func (std::move(el), lh);
              }
	    
            catch (const Exception & e)