v0.14.0
Functions
Distributed mesh manager
MoFEM

Implementation of PETSc DM, managing interactions between mesh data structures and vectors and matrices. More...

Collaboration diagram for Distributed mesh manager:

Functions

PetscErrorCode MoFEM::DMRegister_MoFEM (const char sname[])
 Register MoFEM problem. More...
 
PetscErrorCode MoFEM::DMMoFEMCreateMoFEM (DM dm, MoFEM::Interface *m_field_ptr, const char problem_name[], const MoFEM::BitRefLevel bit_level, const MoFEM::BitRefLevel bit_mask=MoFEM::BitRefLevel().set())
 Must be called by user to set MoFEM data structures. More...
 
PetscErrorCode MoFEM::DMMoFEMCreateSubDM (DM subdm, DM dm, const char problem_name[])
 Must be called by user to set Sub DM MoFEM data structures. More...
 
PetscErrorCode MoFEM::DMoFEMGetInterfacePtr (DM dm, MoFEM::Interface **m_field_ptr)
 Get pointer to MoFEM::Interface. More...
 
PetscErrorCode MoFEM::DMMoFEMGetProblemPtr (DM dm, const MoFEM::Problem **problem_ptr)
 Get pointer to problem data structure. More...
 
PetscErrorCode MoFEM::DMMoFEMSetSquareProblem (DM dm, PetscBool square_problem)
 set squared problem More...
 
PetscErrorCode MoFEM::DMMoFEMGetSquareProblem (DM dm, PetscBool *square_problem)
 get squared problem More...
 
PetscErrorCode MoFEM::DMMoFEMResolveSharedFiniteElements (DM dm, std::string fe_name)
 Resolve shared entities. More...
 
PetscErrorCode MoFEM::DMMoFEMGetProblemFiniteElementLayout (DM dm, std::string fe_name, PetscLayout *layout)
 Get finite elements layout in the problem. More...
 
PetscErrorCode MoFEM::DMMoFEMAddElement (DM dm, std::string fe_name)
 add element to dm More...
 
PetscErrorCode MoFEM::DMMoFEMAddElement (DM dm, std::vector< std::string > fe_name)
 add element to dm More...
 
PetscErrorCode MoFEM::DMMoFEMUnSetElement (DM dm, std::string fe_name)
 unset element from dm More...
 
PetscErrorCode MoFEM::DMoFEMMeshToLocalVector (DM dm, Vec l, InsertMode mode, ScatterMode scatter_mode)
 set local (or ghosted) vector values on mesh for partition only More...
 
PetscErrorCode MoFEM::DMoFEMMeshToGlobalVector (DM dm, Vec g, InsertMode mode, ScatterMode scatter_mode)
 set ghosted vector values on all existing mesh entities More...
 
PetscErrorCode MoFEM::DMoFEMPreProcessFiniteElements (DM dm, MoFEM::FEMethod *method)
 execute finite element method for each element in dm (problem) More...
 
PetscErrorCode MoFEM::DMoFEMPostProcessFiniteElements (DM dm, MoFEM::FEMethod *method)
 execute finite element method for each element in dm (problem) More...
 
