Poisson Struct Reference

Collaboration diagram for Poisson:

Public Member Functions
	Poisson (moab::Core &mb_instance, MoFEM::Core &core, int order)
MoFEMErrorCode	runAnalysis ()

Private Member Functions
MoFEMErrorCode	readMesh ()
MoFEMErrorCode	setupSystem ()
MoFEMErrorCode	setIntegrationRule ()
MoFEMErrorCode	assembleSystem ()
MoFEMErrorCode	applyBc ()
MoFEMErrorCode	solveSystem ()
MoFEMErrorCode	setOutput ()

Private Attributes
boost::shared_ptr< VolEle >	volPipelineLhs
boost::shared_ptr< VolEle >	volPipelineRhs
boost::shared_ptr< VolEle >	postFaceVolume
MoFEM::Interface &	mField
moab::Interface &	moab
Simple *	simpleInterface
DM	dm
KSP	ksp
MatrixDouble	invJac
MPI_Comm	mpiComm
const int	mpiRank
int	order
std::string	domainField

Detailed Description

Definition at line 10 of file poisson_prob.cpp.

Constructor & Destructor Documentation

Poisson()

Poisson::Poisson	(	moab::Core &	mb_instance,
		MoFEM::Core &	core,
		int	order
	)

Definition at line 46 of file poisson_prob.cpp.

    : domainField("U"), moab(mb_instance), mField(core),
      mpiComm(mField.get_comm()), mpiRank(mField.get_comm_rank()),
      order(order) {
  volPipelineLhs = boost::shared_ptr<VolEle>(new VolEle(mField));
  volPipelineRhs = boost::shared_ptr<VolEle>(new VolEle(mField));
}

Member Function Documentation

applyBc()

MoFEMErrorCode Poisson::applyBc ( )

private

Definition at line 171 of file poisson_prob.cpp.

                                {
  MoFEMFunctionBegin;
  Range bdry_entities;
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(mField, BLOCKSET, it)) {
    string name = it->getName();
    if (name.compare(0, 14, "ESSENTIAL") == 0) {
      CHKERR it->getMeshsetIdEntitiesByDimension(mField.get_moab(), 1,
                                                 bdry_entities, true);
    }
  }
 
  cerr << bdry_entities;
  Range bdr_verts;
  CHKERR mField.get_moab().get_connectivity(bdry_entities, bdr_verts, true);
  bdry_entities.merge(bdr_verts);
 
  CHKERR mField.getInterface<ProblemsManager>()->removeDofsOnEntities(
      simpleInterface->getProblemName(), domainField, bdry_entities);
 
  
  MoFEMFunctionReturn(0);
}

assembleSystem()

MoFEMErrorCode Poisson::assembleSystem ( )

private

