ApproxSphere Struct Reference

Collaboration diagram for ApproxSphere:

Public Member Functions
	ApproxSphere (MoFEM::Interface &m_field)
MoFEMErrorCode	runProblem ()
	[Run programme] More...

Private Member Functions
MoFEMErrorCode	getOptions ()
	[Run programme] More...
MoFEMErrorCode	readMesh ()
	[Read mesh] More...
MoFEMErrorCode	setupProblem ()
	[Read mesh] More...
MoFEMErrorCode	setOPs ()
	[Set up problem] More...
MoFEMErrorCode	solveSystem ()
	[Push operators to pipeline] More...
MoFEMErrorCode	outputResults ()
	[Solve] More...

Private Attributes
MoFEM::Interface &	mField

Detailed Description

Examples

approx_sphere.cpp.

Definition at line 262 of file approx_sphere.cpp.

Constructor & Destructor Documentation

ApproxSphere()

ApproxSphere::ApproxSphere ( MoFEM::Interface &  m_field )

inline

Definition at line 264 of file approx_sphere.cpp.

: mField(m_field) {}

Member Function Documentation

getOptions()

MoFEMErrorCode ApproxSphere::getOptions ( )

private

[Run programme]

Examples

approx_sphere.cpp.

Definition at line 292 of file approx_sphere.cpp.

                                        {
  MoFEMFunctionBeginHot;
  MoFEMFunctionReturnHot(0);
}

outputResults()

MoFEMErrorCode ApproxSphere::outputResults ( )

private

[Solve]

Examples

approx_sphere.cpp.

Definition at line 394 of file approx_sphere.cpp.

                                           {
  MoFEMFunctionBegin;
 
  auto x_ptr = boost::make_shared<MatrixDouble>();
  auto det_ptr = boost::make_shared<VectorDouble>();
  auto jac_ptr = boost::make_shared<MatrixDouble>();
  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
 
  auto simple = mField.getInterface<Simple>();
  auto dm = simple->getDM();
 
  auto post_proc_fe =
      boost::make_shared<PostProcBrokenMeshInMoab<DomainEle>>(mField);
 
  post_proc_fe->getOpPtrVector().push_back(
      new OpCalculateVectorFieldValues<3>("HO_POSITIONS", x_ptr));
 
  using OpPPMap = OpPostProcMapInMoab<3, 3>;
 
  post_proc_fe->getOpPtrVector().push_back(
 
      new OpPPMap(post_proc_fe->getPostProcMesh(),
                  post_proc_fe->getMapGaussPts(),
 
                  {},
 
                  {{"HO_POSITIONS", x_ptr}},
 
                  {}, {}
 
                  )
 
  );
 
  CHKERR DMoFEMLoopFiniteElements(dm, "dFE", post_proc_fe);
  CHKERR post_proc_fe->writeFile("out_approx.h5m");
 
  auto error_fe = boost::make_shared<DomainEle>(mField);
 
  error_fe->getOpPtrVector().push_back(
      new OpGetHONormalsOnFace("HO_POSITIONS"));
  error_fe->getOpPtrVector().push_back(
      new OpCalculateVectorFieldValues<3>("HO_POSITIONS", x_ptr));
  error_fe->getOpPtrVector().push_back(new OpError("HO_POSITIONS", x_ptr));
 
  error_fe->preProcessHook = [&]() {
    MoFEMFunctionBegin;
    MOFEM_LOG("EXAMPLE", Sev::inform) << "Create vec ";
    OpError::errorVec = createVectorMPI(
        mField.get_comm(), (!mField.get_comm_rank()) ? 1 : 0, 1);
    VecZeroEntries(OpError::errorVec);
    MoFEMFunctionReturn(0);
  };
 
  error_fe->postProcessHook = [&]() {
    MoFEMFunctionBegin;
    CHKERR VecAssemblyBegin(OpError::errorVec);
    CHKERR VecAssemblyEnd(OpError::errorVec);
    double error2;
    CHKERR VecSum(OpError::errorVec, &error2);
    MOFEM_LOG("EXAMPLE", Sev::inform)
        << "Error " << std::sqrt(error2 / (4 * M_PI * A * A));
    OpError::errorVec.reset();
    MoFEMFunctionReturn(0);
  };
 
  CHKERR DMoFEMLoopFiniteElements(dm, "dFE", error_fe);
 
  CHKERR simple->deleteDM();
  CHKERR simple->deleteFiniteElements();
  if (mField.get_comm_size() > 1)
    CHKERR mField.get_moab().write_file("out_ho_mesh.h5m", "MOAB",
                                        "PARALLEL=WRITE_PART");
  else
    CHKERR mField.get_moab().write_file("out_ho_mesh.h5m");
  MoFEMFunctionReturn(0);
}

readMesh()

MoFEMErrorCode ApproxSphere::readMesh ( )

private

[Read mesh]

Examples

approx_sphere.cpp.