PetscErrorCode MoFEM::DMoFEMLoopFiniteElementsUpAndLowRank (DM dm, const char fe_name[], MoFEM::FEMethod *method, int low_rank, int up_rank, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
 Executes FEMethod for finite elements in DM. More...
 
PetscErrorCode MoFEM::DMoFEMLoopFiniteElementsUpAndLowRank (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, int low_rank, int up_rank, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
 Executes FEMethod for finite elements in DM. More...
 
PetscErrorCode MoFEM::DMoFEMLoopFiniteElements (DM dm, const char fe_name[], MoFEM::FEMethod *method, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
 Executes FEMethod for finite elements in DM. More...
 
PetscErrorCode MoFEM::DMoFEMLoopFiniteElements (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
 Executes FEMethod for finite elements in DM. More...
 
PetscErrorCode MoFEM::DMoFEMLoopDofs (DM dm, const char field_name[], MoFEM::DofMethod *method)
 execute method for dofs on field in problem More...
 
PetscErrorCode MoFEM::DMMoFEMKSPSetComputeRHS (DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
 set KSP right hand side evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMKSPSetComputeRHS (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
 set KSP right hand side evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMKSPSetComputeOperators (DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
 Set KSP operators and push mofem finite element methods. More...
 
PetscErrorCode MoFEM::DMMoFEMKSPSetComputeOperators (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
 Set KSP operators and push mofem finite element methods. More...
 
PetscErrorCode MoFEM::DMMoFEMSNESSetFunction (DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
 set SNES residual evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMSNESSetFunction (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
 set SNES residual evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMSNESSetJacobian (DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
 set SNES Jacobian evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMSNESSetJacobian (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
 set SNES Jacobian evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMTSSetIFunction (DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
 set TS implicit function evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMTSSetIFunction (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
 set TS implicit function evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMTSSetIJacobian (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
 set TS Jacobian evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMTSSetIJacobian (DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
 set TS Jacobian evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMTSSetI2Function (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
 set TS implicit function evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMTSSetI2Function (DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
 set TS implicit function evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMTSSetI2Jacobian (DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
 set TS Jacobian evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMTSSetI2Jacobian (DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
 set TS Jacobian evaluation function More...
 
PetscErrorCode MoFEM::DMMoFEMGetKspCtx (DM dm, MoFEM::KspCtx **ksp_ctx)
 get MoFEM::KspCtx data structure More...
 
PetscErrorCode MoFEM::DMMoFEMGetKspCtx (DM dm, const boost::shared_ptr< MoFEM::KspCtx > &ksp_ctx)
 get MoFEM::KspCtx data structure More...
 
PetscErrorCode MoFEM::DMMoFEMSetKspCtx (DM dm, boost::shared_ptr< MoFEM::KspCtx > ksp_ctx)
 set MoFEM::KspCtx data structure More...
 
PetscErrorCode MoFEM::DMMoFEMGetSnesCtx (DM dm, MoFEM::SnesCtx **snes_ctx)
 get MoFEM::SnesCtx data structure More...
 
PetscErrorCode MoFEM::DMMoFEMGetSnesCtx (DM dm, const boost::shared_ptr< MoFEM::SnesCtx > &snes_ctx)
 get MoFEM::SnesCtx data structure More...
 
PetscErrorCode MoFEM::DMMoFEMSetSnesCtx (DM dm, boost::shared_ptr< MoFEM::SnesCtx > snes_ctx)
 Set MoFEM::SnesCtx data structure. More...
 
PetscErrorCode MoFEM::DMMoFEMGetTsCtx (DM dm, MoFEM::TsCtx **ts_ctx)
 get MoFEM::TsCtx data structure More...
 
PetscErrorCode MoFEM::DMMoFEMGetTsCtx (DM dm, const boost::shared_ptr< MoFEM::TsCtx > &ts_ctx)
 get MoFEM::TsCtx data structure More...
 
PetscErrorCode MoFEM::DMMoFEMSetTsCtx (DM dm, boost::shared_ptr< MoFEM::TsCtx > ts_ctx)
 Set MoFEM::TsCtx data structure. More...
 
PetscErrorCode MoFEM::DMMoFEMSetIsPartitioned (DM dm, PetscBool is_partitioned)
 
PetscErrorCode MoFEM::DMMoFEMGetIsPartitioned (DM dm, PetscBool *is_partitioned)
 
PetscErrorCode MoFEM::DMSetOperators_MoFEM (DM dm)
 Set operators for MoFEM dm. More...
 
PetscErrorCode MoFEM::DMCreate_MoFEM (DM dm)
 Create dm data structure with MoFEM data structure. More...
 
PetscErrorCode MoFEM::DMDestroy_MoFEM (DM dm)
 Destroys dm with MoFEM data structure. More...
 
PetscErrorCode MoFEM::DMCreateGlobalVector_MoFEM (DM dm, Vec *g)
 DMShellSetCreateGlobalVector. More...
 
PetscErrorCode MoFEM::DMCreateGlobalVector_MoFEM (DM dm, SmartPetscObj< Vec > &g_ptr)
 DMShellSetCreateGlobalVector. More...
 
PetscErrorCode MoFEM::DMCreateLocalVector_MoFEM (DM dm, Vec *l)
 DMShellSetCreateLocalVector. More...
 
PetscErrorCode MoFEM::DMCreateMatrix_MoFEM (DM dm, Mat *M)
 
PetscErrorCode MoFEM::DMCreateMatrix_MoFEM (DM dm, SmartPetscObj< Mat > &M)
 
PetscErrorCode MoFEM::DMSetFromOptions_MoFEM (DM dm)
 
PetscErrorCode MoFEM::DMSetUp_MoFEM (DM dm)
 
PetscErrorCode MoFEM::DMSubDMSetUp_MoFEM (DM subdm)
 
PetscErrorCode MoFEM::DMMoFEMAddSubFieldRow (DM dm, const char field_name[])
 
PetscErrorCode MoFEM::DMMoFEMAddSubFieldRow (DM dm, const char field_name[], boost::shared_ptr< Range > r_ptr)
 Add field to sub dm problem on rows. More...
 
PetscErrorCode MoFEM::DMMoFEMAddSubFieldCol (DM dm, const char field_name[])
 
PetscErrorCode MoFEM::DMMoFEMGetIsSubDM (DM dm, PetscBool *is_sub_dm)
 
PetscErrorCode MoFEM::DMMoFEMAddRowCompositeProblem (DM dm, const char prb_name[])
 Add problem to composite DM on row. More...
 
PetscErrorCode MoFEM::DMMoFEMAddColCompositeProblem (DM dm, const char prb_name[])
 Add problem to composite DM on col. More...
 
PetscErrorCode MoFEM::DMMoFEMGetIsCompDM (DM dm, PetscBool *is_comp_dm)
 Get if this DM is composite DM. More...
 
PetscErrorCode MoFEM::DMGlobalToLocalBegin_MoFEM (DM dm, Vec, InsertMode, Vec)
 
PetscErrorCode MoFEM::DMGlobalToLocalEnd_MoFEM (DM dm, Vec, InsertMode, Vec)
 
PetscErrorCode MoFEM::DMLocalToGlobalBegin_MoFEM (DM, Vec, InsertMode, Vec)
 
PetscErrorCode MoFEM::DMLocalToGlobalEnd_MoFEM (DM, Vec, InsertMode, Vec)
 
PetscErrorCode MoFEM::DMCreateFieldIS_MoFEM (DM dm, PetscInt *numFields, char ***fieldNames, IS **fields)
 
PetscErrorCode MoFEM::DMMoFEMGetFieldIS (DM dm, RowColData rc, const char field_name[], IS *is)
 get field is in the problem More...
 
MoFEMErrorCode MoFEM::DMMoFEMCreateBlockMat (DM dm, Mat *mat)
 Create block matrix. More...
 
MoFEMErrorCode MoFEM::DMMoFEMCreateBlockMat (DM dm, SmartPetscObj< Mat > &mat)
 Create block matrix. More...
 
auto MoFEM::createDMMatrix (DM dm)
 Get smart matrix from DM. More...
 
auto MoFEM::createDMVector (DM dm)
 Get smart vector from DM. More...
 
MoFEMErrorCode MoFEM::DMMGViaApproxOrdersPushBackCoarseningIS (DM, IS is, Mat A, bool create_sub_matrix, bool shell_sub_a)
 Push back coarsening level to MG via approximation orders. More...
 
MoFEMErrorCode MoFEM::DMMGViaApproxOrdersPopBackCoarseningIS (DM)
 Pop IS form MG via approximation orders. More...
 
MoFEMErrorCode MoFEM::DMMGViaApproxOrdersClearCoarseningIS (DM)
 Clear approximation orders. More...
 
MoFEMErrorCode MoFEM::DMRegister_MGViaApproxOrders (const char sname[])
 Register DM for Multi-Grid via approximation orders. More...
 
MoFEMErrorCode MoFEM::DMCreateMatrix_MGViaApproxOrders (DM dm, Mat *M)
 Create matrix for Multi-Grid via approximation orders. More...
 
MoFEMErrorCode MoFEM::DMCoarsen_MGViaApproxOrders (DM dm, MPI_Comm comm, DM *dmc)
 Coarsen DM. More...
 
MoFEMErrorCode DMMGViaApproxOrdersPushBackCoarseningIS (DM, IS is, Mat A, Mat *subA, bool create_sub_matrix, bool shell_sub_a)
 Push back coarsening level to MG via approximation orders. More...
 
MoFEMErrorCode DMMGViaApproxOrdersPopBackCoarseningIS (DM)
 Pop is form MG via approximation orders. More...
 
MoFEMErrorCode DMMGViaApproxOrdersClearCoarseningIS (DM)
 Clear approximation orders. More...
 
MoFEMErrorCode DMRegister_MGViaApproxOrders (const char sname[])
 Register DM for Multi-Grid via approximation orders. More...
 
MoFEMErrorCode DMCreateMatrix_MGViaApproxOrders (DM dm, Mat *M)
 Create matrix for Multi-Grid via approximation orders. More...
 
MoFEMErrorCode DMCoarsen_MGViaApproxOrders (DM dm, MPI_Comm comm, DM *dmc)
 Coarsen DM. More...
 

Detailed Description

Implementation of PETSc DM, managing interactions between mesh data structures and vectors and matrices.

DM objects are used to manage communication between the algebraic structures in PETSc (Vec and Mat) and mesh data structures in PDE-based (or other) simulations.

DM is abstract interface, here is it particular implementation for MoFEM code.

Function Documentation

◆ createDMMatrix()

auto MoFEM::createDMMatrix ( DM  dm)
inline

Get smart matrix from DM.

Examples
free_surface.cpp, level_set.cpp, navier_stokes.cpp, plastic.cpp, schur_test_diag_mat.cpp, and simple_interface.cpp.

Definition at line 1056 of file DMMoFEM.hpp.

1056  {
1057  SmartPetscObj<Mat> a;
1058  ierr = DMCreateMatrix_MoFEM(dm, a);
1059  CHKERRABORT(getCommFromPetscObject(reinterpret_cast<PetscObject>(dm)), ierr);
1060  return a;
1061 };

◆ createDMVector()

auto MoFEM::createDMVector ( DM  dm)
inline

Get smart vector from DM.

Examples
adolc_plasticity.cpp, approx_sphere.cpp, child_and_parent.cpp, dg_projection.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, free_surface.cpp, hanging_node_approx.cpp, hcurl_check_approx_in_2d.cpp, heat_equation.cpp, heat_method.cpp, helmholtz.cpp, higher_derivatives.cpp, level_set.cpp, mixed_poisson.cpp, navier_stokes.cpp, nonlinear_elastic.cpp, phase.cpp, photon_diffusion.cpp, plastic.cpp, PlasticOpsMonitor.hpp, plate.cpp, poisson_2d_dis_galerkin.cpp, poisson_2d_homogeneous.cpp, scalar_check_approximation.cpp, schur_test_diag_mat.cpp, seepage.cpp, shallow_wave.cpp, simple_interface.cpp, tensor_divergence_operator.cpp, test_broken_space.cpp, thermo_elastic.cpp, and wave_equation.cpp.

Definition at line 1099 of file DMMoFEM.hpp.

1099  {
1100  SmartPetscObj<Vec> v;
1101  ierr = DMCreateGlobalVector_MoFEM(dm, v);
1102  CHKERRABORT(getCommFromPetscObject(reinterpret_cast<PetscObject>(dm)), ierr);
1103  return v;
1104 };

◆ DMCoarsen_MGViaApproxOrders() [1/2]

MoFEMErrorCode DMCoarsen_MGViaApproxOrders ( DM  dm,
MPI_Comm  comm,
DM *  dmc 
)

Coarsen DM.

Not used directly by user

Parameters
dmDistributed mesh data structure
commCommunicator
dmcCoarse distributed mesh data structure
Returns
Error code

Definition at line 376 of file PCMGSetUpViaApproxOrders.cpp.

376  {
377  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
378  MoFEMFunctionBegin;
379  GET_DM_FIELD(dm);
380  CHKERR PetscObjectGetComm((PetscObject)dm, &comm);
381  CHKERR DMCreate(comm, dmc);
382  (*dmc)->data = dm->data;
383  DMType type;
384  CHKERR DMGetType(dm, &type);
385  CHKERR DMSetType(*dmc, type);
386  dm_field->coarsenDM = SmartPetscObj<DM>(*dmc, false);
387  PetscInfo1(dm, "Coarsen DMMGViaApproxOrders leveldown = %d\n", dm->leveldown);
388  MoFEMFunctionReturn(0);
389 }

◆ DMCoarsen_MGViaApproxOrders() [2/2]

MoFEMErrorCode MoFEM::DMCoarsen_MGViaApproxOrders ( DM  dm,
MPI_Comm  comm,
DM *  dmc 
)

Coarsen DM.

Not used directly by user

Parameters
dmDistributed mesh data structure
commCommunicator
dmcCoarse distributed mesh data structure
Returns
Error code

Definition at line 376 of file PCMGSetUpViaApproxOrders.cpp.

376  {
377  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
378  MoFEMFunctionBegin;
379  GET_DM_FIELD(dm);
380  CHKERR PetscObjectGetComm((PetscObject)dm, &comm);
381  CHKERR DMCreate(comm, dmc);
382  (*dmc)->data = dm->data;
383  DMType type;
384  CHKERR DMGetType(dm, &type);
385  CHKERR DMSetType(*dmc, type);
386  dm_field->coarsenDM = SmartPetscObj<DM>(*dmc, false);
387  PetscInfo1(dm, "Coarsen DMMGViaApproxOrders leveldown = %d\n", dm->leveldown);
388  MoFEMFunctionReturn(0);
389 }

◆ DMCreate_MoFEM()

PetscErrorCode MoFEM::DMCreate_MoFEM ( DM  dm)

Create dm data structure with MoFEM data structure.

Definition at line 71 of file DMMoFEM.cpp.

71  {
72  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
73  MoFEMFunctionBegin;
74  dm->data = new DMCtxImpl();
75  CHKERR DMSetOperators_MoFEM(dm);
76  MoFEMFunctionReturn(0);
77 }

◆ DMCreateFieldIS_MoFEM()

PetscErrorCode MoFEM::DMCreateFieldIS_MoFEM ( DM  dm,
PetscInt *  numFields,
char ***  fieldNames,
IS **  fields 
)

Creates a set of IS objects with the global indices of dofs for each field

Parameters
dmThe number of fields (or NULL if not requested)

Output:

Parameters
numFieldsThe number of fields (or NULL if not requested)
fieldNamesThe name for each field (or NULL if not requested)
fieldsThe global indices for each field (or NULL if not requested)
Returns
error code
Note
The user is responsible for freeing all requested arrays. In particular, every entry of names should be freed with PetscFree(), every entry of fields should be destroyed with ISDestroy(), and both arrays should be freed with PetscFree().

Definition at line 1469 of file DMMoFEM.cpp.

1470  {
1471  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1472  MoFEMFunctionBegin;
1473 
1474  if (numFields) {
1475  *numFields = 0;
1476  }
1477  if (fieldNames) {
1478  *fieldNames = NULL;
1479  }
1480  if (fields) {
1481  *fields = NULL;
1482  }
1483 
1484  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1485  auto fields_ptr = dm_field->mField_ptr->get_fields();
1486  Field_multiIndex::iterator fit, hi_fit;
1487  fit = fields_ptr->begin();
1488  hi_fit = fields_ptr->end();
1489  *numFields = std::distance(fit, hi_fit);
1490 
1491  if (fieldNames) {
1492  CHKERR PetscMalloc1(*numFields, fieldNames);
1493  }
1494  if (fields) {
1495  CHKERR PetscMalloc1(*numFields, fields);
1496  }
1497 
1498  for (int f = 0; fit != hi_fit; fit++, f++) {
1499  if (fieldNames) {
1500  CHKERR PetscStrallocpy(fit->get()->getName().c_str(),
1501  (char **)&(*fieldNames)[f]);
1502  }
1503  if (fields) {
1504  CHKERR dm_field->mField_ptr->getInterface<ISManager>()
1505  ->isCreateProblemFieldAndRank(
1506  dm_field->problemPtr->getName(), ROW, fit->get()->getName(), 0,
1507  fit->get()->getNbOfCoeffs(), &(*fields)[f]);
1508  }
1509  }
1510 
1511  MoFEMFunctionReturn(0);
1512 }

◆ DMCreateGlobalVector_MoFEM() [1/2]

PetscErrorCode MoFEM::DMCreateGlobalVector_MoFEM ( DM  dm,
SmartPetscObj< Vec > &  g_ptr 
)

DMShellSetCreateGlobalVector.

sets the routine to create a global vector associated with the shell DM

Definition at line 1177 of file DMMoFEM.cpp.

1177  {
1178  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1179  MoFEMFunctionBegin;
1180  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1181  CHKERR dm_field->mField_ptr->getInterface<VecManager>()->vecCreateGhost(
1182  dm_field->problemName, COL, g_ptr);
1183  CHKERR VecSetDM(g_ptr, dm);
1184  MoFEMFunctionReturn(0);
1185 }

◆ DMCreateGlobalVector_MoFEM() [2/2]

PetscErrorCode MoFEM::DMCreateGlobalVector_MoFEM ( DM  dm,
Vec *  g 
)

DMShellSetCreateGlobalVector.

sets the routine to create a global vector associated with the shell DM

Examples
elasticity.cpp, elasticity_mixed_formulation.cpp, ep.cpp, reaction_diffusion.cpp, and test_cache_on_entities.cpp.

Definition at line 1167 of file DMMoFEM.cpp.

1167  {
1168  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1169  MoFEMFunctionBegin;
1170  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1171  CHKERR dm_field->mField_ptr->getInterface<VecManager>()->vecCreateGhost(
1172  dm_field->problemName, COL, g);
1173  CHKERR VecSetDM(*g, dm);
1174  MoFEMFunctionReturn(0);
1175 }

◆ DMCreateLocalVector_MoFEM()

PetscErrorCode MoFEM::DMCreateLocalVector_MoFEM ( DM  dm,
Vec *  l 
)

DMShellSetCreateLocalVector.

sets the routine to create a local vector associated with the shell DM

Definition at line 1187 of file DMMoFEM.cpp.

1187  {
1188  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1189  MoFEMFunctionBegin;
1190  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1191  CHKERR dm_field->mField_ptr->getInterface<VecManager>()->vecCreateSeq(
1192  dm_field->problemName, COL, l);
1193  CHKERR VecSetDM(*l, dm);
1194  MoFEMFunctionReturn(0);
1195 }

◆ DMCreateMatrix_MGViaApproxOrders() [1/2]

MoFEMErrorCode DMCreateMatrix_MGViaApproxOrders ( DM  dm,
Mat *  M 
)

Create matrix for Multi-Grid via approximation orders.

Not used directly by user

Parameters
dmDistributed mesh data structure
MMatrix
Returns
Error code

Definition at line 344 of file PCMGSetUpViaApproxOrders.cpp.

344  {
345 
346  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
347  MoFEMFunctionBegin;
348  GET_DM_FIELD(dm);
349 
350  int leveldown = dm->leveldown;
351 
352  if (dm_field->kspOperators.empty()) {
353  CHKERR DMCreateMatrix_MoFEM(dm, M);
354  } else {
355  MPI_Comm comm;
356  CHKERR PetscObjectGetComm((PetscObject)dm, &comm);
357  if (dm_field->kspOperators.empty()) {
358  SETERRQ(comm, MOFEM_DATA_INCONSISTENCY,
359  "data inconsistency, operator can not be set");
360  }
361  if (static_cast<int>(dm_field->kspOperators.size()) < leveldown) {
362  SETERRQ(comm, MOFEM_DATA_INCONSISTENCY,
363  "data inconsistency, no IS for that level");
364  }
365  *M = dm_field->kspOperators[dm_field->kspOperators.size() - 1 - leveldown];
366  CHKERR PetscObjectReference((PetscObject)*M);
367  CHKERR MatSetDM(*M, dm);
368  }
369 
370  PetscInfo1(dm, "Create Matrix DMMGViaApproxOrders leveldown = %d\n",
371  leveldown);
372 
373  MoFEMFunctionReturn(0);
374 }

◆ DMCreateMatrix_MGViaApproxOrders() [2/2]

MoFEMErrorCode MoFEM::DMCreateMatrix_MGViaApproxOrders ( DM  dm,
Mat *  M 
)

Create matrix for Multi-Grid via approximation orders.

Not used directly by user

Parameters
dmDistributed mesh data structure
MMatrix
Returns
Error code

Definition at line 344 of file PCMGSetUpViaApproxOrders.cpp.

344  {
345 
346  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
347  MoFEMFunctionBegin;
348  GET_DM_FIELD(dm);
349 
350  int leveldown = dm->leveldown;
351 
352  if (dm_field->kspOperators.empty()) {
353  CHKERR DMCreateMatrix_MoFEM(dm, M);
354  } else {
355  MPI_Comm comm;
356  CHKERR PetscObjectGetComm((PetscObject)dm, &comm);
357  if (dm_field->kspOperators.empty()) {
358  SETERRQ(comm, MOFEM_DATA_INCONSISTENCY,
359  "data inconsistency, operator can not be set");
360  }
361  if (static_cast<int>(dm_field->kspOperators.size()) < leveldown) {
362  SETERRQ(comm, MOFEM_DATA_INCONSISTENCY,
363  "data inconsistency, no IS for that level");
364  }
365  *M = dm_field->kspOperators[dm_field->kspOperators.size() - 1 - leveldown];
366  CHKERR PetscObjectReference((PetscObject)*M);
367  CHKERR MatSetDM(*M, dm);
368  }
369 
370  PetscInfo1(dm, "Create Matrix DMMGViaApproxOrders leveldown = %d\n",
371  leveldown);
372 
373  MoFEMFunctionReturn(0);
374 }

◆ DMCreateMatrix_MoFEM() [1/2]

PetscErrorCode MoFEM::DMCreateMatrix_MoFEM ( DM  dm,
Mat *  M 
)

DMShellSetCreateMatrix

sets the routine to create a matrix associated with the shell DM

Examples
dm_build_partitioned_mesh.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, ep.cpp, and reaction_diffusion.cpp.

Definition at line 1197 of file DMMoFEM.cpp.

1197  {
1198  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1199  MoFEMFunctionBegin;
1200  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1201 
1202  if (strcmp(dm->mattype, MATSHELL) == 0) {
1203 
1204  if (dm_field->blocMatDataPtr) {
1205  CHKERR DMMoFEMCreateBlockMat(dm, M);
1206  MoFEMFunctionReturnHot(0);
1207  } else {
1208  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
1209  "Matrix shell data not set, or matrix type not implemented");
1210  }
1211 
1212  } else if (strcmp(dm->mattype, MATMPIAIJ) == 0) {
1213  CHKERR dm_field->mField_ptr->getInterface<MatrixManager>()
1214  ->createMPIAIJWithArrays<PetscGlobalIdx_mi_tag>(dm_field->problemName,
1215  M);
1216  } else if (strcmp(dm->mattype, MATAIJ) == 0) {
1217  CHKERR dm_field->mField_ptr->getInterface<MatrixManager>()
1218  ->createSeqAIJWithArrays<PetscLocalIdx_mi_tag>(dm_field->problemName,
1219  M);
1220  } else if (strcmp(dm->mattype, MATAIJCUSPARSE) == 0) {
1221  CHKERR dm_field->mField_ptr->getInterface<MatrixManager>()
1222  ->createMPIAIJCUSPARSEWithArrays<PetscGlobalIdx_mi_tag>(
1223  dm_field->problemName, M);
1224  } else if (strcmp(dm->mattype, MATSEQAIJCUSPARSE) == 0) {
1225  CHKERR dm_field->mField_ptr->getInterface<MatrixManager>()
1226  ->createSeqAIJCUSPARSEWithArrays<PetscLocalIdx_mi_tag>(
1227  dm_field->problemName, M);
1228  } else {
1229  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
1230  "Matrix type not implemented");
1231  }
1232  CHKERR MatSetDM(*M, dm);
1233  MoFEMFunctionReturn(0);
1234 }

◆ DMCreateMatrix_MoFEM() [2/2]

PetscErrorCode MoFEM::DMCreateMatrix_MoFEM ( DM  dm,
SmartPetscObj< Mat > &  M 
)

DMShellSetCreateMatrix

sets the routine to create a matrix associated with the shell DM

Definition at line 1236 of file DMMoFEM.cpp.

1236  {
1237  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1238  MoFEMFunctionBegin;
1239  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1240 
1241  if (strcmp(dm->mattype, MATSHELL) == 0) {
1242  if (dm_field->blocMatDataPtr) {
1243  Mat mat_raw;
1244  CHKERR DMMoFEMCreateBlockMat(dm, &mat_raw);
1245  M = SmartPetscObj<Mat>(mat_raw);
1246  MoFEMFunctionReturnHot(0);
1247  } else {
1248  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
1249  "Matrix shell data not set, or matrix type not implemented");
1250  }
1251  } else if (strcmp(dm->mattype, MATMPIAIJ) == 0) {
1252  CHKERR dm_field->mField_ptr->getInterface<MatrixManager>()
1253  ->createMPIAIJWithArrays<PetscGlobalIdx_mi_tag>(dm_field->problemName,
1254  M);
1255  } else if (strcmp(dm->mattype, MATAIJ) == 0) {
1256  CHKERR dm_field->mField_ptr->getInterface<MatrixManager>()
1257  ->createSeqAIJWithArrays<PetscLocalIdx_mi_tag>(dm_field->problemName,
1258  M);
1259  } else if (strcmp(dm->mattype, MATAIJCUSPARSE) == 0) {
1260  CHKERR dm_field->mField_ptr->getInterface<MatrixManager>()
1261  ->createMPIAIJCUSPARSEWithArrays<PetscGlobalIdx_mi_tag>(
1262  dm_field->problemName, M);
1263  } else if (strcmp(dm->mattype, MATSEQAIJCUSPARSE) == 0) {
1264  CHKERR dm_field->mField_ptr->getInterface<MatrixManager>()
1265  ->createSeqAIJCUSPARSEWithArrays<PetscLocalIdx_mi_tag>(
1266  dm_field->problemName, M);
1267  } else {
1268  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
1269  "Matrix type not implemented");
1270  }
1271  CHKERR MatSetDM(M, dm);
1272  MoFEMFunctionReturn(0);
1273 }

◆ DMDestroy_MoFEM()

PetscErrorCode MoFEM::DMDestroy_MoFEM ( DM  dm)

Destroys dm with MoFEM data structure.

destroy the MoFEM structure

Definition at line 79 of file DMMoFEM.cpp.

79  {
80  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
81  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
82  MoFEMFunctionBeginHot;
83 
84  MPI_Comm comm;
85  CHKERR PetscObjectGetComm((PetscObject)dm, &comm);
86 
87  int result;
88  MPI_Comm_compare(comm, PETSC_COMM_SELF, &result);
89  if (result == MPI_IDENT)
90  MOFEM_LOG("DMSELF", Sev::noisy)
91  << "MoFEM DM destroy for problem " << dm_field->problemName
92  << " referenceNumber " << dm_field->referenceNumber;
93  else
94  MOFEM_LOG("DMWORLD", Sev::noisy)
95  << "MoFEM DM destroy for problem " << dm_field->problemName
96  << " referenceNumber " << dm_field->referenceNumber;
97 
98  if (dm_field->referenceNumber == 0) {
99  if (dm_field->destroyProblem) {
100 
101  if (dm_field->mField_ptr->check_problem(dm_field->problemName)) {
102  dm_field->mField_ptr->delete_problem(dm_field->problemName);
103  } // else problem has to be deleted by the user
104  }
105 
106  delete static_cast<DMCtxImpl *>(dm->data);
107 
108  } else
109  --dm_field->referenceNumber;
110 
111  MoFEMFunctionReturnHot(0);
112 }

◆ DMGlobalToLocalBegin_MoFEM()

PetscErrorCode MoFEM::DMGlobalToLocalBegin_MoFEM ( DM  dm,
Vec  g,
InsertMode  mode,
Vec  l 
)

DMShellSetGlobalToLocal

the routine that begins the global to local scatter

Definition at line 1397 of file DMMoFEM.cpp.

1398  {
1399  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1400  MoFEMFunctionBegin;
1401  CHKERR VecGhostUpdateBegin(g, INSERT_VALUES, SCATTER_FORWARD);
1402  MoFEMFunctionReturn(0);
1403 }

◆ DMGlobalToLocalEnd_MoFEM()

PetscErrorCode MoFEM::DMGlobalToLocalEnd_MoFEM ( DM  dm,
Vec  g,
InsertMode  mode,
Vec  l 
)

DMShellSetGlobalToLocal

the routine that begins the global to local scatter

Definition at line 1405 of file DMMoFEM.cpp.

1405  {
1406  MoFEMFunctionBeginHot;
1407  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1408  MoFEMFunctionBegin;
1409 
1410  CHKERR VecGhostUpdateEnd(g, INSERT_VALUES, SCATTER_FORWARD);
1411 
1412  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1413  int nb_dofs = dm_field->problemPtr->getNbLocalDofsRow();
1414  int nb_ghost = dm_field->problemPtr->getNbGhostDofsRow();
1415 
1416  double *array_loc, *array_glob;
1417  CHKERR VecGetArray(l, &array_loc);
1418  CHKERR VecGetArray(g, &array_glob);
1419  switch (mode) {
1420  case INSERT_VALUES:
1421  cblas_dcopy(nb_dofs + nb_ghost, array_glob, 1, array_loc, 1);
1422  break;
1423  case ADD_VALUES:
1424  cblas_daxpy(nb_dofs + nb_ghost, 1, array_glob, 1, array_loc, 1);
1425  break;
1426  default:
1427  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "not implemented");
1428  }
1429  CHKERR VecRestoreArray(l, &array_loc);
1430  CHKERR VecRestoreArray(g, &array_glob);
1431  MoFEMFunctionReturn(0);
1432 }

◆ DMLocalToGlobalBegin_MoFEM()

PetscErrorCode MoFEM::DMLocalToGlobalBegin_MoFEM ( DM  dm,
Vec  l,
InsertMode  mode,
Vec  g 
)

DMShellSetLocalToGlobal

the routine that begins the local to global scatter

DMShellSetLocalToGlobal

the routine that ends the local to global scatter

Definition at line 1434 of file DMMoFEM.cpp.

1435  {
1436 
1437  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1438  MoFEMFunctionBegin;
1439 
1440  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1441  int nb_dofs = dm_field->problemPtr->getNbLocalDofsRow();
1442  int nb_ghost = dm_field->problemPtr->getNbGhostDofsRow();
1443 
1444  double *array_loc, *array_glob;
1445  CHKERR VecGetArray(l, &array_loc);
1446  CHKERR VecGetArray(g, &array_glob);
1447  switch (mode) {
1448  case INSERT_VALUES:
1449  cblas_dcopy(nb_dofs + nb_ghost, array_loc, 1, array_glob, 1);
1450  break;
1451  case ADD_VALUES:
1452  cblas_daxpy(nb_dofs + nb_ghost, 1, array_loc, 1, array_glob, 1);
1453  break;
1454  default:
1455  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "not implemented");
1456  }
1457  CHKERR VecRestoreArray(l, &array_loc);
1458  CHKERR VecRestoreArray(g, &array_glob);
1459 
1460  MoFEMFunctionReturn(0);
1461 }

◆ DMLocalToGlobalEnd_MoFEM()

PetscErrorCode MoFEM::DMLocalToGlobalEnd_MoFEM ( DM  ,
Vec  l,
InsertMode  mode,
Vec  g 
)

DMShellSetLocalToGlobal

the routine that ends the local to global scatter

Definition at line 1463 of file DMMoFEM.cpp.

1463  {
1464  //
1465  MoFEMFunctionBeginHot;
1466  MoFEMFunctionReturnHot(0);
1467 }

◆ DMMGViaApproxOrdersClearCoarseningIS() [1/2]

MoFEMErrorCode DMMGViaApproxOrdersClearCoarseningIS ( DM  )

Clear approximation orders.

Parameters
DMdm
Returns
Error code

Definition at line 293 of file PCMGSetUpViaApproxOrders.cpp.

293  {
294  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
295  MoFEMFunctionBegin;
296  GET_DM_FIELD(dm);
297  CHKERR dm_field->destroyCoarseningIS();
298  PetscInfo(dm, "Clear DMs data structures\n");
299  MoFEMFunctionReturn(0);
300 }

◆ DMMGViaApproxOrdersClearCoarseningIS() [2/2]

MoFEMErrorCode MoFEM::DMMGViaApproxOrdersClearCoarseningIS ( DM  )

Clear approximation orders.

Parameters
DMdm
Returns
Error code

Definition at line 293 of file PCMGSetUpViaApproxOrders.cpp.

293  {
294  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
295  MoFEMFunctionBegin;
296  GET_DM_FIELD(dm);
297  CHKERR dm_field->destroyCoarseningIS();
298  PetscInfo(dm, "Clear DMs data structures\n");
299  MoFEMFunctionReturn(0);
300 }

◆ DMMGViaApproxOrdersPopBackCoarseningIS() [1/2]

MoFEMErrorCode DMMGViaApproxOrdersPopBackCoarseningIS ( DM  )

Pop is form MG via approximation orders.

Parameters
DMdm
ispop back IS
Returns
error code

Pop is form MG via approximation orders.

Parameters
DMdm
Returns
error code

Definition at line 281 of file PCMGSetUpViaApproxOrders.cpp.

281  {
282  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
283  MoFEMFunctionBegin;
284  GET_DM_FIELD(dm);
285  if (dm_field->coarseningIS.back()) {
286  dm_field->coarseningIS.pop_back();
287  }
288  dm_field->kspOperators.pop_back();
289  PetscInfo(dm, "Pop back IS to DMMGViaApproxOrders\n");
290  MoFEMFunctionReturn(0);
291 }

◆ DMMGViaApproxOrdersPopBackCoarseningIS() [2/2]

MoFEMErrorCode MoFEM::DMMGViaApproxOrdersPopBackCoarseningIS ( DM  )

Pop IS form MG via approximation orders.

Pop is form MG via approximation orders.

Parameters
DMdm
Returns
error code

Definition at line 281 of file PCMGSetUpViaApproxOrders.cpp.

281  {
282  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
283  MoFEMFunctionBegin;
284  GET_DM_FIELD(dm);
285  if (dm_field->coarseningIS.back()) {
286  dm_field->coarseningIS.pop_back();
287  }
288  dm_field->kspOperators.pop_back();
289  PetscInfo(dm, "Pop back IS to DMMGViaApproxOrders\n");
290  MoFEMFunctionReturn(0);
291 }

◆ DMMGViaApproxOrdersPushBackCoarseningIS() [1/2]

MoFEMErrorCode MoFEM::DMMGViaApproxOrdersPushBackCoarseningIS ( DM  ,
IS  is,
Mat  A,
bool  create_sub_matrix,
bool  shell_sub_a 
)

Push back coarsening level to MG via approximation orders.

Parameters
DMdiscrete manager
isPush back IS used for coarsening
AGet sub-matrix of A using is (that sets operators for coarsening levels)
subAReturning pointer to created sub matrix
subAIf true create sub matrix, otherwise in subA has to be valid pointer to subA
Returns
Error code

Definition at line 227 of file PCMGSetUpViaApproxOrders.cpp.

229  {
230  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
231  MoFEMFunctionBegin;
232  GET_DM_FIELD(dm);
233  dm_field->coarseningIS.emplace_back(SmartPetscObj<IS>(is, true));
234  dm_field->shellMatrixCtxPtr.push_back(new PCMGSubMatrixCtx(A, is));
235  if (is) {
236  auto is2 = SmartPetscObj<IS>(is, true);
237  if (dm_field->aO) {
238  is2 = isDuplicate(is);
239  CHKERR ISCopy(is, is2);
240  CHKERR AOApplicationToPetscIS(dm_field->aO, is2);
241  }
242  if (create_sub_matrix) {
243  if (shell_sub_a) {
244  int n, N;
245  CHKERR ISGetSize(is, &N);
246  CHKERR ISGetLocalSize(is, &n);
247  MPI_Comm comm;
248  CHKERR PetscObjectGetComm((PetscObject)A, &comm);
249  Mat sub_a_raw;
250  CHKERR MatCreateShell(comm, n, n, N, N,
251  &(dm_field->shellMatrixCtxPtr.back()),
252  &sub_a_raw);
253  CHKERR MatShellSetOperation(sub_a_raw, MATOP_MULT,
254  (void (*)(void))sub_mat_mult);
255  CHKERR MatShellSetOperation(sub_a_raw, MATOP_MULT_ADD,
256  (void (*)(void))sub_mat_mult_add);
257  CHKERR MatShellSetOperation(sub_a_raw, MATOP_SOR,
258  (void (*)(void))sub_mat_sor);
259  dm_field->kspOperators.emplace_back(
260  SmartPetscObj<Mat>(sub_a_raw, false));
261  } else {
262  Mat sub_a_raw;
263 #if PETSC_VERSION_GE(3, 8, 0)
264  CHKERR MatCreateSubMatrix(A, is2, is2, MAT_INITIAL_MATRIX, &sub_a_raw);
265 #else
266  CHKERR MatGetSubMatrix(A, is2, is2, MAT_INITIAL_MATRIX, &sub_a_raw);
267 #endif
268  dm_field->kspOperators.emplace_back(
269  SmartPetscObj<Mat>(sub_a_raw, false));
270  }
271  } else {
272  dm_field->kspOperators.emplace_back(SmartPetscObj<Mat>(A, true));
273  }
274  } else {
275  SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
276  }
277  PetscInfo(dm, "Push back IS to DMMGViaApproxOrders\n");
278  MoFEMFunctionReturn(0);
279 }

◆ DMMGViaApproxOrdersPushBackCoarseningIS() [2/2]

MoFEMErrorCode DMMGViaApproxOrdersPushBackCoarseningIS ( DM  ,
IS  is,
Mat  A,
Mat *  subA,
bool  create_sub_matrix,
bool  shell_sub_a 
)

Push back coarsening level to MG via approximation orders.

Parameters
DMdiscrete manager
isPush back IS used for coarsening
AGet sub-matrix of A using is (that sets operators for coarsening levels)
subAReturning pointer to created sub matrix
subAIf true create sub matrix, otherwise in subA has to be valid pointer to subA
Returns
Error code

Definition at line 205 of file PCMGSetUpViaApproxOrders.cpp.

208  {
209  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
210  MoFEMFunctionBegin;
211  GET_DM_FIELD(dm);
212  dm_field->coarseningIS.push_back(is);
213  dm_field->shellMatrixCtxPtr.push_back(new PCMGSubMatrixCtx_private(A, is));
214  if (is) {
215  CHKERR PetscObjectReference((PetscObject)is);
216  }
217  if (is) {
218  IS is2 = is;
219  if (dm_field->aO) {
220  CHKERR ISDuplicate(is, &is2);
221  CHKERR ISCopy(is, is2);
222  CHKERR AOApplicationToPetscIS(dm_field->aO, is2);
223  }
224  if (create_sub_matrix) {
225  if (shell_sub_a) {
226  int n, N;
227  CHKERR ISGetSize(is, &N);
228  CHKERR ISGetLocalSize(is, &n);
229  MPI_Comm comm;
230  CHKERR PetscObjectGetComm((PetscObject)A, &comm);
231  CHKERR MatCreateShell(comm, n, n, N, N,
232  &(dm_field->shellMatrixCtxPtr.back()), subA);
233  CHKERR MatShellSetOperation(*subA, MATOP_MULT,
234  (void (*)(void))sub_mat_mult);
235  CHKERR MatShellSetOperation(*subA, MATOP_MULT_ADD,
236  (void (*)(void))sub_mat_mult_add);
237  CHKERR MatShellSetOperation(*subA, MATOP_SOR,
238  (void (*)(void))sub_mat_sor);
239  } else {
240 #if PETSC_VERSION_GE(3, 8, 0)
241  CHKERR MatCreateSubMatrix(A, is2, is2, MAT_INITIAL_MATRIX, subA);
242 #else
243  CHKERR MatGetSubMatrix(A, is2, is2, MAT_INITIAL_MATRIX, subA);
244 #endif
245  }
246  }
247  if (dm_field->aO) {
248  CHKERR ISDestroy(&is2);
249  }
250  dm_field->kspOperators.push_back(*subA);
251  CHKERR PetscObjectReference((PetscObject)(*subA));
252  } else {
253  SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
254  }
255  PetscInfo(dm, "Push back IS to DMMGViaApproxOrders\n");
256  MoFEMFunctionReturn(0);
257 }

◆ DMMoFEMAddColCompositeProblem()

PetscErrorCode MoFEM::DMMoFEMAddColCompositeProblem ( DM  dm,
const char  prb_name[] 
)

Add problem to composite DM on col.

This create block on col with DOFs from problem of given name

Parameters
dmthe DM object
prb_nameadd problem name
Returns
error code

Definition at line 385 of file DMMoFEM.cpp.

385  {
386  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
387  MoFEMFunctionBegin;
388  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
389  if (!dm->data) {
390  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
391  "data structure for MoFEM not yet created");
392  }
393  if (!dm_field->isCompDM) {
394  dm_field->isCompDM = PETSC_TRUE;
395  }
396  if (dm_field->isSquareMatrix) {
397  SETERRQ(PETSC_COMM_SELF, MOFEM_INVALID_DATA,
398  "No need to add problem on column when problem block structurally "
399  "symmetric");
400  }
401  dm_field->colCompPrb.push_back(prb_name);
402  MoFEMFunctionReturn(0);
403 }

◆ DMMoFEMAddElement() [1/2]

PetscErrorCode MoFEM::DMMoFEMAddElement ( DM  dm,
std::string  fe_name 
)

add element to dm

Note
add_file is a collective, should be executed on all processors. Otherwise could lead to deadlock.
Examples
analytical_poisson_field_split.cpp, cell_forces.cpp, dm_build_partitioned_mesh.cpp, dm_create_subdm.cpp, dm_partitioned_no_field.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, EshelbianPlasticity.cpp, free_surface.cpp, level_set.cpp, lorentz_force.cpp, MagneticElement.hpp, minimal_surface_area.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, NonlinearElasticElementInterface.hpp, photon_diffusion.cpp, plastic.cpp, Remodeling.cpp, schur_test_diag_mat.cpp, simple_contact.cpp, simple_contact_thermal.cpp, simple_elasticity.cpp, test_broken_space.cpp, test_jacobian_of_simple_contact_element.cpp, and UnsaturatedFlow.hpp.

Definition at line 497 of file DMMoFEM.cpp.

497  {
498  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
499  MoFEMFunctionBegin;
500  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
501  CHKERR dm_field->mField_ptr->modify_problem_add_finite_element(
502  dm_field->problemName, fe_name);
503  MoFEMFunctionReturn(0);
504 }

◆ DMMoFEMAddElement() [2/2]

PetscErrorCode MoFEM::DMMoFEMAddElement ( DM  dm,
std::vector< std::string >  fe_name 
)

add element to dm

Note
add_file is a collective, should be executed on all processors. Otherwise could lead to deadlock.

Definition at line 506 of file DMMoFEM.cpp.

506  {
507  MoFEMFunctionBegin;
508  for (auto fe : fe_name) {
509  CHKERR DMMoFEMAddElement(dm, fe);
510  }
511  MoFEMFunctionReturn(0);
512 }

◆ DMMoFEMAddRowCompositeProblem()

PetscErrorCode MoFEM::DMMoFEMAddRowCompositeProblem ( DM  dm,
const char  prb_name[] 
)

Add problem to composite DM on row.

This create block on row with DOFs from problem of given name

Parameters
dmthe DM object
prb_nameadd problem name
Returns
error code

Definition at line 367 of file DMMoFEM.cpp.

367  {
368  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
369  MoFEMFunctionBegin;
370  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
371  if (!dm->data) {
372  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
373  "data structure for MoFEM not yet created");
374  }
375  if (!dm_field->isCompDM) {
376  dm_field->isCompDM = PETSC_TRUE;
377  }
378  dm_field->rowCompPrb.push_back(prb_name);
379  if (dm_field->isSquareMatrix) {
380  dm_field->colCompPrb.push_back(prb_name);
381  }
382  MoFEMFunctionReturn(0);
383 }

◆ DMMoFEMAddSubFieldCol()

PetscErrorCode MoFEM::DMMoFEMAddSubFieldCol ( DM  dm,
const char  field_name[] 
)

Add field to sub dm problem on columns

Examples
analytical_poisson_field_split.cpp, cell_forces.cpp, dm_create_subdm.cpp, free_surface.cpp, plastic.cpp, schur_test_diag_mat.cpp, and test_broken_space.cpp.

Definition at line 280 of file DMMoFEM.cpp.

280  {
281  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
282  MoFEMFunctionBeginHot;
283  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
284  if (!dm->data) {
285  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
286  "data structure for MoFEM not yet created");
287  }
288  if (!dm_field->isSubDM) {
289  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "this is not sub-dm");
290  }
291  dm_field->colSubFields.push_back(field_name);
292  dm_field->mapTypeCol.erase(field_name);
293  MoFEMFunctionReturnHot(0);
294 }

◆ DMMoFEMAddSubFieldRow() [1/2]

PetscErrorCode MoFEM::DMMoFEMAddSubFieldRow ( DM  dm,
const char  field_name[] 
)

Add field to sub dm problem on rows

Examples
analytical_poisson_field_split.cpp, cell_forces.cpp, dm_create_subdm.cpp, EshelbianPlasticity.cpp, free_surface.cpp, level_set.cpp, plastic.cpp, schur_test_diag_mat.cpp, and test_broken_space.cpp.

Definition at line 238 of file DMMoFEM.cpp.

238  {
239  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
240  MoFEMFunctionBeginHot;
241  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
242  if (!dm->data) {
243  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
244  "data structure for MoFEM not yet created");
245  }
246  if (!dm_field->isSubDM) {
247  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "this is not sub-dm");
248  }
249  dm_field->rowSubFields.push_back(field_name);
250  dm_field->mapTypeRow.erase(field_name);
251  MoFEMFunctionReturnHot(0);
252 }

◆ DMMoFEMAddSubFieldRow() [2/2]

PetscErrorCode MoFEM::DMMoFEMAddSubFieldRow ( DM  dm,
const char  field_name[],
boost::shared_ptr< Range >  r_ptr 
)

Add field to sub dm problem on rows.

Parameters
dm
field_name
m
Returns
PetscErrorCode

Definition at line 258 of file DMMoFEM.cpp.

259  {
260  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
261  MoFEMFunctionBeginHot;
262  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
263  if (!dm->data) {
264  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
265  "data structure for MoFEM not yet created");
266  }
267  if (!dm_field->isSubDM) {
268  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "this is not sub-dm");
269  }
270  dm_field->rowSubFields.push_back(field_name);
271  dm_field->mapTypeRow[field_name] = r_ptr;
272  MoFEMFunctionReturnHot(0);
273 }

◆ DMMoFEMCreateBlockMat() [1/2]

MoFEMErrorCode MoFEM::DMMoFEMCreateBlockMat ( DM  dm,
Mat *  mat 
)

Create block matrix.

Parameters
dm
mat
Returns
MoFEMErrorCode

Definition at line 1551 of file DMMoFEM.cpp.

1551  {
1552  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1553  MoFEMFunctionBegin;
1554  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1555  auto mat_data = createBlockMat(dm, dm_field->blocMatDataPtr);
1556  *mat = mat_data.first;
1557  CHKERR PetscObjectReference((PetscObject)(*mat));
1558  MoFEMFunctionReturn(0);
1559 }

◆ DMMoFEMCreateBlockMat() [2/2]

MoFEMErrorCode MoFEM::DMMoFEMCreateBlockMat ( DM  dm,
SmartPetscObj< Mat > &  mat 
)

Create block matrix.

Parameters
dm
matsmart pointer
Returns
MoFEMErrorCode

◆ DMMoFEMCreateMoFEM()

PetscErrorCode MoFEM::DMMoFEMCreateMoFEM ( DM  dm,
MoFEM::Interface *  m_field_ptr,
const char  problem_name[],
const MoFEM::BitRefLevel  bit_level,
const MoFEM::BitRefLevel  bit_mask = MoFEM::BitRefLevel().set() 
)

Must be called by user to set MoFEM data structures.

Note
If problem exist function create DM for it. If you set bit levels, those bits are to existing bits. Thus if you do not like to change bit ref level for existing problem, set bits to zero.
Examples
cell_forces.cpp, dm_build_partitioned_mesh.cpp, dm_create_subdm.cpp, dm_partitioned_no_field.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, EshelbianPlasticity.cpp, free_surface.cpp, lorentz_force.cpp, MagneticElement.hpp, minimal_surface_area.cpp, mortar_contact.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, Remodeling.cpp, simple_contact.cpp, simple_contact_thermal.cpp, test_jacobian_of_simple_contact_element.cpp, and UnsaturatedFlow.hpp.

Definition at line 114 of file DMMoFEM.cpp.

117  {
118  MoFEMFunctionBegin;
119 
120  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
121  if (!dm->data) {
122  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
123  "data structure for MoFEM not yet created");
124  }
125  if (!m_field_ptr) {
126  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
127  "DM function not implemented into MoFEM");
128  }
129  dm_field->mField_ptr = m_field_ptr;
130  dm_field->problemName = problem_name;
131  if (!m_field_ptr->check_problem(dm_field->problemName)) {
132  // problem is not defined, declare problem internally set bool to
133  // destroyProblem problem with DM
134  dm_field->destroyProblem = PETSC_TRUE;
135  CHKERR dm_field->mField_ptr->add_problem(dm_field->problemName, MF_EXCL,
136  dm_field->verbosity);
137  } else {
138  dm_field->destroyProblem = PETSC_FALSE;
139  }
140  CHKERR dm_field->mField_ptr->modify_problem_ref_level_add_bit(
141  dm_field->problemName, bit_level);
142  CHKERR dm_field->mField_ptr->modify_problem_mask_ref_level_add_bit(
143  dm_field->problemName, bit_mask);
144  dm_field->kspCtx =
145  boost::shared_ptr<KspCtx>(new KspCtx(*m_field_ptr, problem_name));
146  dm_field->snesCtx =
147  boost::shared_ptr<SnesCtx>(new SnesCtx(*m_field_ptr, problem_name));
148  dm_field->tsCtx =
149  boost::shared_ptr<TsCtx>(new TsCtx(*m_field_ptr, problem_name));
150 
151  MPI_Comm comm;
152  CHKERR PetscObjectGetComm((PetscObject)dm, &comm);
153  int result = 0;
154  MPI_Comm_compare(comm, m_field_ptr->get_comm(), &result);
155  if (result > MPI_CONGRUENT) {
156  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
157  "MoFEM and DM using different communicators");
158  }
159  MPI_Comm_size(comm, &dm_field->sIze);
160  MPI_Comm_rank(comm, &dm_field->rAnk);
161 
162  // problem structure
163  CHKERR dm_field->mField_ptr->get_problem(dm_field->problemName,
164  &dm_field->problemPtr);
165 
166  MPI_Comm_compare(comm, PETSC_COMM_SELF, &result);
167  if (result == MPI_IDENT) {
168  MOFEM_LOG("DMSELF", Sev::noisy)
169  << "MoFEM DM created for problem " << dm_field->problemName;
170  MOFEM_LOG("DMSELF", Sev::noisy) << *dm_field->problemPtr;
171  } else {
172  MOFEM_LOG("DMWORLD", Sev::noisy)
173  << "MoFEM DM created for problem " << dm_field->problemName;
174  MOFEM_LOG("DMWORLD", Sev::noisy) << *dm_field->problemPtr;
175  }
176 
177  MoFEMFunctionReturn(0);
178 }

◆ DMMoFEMCreateSubDM()

PetscErrorCode MoFEM::DMMoFEMCreateSubDM ( DM  subdm,
DM  dm,
const char  problem_name[] 
)

Must be called by user to set Sub DM MoFEM data structures.

Examples
analytical_poisson_field_split.cpp, cell_forces.cpp, dm_create_subdm.cpp, EshelbianPlasticity.cpp, free_surface.cpp, level_set.cpp, plastic.cpp, schur_test_diag_mat.cpp, and test_broken_space.cpp.

Definition at line 215 of file DMMoFEM.cpp.

215  {
216  MoFEMFunctionBegin;
217 
218  if (!dm->data) {
219  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
220  "data structure for MoFEM not yet created");
221  }
222  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
223 
224  CHKERR DMMoFEMCreateMoFEM(subdm, dm_field->mField_ptr, problem_name,
225  dm_field->problemPtr->getBitRefLevel(),
226  dm_field->problemPtr->getBitRefLevelMask());
227 
228  DMCtxImpl *subdm_field = (DMCtxImpl *)subdm->data;
229  subdm_field->isSubDM = PETSC_TRUE;
230  subdm_field->problemMainOfSubPtr = dm_field->problemPtr;
231  subdm_field->isPartitioned = dm_field->isPartitioned;
232  subdm_field->isSquareMatrix = PETSC_FALSE;
233  subdm->ops->setup = DMSubDMSetUp_MoFEM;
234 
235  MoFEMFunctionReturn(0);
236 }

◆ DMMoFEMGetFieldIS()

PetscErrorCode MoFEM::DMMoFEMGetFieldIS ( DM  dm,
RowColData  rc,
const char  field_name[],
IS *  is 
)

get field is in the problem

Parameters
dmthe DM objects
rcROW or COL (no difference is problem is squared)
field_namename of the field
isreturned the IS object
Returns
error code
IS is;
ierr = DMMoFEMGetFieldIS(dm,ROW,"DISP",&is_disp); CHKERRG(ierr);

Definition at line 1514 of file DMMoFEM.cpp.

1515  {
1516  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1517  MoFEMFunctionBegin;
1518  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1519  CHKERR dm_field->mField_ptr->getInterface<ISManager>()
1520  ->isCreateProblemFieldAndRank(dm_field->problemPtr->getName(), ROW,
1521  field_name, 0, 1000, is);
1522  MoFEMFunctionReturn(0);
1523 }

◆ DMMoFEMGetIsCompDM()

PetscErrorCode MoFEM::DMMoFEMGetIsCompDM ( DM  dm,
PetscBool *  is_comp_dm 
)

Get if this DM is composite DM.

Parameters
dmthe DM object
is_comp_dmreturn true if composite problem here
Returns
error code

Definition at line 405 of file DMMoFEM.cpp.

405  {
406  MoFEMFunctionBeginHot;
407  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
408  MoFEMFunctionBeginHot;
409  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
410  *is_comp_dm = dm_field->isCompDM;
411  MoFEMFunctionReturnHot(0);
412 }

◆ DMMoFEMGetIsPartitioned()

PetscErrorCode MoFEM::DMMoFEMGetIsPartitioned ( DM  dm,
PetscBool *  is_partitioned 
)

get if read mesh is partitioned

Definition at line 1134 of file DMMoFEM.cpp.

1134  {
1135  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1136  MoFEMFunctionBeginHot;
1137  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1138  *is_partitioned = dm_field->isPartitioned;
1139  MoFEMFunctionReturnHot(0);
1140 }

◆ DMMoFEMGetIsSubDM()

PetscErrorCode MoFEM::DMMoFEMGetIsSubDM ( DM  dm,
PetscBool *  is_sub_dm 
)

Return true if this DM is sub problem

Parameters
dmthe DM object
is_subproblemtrue if subproblem
Returns
error code

Definition at line 322 of file DMMoFEM.cpp.

322  {
323  MoFEMFunctionBeginHot;
324  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
325  MoFEMFunctionBeginHot;
326  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
327  *is_sub_dm = dm_field->isSubDM;
328  MoFEMFunctionReturnHot(0);
329 }

◆ DMMoFEMGetKspCtx() [1/2]

PetscErrorCode MoFEM::DMMoFEMGetKspCtx ( DM  dm,
const boost::shared_ptr< MoFEM::KspCtx > &  ksp_ctx 
)

get MoFEM::KspCtx data structure

Definition at line 1077 of file DMMoFEM.cpp.

1077  {
1078  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1079  MoFEMFunctionBeginHot;
1080  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1081  const_cast<boost::shared_ptr<MoFEM::KspCtx> &>(ksp_ctx) = dm_field->kspCtx;
1082  MoFEMFunctionReturnHot(0);
1083 }

◆ DMMoFEMGetKspCtx() [2/2]

PetscErrorCode MoFEM::DMMoFEMGetKspCtx ( DM  dm,
MoFEM::KspCtx **  ksp_ctx 
)

get MoFEM::KspCtx data structure

Definition at line 1068 of file DMMoFEM.cpp.

1068  {
1069  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1070  MoFEMFunctionBeginHot;
1071  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1072  *ksp_ctx = dm_field->kspCtx.get();
1073  MoFEMFunctionReturnHot(0);
1074 }

◆ DMMoFEMGetProblemFiniteElementLayout()

PetscErrorCode MoFEM::DMMoFEMGetProblemFiniteElementLayout ( DM  dm,
std::string  fe_name,
PetscLayout *  layout 
)

Get finite elements layout in the problem.

In layout is stored information how many elements is on each processor, for more information look int petsc documentation http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/IS/PetscLayoutCreate.html#PetscLayoutCreate

Parameters
dmdiscrete manager for this problem
fe_namefinite element name
layoutpointer to layout, for created layout user takes responsibility for destroying it.
Returns
error code

Definition at line 473 of file DMMoFEM.cpp.

474  {
475  MoFEMFunctionBeginHot;
476  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
477  MoFEMFunctionBegin;
478  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
479 
480  MPI_Comm comm;
481  CHKERR PetscObjectGetComm((PetscObject)dm, &comm);
482  CHKERR dm_field->problemPtr->getNumberOfElementsByNameAndPart(comm, fe_name,
483  layout);
484  MoFEMFunctionReturn(0);
485 }

◆ DMMoFEMGetProblemPtr()

PetscErrorCode MoFEM::DMMoFEMGetProblemPtr ( DM  dm,
const MoFEM::Problem **  problem_ptr 
)

Get pointer to problem data structure.

Examples
field_evaluator.cpp, loop_entities.cpp, lorentz_force.cpp, MagneticElement.hpp, partition_mesh.cpp, poisson_2d_dis_galerkin.cpp, Remodeling.cpp, and test_cache_on_entities.cpp.

Definition at line 426 of file DMMoFEM.cpp.

426  {
427  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
428  MoFEMFunctionBeginHot;
429  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
430  if (!dm->data) {
431  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
432  "data structure for MoFEM not yet created");
433  }
434  *problem_ptr = dm_field->problemPtr;
435  MoFEMFunctionReturnHot(0);
436 }

◆ DMMoFEMGetSnesCtx() [1/2]

PetscErrorCode MoFEM::DMMoFEMGetSnesCtx ( DM  dm,
const boost::shared_ptr< MoFEM::SnesCtx > &  snes_ctx 
)

get MoFEM::SnesCtx data structure

Definition at line 1103 of file DMMoFEM.cpp.

1103  {
1104  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1105  MoFEMFunctionBeginHot;
1106  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1107  const_cast<boost::shared_ptr<MoFEM::SnesCtx> &>(snes_ctx) = dm_field->snesCtx;
1108  MoFEMFunctionReturnHot(0);
1109 }

◆ DMMoFEMGetSnesCtx() [2/2]

PetscErrorCode MoFEM::DMMoFEMGetSnesCtx ( DM  dm,
MoFEM::SnesCtx **  snes_ctx 
)

get MoFEM::SnesCtx data structure

Examples
minimal_surface_area.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, simple_contact.cpp, simple_contact_thermal.cpp, test_jacobian_of_simple_contact_element.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 1094 of file DMMoFEM.cpp.

1094  {
1095  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1096  MoFEMFunctionBeginHot;
1097  DMCtxImpl *dm_field = (DMCtxImpl *)dm->data;
1098  *snes_ctx = dm_field->snesCtx.get();
1099  MoFEMFunctionReturnHot(0);
1100 }

◆ DMMoFEMGetSquareProblem()

PetscErrorCode MoFEM::DMMoFEMGetSquareProblem ( DM  dm,
PetscBool *  square_problem 
)

get squared problem

It if true is assumed that matrix has the same indexing on rows and columns. This reduces interprocessor communication.

Definition at line 487 of file DMMoFEM.cpp.

487  {
488  MoFEMFunctionBeginHot;
489  MoFEMFunctionBeginHot;
490  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
491  MoFEMFunctionBeginHot;
492  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
493  *square_problem = dm_field->isSquareMatrix;
494  MoFEMFunctionReturnHot(0);
495 }

◆ DMMoFEMGetTsCtx() [1/2]

PetscErrorCode MoFEM::DMMoFEMGetTsCtx ( DM  dm,
const boost::shared_ptr< MoFEM::TsCtx > &  ts_ctx 
)

get MoFEM::TsCtx data structure

Definition at line 1150 of file DMMoFEM.cpp.

1151  {
1152  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1153  MoFEMFunctionBeginHot;
1154  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1155  const_cast<boost::shared_ptr<MoFEM::TsCtx> &>(ts_ctx) = dm_field->tsCtx;
1156  MoFEMFunctionReturnHot(0);
1157 }

◆ DMMoFEMGetTsCtx() [2/2]

PetscErrorCode MoFEM::DMMoFEMGetTsCtx ( DM  dm,
MoFEM::TsCtx **  ts_ctx 
)

get MoFEM::TsCtx data structure

Examples
EshelbianPlasticity.cpp, Remodeling.cpp, and UnsaturatedFlow.hpp.

Definition at line 1142 of file DMMoFEM.cpp.

1142  {
1143  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1144  MoFEMFunctionBeginHot;
1145  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1146  *ts_ctx = dm_field->tsCtx.get();
1147  MoFEMFunctionReturnHot(0);
1148 }

◆ DMMoFEMKSPSetComputeOperators() [1/2]

PetscErrorCode MoFEM::DMMoFEMKSPSetComputeOperators ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
MoFEM::BasicMethod *  pre_only,
MoFEM::BasicMethod *  post_only 
)