Definition at line 84 of file poisson_prob.cpp.

                                       {
  MoFEMFunctionBegin;
  auto det_ptr = boost::make_shared<VectorDouble>();
  auto jac_ptr = boost::make_shared<MatrixDouble>();
  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
  volPipelineLhs->getOpPtrVector().push_back(
      new OpCalculateHOJacForFace(jac_ptr));
  volPipelineLhs->getOpPtrVector().push_back(
      new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
  volPipelineLhs->getOpPtrVector().push_back(
      new OpLhsUU(domainField, domainField));
  volPipelineRhs->getOpPtrVector().push_back(new OpRhsU(domainField));
 
  // dm = simpleInterface->getDM();
  CHKERR simpleInterface->getDM(&dm);
  boost::shared_ptr<VolEle> null;
 
  CHKERR setIntegrationRule();
 
  CHKERR DMMoFEMKSPSetComputeOperators(dm, simpleInterface->getDomainFEName(), volPipelineLhs, null, null);
  CHKERR DMMoFEMKSPSetComputeRHS(dm, simpleInterface->getDomainFEName(), volPipelineRhs, null, null);
  MoFEMFunctionReturn(0);
}

readMesh()

MoFEMErrorCode Poisson::readMesh ( )

private

Definition at line 54 of file poisson_prob.cpp.

                                 {
  MoFEMFunctionBegin;
 
  CHKERR mField.getInterface(simpleInterface);
  CHKERR simpleInterface->getOptions();
  CHKERR simpleInterface->loadFile();
 
  MoFEMFunctionReturn(0);
}

runAnalysis()

MoFEMErrorCode Poisson::runAnalysis

(

)

Definition at line 194 of file poisson_prob.cpp.

                                    {
  MoFEMFunctionBegin;
 
  CHKERR readMesh();
  CHKERR setupSystem();
  CHKERR applyBc();
  
  CHKERR assembleSystem();
  
  CHKERR setOutput();
  
  CHKERR solveSystem();
 
  MoFEMFunctionReturn(0);
}

setIntegrationRule()

MoFEMErrorCode Poisson::setIntegrationRule ( )

private

Definition at line 75 of file poisson_prob.cpp.

                                           {
  MoFEMFunctionBegin;
  auto vol_ruleLhs = [](int, int, int p) -> int { return 2 * (p - 1); };
  auto vol_ruleRhs = [](int, int, int p) -> int { return 2 * (p - 1); };
  volPipelineLhs->getRuleHook = vol_ruleLhs;
  volPipelineRhs->getRuleHook = vol_ruleRhs;
  MoFEMFunctionReturn(0);
}

setOutput()

MoFEMErrorCode Poisson::setOutput ( )

private

Definition at line 152 of file poisson_prob.cpp.

                                  {
  MoFEMFunctionBegin;
 
  postFaceVolume =
      boost::shared_ptr<VolEle>(new PostProcFaceOnRefinedMesh(mField));
 
  CHKERR
  boost::static_pointer_cast<PostProcFaceOnRefinedMesh>(postFaceVolume)
      ->generateReferenceElementMesh();
  CHKERR
  boost::static_pointer_cast<PostProcFaceOnRefinedMesh>(postFaceVolume)
      ->addFieldValuesPostProc(domainField);
  // CHKERR
  // boost::static_pointer_cast<PostProcFaceOnRefinedMesh>(postFaceVolume)
  //     ->addFieldValuesGradientPostProc(domainField);
 
  MoFEMFunctionReturn(0);
}

setupSystem()

MoFEMErrorCode Poisson::setupSystem ( )

private

Definition at line 64 of file poisson_prob.cpp.

                                    {
  MoFEMFunctionBegin;
 
  CHKERR simpleInterface->addDomainField(domainField, H1,
                                         AINSWORTH_LEGENDRE_BASE, 1);
  CHKERR simpleInterface->setFieldOrder(domainField, order);
  CHKERR simpleInterface->setUp();
 
  MoFEMFunctionReturn(0);
}

solveSystem()

MoFEMErrorCode Poisson::solveSystem ( )

private

Definition at line 108 of file poisson_prob.cpp.

                                    {
  MoFEMFunctionBegin;
 
  Vec global_rhs, global_solution;
  CHKERR DMCreateGlobalVector(dm, &global_rhs);
  CHKERR VecDuplicate(global_rhs, &global_solution);
 
 
  // ksp = createKSP(mField.get_comm());
 
  CHKERR KSPCreate(PETSC_COMM_WORLD, &ksp);
 
  // CHKERR KSPSetOperators(ksp, volPipelineLhs->getFEMethod()->ksp_B,
  //                        volPipelineLhs->getFEMethod()->ksp_B);
 
  CHKERR KSPSetFromOptions(ksp);
 
  CHKERR KSPSetDM(ksp, dm);
 
  CHKERR KSPSetUp(ksp);
 
  CHKERR KSPSolve(ksp, global_rhs, global_solution);
 
  CHKERR DMoFEMMeshToGlobalVector(dm, global_solution, INSERT_VALUES,
                                  SCATTER_REVERSE);
 
  
  
  CHKERR KSPDestroy(&ksp);
  CHKERR VecDestroy(&global_solution);
  CHKERR VecDestroy(&global_rhs);
 
  CHKERR DMoFEMLoopFiniteElements(dm, simpleInterface->getDomainFEName(),
                                  postFaceVolume);
 
  CHKERR
  boost::static_pointer_cast<PostProcFaceOnRefinedMesh>(postFaceVolume)
      ->writeFile("out_vol.h5m");
 
  CHKERR DMDestroy(&dm);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dm

DM Poisson::dm

private

Definition at line 33 of file poisson_prob.cpp.

domainField

std::string Poisson::domainField

private

Definition at line 43 of file poisson_prob.cpp.

invJac

MatrixDouble Poisson::invJac

private

Definition at line 36 of file poisson_prob.cpp.

ksp

KSP Poisson::ksp

private

Definition at line 34 of file poisson_prob.cpp.

mField

MoFEM::Interface& Poisson::mField

private

Definition at line 29 of file poisson_prob.cpp.

moab

moab::Interface& Poisson::moab

private

Definition at line 30 of file poisson_prob.cpp.

mpiComm

MPI_Comm Poisson::mpiComm

private

Definition at line 38 of file poisson_prob.cpp.

mpiRank

const int Poisson::mpiRank

private

Definition at line 39 of file poisson_prob.cpp.

order

int Poisson::order

private

Definition at line 41 of file poisson_prob.cpp.

postFaceVolume

boost::shared_ptr<VolEle> Poisson::postFaceVolume

private

Definition at line 27 of file poisson_prob.cpp.

simpleInterface

Simple* Poisson::simpleInterface

private

Definition at line 31 of file poisson_prob.cpp.

volPipelineLhs

boost::shared_ptr<VolEle> Poisson::volPipelineLhs

private

Definition at line 25 of file poisson_prob.cpp.

volPipelineRhs

boost::shared_ptr<VolEle> Poisson::volPipelineRhs

private

Definition at line 26 of file poisson_prob.cpp.

---

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

• users_modules/softmech/poisson_problem/poisson_prob.cpp