Method FinalizeLevel

comp/preconditioner.cpp:705–793

  }
  
  void LocalPreconditioner :: FinalizeLevel (const BaseMatrix * mat)
  {
      cout << IM(3) << "Update Local Preconditioner" << flush;
      timestamp = bfa->GetTimeStamp();

      if (flags.StringFlagDefined("blocktype") || flags.StringListFlagDefined("blocktype"))
        {
          if (additional_dirichlet_constraints)
            flags.SetFlag ("additional_dirichlet_constraints", std::any(*additional_dirichlet_constraints));
          
          auto blocks = bfa->GetFESpace()->CreateSmoothingBlocks(flags);
          shared_ptr<BaseMatrix> mat = bfa->GetMatrixPtr();          
          auto spmat = dynamic_pointer_cast<BaseSparseMatrix> (mat);
          
          if (GaussSeidel)
            jacobi = make_shared<SymmetricBlockGaussSeidelPrecond>(spmat, blocks);
          else
            jacobi = spmat -> CreateBlockJacobiPrecond(blocks);

            /*
            jacobi = dynamic_cast<const BaseSparseMatrix&> (bfa->GetMatrix())
              .CreateBlockJacobiPrecond(blocks);
            */
          return;
        }


      
      int blocktype = int (flags.GetNumFlag ( "blocktype", -1));
      
      // if (MyMPI_GetNTasks() != 1) return;
      bool parallel = (this->on_proc == -1);
      /*
        if ( !parallel && id != this->on_proc )
        {
	jacobi = 0; 
	return;
        }
      */

//     BaseBlockJacobiPrecond::COARSE_TYPE bbct;
//     switch ( ct  )
//       {
//       case NO_COARSE:
// 	bbct = NO_COARSE; break;
//       case DIRECT_COARSE:
// 	bbct = DIRECT_COARSE; break;
//       }

      if (blockcreator)
        {
          shared_ptr<Table<int>> blocks = blockcreator(*bfa->GetFESpace());
          // cout << "created blocks: " << *blocks << endl;
          jacobi = dynamic_cast<const BaseSparseMatrix&> (bfa->GetMatrix())
            .CreateBlockJacobiPrecond(blocks); // , 0, parallel, bfa->GetFESpace()->GetFreeDofs());
          return;
        }