Set KSP operators and push mofem finite element methods.

Parameters
dmDM
fe_namefinite element name
methodmethod on the element (executed for each element in the problem which given name)
pre_onlymethod for pre-process before element method
post_onlymethod for post-process after element method
Returns
error code
Examples
analytical_poisson.cpp, elasticity_mixed_formulation.cpp, level_set.cpp, and simple_elasticity.cpp.

Definition at line 678 of file DMMoFEM.cpp.

681  {
682  return DMMoFEMKSPSetComputeOperators<const char *, MoFEM::FEMethod *,
683  MoFEM::BasicMethod *,
684  MoFEM::BasicMethod *>(
685  dm, fe_name, method, pre_only, post_only);
686 }

◆ DMMoFEMKSPSetComputeOperators() [2/2]

PetscErrorCode MoFEM::DMMoFEMKSPSetComputeOperators ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
boost::shared_ptr< MoFEM::BasicMethod >  pre_only,
boost::shared_ptr< MoFEM::BasicMethod >  post_only 
)

Set KSP operators and push mofem finite element methods.

Parameters
dmDM
fe_namefinite element name
methodmethod on the element (executed for each element in the problem which given name)
pre_onlymethod for pre-process before element method
post_onlymethod for post-process after element method
Returns
error code

Definition at line 689 of file DMMoFEM.cpp.

