Method FinalizeUpdate

comp/fespace.cpp:518–762

  void FESpace :: FinalizeUpdate()
  {
    static Timer timer ("FESpace::FinalizeUpdate");
    /*
    static Timer timer1 ("FESpace::FinalizeUpdate 1");
    static Timer timer2 ("FESpace::FinalizeUpdate 2");
    static Timer timer3 ("FESpace::FinalizeUpdate 3");
    static Timer tcol ("FESpace::FinalizeUpdate - coloring");
    static Timer tcolbits ("FESpace::FinalizeUpdate - bitarrays");
    static Timer tcolmutex ("FESpace::FinalizeUpdate - coloring, init mutex");
    */
    if (low_order_space) low_order_space -> FinalizeUpdate();

    UpdateFreeDofs();
    
    
    UpdateParallelDofs();

    if (print)
      *testout << "coloring ... " << flush;

    if (low_order_space)
      {
	for(auto vb : {VOL, BND, BBND, BBBND})
	  element_coloring[vb] = Table<int>(low_order_space->element_coloring[vb]);
      }
    else
      {
        // tcolmutex.Start();
      Array<MyMutex> locks(GetNDof());
      // tcolmutex.Stop();
      
      for (auto vb : { VOL, BND, BBND, BBBND })
      {
        /*
        tcol.Start();
        Array<int> col(ma->GetNE(vb));
        col = -1;

        int maxcolor = 0;
        
        int basecol = 0;
        Array<unsigned int> mask(GetNDof());

        size_t cnt = 0, found = 0;
        for (ElementId el : Elements(vb)) { cnt++; (void)el; } // no warning 

        do
          {
            mask = 0;

            Array<DofId> dofs;
            // for (auto el : Elements(vb))
            for (auto el : ma->Elements(vb))
              {
                if (!DefinedOn(el)) continue;
                if (col[el.Nr()] >= 0) continue;