SetUpSchurImpl Struct Reference

Inheritance diagram for SetUpSchurImpl:

Collaboration diagram for SetUpSchurImpl:

Public Member Functions
	SetUpSchurImpl (MoFEM::Interface &m_field, SmartPetscObj< DM > sub_dm, SmartPetscObj< IS > field_split_is, SmartPetscObj< AO > ao_up)
virtual	~SetUpSchurImpl ()
MoFEMErrorCode	setUp (KSP solver)
MoFEMErrorCode	preProc ()
MoFEMErrorCode	postProc ()
	SetUpSchurImpl (MoFEM::Interface &m_field)
virtual	~SetUpSchurImpl ()
MoFEMErrorCode	setUp (SmartPetscObj< KSP > solver)
MoFEMErrorCode	preProc ()
MoFEMErrorCode	postProc ()
virtual MoFEMErrorCode	setUp (KSP solver)=0
virtual MoFEMErrorCode	preProc ()=0
virtual MoFEMErrorCode	postProc ()=0
virtual MoFEMErrorCode	setUp (SmartPetscObj< KSP > solver)=0
virtual MoFEMErrorCode	preProc ()=0
virtual MoFEMErrorCode	postProc ()=0

Private Member Functions
MoFEMErrorCode	setEntities ()
MoFEMErrorCode	setUpSubDM ()
MoFEMErrorCode	setOperator ()
MoFEMErrorCode	setPC (PC pc)

Private Attributes
SmartPetscObj< Mat >	S
MoFEM::Interface &	mField
SmartPetscObj< DM >	subDM
	field split sub dm More...
SmartPetscObj< IS >	fieldSplitIS
	IS for split Schur block. More...
SmartPetscObj< AO >	aoUp
Range	volEnts
Range	subEnts

Additional Inherited Members
Static Public Member Functions inherited from SetUpSchur
static boost::shared_ptr< SetUpSchur >	createSetUpSchur (MoFEM::Interface &m_field, SmartPetscObj< DM > sub_dm, SmartPetscObj< IS > field_split_it, SmartPetscObj< AO > ao_map)
	Create data structure for handling Schur complement. More...
static boost::shared_ptr< SetUpSchur >	createSetUpSchur (MoFEM::Interface &m_field)
Protected Member Functions inherited from SetUpSchur
	SetUpSchur ()=default
	SetUpSchur ()=default

Detailed Description

Examples

EshelbianPlasticity.cpp, and plastic.cpp.

Definition at line 1231 of file plastic.cpp.

Constructor & Destructor Documentation

SetUpSchurImpl() [1/2]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field, SmartPetscObj< DM >  sub_dm, SmartPetscObj< IS >  field_split_is, SmartPetscObj< AO >  ao_up  )

inline

Definition at line 1233 of file plastic.cpp.

      : SetUpSchur(), mField(m_field), subDM(sub_dm),
        fieldSplitIS(field_split_is), aoUp(ao_up) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
  }

~SetUpSchurImpl() [1/2]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

inlinevirtual

Examples

plastic.cpp.

Definition at line 1244 of file plastic.cpp.

{ S.reset(); }

SetUpSchurImpl() [2/2]

SetUpSchurImpl::SetUpSchurImpl ( MoFEM::Interface &  m_field )

inline

Definition at line 733 of file elastic.cpp.

                                          : SetUpSchur(), mField(m_field) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
  }

~SetUpSchurImpl() [2/2]

virtual SetUpSchurImpl::~SetUpSchurImpl ( )

inlinevirtual

Definition at line 741 of file elastic.cpp.

{ S.reset(); }

Member Function Documentation

postProc() [1/2]

MoFEMErrorCode SetUpSchurImpl::postProc ( )

virtual

Implements SetUpSchur.

Examples

plastic.cpp.

Definition at line 1335 of file plastic.cpp.

                                        {
  MoFEMFunctionBegin;
  MOFEM_LOG("TIMER", Sev::verbose) << "Lhs Assemble End";
  if (S) {
    CHKERR MatAssemblyBegin(S, MAT_FINAL_ASSEMBLY);
    CHKERR MatAssemblyEnd(S, MAT_FINAL_ASSEMBLY);
  }
  MoFEMFunctionReturn(0);
}

postProc() [2/2]

MoFEMErrorCode SetUpSchurImpl::postProc ( )

virtual

Implements SetUpSchur.

preProc() [1/2]

MoFEMErrorCode SetUpSchurImpl::preProc ( )

virtual

Implements SetUpSchur.

Examples

plastic.cpp.