692  {
693  return DMMoFEMKSPSetComputeOperators<const std::string,
694  boost::shared_ptr<MoFEM::FEMethod>>(
695  dm, fe_name, method, pre_only, post_only);
696 }

◆ DMMoFEMKSPSetComputeRHS() [1/2]

PetscErrorCode MoFEM::DMMoFEMKSPSetComputeRHS ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
MoFEM::BasicMethod *  pre_only,
MoFEM::BasicMethod *  post_only 
)

set KSP right hand side evaluation function

Examples
analytical_poisson.cpp, level_set.cpp, and simple_elasticity.cpp.

Definition at line 637 of file DMMoFEM.cpp.

640  {
641  return DMMoFEMKSPSetComputeRHS<const char *, MoFEM::FEMethod *,
642  MoFEM::BasicMethod *, MoFEM::BasicMethod *>(
643  dm, fe_name, method, pre_only, post_only);
644 }

◆ DMMoFEMKSPSetComputeRHS() [2/2]

PetscErrorCode MoFEM::DMMoFEMKSPSetComputeRHS ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
boost::shared_ptr< MoFEM::BasicMethod >  pre_only,
boost::shared_ptr< MoFEM::BasicMethod >  post_only 
)

set KSP right hand side evaluation function

Definition at line 647 of file DMMoFEM.cpp.