Definition at line 298 of file approx_sphere.cpp.

                                      {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
  CHKERR simple->getOptions();
  CHKERR simple->loadFile();
  MoFEMFunctionReturn(0);
}

runProblem()

MoFEMErrorCode ApproxSphere::runProblem

(

)

[Run programme]

Examples

approx_sphere.cpp.

Definition at line 280 of file approx_sphere.cpp.

                                        {
  MoFEMFunctionBegin;
  CHKERR getOptions();
  CHKERR readMesh();
  CHKERR setupProblem();
  CHKERR setOPs();
  CHKERR solveSystem();
  CHKERR outputResults();
  MoFEMFunctionReturn(0);
}

setOPs()

MoFEMErrorCode ApproxSphere::setOPs ( )

private

[Set up problem]

[Push operators to pipeline]

Examples

approx_sphere.cpp.

Definition at line 325 of file approx_sphere.cpp.

                                    {
  MoFEMFunctionBegin;
  PipelineManager *pipeline_mng = mField.getInterface<PipelineManager>();
 
  auto integration_rule = [](int, int, int approx_order) {
    return 3 * approx_order;
  };
  CHKERR pipeline_mng->setDomainRhsIntegrationRule(integration_rule);
  CHKERR pipeline_mng->setDomainLhsIntegrationRule(integration_rule);
 
  auto x_ptr = boost::make_shared<MatrixDouble>();
  auto dot_x_ptr = boost::make_shared<MatrixDouble>();
  auto det_ptr = boost::make_shared<VectorDouble>();
  auto jac_ptr = boost::make_shared<MatrixDouble>();
  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
 
  auto def_ops = [&](auto &pipeline) {
    pipeline.push_back(
        new OpCalculateVectorFieldValues<3>("HO_POSITIONS", x_ptr));
    pipeline.push_back(
        new OpCalculateVectorFieldValuesDot<3>("HO_POSITIONS", dot_x_ptr));
  };
 
  def_ops(pipeline_mng->getOpDomainRhsPipeline());
  def_ops(pipeline_mng->getOpDomainLhsPipeline());
 
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpRhs("HO_POSITIONS", x_ptr, dot_x_ptr));
  pipeline_mng->getOpDomainLhsPipeline().push_back(
      new OpLhs("HO_POSITIONS", x_ptr, dot_x_ptr));
 
  MoFEMFunctionReturn(0);
}

setupProblem()

MoFEMErrorCode ApproxSphere::setupProblem ( )

private

[Read mesh]

[Set up problem]

Examples

approx_sphere.cpp.

Definition at line 308 of file approx_sphere.cpp.

                                          {
  MoFEMFunctionBegin;
 
  auto simple = mField.getInterface<Simple>();
  CHKERR simple->addDomainField("HO_POSITIONS", H1, AINSWORTH_LEGENDRE_BASE, 3);
  CHKERR simple->addDataField("HO_POSITIONS", H1, AINSWORTH_LEGENDRE_BASE, 3);
 
  int order = 3;
  CHKERR PetscOptionsGetInt(PETSC_NULL, "", "-order", &order, PETSC_NULL);
  CHKERR simple->setFieldOrder("HO_POSITIONS", order);
  CHKERR simple->setUp();
 
  MoFEMFunctionReturn(0);
}

solveSystem()

MoFEMErrorCode ApproxSphere::solveSystem ( )

private

[Push operators to pipeline]

[Solve]

Examples

approx_sphere.cpp.

Definition at line 361 of file approx_sphere.cpp.

                                         {
  MoFEMFunctionBegin;
 
  // Project HO geometry from mesh
  Projection10NodeCoordsOnField ent_method_material(mField, "HO_POSITIONS");
  CHKERR mField.loop_dofs("HO_POSITIONS", ent_method_material);
 
  auto *simple = mField.getInterface<Simple>();
  auto *pipeline_mng = mField.getInterface<PipelineManager>();
 
  auto dm = simple->getDM();
  MoFEM::SmartPetscObj<TS> ts;
  ts = pipeline_mng->createTSIM();
 
  double ftime = 1;
  CHKERR TSSetMaxSteps(ts, 1);
  CHKERR TSSetExactFinalTime(ts, TS_EXACTFINALTIME_MATCHSTEP);
 
  auto T = createDMVector(simple->getDM());
  CHKERR DMoFEMMeshToLocalVector(simple->getDM(), T, INSERT_VALUES,
                                 SCATTER_FORWARD);
  CHKERR TSSetSolution(ts, T);
  CHKERR TSSetFromOptions(ts);
 
  CHKERR TSSolve(ts, NULL);
  CHKERR TSGetTime(ts, &ftime);
 
  CHKERR mField.getInterface<FieldBlas>()->fieldScale(A, "HO_POSITIONS");
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

mField

MoFEM::Interface& ApproxSphere::mField

private

Definition at line 269 of file approx_sphere.cpp.

---

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

• tutorials/vec-4/approx_sphere.cpp