Definition at line 1326 of file plastic.cpp.

                                       {
  MoFEMFunctionBegin;
  if (SetUpSchurImpl::S) {
    CHKERR MatZeroEntries(S);
  }
  MOFEM_LOG("TIMER", Sev::verbose) << "Lhs Assemble Begin";
  MoFEMFunctionReturn(0);
}

preProc() [2/2]

MoFEMErrorCode SetUpSchurImpl::preProc ( )

virtual

Implements SetUpSchur.

setEntities()

MoFEMErrorCode SetUpSchurImpl::setEntities ( )

private

Definition at line 795 of file elastic.cpp.

                                           {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  CHKERR mField.get_moab().get_entities_by_dimension(simple->getMeshset(),
                                                     SPACE_DIM, volEnts);
  CHKERR mField.get_moab().get_entities_by_handle(simple->getMeshset(),
                                                  subEnts);
  subEnts = subtract(subEnts, volEnts);
  MoFEMFunctionReturn(0);
};

setOperator()

MoFEMErrorCode SetUpSchurImpl::setOperator ( )

private

Definition at line 820 of file elastic.cpp.

                                           {
  MoFEMFunctionBegin;
  auto pip = mField.getInterface<PipelineManager>();
  // Boundary
  auto dm_is = getDMSubData(subDM)->getSmartRowIs();
  auto ao_up = createAOMappingIS(dm_is, PETSC_NULL);
  pip->getOpBoundaryLhsPipeline().push_front(new OpSchurAssembleBegin());
  pip->getOpBoundaryLhsPipeline().push_back(new OpSchurAssembleEnd<SCHUR_DSYSV>(
      {"U"}, {boost::make_shared<Range>(volEnts)}, {ao_up}, {S}, {true}));
  // Domain
  pip->getOpDomainLhsPipeline().push_front(new OpSchurAssembleBegin());
  pip->getOpDomainLhsPipeline().push_back(new OpSchurAssembleEnd<SCHUR_DSYSV>(
      {"U"}, {boost::make_shared<Range>(volEnts)}, {ao_up}, {S}, {true}));
  MoFEMFunctionReturn(0);
}

setPC()

MoFEMErrorCode SetUpSchurImpl::setPC ( PC  pc )

private

Definition at line 836 of file elastic.cpp.

                                          {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  SmartPetscObj<IS> vol_is;
  mField.getInterface<ISManager>()->isCreateProblemFieldAndRank(
      simple->getProblemName(), ROW, "U", 0, SPACE_DIM, vol_is, &volEnts);
  CHKERR PCFieldSplitSetIS(pc, NULL, vol_is);
  CHKERR PCFieldSplitSetSchurPre(pc, PC_FIELDSPLIT_SCHUR_PRE_USER, S);
  MoFEMFunctionReturn(0);
}

setUp() [1/2]

MoFEMErrorCode SetUpSchurImpl::setUp ( KSP  solver )

virtual

Implements SetUpSchur.

Examples

plastic.cpp.

Definition at line 1259 of file plastic.cpp.

                                               {
  MoFEMFunctionBegin;
  auto pip = mField.getInterface<PipelineManager>();
  PC pc;
  CHKERR KSPSetFromOptions(solver);
  CHKERR KSPGetPC(solver, &pc);
  PetscBool is_pcfs = PETSC_FALSE;
  PetscObjectTypeCompare((PetscObject)pc, PCFIELDSPLIT, &is_pcfs);
  if (is_pcfs) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
    S = createDMMatrix(subDM);
 
    auto set_ops = [&]() {
      MoFEMFunctionBegin;
      auto pip = mField.getInterface<PipelineManager>();
      // Boundary
      pip->getOpBoundaryLhsPipeline().push_front(new OpSchurAssembleBegin());
      pip->getOpBoundaryLhsPipeline().push_back(
          new OpSchurAssembleEnd<SCHUR_DGESV>(
 
              {"EP", "TAU"}, {nullptr, nullptr}, {SmartPetscObj<AO>(), aoUp},
              {SmartPetscObj<Mat>(), S}, {false, false}
 
              ));
      // Domain
      pip->getOpDomainLhsPipeline().push_front(new OpSchurAssembleBegin());
      pip->getOpDomainLhsPipeline().push_back(
          new OpSchurAssembleEnd<SCHUR_DGESV>(
 
              {"EP", "TAU"}, {nullptr, nullptr}, {SmartPetscObj<AO>(), aoUp},
              {SmartPetscObj<Mat>(), S}, {false, false}
 
              ));
      MoFEMFunctionReturn(0);
    };
 
    auto set_pc = [&]() {
      MoFEMFunctionBegin;
      CHKERR PCFieldSplitSetIS(pc, NULL, fieldSplitIS);
      CHKERR PCFieldSplitSetSchurPre(pc, PC_FIELDSPLIT_SCHUR_PRE_USER, S);
      MoFEMFunctionReturn(0);
    };
 
    CHKERR set_ops();
    CHKERR set_pc();
 
  } else {
    pip->getOpBoundaryLhsPipeline().push_front(new OpSchurAssembleBegin());
    pip->getOpBoundaryLhsPipeline().push_back(
        new OpSchurAssembleEnd<SCHUR_DGESV>({}, {}, {}, {}, {}));
    pip->getOpDomainLhsPipeline().push_front(new OpSchurAssembleBegin());
    pip->getOpDomainLhsPipeline().push_back(
        new OpSchurAssembleEnd<SCHUR_DGESV>({}, {}, {}, {}, {}));
  }
 
  // we do not those anymore
  subDM.reset();
  fieldSplitIS.reset();
  aoUp.reset();
  MoFEMFunctionReturn(0);
}