650  {
651  return DMMoFEMKSPSetComputeRHS<const std::string,
652  boost::shared_ptr<MoFEM::FEMethod>,
653  boost::shared_ptr<MoFEM::BasicMethod>,
654  boost::shared_ptr<MoFEM::BasicMethod>>(
655  dm, fe_name, method, pre_only, post_only);
656 }

◆ DMMoFEMResolveSharedFiniteElements()

PetscErrorCode MoFEM::DMMoFEMResolveSharedFiniteElements ( DM  dm,
std::string  fe_name 
)

Resolve shared entities.

Parameters
dmdm
fe_namefinite element for which shared entities are resolved
Returns
error code
Note
This function is valid for parallel algebra and serial mesh. It should be run collectively, i.e. on all processors.

This allows for tag reduction or tag exchange, f.e.

CHKERR DMMoFEMResolveSharedFiniteElements(dm,"SHELL_ELEMENT");
Tag th;
CHKERR mField.get_moab().tag_get_handle("ADAPT_ORDER",th);
ParallelComm* pcomm =
ParallelComm::get_pcomm(&mField.get_moab(),MYPCOMM_INDEX);
// CHKERR pcomm->reduce_tags(th,MPI_SUM,prisms);
CHKERR pcomm->exchange_tags(th,prisms);

Definition at line 464 of file DMMoFEM.cpp.

464  {
465  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
466  MoFEMFunctionBegin;
467  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
468  CHKERR dm_field->mField_ptr->getInterface<CommInterface>()
469  ->resolveSharedFiniteElements(dm_field->problemPtr, fe_name);
470  MoFEMFunctionReturn(0);
471 }

◆ DMMoFEMSetIsPartitioned()

PetscErrorCode MoFEM::DMMoFEMSetIsPartitioned ( DM  dm,
PetscBool  is_partitioned 
)

sets if read mesh is partitioned

get if read mesh is partitioned

Examples
cell_forces.cpp, dm_build_partitioned_mesh.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, EshelbianPlasticity.cpp, MagneticElement.hpp, mortar_contact.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, Remodeling.cpp, simple_contact.cpp, simple_contact_thermal.cpp, simple_elasticity.cpp, test_jacobian_of_simple_contact_element.cpp, and UnsaturatedFlow.hpp.

Definition at line 1123 of file DMMoFEM.cpp.

1123  {
1124  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1125  MoFEMFunctionBeginHot;
1126  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1127  dm_field->isPartitioned = is_partitioned;
1128  MoFEMFunctionReturnHot(0);
1129 }

◆ DMMoFEMSetKspCtx()

PetscErrorCode MoFEM::DMMoFEMSetKspCtx ( DM  dm,
boost::shared_ptr< MoFEM::KspCtx >  ksp_ctx 
)

set MoFEM::KspCtx data structure

Definition at line 1085 of file DMMoFEM.cpp.

1086  {
1087  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1088  MoFEMFunctionBeginHot;
1089  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1090  dm_field->kspCtx = ksp_ctx;
1091  MoFEMFunctionReturnHot(0);
1092 }

◆ DMMoFEMSetSnesCtx()

PetscErrorCode MoFEM::DMMoFEMSetSnesCtx ( DM  dm,
boost::shared_ptr< MoFEM::SnesCtx >  snes_ctx 
)

Set MoFEM::SnesCtx data structure.

Definition at line 1111 of file DMMoFEM.cpp.

1112  {
1113  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1114  MoFEMFunctionBeginHot;
1115  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1116  dm_field->snesCtx = snes_ctx;
1117  MoFEMFunctionReturnHot(0);
1118 }

◆ DMMoFEMSetSquareProblem()

PetscErrorCode MoFEM::DMMoFEMSetSquareProblem ( DM  dm,
PetscBool  square_problem 
)

set squared problem

It if true is assumed that matrix has the same indexing on rows and columns. This reduces interprocessor communication.

Examples
analytical_poisson_field_split.cpp, cell_forces.cpp, dm_build_partitioned_mesh.cpp, dm_create_subdm.cpp, dm_partitioned_no_field.cpp, EshelbianPlasticity.cpp, free_surface.cpp, level_set.cpp, plastic.cpp, schur_test_diag_mat.cpp, and test_broken_space.cpp.

