MoFEM::PCMGSetUpViaApproxOrdersCtx Struct Reference

Collaboration diagram for MoFEM::PCMGSetUpViaApproxOrdersCtx:

Public Member Functions
	PCMGSetUpViaApproxOrdersCtx (DM dm, Mat a, bool shell_sub_a)
virtual	~PCMGSetUpViaApproxOrdersCtx ()=default
virtual MoFEMErrorCode	getOptions ()
	get options from line command
virtual MoFEMErrorCode	createIsAtLevel (int kk, IS *is)
	Set IS for levels.
virtual MoFEMErrorCode	destroyIsAtLevel (int kk, IS *is)
	Destroy IS if internally created.
virtual MoFEMErrorCode	buildProlongationOperator (bool use_mat_a, int verb=0)
	Set up data structures for MG.

Public Attributes
DM	dM
	Distributed mesh manager.
Mat	A
	Matrix at fine level.
int	nbLevels
	number of multi-grid levels
int	coarseOrder
	approximation order of coarse level
int	orderAtLastLevel
	set maximal evaluated order
bool	shellSubA
int	verboseLevel

Detailed Description

Definition at line 92 of file PCMGSetUpViaApproxOrders.cpp.

Constructor & Destructor Documentation

PCMGSetUpViaApproxOrdersCtx()

MoFEM::PCMGSetUpViaApproxOrdersCtx::PCMGSetUpViaApproxOrdersCtx ( DM  dm, Mat  a, bool  shell_sub_a  )

inline

Definition at line 94 of file PCMGSetUpViaApproxOrders.cpp.

      : dM(dm), A(a), nbLevels(2), coarseOrder(2), orderAtLastLevel(1000),
        shellSubA(shell_sub_a), verboseLevel(0) {}

~PCMGSetUpViaApproxOrdersCtx()

virtual MoFEM::PCMGSetUpViaApproxOrdersCtx::~PCMGSetUpViaApproxOrdersCtx ( )

virtualdefault

Member Function Documentation

buildProlongationOperator()

MoFEMErrorCode MoFEM::PCMGSetUpViaApproxOrdersCtx::buildProlongationOperator ( bool  use_mat_a, int  verb = 0  )

virtual

Set up data structures for MG.

Parameters

pc	MG pre-conditioner http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCMG.html
verb	verbosity level

Returns

error code

Definition at line 590 of file PCMGSetUpViaApproxOrders.cpp.

                                                                 {
  MoFEMFunctionBegin;
  verb = verb > verboseLevel ? verb : verboseLevel;
 
  MPI_Comm comm;
  CHKERR PetscObjectGetComm((PetscObject)dM, &comm);
 
  MOFEM_LOG_C("PCMGViaApproximationOrders", Sev::inform, "set MG levels %u\n",
              nbLevels);
 
  MOFEM_LOG_CHANNEL("SYNC");
  MOFEM_LOG_TAG("SYNC", "PCMGViaApproximationOrders")
 
  std::vector<IS> is_vec(nbLevels + 1);
  std::vector<int> is_glob_size(nbLevels + 1), is_loc_size(nbLevels + 1);
 
  for (int kk = 0; kk < nbLevels; kk++) {
 
    // get indices up to up to give approximation order
    CHKERR createIsAtLevel(kk, &is_vec[kk]);
    CHKERR ISGetSize(is_vec[kk], &is_glob_size[kk]);
    CHKERR ISGetLocalSize(is_vec[kk], &is_loc_size[kk]);
 
    MOFEM_LOG_C("SYNC", Sev::inform,
                "Nb. dofs at level [ %d ] global %u local %d\n", kk,
                is_glob_size[kk], is_loc_size[kk]);
 
    // if no dofs on level kk finish here
    if (is_glob_size[kk] == 0) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "no dofs at level");
    }
  }
 
  for (int kk = 0; kk != nbLevels; kk++) {
    if (kk == nbLevels - 1 && use_mat_a) {
      CHKERR DMMGViaApproxOrdersPushBackCoarseningIS(dM, is_vec[kk], A, false,
                                                     false);
    } else {
      if (kk > 0) {
        // Not coarse level
        CHKERR DMMGViaApproxOrdersPushBackCoarseningIS(dM, is_vec[kk], A, true,
                                                       shellSubA);
      } else {
        // Coarse lave is compressed matrix allowing for factorization when
        // needed
        CHKERR DMMGViaApproxOrdersPushBackCoarseningIS(dM, is_vec[kk], A, true,
                                                       false);
      }
    }
  }
 
  for (unsigned int kk = 0; kk < is_vec.size(); kk++) {
    CHKERR destroyIsAtLevel(kk, &is_vec[kk]);
  }
 
  MOFEM_LOG_SEVERITY_SYNC(comm, Sev::inform);
  MOFEM_LOG_CHANNEL("SYNC");
 
  MoFEMFunctionReturn(0);
}

createIsAtLevel()

MoFEMErrorCode MoFEM::PCMGSetUpViaApproxOrdersCtx::createIsAtLevel ( int  kk, IS *  is  )

virtual

Set IS for levels.

Parameters

