Method Update

comp/facetfespace.cpp:401–439

    }

    virtual void Update (const FESpace & bfes) override
    {
      auto & fes = dynamic_cast<const FacetFESpace&> (bfes);
      
      size_t oldnfacet = prol_class.Size();
      size_t nfacet = ma->GetNFacets();
      cout << IM(3) << "update prol, level = " << ma->GetNLevels() <<  ", nfacet = " << nfacet << endl;

      while (facets_on_level.Size() < ma->GetNLevels())
        facets_on_level.Append(oldnfacet);
      facets_on_level[ma->GetNLevels()-1] = nfacet;

      // *testout << "update first_dofs, first_dofs.Size() == " << first_dofs.Size() << ", levels = " << ma->GetNLevels() << endl;
      if (first_dofs.Size() < ma->GetNLevels())
        first_dofs += fes.GetFirstFacetDof();
      
      // *testout << /* IM(3) << */ "facets_on_level = " << endl << facets_on_level << endl;

      if (ma->GetNLevels() == 1)
        {
          prol_class.SetSize(nfacet);
          prol_class = 0;
          return;
        }
      
      CalcMatrices();

      prol_class.SetSize(nfacet);

      // Array<size_t> verts(4);
      for (size_t i = oldnfacet; i < nfacet; i++)
        if (auto parents = get<1>(ma->GetParentEdges(i)); parents[1] == -1)
          {
            auto myverts = ma->GetEdgePNums(i);
            auto paverts = ma->GetEdgePNums(parents[0]);
            auto minpa = Min(paverts);
            prol_class[i] = myverts.Contains(minpa) ? 0 : 1;
          }
    }
  };