Definition at line 456 of file DMMoFEM.cpp.

456  {
457  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
458  MoFEMFunctionBeginHot;
459  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
460  dm_field->isSquareMatrix = square_problem;
461  MoFEMFunctionReturnHot(0);
462 }

◆ DMMoFEMSetTsCtx()

PetscErrorCode MoFEM::DMMoFEMSetTsCtx ( DM  dm,
boost::shared_ptr< MoFEM::TsCtx >  ts_ctx 
)

Set MoFEM::TsCtx data structure.

It take over pointer, do not delete it, DM will destroy pointer when is destroyed.

Definition at line 1159 of file DMMoFEM.cpp.

1159  {
1160  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1161  MoFEMFunctionBeginHot;
1162  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1163  dm_field->tsCtx = ts_ctx;
1164  MoFEMFunctionReturnHot(0);
1165 }

◆ DMMoFEMSNESSetFunction() [1/2]

PetscErrorCode MoFEM::DMMoFEMSNESSetFunction ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
MoFEM::BasicMethod *  pre_only,
MoFEM::BasicMethod *  post_only 
)

set SNES residual evaluation function

Examples
analytical_nonlinear_poisson.cpp, mesh_smoothing.cpp, minimal_surface_area.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, NonlinearElasticElementInterface.hpp, simple_contact.cpp, simple_contact_thermal.cpp, test_jacobian_of_simple_contact_element.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 718 of file DMMoFEM.cpp.

721  {
722  return DMMoFEMSNESSetFunction<const char *, MoFEM::FEMethod *,
723  MoFEM::BasicMethod *, MoFEM::BasicMethod *>(
724  dm, fe_name, method, pre_only, post_only);
725 }

◆ DMMoFEMSNESSetFunction() [2/2]

PetscErrorCode MoFEM::DMMoFEMSNESSetFunction ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
boost::shared_ptr< MoFEM::BasicMethod >  pre_only,
boost::shared_ptr< MoFEM::BasicMethod >  post_only 
)

set SNES residual evaluation function

Definition at line 728 of file DMMoFEM.cpp.

731  {
732  return DMMoFEMSNESSetFunction<const std::string,
733  boost::shared_ptr<MoFEM::FEMethod>,
734  boost::shared_ptr<MoFEM::BasicMethod>,
735  boost::shared_ptr<MoFEM::BasicMethod>>(
736  dm, fe_name, method, pre_only, post_only);
737 }

◆ DMMoFEMSNESSetJacobian() [1/2]

PetscErrorCode MoFEM::DMMoFEMSNESSetJacobian ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
MoFEM::BasicMethod *  pre_only,
MoFEM::BasicMethod *  post_only 
)

set SNES Jacobian evaluation function

Examples
analytical_nonlinear_poisson.cpp, mesh_smoothing.cpp, minimal_surface_area.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, NonlinearElasticElementInterface.hpp, simple_contact.cpp, simple_contact_thermal.cpp, test_jacobian_of_simple_contact_element.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 759 of file DMMoFEM.cpp.

762  {
763  return DMMoFEMSNESSetJacobian<const char *, MoFEM::FEMethod *,
764  MoFEM::BasicMethod *, MoFEM::BasicMethod *>(
765  dm, fe_name, method, pre_only, post_only);
766 }

◆ DMMoFEMSNESSetJacobian() [2/2]

PetscErrorCode MoFEM::DMMoFEMSNESSetJacobian ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
boost::shared_ptr< MoFEM::BasicMethod >  pre_only,
boost::shared_ptr< MoFEM::BasicMethod >  post_only 
)

set SNES Jacobian evaluation function

Definition at line 769 of file DMMoFEM.cpp.

772  {
773  return DMMoFEMSNESSetJacobian<const std::string,
774  boost::shared_ptr<MoFEM::FEMethod>,
775  boost::shared_ptr<MoFEM::BasicMethod>,
776  boost::shared_ptr<MoFEM::BasicMethod>>(
777  dm, fe_name, method, pre_only, post_only);
778 }

◆ DMMoFEMTSSetI2Function() [1/2]

PetscErrorCode MoFEM::DMMoFEMTSSetI2Function ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
MoFEM::BasicMethod *  pre_only,
MoFEM::BasicMethod *  post_only 
)

set TS implicit function evaluation function

Definition at line 964 of file DMMoFEM.cpp.

967  {
968  return DMMoFEMTSSetI2Function<const char *, MoFEM::FEMethod *,
969  MoFEM::BasicMethod *, MoFEM::BasicMethod *>(
970  dm, fe_name, method, pre_only, post_only);
971  MoFEMFunctionReturnHot(0);
972 }

◆ DMMoFEMTSSetI2Function() [2/2]

PetscErrorCode MoFEM::DMMoFEMTSSetI2Function ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
boost::shared_ptr< MoFEM::BasicMethod >  pre_only,
boost::shared_ptr< MoFEM::BasicMethod >  post_only 
)

set TS implicit function evaluation function

Examples
EshelbianPlasticity.cpp, and NonlinearElasticElementInterface.hpp.

Definition at line 975 of file DMMoFEM.cpp.

978  {
979  return DMMoFEMTSSetI2Function<const std::string,
980  boost::shared_ptr<MoFEM::FEMethod>,
981  boost::shared_ptr<MoFEM::BasicMethod>,
982  boost::shared_ptr<MoFEM::BasicMethod>>(
983  dm, fe_name, method, pre_only, post_only);
984  MoFEMFunctionReturnHot(0);
985 }

◆ DMMoFEMTSSetI2Jacobian() [1/2]

PetscErrorCode MoFEM::DMMoFEMTSSetI2Jacobian ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
MoFEM::BasicMethod *  pre_only,
MoFEM::BasicMethod *  post_only 
)

set TS Jacobian evaluation function

Definition at line 1007 of file DMMoFEM.cpp.

1010  {
1011  return DMMoFEMTSSetI2Jacobian<const char *, FEMethod *, MoFEM::BasicMethod *,
1012  MoFEM::BasicMethod *>(dm, fe_name, method,
1013  pre_only, post_only);
1014 }

◆ DMMoFEMTSSetI2Jacobian() [2/2]

PetscErrorCode MoFEM::DMMoFEMTSSetI2Jacobian ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
boost::shared_ptr< MoFEM::BasicMethod >  pre_only,
boost::shared_ptr< MoFEM::BasicMethod >  post_only 
)

set TS Jacobian evaluation function

Examples
EshelbianPlasticity.cpp, and NonlinearElasticElementInterface.hpp.

Definition at line 1017 of file DMMoFEM.cpp.

1020  {
1021  return DMMoFEMTSSetI2Jacobian<const std::string,
1022  boost::shared_ptr<MoFEM::FEMethod>,
1023  boost::shared_ptr<MoFEM::BasicMethod>,
1024  boost::shared_ptr<MoFEM::BasicMethod>>(
1025  dm, fe_name, method, pre_only, post_only);
1026 }

◆ DMMoFEMTSSetIFunction() [1/2]

PetscErrorCode MoFEM::DMMoFEMTSSetIFunction ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
MoFEM::BasicMethod *  pre_only,
MoFEM::BasicMethod *  post_only 
)

set TS implicit function evaluation function

Examples
bone_adaptation.cpp, EshelbianPlasticity.cpp, NonlinearElasticElementInterface.hpp, reaction_diffusion.cpp, and UnsaturatedFlow.hpp.

Definition at line 800 of file DMMoFEM.cpp.

803  {
804  return DMMoFEMTSSetIFunction<const char *, MoFEM::FEMethod *,
805  MoFEM::BasicMethod *, MoFEM::BasicMethod *>(
806  dm, fe_name, method, pre_only, post_only);
807  MoFEMFunctionReturnHot(0);
808 }

◆ DMMoFEMTSSetIFunction() [2/2]

PetscErrorCode MoFEM::DMMoFEMTSSetIFunction ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
boost::shared_ptr< MoFEM::BasicMethod >  pre_only,
boost::shared_ptr< MoFEM::BasicMethod >  post_only 
)

set TS implicit function evaluation function

Definition at line 811 of file DMMoFEM.cpp.

814  {
815  return DMMoFEMTSSetIFunction<const std::string,
816  boost::shared_ptr<MoFEM::FEMethod>,
817  boost::shared_ptr<MoFEM::BasicMethod>,
818  boost::shared_ptr<MoFEM::BasicMethod>>(
819  dm, fe_name, method, pre_only, post_only);
820  MoFEMFunctionReturnHot(0);
821 }

◆ DMMoFEMTSSetIJacobian() [1/2]

PetscErrorCode MoFEM::DMMoFEMTSSetIJacobian ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
MoFEM::BasicMethod *  pre_only,
MoFEM::BasicMethod *  post_only 
)

set TS Jacobian evaluation function

Definition at line 843 of file DMMoFEM.cpp.

846  {
847  return DMMoFEMTSSetIJacobian<const char *, FEMethod *, MoFEM::BasicMethod *,
848  MoFEM::BasicMethod *>(dm, fe_name, method,
849  pre_only, post_only);
850 }

◆ DMMoFEMTSSetIJacobian() [2/2]

PetscErrorCode MoFEM::DMMoFEMTSSetIJacobian ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
boost::shared_ptr< MoFEM::BasicMethod >  pre_only,
boost::shared_ptr< MoFEM::BasicMethod >  post_only 
)

set TS Jacobian evaluation function

Examples
bone_adaptation.cpp, EshelbianPlasticity.cpp, NonlinearElasticElementInterface.hpp, reaction_diffusion.cpp, and UnsaturatedFlow.hpp.

Definition at line 853 of file DMMoFEM.cpp.

856  {
857  return DMMoFEMTSSetIJacobian<const std::string,
858  boost::shared_ptr<MoFEM::FEMethod>,
859  boost::shared_ptr<MoFEM::BasicMethod>,
860  boost::shared_ptr<MoFEM::BasicMethod>>(
861  dm, fe_name, method, pre_only, post_only);
862 }

◆ DMMoFEMUnSetElement()

PetscErrorCode MoFEM::DMMoFEMUnSetElement ( DM  dm,
std::string  fe_name 
)

unset element from dm

Definition at line 514 of file DMMoFEM.cpp.

514  {
515  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
516  MoFEMFunctionBegin;
517  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
518  CHKERR dm_field->mField_ptr->modify_problem_unset_finite_element(
519  dm_field->problemName, fe_name);
520  MoFEMFunctionReturn(0);
521 }

◆ DMoFEMGetInterfacePtr()

PetscErrorCode MoFEM::DMoFEMGetInterfacePtr ( DM  dm,
MoFEM::Interface **  m_field_ptr 
)

Get pointer to MoFEM::Interface.

Parameters
dmDistributed mesh manager
m_field_ptrPointer to pointer of field interface
Returns
Error code
Examples
adolc_plasticity.cpp, free_surface.cpp, HookeElement.cpp, mesh_smoothing.cpp, and PlasticOpsMonitor.hpp.

Definition at line 414 of file DMMoFEM.cpp.

414  {
415  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
416  MoFEMFunctionBeginHot;
417  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
418  if (!dm->data) {
419  SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
420  "data structure for MoFEM not yet created");
421  }
422  *m_field_ptr = dm_field->mField_ptr;
423  MoFEMFunctionReturnHot(0);
424 }

◆ DMoFEMLoopDofs()

PetscErrorCode MoFEM::DMoFEMLoopDofs ( DM  dm,
const char  field_name[],
MoFEM::DofMethod *  method 
)

execute method for dofs on field in problem

Definition at line 605 of file DMMoFEM.cpp.

606  {
607  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
608  MoFEMFunctionBeginHot;
609  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
610  ierr =
611  dm_field->mField_ptr->loop_dofs(dm_field->problemPtr, field_name, COL,
612  *method, dm_field->rAnk, dm_field->rAnk);
613  CHKERRG(ierr);
614  MoFEMFunctionReturnHot(0);
615 }

◆ DMoFEMLoopFiniteElements() [1/2]

PetscErrorCode MoFEM::DMoFEMLoopFiniteElements ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
CacheTupleWeakPtr  cache_ptr = CacheTupleSharedPtr() 
)

Executes FEMethod for finite elements in DM.

Parameters
dmMoFEM discrete manager
fe_namename of element
methodpointer to MOFEM::FEMethod
Returns
Error code
Examples
adolc_plasticity.cpp, analytical_nonlinear_poisson.cpp, analytical_poisson.cpp, analytical_poisson_field_split.cpp, approx_sphere.cpp, boundary_marker.cpp, cell_forces.cpp, continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, dm_build_partitioned_mesh.cpp, dynamic_first_order_con_law.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, EshelbianPlasticity.cpp, free_surface.cpp, heat_equation.cpp, HookeElement.cpp, level_set.cpp, MagneticElement.hpp, mesh_smoothing.cpp, minimal_surface_area.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, nonlinear_elastic.cpp, NonlinearElasticElementInterface.hpp, operators_tests.cpp, photon_diffusion.cpp, PlasticOpsMonitor.hpp, reaction_diffusion.cpp, Remodeling.cpp, scalar_check_approximation.cpp, schur_test_diag_mat.cpp, seepage.cpp, shallow_wave.cpp, simple_contact.cpp, simple_contact_thermal.cpp, simple_elasticity.cpp, simple_interface.cpp, tensor_divergence_operator.cpp, test_broken_space.cpp, test_cache_on_entities.cpp, thermo_elastic.cpp, UnsaturatedFlow.hpp, and wave_equation.cpp.

Definition at line 586 of file DMMoFEM.cpp.

588  {
589  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
590  MoFEMFunctionBeginHot;
591  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
592  ierr = DMoFEMLoopFiniteElementsUpAndLowRank(
593  dm, fe_name, method, dm_field->rAnk, dm_field->rAnk, cache_ptr);
594  CHKERRG(ierr);
595  MoFEMFunctionReturnHot(0);
596 }

◆ DMoFEMLoopFiniteElements() [2/2]

PetscErrorCode MoFEM::DMoFEMLoopFiniteElements ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
CacheTupleWeakPtr  cache_ptr = CacheTupleSharedPtr() 
)

Executes FEMethod for finite elements in DM.

Parameters
dmMoFEM discrete manager
fe_namename of element
methodpointer to MOFEM::FEMethod
Returns
Error code

Definition at line 599 of file DMMoFEM.cpp.

601  {
602  return DMoFEMLoopFiniteElements(dm, fe_name.c_str(), method.get(), cache_ptr);
603 }

◆ DMoFEMLoopFiniteElementsUpAndLowRank() [1/2]

PetscErrorCode MoFEM::DMoFEMLoopFiniteElementsUpAndLowRank ( DM  dm,
const char  fe_name[],
MoFEM::FEMethod *  method,
int  low_rank,
int  up_rank,
CacheTupleWeakPtr  cache_ptr = CacheTupleSharedPtr() 
)

