Method Update

comp/preconditioner.cpp:426–472

  
  
  void MGPreconditioner :: Update ()
  {
    static Timer t("MGPreconditioner::Update"); RegionTimer reg(t);
    
    shared_ptr<BilinearForm> lo_bfa = bfa->GetLowOrderBilinearForm();

    string invtype, loinvtype = default_inversetype;
    invtype = dynamic_cast<const BaseSparseMatrix & > (bfa->GetMatrix()).SetInverseType (inversetype);
    if (lo_bfa)
      loinvtype = dynamic_cast<const BaseSparseMatrix & > (lo_bfa->GetMatrix()) .SetInverseType (inversetype);


    mgp->Update();

    if (coarse_pre)
      {
	mgp->SetCoarseGridPreconditioner (coarse_pre);
      }

    if (bfa->GetLowOrderBilinearForm()) //  || ntasks > 1) not supported anymore
      {
        static Timer t("MGPreconditioner::Update - fine precond"); RegionTimer reg(t);
        if (additional_dirichlet_constraints)
          flags.SetFlag ("additional_dirichlet_constraints", std::any(*additional_dirichlet_constraints));
        auto fine_smoother = make_shared<BlockSmoother> (*bfa->GetMeshAccess(), *bfa, flags);
        fine_smoother -> SetAdditionalDirichletConstraints(additional_dirichlet_constraints);
        GetMemoryTracer().Track(*fine_smoother, "FineSmoother");
        tlp = make_shared<TwoLevelMatrix> (&bfa->GetMatrix(),
                                           &*mgp,
                                           fine_smoother,
                                           bfa->GetMeshAccess()->GetNLevels()-1);
        tlp -> SetSmoothingSteps (finesmoothingsteps);
        if (bfa->GetFESpace()->LowOrderEmbedding())
          tlp->SetEmbedding (bfa->GetFESpace()->LowOrderEmbedding());
	tlp -> Update();
      }
    else
      tlp = nullptr;

    if (timing) Timing();
    if (test) Test();
    if (mgtest) MgTest();

    dynamic_cast<const BaseSparseMatrix & > (bfa->GetMatrix()).SetInverseType ( invtype );
    if (lo_bfa)
      dynamic_cast<const BaseSparseMatrix & > (lo_bfa->GetMatrix()) .SetInverseType ( loinvtype );
  }
  
  
  void MGPreconditioner :: CleanUpLevel ()