v0.13.1
Loading...
Searching...
No Matches
Public Member Functions | Private Member Functions | Private Attributes | List of all members
ApproxSphere Struct Reference
Collaboration diagram for ApproxSphere:
[legend]

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::InterfacemField
 

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.

264: mField(m_field) {}
MoFEM::Interface & mField

Member Function Documentation

◆ getOptions()

MoFEMErrorCode ApproxSphere::getOptions ( )
private

[Run programme]

Examples
approx_sphere.cpp.

Definition at line 292 of file approx_sphere.cpp.

292 {
295}
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:447
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:440

◆ outputResults()

MoFEMErrorCode ApproxSphere::outputResults ( )
private

[Solve]

Examples
approx_sphere.cpp.

Definition at line 394 of file approx_sphere.cpp.

394 {
396
397 auto x_ptr = boost::make_shared<MatrixDouble>();
398 auto det_ptr = boost::make_shared<VectorDouble>();
399 auto jac_ptr = boost::make_shared<MatrixDouble>();
400 auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
401
403 auto dm = simple->getDM();
404
405 auto post_proc_fe =
406 boost::make_shared<PostProcBrokenMeshInMoab<DomainEle>>(mField);
407
408 post_proc_fe->getOpPtrVector().push_back(
409 new OpCalculateVectorFieldValues<3>("HO_POSITIONS", x_ptr));
410
412
413 post_proc_fe->getOpPtrVector().push_back(
414
415 new OpPPMap(post_proc_fe->getPostProcMesh(),
416 post_proc_fe->getMapGaussPts(),
417
418 {},
419
420 {{"HO_POSITIONS", x_ptr}},
421
422 {}, {}
423
424 )
425
426 );
427
428 CHKERR DMoFEMLoopFiniteElements(dm, "dFE", post_proc_fe);
429 CHKERR post_proc_fe->writeFile("out_approx.h5m");
430
431 auto error_fe = boost::make_shared<DomainEle>(mField);
432
433 error_fe->getOpPtrVector().push_back(
434 new OpGetHONormalsOnFace("HO_POSITIONS"));
435 error_fe->getOpPtrVector().push_back(
436 new OpCalculateVectorFieldValues<3>("HO_POSITIONS", x_ptr));
437 error_fe->getOpPtrVector().push_back(new OpError("HO_POSITIONS", x_ptr));
438
439 error_fe->preProcessHook = [&]() {
441 MOFEM_LOG("EXAMPLE", Sev::inform) << "Create vec ";
443 mField.get_comm(), (!mField.get_comm_rank()) ? 1 : 0, 1);
444 VecZeroEntries(OpError::errorVec);
446 };
447
448 error_fe->postProcessHook = [&]() {
450 CHKERR VecAssemblyBegin(OpError::errorVec);
451 CHKERR VecAssemblyEnd(OpError::errorVec);
452 double error2;
453 CHKERR VecSum(OpError::errorVec, &error2);
454 MOFEM_LOG("EXAMPLE", Sev::inform)
455 << "Error " << std::sqrt(error2 / (4 * M_PI * A * A));
456 OpError::errorVec.reset();
458 };
459
460 CHKERR DMoFEMLoopFiniteElements(dm, "dFE", error_fe);
461
462 CHKERR simple->deleteDM();
463 CHKERR simple->deleteFiniteElements();
464 if (mField.get_comm_size() > 1)
465 CHKERR mField.get_moab().write_file("out_ho_mesh.h5m", "MOAB",
466 "PARALLEL=WRITE_PART");
467 else
468 CHKERR mField.get_moab().write_file("out_ho_mesh.h5m");
470}
void simple(double P1[], double P2[], double P3[], double c[], const int N)
Definition: acoustic.cpp:69
constexpr double A
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:346
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:416
#define CHKERR
Inline error check.
Definition: definitions.h:535
PetscErrorCode DMoFEMLoopFiniteElements(DM dm, const char fe_name[], MoFEM::FEMethod *method, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
Executes FEMethod for finite elements in DM.
Definition: DMMMoFEM.cpp:533
#define MOFEM_LOG(channel, severity)
Log.
Definition: LogManager.hpp:301
auto createSmartVectorMPI(MPI_Comm comm, PetscInt n, PetscInt N)
Create MPI Vector.
OpPostProcMapInMoab< SPACE_DIM, SPACE_DIM > OpPPMap
virtual int get_comm_size() const =0
virtual moab::Interface & get_moab()=0
virtual MPI_Comm & get_comm() const =0
virtual int get_comm_rank() const =0
Get values at integration pts for tensor filed rank 1, i.e. vector field.
Calculate normals at Gauss points of triangle element.
Post post-proc data at points from hash maps.
Simple interface for fast problem set-up.
Definition: Simple.hpp:26
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.
static SmartPetscObj< Vec > errorVec

◆ readMesh()

MoFEMErrorCode ApproxSphere::readMesh ( )
private

[Read mesh]

Examples
approx_sphere.cpp.

Definition at line 298 of file approx_sphere.cpp.

298 {
301 CHKERR simple->getOptions();
302 CHKERR simple->loadFile();
304}

◆ runProblem()

MoFEMErrorCode ApproxSphere::runProblem ( )

[Run programme]

Examples
approx_sphere.cpp.

Definition at line 280 of file approx_sphere.cpp.

280 {
285 CHKERR setOPs();
289}
MoFEMErrorCode readMesh()
[Read mesh]
MoFEMErrorCode solveSystem()
[Push operators to pipeline]
MoFEMErrorCode setupProblem()
[Read mesh]
MoFEMErrorCode setOPs()
[Set up problem]
MoFEMErrorCode getOptions()
[Run programme]
MoFEMErrorCode outputResults()
[Solve]

◆ 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.

325 {
328
329 auto integration_rule = [](int, int, int approx_order) {
330 return 3 * approx_order;
331 };
334
335 auto x_ptr = boost::make_shared<MatrixDouble>();
336 auto dot_x_ptr = boost::make_shared<MatrixDouble>();
337 auto det_ptr = boost::make_shared<VectorDouble>();
338 auto jac_ptr = boost::make_shared<MatrixDouble>();
339 auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
340
341 auto def_ops = [&](auto &pipeline) {
342 pipeline.push_back(
343 new OpCalculateVectorFieldValues<3>("HO_POSITIONS", x_ptr));
344 pipeline.push_back(
345 new OpCalculateVectorFieldValuesDot<3>("HO_POSITIONS", dot_x_ptr));
346 };
347
348 def_ops(pipeline_mng->getOpDomainRhsPipeline());
349 def_ops(pipeline_mng->getOpDomainLhsPipeline());
350
351 pipeline_mng->getOpDomainRhsPipeline().push_back(
352 new OpRhs("HO_POSITIONS", x_ptr, dot_x_ptr));
353 pipeline_mng->getOpDomainLhsPipeline().push_back(
354 new OpLhs("HO_POSITIONS", x_ptr, dot_x_ptr));
355
357}
auto integration_rule
boost::ptr_vector< UserDataOperator > & getOpDomainLhsPipeline()
Get the Op Domain Lhs Pipeline object.
boost::ptr_vector< UserDataOperator > & getOpDomainRhsPipeline()
Get the Op Domain Rhs Pipeline object.
static constexpr int approx_order
Approximate field valuse for given petsc vector.
PipelineManager interface.
MoFEMErrorCode setDomainRhsIntegrationRule(RuleHookFun rule)
MoFEMErrorCode setDomainLhsIntegrationRule(RuleHookFun rule)

◆ setupProblem()

MoFEMErrorCode ApproxSphere::setupProblem ( )
private

[Read mesh]

[Set up problem]

Examples
approx_sphere.cpp.

Definition at line 308 of file approx_sphere.cpp.

308 {
310
312 CHKERR simple->addDomainField("HO_POSITIONS", H1, AINSWORTH_LEGENDRE_BASE, 3);
313 CHKERR simple->addDataField("HO_POSITIONS", H1, AINSWORTH_LEGENDRE_BASE, 3);
314
315 int order = 3;
316 CHKERR PetscOptionsGetInt(PETSC_NULL, "", "-order", &order, PETSC_NULL);
317 CHKERR simple->setFieldOrder("HO_POSITIONS", order);
318 CHKERR simple->setUp();
319
321}
@ AINSWORTH_LEGENDRE_BASE
Ainsworth Cole (Legendre) approx. base .
Definition: definitions.h:60
@ H1
continuous field
Definition: definitions.h:85
PetscErrorCode PetscOptionsGetInt(PetscOptions *, const char pre[], const char name[], PetscInt *ivalue, PetscBool *set)

◆ solveSystem()

MoFEMErrorCode ApproxSphere::solveSystem ( )
private

[Push operators to pipeline]

[Solve]

Examples
approx_sphere.cpp.

Definition at line 361 of file approx_sphere.cpp.

361 {
363
364 // Project HO geometry from mesh
365 Projection10NodeCoordsOnField ent_method_material(mField, "HO_POSITIONS");
366 CHKERR mField.loop_dofs("HO_POSITIONS", ent_method_material);
367
368 auto *simple = mField.getInterface<Simple>();
369 auto *pipeline_mng = mField.getInterface<PipelineManager>();
370
371 auto dm = simple->getDM();
373 ts = pipeline_mng->createTSIM();
374
375 double ftime = 1;
376 CHKERR TSSetMaxSteps(ts, 1);
377 CHKERR TSSetExactFinalTime(ts, TS_EXACTFINALTIME_MATCHSTEP);
378
379 auto T = smartCreateDMVector(simple->getDM());
380 CHKERR DMoFEMMeshToLocalVector(simple->getDM(), T, INSERT_VALUES,
381 SCATTER_FORWARD);
382 CHKERR TSSetSolution(ts, T);
383 CHKERR TSSetFromOptions(ts);
384
385 CHKERR TSSolve(ts, NULL);
386 CHKERR TSGetTime(ts, &ftime);
387
388 CHKERR mField.getInterface<FieldBlas>()->fieldScale(A, "HO_POSITIONS");
389
391}
PetscErrorCode DMoFEMMeshToLocalVector(DM dm, Vec l, InsertMode mode, ScatterMode scatter_mode)
set local (or ghosted) vector values on mesh for partition only
Definition: DMMMoFEM.cpp:470
auto smartCreateDMVector(DM dm)
Get smart vector from DM.
Definition: DMMoFEM.hpp:965
virtual MoFEMErrorCode loop_dofs(const Problem *problem_ptr, const std::string &field_name, RowColData rc, DofMethod &method, int lower_rank, int upper_rank, int verb=DEFAULT_VERBOSITY)=0
Make a loop over dofs.
const double T
Basic algebra on fields.
Definition: FieldBlas.hpp:21
Projection of edge entities with one mid-node on hierarchical basis.
intrusive_ptr for managing petsc objects

Member Data Documentation

◆ mField

MoFEM::Interface& ApproxSphere::mField
private
Examples
approx_sphere.cpp.

Definition at line 269 of file approx_sphere.cpp.


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