Executes FEMethod for finite elements in DM.

Parameters
dmMoFEM discrete manager
fe_namename of finite element
methodpointer to MoFEM::FEMethod
low_ranklowest rank of processor
up_rankupper run of processor
Returns
Error code
Examples
mortar_contact_thermal.cpp, and simple_contact_thermal.cpp.

Definition at line 567 of file DMMoFEM.cpp.

569  {
570  MoFEMFunctionBeginHot;
571  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
572  ierr = dm_field->mField_ptr->loop_finite_elements(
573  dm_field->problemPtr, fe_name, *method, low_rank, up_rank, nullptr,
574  MF_EXIST, cache_ptr);
575  CHKERRG(ierr);
576  MoFEMFunctionReturnHot(0);
577 }

◆ DMoFEMLoopFiniteElementsUpAndLowRank() [2/2]

PetscErrorCode MoFEM::DMoFEMLoopFiniteElementsUpAndLowRank ( DM  dm,
const std::string  fe_name,
boost::shared_ptr< MoFEM::FEMethod >  method,
int  low_rank,
int  up_rank,
CacheTupleWeakPtr  cache_ptr = CacheTupleSharedPtr() 
)

Executes FEMethod for finite elements in DM.

Parameters
dmMoFEM discrete manager
fe_namename of finite element
methodpointer to MoFEM::FEMethod
low_ranklowest rank of processor
up_rankupper run of processor
Returns
Error code

Definition at line 579 of file DMMoFEM.cpp.

581  {
582  return DMoFEMLoopFiniteElementsUpAndLowRank(dm, fe_name.c_str(), method.get(),
583  low_rank, up_rank, cache_ptr);
584 }

◆ DMoFEMMeshToGlobalVector()

PetscErrorCode MoFEM::DMoFEMMeshToGlobalVector ( DM  dm,
Vec  g,
InsertMode  mode,
ScatterMode  scatter_mode 
)

set ghosted vector values on all existing mesh entities

Parameters
gvector
modesee petsc manual for VecSetValue (ADD_VALUES or INSERT_VALUES)
scatter_modesee petsc manual for ScatterMode (The available modes are: SCATTER_FORWARD or SCATTER_REVERSE)

SCATTER_REVERSE set data to field entities from V vector.

SCATTER_FORWARD set vector V from data field entities

Examples
analytical_nonlinear_poisson.cpp, analytical_poisson.cpp, analytical_poisson_field_split.cpp, elasticity_mixed_formulation.cpp, mesh_smoothing.cpp, minimal_surface_area.cpp, mortar_contact.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, simple_contact.cpp, simple_contact_thermal.cpp, and simple_elasticity.cpp.

Definition at line 535 of file DMMoFEM.cpp.

536  {
537  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
538  MoFEMFunctionBeginHot;
539  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
540  ierr = dm_field->mField_ptr->getInterface<VecManager>()->setGlobalGhostVector(
541  dm_field->problemPtr, COL, g, mode, scatter_mode);
542  CHKERRG(ierr);
543  MoFEMFunctionReturnHot(0);
544 }

◆ DMoFEMMeshToLocalVector()

PetscErrorCode MoFEM::DMoFEMMeshToLocalVector ( DM  dm,
Vec  l,
InsertMode  mode,
ScatterMode  scatter_mode 
)

set local (or ghosted) vector values on mesh for partition only

Parameters
lvector
modesee petsc manual for VecSetValue (ADD_VALUES or INSERT_VALUES)
scatter_modesee petsc manual for ScatterMode (The available modes are: SCATTER_FORWARD or SCATTER_REVERSE)

SCATTER_REVERSE set data to field entities from V vector.

SCATTER_FORWARD set vector V from data field entities

Examples
adolc_plasticity.cpp, approx_sphere.cpp, bone_adaptation.cpp, cell_forces.cpp, child_and_parent.cpp, dg_projection.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, ep.cpp, EshelbianPlasticity.cpp, free_surface.cpp, hanging_node_approx.cpp, hcurl_check_approx_in_2d.cpp, heat_equation.cpp, heat_method.cpp, helmholtz.cpp, higher_derivatives.cpp, level_set.cpp, MagneticElement.hpp, mesh_smoothing.cpp, minimal_surface_area.cpp, mixed_poisson.cpp, mortar_contact_thermal.cpp, nonlinear_elastic.cpp, phase.cpp, photon_diffusion.cpp, plate.cpp, poisson_2d_dis_galerkin.cpp, poisson_2d_homogeneous.cpp, reaction_diffusion.cpp, scalar_check_approximation.cpp, shallow_wave.cpp, simple_contact.cpp, simple_contact_thermal.cpp, tensor_divergence_operator.cpp, test_broken_space.cpp, test_cache_on_entities.cpp, test_jacobian_of_simple_contact_element.cpp, testing_jacobian_of_hook_element.cpp, testing_jacobian_of_hook_scaled_with_density_element.cpp, thermo_elastic.cpp, UnsaturatedFlow.hpp, and wave_equation.cpp.

Definition at line 523 of file DMMoFEM.cpp.

524  {
525  MoFEMFunctionBeginHot;
526  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
527  MoFEMFunctionBeginHot;
528  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
529  ierr = dm_field->mField_ptr->getInterface<VecManager>()->setLocalGhostVector(
530  dm_field->problemPtr, COL, l, mode, scatter_mode);
531  CHKERRG(ierr);
532  MoFEMFunctionReturnHot(0);
533 }

◆ DMoFEMPostProcessFiniteElements()

PetscErrorCode MoFEM::DMoFEMPostProcessFiniteElements ( DM  dm,
MoFEM::FEMethod *  method 
)

execute finite element method for each element in dm (problem)

Examples
cell_forces.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, free_surface.cpp, scalar_check_approximation.cpp, and simple_elasticity.cpp.

Definition at line 556 of file DMMoFEM.cpp.

556  {
557  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
558  MoFEMFunctionBeginHot;
559  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
560  ierr = dm_field->mField_ptr->problem_basic_method_postProcess(
561  dm_field->problemPtr, *method);
562  CHKERRG(ierr);
563  MoFEMFunctionReturnHot(0);
564 }

◆ DMoFEMPreProcessFiniteElements()

PetscErrorCode MoFEM::DMoFEMPreProcessFiniteElements ( DM  dm,
MoFEM::FEMethod *  method 
)

execute finite element method for each element in dm (problem)

Examples
cell_forces.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, free_surface.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, scalar_check_approximation.cpp, simple_contact.cpp, and simple_contact_thermal.cpp.

Definition at line 546 of file DMMoFEM.cpp.

546  {
547  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
548  MoFEMFunctionBeginHot;
549  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
550  ierr = dm_field->mField_ptr->problem_basic_method_preProcess(
551  dm_field->problemPtr, *method);
552  CHKERRG(ierr);
553  MoFEMFunctionReturnHot(0);
554 }

◆ DMRegister_MGViaApproxOrders() [1/2]

MoFEMErrorCode DMRegister_MGViaApproxOrders ( const char  sname[])

Register DM for Multi-Grid via approximation orders.

Parameters
snameproblem/dm registered name
Returns
error code

Definition at line 302 of file PCMGSetUpViaApproxOrders.cpp.

302  {
303  MoFEMFunctionBegin;
304  CHKERR DMRegister(sname, DMCreate_MGViaApproxOrders);
305  MoFEMFunctionReturn(0);
306 }

◆ DMRegister_MGViaApproxOrders() [2/2]

MoFEMErrorCode MoFEM::DMRegister_MGViaApproxOrders ( const char  sname[])

Register DM for Multi-Grid via approximation orders.

Parameters
snameproblem/dm registered name
Returns
error code
Examples
elasticity.cpp, and test_broken_space.cpp.

Definition at line 302 of file PCMGSetUpViaApproxOrders.cpp.

302  {
303  MoFEMFunctionBegin;
304  CHKERR DMRegister(sname, DMCreate_MGViaApproxOrders);
305  MoFEMFunctionReturn(0);
306 }

◆ DMRegister_MoFEM()

PetscErrorCode MoFEM::DMRegister_MoFEM ( const char  sname[])

Register MoFEM problem.

Examples
adolc_plasticity.cpp, analytical_nonlinear_poisson.cpp, analytical_poisson.cpp, analytical_poisson_field_split.cpp, approx_sphere.cpp, boundary_marker.cpp, cell_forces.cpp, child_and_parent.cpp, continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, dg_projection.cpp, dm_build_partitioned_mesh.cpp, dm_create_subdm.cpp, dm_partitioned_no_field.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, elasticity_mixed_formulation.cpp, ep.cpp, field_evaluator.cpp, free_surface.cpp, hanging_node_approx.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, heat_equation.cpp, heat_method.cpp, hello_world.cpp, helmholtz.cpp, higher_derivatives.cpp, level_set.cpp, loop_entities.cpp, lorentz_force.cpp, MagneticElement.hpp, mesh_smoothing.cpp, minimal_surface_area.cpp, mixed_poisson.cpp, mortar_contact.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, nonlinear_elastic.cpp, operators_tests.cpp, partition_mesh.cpp, phase.cpp, photon_diffusion.cpp, plastic.cpp, plate.cpp, plot_base.cpp, poisson_2d_dis_galerkin.cpp, poisson_2d_homogeneous.cpp, reaction_diffusion.cpp, Remodeling.cpp, scalar_check_approximation.cpp, schur_test_diag_mat.cpp, seepage.cpp, shallow_wave.cpp, simple_contact.cpp, simple_contact_thermal.cpp, simple_elasticity.cpp, simple_interface.cpp, simple_l2_only.cpp, tensor_divergence_operator.cpp, test_broken_space.cpp, test_cache_on_entities.cpp, test_jacobian_of_simple_contact_element.cpp, testing_jacobian_of_hook_element.cpp, testing_jacobian_of_hook_scaled_with_density_element.cpp, thermo_elastic.cpp, unsaturated_transport.cpp, and wave_equation.cpp.

Definition at line 43 of file DMMoFEM.cpp.

43  {
44  MoFEMFunctionBegin;
45  CHKERR DMRegister(sname, DMCreate_MoFEM);
46  MoFEMFunctionReturn(0);
47 }

◆ DMSetFromOptions_MoFEM()

PetscErrorCode MoFEM::DMSetFromOptions_MoFEM ( DM  dm)

Set options for MoFEM DM

Definition at line 1281 of file DMMoFEM.cpp.

1281  {
1282 #endif
1283 
1284  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1285  MoFEMFunctionBeginHot;
1286  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1287 #if PETSC_VERSION_GE(3, 5, 3)
1288  ierr = PetscOptionsHead(PetscOptionsObject, "DMMoFEM Options");
1289  CHKERRG(ierr);
1290 #else
1291  ierr = PetscOptionsHead("DMMoFEM Options");
1292  CHKERRG(ierr);
1293 #endif
1294  ierr = PetscOptionsBool("-dm_is_partitioned",
1295  "set if mesh is partitioned (works which native MOAB "
1296  "file format, i.e. h5m",
1297  "DMSetUp", dm_field->isPartitioned,
1298  &dm_field->isPartitioned, NULL);
1299  CHKERRG(ierr);
1300  MoFEMFunctionReturnHot(0);
1301 }

◆ DMSetOperators_MoFEM()

PetscErrorCode MoFEM::DMSetOperators_MoFEM ( DM  dm)

Set operators for MoFEM dm.

Parameters
dm
Returns
error code

Definition at line 49 of file DMMoFEM.cpp.

49  {
50  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
51  MoFEMFunctionBegin;
52 
53  dm->ops->createglobalvector = DMCreateGlobalVector_MoFEM;
54  dm->ops->createlocalvector = DMCreateLocalVector_MoFEM;
55  dm->ops->creatematrix = DMCreateMatrix_MoFEM;
56  dm->ops->setup = DMSetUp_MoFEM;
57  dm->ops->destroy = DMDestroy_MoFEM;
58  dm->ops->setfromoptions = DMSetFromOptions_MoFEM;
59  dm->ops->globaltolocalbegin = DMGlobalToLocalBegin_MoFEM;
60  dm->ops->globaltolocalend = DMGlobalToLocalEnd_MoFEM;
61  dm->ops->localtoglobalbegin = DMLocalToGlobalBegin_MoFEM;
62  dm->ops->localtoglobalend = DMLocalToGlobalEnd_MoFEM;
63  dm->ops->createfieldis = DMCreateFieldIS_MoFEM;
64 
65  // Default matrix type
66  CHKERR DMSetMatType(dm, MATMPIAIJ);
67 
68  MoFEMFunctionReturn(0);
69 }

◆ DMSetUp_MoFEM()

PetscErrorCode MoFEM::DMSetUp_MoFEM ( DM  dm)

sets up the MoFEM structures inside a DM object

Examples
mixed_poisson.cpp.

Definition at line 1303 of file DMMoFEM.cpp.

1303  {
1304  PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1305  ProblemsManager *prb_mng_ptr;
1306  MoFEMFunctionBegin;
1307  DMCtxImpl *dm_field = static_cast<DMCtxImpl *>(dm->data);
1308  CHKERR dm_field->mField_ptr->getInterface(prb_mng_ptr);
1309 
1310  if (dm_field->isCompDM) {
1311  // It is composite probelm
1312  CHKERR prb_mng_ptr->buildComposedProblem(
1313  dm_field->problemName, dm_field->rowCompPrb, dm_field->colCompPrb,
1314  dm_field->isSquareMatrix == PETSC_TRUE, dm_field->verbosity);
1315  } else {
1316  if (dm_field->isPartitioned) {
1317  CHKERR prb_mng_ptr->buildProblemOnDistributedMesh(
1318  dm_field->problemName, dm_field->isSquareMatrix == PETSC_TRUE,
1319  dm_field->verbosity);
1320  } else {
1321  CHKERR prb_mng_ptr->buildProblem(dm_field->problemName,
1322  dm_field->isSquareMatrix == PETSC_TRUE,
1323  dm_field->verbosity);
1324  CHKERR prb_mng_ptr->partitionProblem(dm_field->problemName,
1325  dm_field->verbosity);
1326  }
1327  }
1328 
1329  // Partition finite elements
1330  if (dm_field->isPartitioned) {
1331  CHKERR prb_mng_ptr->partitionFiniteElements(
1332  dm_field->problemName, true, 0, dm_field->sIze, dm_field->verbosity);
1333  CHKERR prb_mng_ptr->partitionGhostDofsOnDistributedMesh(
1334  dm_field->problemName, dm_field->verbosity);
1335  } else {
1336  // partition finite elemnets
1337  CHKERR prb_mng_ptr->partitionFiniteElements(dm_field->problemName, false,
1338  -1, -1, dm_field->verbosity);
1339  // Get ghost DOFs
1340  CHKERR prb_mng_ptr->partitionGhostDofs(dm_field->problemName,
1341  dm_field->verbosity);
1342  }
1343 
1344  // Set flag that problem is build and partitioned
1345  dm_field->isProblemBuild = PETSC_TRUE;
1346 
1347  MoFEMFunctionReturn(0);
1348 }