kk	level
is	pointer to IS

Returns

error code

Definition at line 556 of file PCMGSetUpViaApproxOrders.cpp.

                                                                          {
  MoFEM::Interface *m_field_ptr;
  MoFEM::ISManager *is_manager_ptr;
  MoFEMFunctionBegin;
  // if is last level, take all remaining orders dofs, if any left
  CHKERR DMoFEMGetInterfacePtr(dM, &m_field_ptr);
  CHKERR m_field_ptr->getInterface(is_manager_ptr);
  const Problem *problem_ptr;
  CHKERR DMMoFEMGetProblemPtr(dM, &problem_ptr);
  int order_at_next_level = kk + coarseOrder;
  if (kk == nbLevels - 1) {
    int first = problem_ptr->getNumeredRowDofsPtr()
                    ->get<PetscGlobalIdx_mi_tag>()
                    .lower_bound(0)
                    ->get()
                    ->getPetscGlobalDofIdx();
    CHKERR ISCreateStride(PETSC_COMM_WORLD, problem_ptr->getNbLocalDofsRow(),
                          first, 1, is);
    MoFEMFunctionReturnHot(0);
    // order_at_next_level = orderAtLastLevel;
  }
  string problem_name = problem_ptr->getName();
  CHKERR is_manager_ptr->isCreateProblemOrder(problem_name, ROW, 0,
                                              order_at_next_level, is);
  MoFEMFunctionReturn(0);
}

destroyIsAtLevel()

MoFEMErrorCode MoFEM::PCMGSetUpViaApproxOrdersCtx::destroyIsAtLevel ( int  kk, IS *  is  )

virtual

Destroy IS if internally created.

Parameters

kk	level
is	pointer to is

Returns

error code

Definition at line 583 of file PCMGSetUpViaApproxOrders.cpp.

                                                                           {
  MoFEMFunctionBegin;
  CHKERR ISDestroy(is);
  MoFEMFunctionReturn(0);
}

getOptions()

MoFEMErrorCode MoFEM::PCMGSetUpViaApproxOrdersCtx::getOptions ( )

virtual

get options from line command

Returns

error code

Definition at line 533 of file PCMGSetUpViaApproxOrders.cpp.

                                                       {
  MoFEMFunctionBeginHot;
  PetscOptionsBegin(PETSC_COMM_WORLD, "",
                           "MOFEM Multi-Grid (Orders) pre-conditioner", "none");
 
  CHKERR PetscOptionsInt("-mofem_mg_levels", "nb levels of multi-grid solver",
                         "", 2, &nbLevels, PETSC_NULLPTR);
  CHKERR PetscOptionsInt("-mofem_mg_coarse_order",
                         "approximation order of coarse level", "", 2,
                         &coarseOrder, PETSC_NULLPTR);
  CHKERR PetscOptionsInt("-mofem_mg_order_at_last_level", "order at last level",
                         "", 100, &orderAtLastLevel, PETSC_NULLPTR);
  CHKERR PetscOptionsInt("-mofem_mg_verbose", "nb levels of multi-grid solver",
                         "", 0, &verboseLevel, PETSC_NULLPTR);
  PetscBool shell_sub_a = shellSubA ? PETSC_TRUE : PETSC_FALSE;
  CHKERR PetscOptionsBool("-mofem_mg_shell_a", "use shell matrix as sub matrix",
                          "", shell_sub_a, &shell_sub_a, NULL);
  shellSubA = (shellSubA == PETSC_TRUE);
 
  PetscOptionsEnd();
  MoFEMFunctionReturnHot(0);
}

Member Data Documentation

A

Mat MoFEM::PCMGSetUpViaApproxOrdersCtx::A

Matrix at fine level.

Definition at line 133 of file PCMGSetUpViaApproxOrders.cpp.

coarseOrder

int MoFEM::PCMGSetUpViaApproxOrdersCtx::coarseOrder

approximation order of coarse level

Definition at line 136 of file PCMGSetUpViaApproxOrders.cpp.

dM

DM MoFEM::PCMGSetUpViaApproxOrdersCtx::dM

Distributed mesh manager.

Definition at line 132 of file PCMGSetUpViaApproxOrders.cpp.

nbLevels

int MoFEM::PCMGSetUpViaApproxOrdersCtx::nbLevels

number of multi-grid levels

Definition at line 135 of file PCMGSetUpViaApproxOrders.cpp.

orderAtLastLevel

int MoFEM::PCMGSetUpViaApproxOrdersCtx::orderAtLastLevel

set maximal evaluated order

Definition at line 137 of file PCMGSetUpViaApproxOrders.cpp.

shellSubA

bool MoFEM::PCMGSetUpViaApproxOrdersCtx::shellSubA

Definition at line 139 of file PCMGSetUpViaApproxOrders.cpp.

verboseLevel

int MoFEM::PCMGSetUpViaApproxOrdersCtx::verboseLevel

Definition at line 140 of file PCMGSetUpViaApproxOrders.cpp.

---

The documentation for this struct was generated from the following file:

• src/petsc/impl/PCMGSetUpViaApproxOrders.cpp