setUp() [2/2]

MoFEMErrorCode SetUpSchurImpl::setUp ( SmartPetscObj< KSP >  solver )

virtual

Implements SetUpSchur.

Definition at line 763 of file elastic.cpp.

                                                              {
  MoFEMFunctionBegin;
  auto pip = mField.getInterface<PipelineManager>();
  PC pc;
  CHKERR KSPGetPC(solver, &pc);
  PetscBool is_pcfs = PETSC_FALSE;
  PetscObjectTypeCompare((PetscObject)pc, PCFIELDSPLIT, &is_pcfs);
  if (is_pcfs) {
    if (S) {
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "Is expected that schur matrix is not allocated. This is "
          "possible only is PC is set up twice");
    }
    CHKERR setEntities();
    CHKERR setUpSubDM();
    S = createDMMatrix(subDM);
    // CHKERR MatSetBlockSize(S, SPACE_DIM);
    // CHKERR MatSetOption(S, MAT_SYMMETRIC, PETSC_TRUE);
    CHKERR setOperator();
    CHKERR setPC(pc);
  } else {
    pip->getOpBoundaryLhsPipeline().push_front(new OpSchurAssembleBegin());
    pip->getOpBoundaryLhsPipeline().push_back(
        new OpSchurAssembleEnd<SCHUR_DSYSV>({}, {}, {}, {}, {}));
    pip->getOpDomainLhsPipeline().push_front(new OpSchurAssembleBegin());
    pip->getOpDomainLhsPipeline().push_back(
        new OpSchurAssembleEnd<SCHUR_DSYSV>({}, {}, {}, {}, {}));
  }
  MoFEMFunctionReturn(0);
}

setUpSubDM()

MoFEMErrorCode SetUpSchurImpl::setUpSubDM ( )

private

Definition at line 806 of file elastic.cpp.

                                          {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  subDM = createDM(mField.get_comm(), "DMMOFEM");
  CHKERR DMMoFEMCreateSubDM(subDM, simple->getDM(), "SUB");
  CHKERR DMMoFEMSetSquareProblem(subDM, PETSC_TRUE);
  CHKERR DMMoFEMAddElement(subDM, simple->getDomainFEName());
  auto sub_ents_ptr = boost::make_shared<Range>(subEnts);
  CHKERR DMMoFEMAddSubFieldRow(subDM, "U", sub_ents_ptr);
  CHKERR DMMoFEMAddSubFieldCol(subDM, "U", sub_ents_ptr);
  CHKERR DMSetUp(subDM);
  MoFEMFunctionReturn(0);
}

Member Data Documentation

aoUp

SmartPetscObj<AO> SetUpSchurImpl::aoUp

private

Examples

plastic.cpp.

Definition at line 1256 of file plastic.cpp.

fieldSplitIS

SmartPetscObj<IS> SetUpSchurImpl::fieldSplitIS

private

IS for split Schur block.

Examples

plastic.cpp.

Definition at line 1255 of file plastic.cpp.

mField

MoFEM::Interface & SetUpSchurImpl::mField

private

Examples

plastic.cpp.

Definition at line 1253 of file plastic.cpp.

S

SmartPetscObj< Mat > SetUpSchurImpl::S

private

Examples

plastic.cpp.

Definition at line 1251 of file plastic.cpp.

subDM

SmartPetscObj< DM > SetUpSchurImpl::subDM

private

field split sub dm

Examples

plastic.cpp.

Definition at line 1254 of file plastic.cpp.

subEnts

Range SetUpSchurImpl::subEnts

private

Definition at line 759 of file elastic.cpp.

volEnts

Range SetUpSchurImpl::volEnts

private

Definition at line 758 of file elastic.cpp.

---

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

• users_modules/tutorials/adv-0/plastic.cpp
• users_modules/tutorials/vec-0/elastic.cpp