◆ DMSubDMSetUp_MoFEM()

PetscErrorCode MoFEM::DMSubDMSetUp_MoFEM ( DM  subdm)

Sets up the MoFEM structures inside a DM object for sub dm

Examples
level_set.cpp.

Definition at line 1350 of file DMMoFEM.cpp.

1350  {
1351  PetscValidHeaderSpecific(subdm, DM_CLASSID, 1);
1352  ProblemsManager *prb_mng_ptr;
1353  MoFEMFunctionBegin;
1354 
1355  DMCtxImpl *subdm_field = static_cast<DMCtxImpl *>(subdm->data);
1356 
1357  // build sub dm problem
1358  CHKERR subdm_field->mField_ptr->getInterface(prb_mng_ptr);
1359 
1360  map<std::string, boost::shared_ptr<Range>> *entity_map_row = nullptr;
1361  map<std::string, boost::shared_ptr<Range>> *entity_map_col = nullptr;
1362 
1363  if (subdm_field->mapTypeRow.size())
1364  entity_map_row = &subdm_field->mapTypeRow;
1365  if (subdm_field->mapTypeCol.size())
1366  entity_map_col = &subdm_field->mapTypeCol;
1367 
1368  CHKERR prb_mng_ptr->buildSubProblem(
1369  subdm_field->problemName, subdm_field->rowSubFields,
1370  subdm_field->colSubFields, subdm_field->problemMainOfSubPtr->getName(),
1371  subdm_field->isSquareMatrix == PETSC_TRUE, entity_map_row, entity_map_col,
1372  subdm_field->verbosity);
1373 
1374  // partition problem
1375  subdm_field->isPartitioned = subdm_field->isPartitioned;
1376  if (subdm_field->isPartitioned) {
1377  CHKERR prb_mng_ptr->partitionFiniteElements(subdm_field->problemName, true,
1378  0, subdm_field->sIze,
1379  subdm_field->verbosity);
1380  // set ghost nodes
1381  CHKERR prb_mng_ptr->partitionGhostDofsOnDistributedMesh(
1382  subdm_field->problemName, subdm_field->verbosity);
1383  } else {
1384  // partition finite elements
1385  CHKERR prb_mng_ptr->partitionFiniteElements(subdm_field->problemName, false,
1386  -1, -1, subdm_field->verbosity);
1387  // set ghost nodes
1388  CHKERR prb_mng_ptr->partitionGhostDofs(subdm_field->problemName,
1389  subdm_field->verbosity);
1390  }
1391 
1392  subdm_field->isProblemBuild = PETSC_TRUE;
1393 
1394  MoFEMFunctionReturn(0);
1395 }
MoFEMFunctionReturnHot
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:460
MoFEM::DMSetFromOptions_MoFEM
PetscErrorCode DMSetFromOptions_MoFEM(DM dm)
Definition: DMMoFEM.cpp:1281
g
constexpr double g
Definition: shallow_wave.cpp:63
MYPCOMM_INDEX
#define MYPCOMM_INDEX
default communicator number PCOMM
Definition: definitions.h:228
MoFEM::DMMoFEMSNESSetJacobian
static PetscErrorCode DMMoFEMSNESSetJacobian(DM dm, S fe_name, T0 method, T1 pre_only, T2 post_only)
Definition: DMMoFEM.cpp:740
MoFEM::FEMethod
structure for User Loop Methods on finite elements
Definition: LoopMethods.hpp:369
MoFEM::CoreInterface::get_comm
virtual MPI_Comm & get_comm() const =0
MoFEM::DMCreateGlobalVector_MoFEM
PetscErrorCode DMCreateGlobalVector_MoFEM(DM dm, SmartPetscObj< Vec > &g_ptr)
DMShellSetCreateGlobalVector.
Definition: DMMoFEM.cpp:1177
MoFEM::DMMoFEMTSSetI2Function
static PetscErrorCode DMMoFEMTSSetI2Function(DM dm, S fe_name, T0 method, T1 pre_only, T2 post_only)
Definition: DMMoFEM.cpp:945
MoFEM::DMGlobalToLocalBegin_MoFEM
PetscErrorCode DMGlobalToLocalBegin_MoFEM(DM dm, Vec, InsertMode, Vec)
Definition: DMMoFEM.cpp:1397
MoFEM::th
Tag th
Definition: Projection10NodeCoordsOnField.cpp:122
MoFEM::DMMoFEMTSSetIJacobian
static PetscErrorCode DMMoFEMTSSetIJacobian(DM dm, S fe_name, T0 method, T1 pre_only, T2 post_only)
Definition: DMMoFEM.cpp:824
MoFEM::DMSetUp_MoFEM
PetscErrorCode DMSetUp_MoFEM(DM dm)
Definition: DMMoFEM.cpp:1303
A
constexpr AssemblyType A
Definition: operators_tests.cpp:30
ts_ctx
MoFEM::TsCtx * ts_ctx
Definition: level_set.cpp:1932
MoFEM::DMLocalToGlobalBegin_MoFEM
PetscErrorCode DMLocalToGlobalBegin_MoFEM(DM, Vec, InsertMode, Vec)
Definition: DMMoFEM.cpp:1434
GET_DM_FIELD
#define GET_DM_FIELD(DM)
Definition: PCMGSetUpViaApproxOrders.cpp:214
DMCreate_MGViaApproxOrders
MoFEMErrorCode DMCreate_MGViaApproxOrders(DM dm)
Create DM data structure for Multi-Grid via approximation orders.
Definition: PCMGSetUpViaApproxOrders.cpp:408
MoFEM::DMMoFEMAddElement
PetscErrorCode DMMoFEMAddElement(DM dm, std::string fe_name)
add element to dm
Definition: DMMoFEM.cpp:497
MoFEM::DMMoFEMTSSetI2Jacobian
static PetscErrorCode DMMoFEMTSSetI2Jacobian(DM dm, S fe_name, T0 method, T1 pre_only, T2 post_only)
Definition: DMMoFEM.cpp:988
MoFEM::CoreInterface::check_problem
virtual bool check_problem(const std::string name)=0
check if problem exist
ROW
@ ROW
Definition: definitions.h:136
MoFEM::DMCreateGlobalVector_MoFEM
PetscErrorCode DMCreateGlobalVector_MoFEM(DM dm, Vec *g)
DMShellSetCreateGlobalVector.
Definition: DMMoFEM.cpp:1167
PlasticOps::M
FTensor::Index< 'M', 3 > M
Definition: PlasticOps.hpp:117
MoFEM::DMCreateMatrix_MoFEM
PetscErrorCode DMCreateMatrix_MoFEM(DM dm, Mat *M)
Definition: DMMoFEM.cpp:1197
FEMethod
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:548
MoFEM::createBlockMat
SchurShellMatData createBlockMat(DM dm, boost::shared_ptr< BlockStructure > data)
Create a Schur Mat object.
Definition: Schur.cpp:1503
MoFEM::DMoFEMLoopFiniteElementsUpAndLowRank
PetscErrorCode DMoFEMLoopFiniteElementsUpAndLowRank(DM dm, const char fe_name[], MoFEM::FEMethod *method, int low_rank, int up_rank, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
Executes FEMethod for finite elements in DM.
Definition: DMMoFEM.cpp:567
a
constexpr double a
Definition: approx_sphere.cpp:30
MoFEM::DMDestroy_MoFEM
PetscErrorCode DMDestroy_MoFEM(DM dm)
Destroys dm with MoFEM data structure.
Definition: DMMoFEM.cpp:79
MoFEM::DMSubDMSetUp_MoFEM
PetscErrorCode DMSubDMSetUp_MoFEM(DM subdm)
Definition: DMMoFEM.cpp:1350
sub_mat_mult
MoFEMErrorCode sub_mat_mult(Mat a, Vec x, Vec f)
Definition: PCMGSetUpViaApproxOrders.cpp:106
convert.type
type
Definition: convert.py:64
COL
@ COL
Definition: definitions.h:136
PCMGSubMatrixCtx_private
Definition: PCMGSetUpViaApproxOrders.cpp:48
PCMGSubMatrixCtx
Definition: PCMGSetUpViaApproxOrders.hpp:9
MoFEM::DMMoFEMTSSetIFunction
static PetscErrorCode DMMoFEMTSSetIFunction(DM dm, S fe_name, T0 method, T1 pre_only, T2 post_only)
Definition: DMMoFEM.cpp:781
MoFEM::DMGlobalToLocalEnd_MoFEM
PetscErrorCode DMGlobalToLocalEnd_MoFEM(DM dm, Vec, InsertMode, Vec)
Definition: DMMoFEM.cpp:1405
MoFEM::BasicMethod
Data structure to exchange data between mofem and User Loop Methods.
Definition: LoopMethods.hpp:183
MoFEM::DMMoFEMCreateMoFEM
PetscErrorCode DMMoFEMCreateMoFEM(DM dm, MoFEM::Interface *m_field_ptr, const char problem_name[], const MoFEM::BitRefLevel bit_level, const MoFEM::BitRefLevel bit_mask=MoFEM::BitRefLevel().set())
Must be called by user to set MoFEM data structures.
Definition: DMMoFEM.cpp:114
MoFEM::isDuplicate
auto isDuplicate(IS is)
Definition: PetscSmartObj.hpp:396
MoFEM::DMMoFEMKSPSetComputeRHS
static PetscErrorCode DMMoFEMKSPSetComputeRHS(DM dm, S fe_name, T0 method, T1 pre_only, T2 post_only)
Definition: DMMoFEM.cpp:618
MoFEM::DMSetOperators_MoFEM
PetscErrorCode DMSetOperators_MoFEM(DM dm)
Set operators for MoFEM dm.
Definition: DMMoFEM.cpp:49
MoFEM::DMCreateFieldIS_MoFEM
PetscErrorCode DMCreateFieldIS_MoFEM(DM dm, PetscInt *numFields, char ***fieldNames, IS **fields)
Definition: DMMoFEM.cpp:1469
field_name
constexpr auto field_name
Definition: poisson_2d_homogeneous.cpp:13
MoFEM::DMCreateMatrix_MoFEM
PetscErrorCode DMCreateMatrix_MoFEM(DM dm, SmartPetscObj< Mat > &M)
Definition: DMMoFEM.cpp:1236
MoFEM::DMMoFEMResolveSharedFiniteElements
PetscErrorCode DMMoFEMResolveSharedFiniteElements(DM dm, std::string fe_name)
Resolve shared entities.
Definition: DMMoFEM.cpp:464
convert.n
n
Definition: convert.py:82
sub_mat_mult_add
MoFEMErrorCode sub_mat_mult_add(Mat a, Vec x, Vec f)
Definition: PCMGSetUpViaApproxOrders.cpp:110
N
const int N
Definition: speed_test.cpp:3
v
const double v
phase velocity of light in medium (cm/ns)
Definition: initial_diffusion.cpp:40
SnesCtx
MoFEM::DMCreate_MoFEM
PetscErrorCode DMCreate_MoFEM(DM dm)
Create dm data structure with MoFEM data structure.
Definition: DMMoFEM.cpp:71
MOFEM_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition: LogManager.hpp:308
MoFEM::DMMoFEMGetFieldIS
PetscErrorCode DMMoFEMGetFieldIS(DM dm, RowColData rc, const char field_name[], IS *is)
get field is in the problem
Definition: DMMoFEM.cpp:1514
HenckyOps::f
auto f
Definition: HenckyOps.hpp:15
MoFEM::Exceptions::ierr
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
Definition: Exceptions.hpp:76
MOFEM_DATA_INCONSISTENCY
@ MOFEM_DATA_INCONSISTENCY
Definition: definitions.h:31
MoFEM::getCommFromPetscObject
MPI_Comm getCommFromPetscObject(PetscObject obj)
Get the Comm From Petsc Object object.
Definition: PetscSmartObj.hpp:160
MoFEM::DMMoFEMCreateBlockMat
MoFEMErrorCode DMMoFEMCreateBlockMat(DM dm, Mat *mat)
Create block matrix.
Definition: DMMoFEM.cpp:1551
MoFEM::DMMoFEMKSPSetComputeOperators
static PetscErrorCode DMMoFEMKSPSetComputeOperators(DM dm, S fe_name, T0 method, T1 pre_only, T2 post_only)
Definition: DMMoFEM.cpp:659
MoFEMFunctionBeginHot
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:453
MF_EXCL
@ MF_EXCL
Definition: definitions.h:112
MoFEM::DMLocalToGlobalEnd_MoFEM
PetscErrorCode DMLocalToGlobalEnd_MoFEM(DM, Vec, InsertMode, Vec)
Definition: DMMoFEM.cpp:1463
MoFEM::DMoFEMLoopFiniteElements
PetscErrorCode DMoFEMLoopFiniteElements(DM dm, const char fe_name[], MoFEM::FEMethod *method, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
Executes FEMethod for finite elements in DM.
Definition: DMMoFEM.cpp:586
MF_EXIST
@ MF_EXIST
Definition: definitions.h:113
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
MOFEM_NOT_IMPLEMENTED
@ MOFEM_NOT_IMPLEMENTED
Definition: definitions.h:32
GET_DM_FIELD
#define GET_DM_FIELD(DM)
Definition: PCMGSetUpViaApproxOrders.cpp:169
MoFEM::DMCreateLocalVector_MoFEM
PetscErrorCode DMCreateLocalVector_MoFEM(DM dm, Vec *l)
DMShellSetCreateLocalVector.
Definition: DMMoFEM.cpp:1187
CHKERRG
#define CHKERRG(n)
Check error code of MoFEM/MOAB/PETSc function.
Definition: definitions.h:496
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
MOFEM_INVALID_DATA
@ MOFEM_INVALID_DATA
Definition: definitions.h:36
sub_mat_sor
MoFEMErrorCode sub_mat_sor(Mat mat, Vec b, PetscReal omega, MatSORType flag, PetscReal shift, PetscInt its, PetscInt lits, Vec x)
Definition: PCMGSetUpViaApproxOrders.cpp:114
MoFEM::DMMoFEMSNESSetFunction
static PetscErrorCode DMMoFEMSNESSetFunction(DM dm, S fe_name, T0 method, T1 pre_only, T2 post_only)
Definition: DMMoFEM.cpp:699
Generated by Doxygen 1.8.17 and hosted at University of Glasgow