v0.14.0
Public Member Functions | Private Member Functions | Private Attributes | List of all members
Poisson2DiscontGalerkin Struct Reference
Collaboration diagram for Poisson2DiscontGalerkin:
[legend]

Public Member Functions

 Poisson2DiscontGalerkin (MoFEM::Interface &m_field)
 
MoFEMErrorCode runProgram ()
 [Output results] More...
 

Private Member Functions

MoFEMErrorCode readMesh ()
 [Read mesh] More...
 
MoFEMErrorCode setupProblem ()
 [Read mesh] More...
 
MoFEMErrorCode boundaryCondition ()
 [Setup problem] More...
 
MoFEMErrorCode assembleSystem ()
 [Boundary condition] More...
 
MoFEMErrorCode setIntegrationRules ()
 [Assemble system] More...
 
MoFEMErrorCode solveSystem ()
 [Set integration rules] More...
 
MoFEMErrorCode checkResults ()
 [Solve system] More...
 
MoFEMErrorCode outputResults ()
 [Output results] More...
 

Private Attributes

MoFEM::Interface & mField
 
Simple * simpleInterface
 
std::string domainField
 
int oRder
 

Detailed Description

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 73 of file poisson_2d_dis_galerkin.cpp.

Constructor & Destructor Documentation

◆ Poisson2DiscontGalerkin()

Poisson2DiscontGalerkin::Poisson2DiscontGalerkin ( MoFEM::Interface &  m_field)
Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 100 of file poisson_2d_dis_galerkin.cpp.

101  : domainField("U"), mField(m_field), oRder(4) {}

Member Function Documentation

◆ assembleSystem()

MoFEMErrorCode Poisson2DiscontGalerkin::assembleSystem ( )
private

[Boundary condition]

[Assemble system]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 173 of file poisson_2d_dis_galerkin.cpp.

173  {
174  MoFEMFunctionBegin;
175 
176  auto pipeline_mng = mField.getInterface<PipelineManager>();
177 
178  auto add_base_ops = [&](auto &pipeline) {
179  auto det_ptr = boost::make_shared<VectorDouble>();
180  auto jac_ptr = boost::make_shared<MatrixDouble>();
181  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
182  pipeline.push_back(new OpCalculateHOJacForFace(jac_ptr));
183  pipeline.push_back(new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
184  pipeline.push_back(new OpSetInvJacL2ForFace(inv_jac_ptr));
185  };
186 
187  add_base_ops(pipeline_mng->getOpDomainLhsPipeline());
188  pipeline_mng->getOpDomainLhsPipeline().push_back(new OpDomainGradGrad(
189  domainField, domainField,
190  [](const double, const double, const double) { return 1; }));
191  pipeline_mng->getOpDomainRhsPipeline().push_back(
192  new OpDomainSource(domainField, source));
193 
194  // Push operators to the Pipeline for Skeleton
195  auto side_fe_ptr = boost::make_shared<FaceSideEle>(mField);
196  add_base_ops(side_fe_ptr->getOpPtrVector());
197  side_fe_ptr->getOpPtrVector().push_back(
198  new OpCalculateSideData(domainField, domainField));
199 
200  // Push operators to the Pipeline for Skeleton
201  pipeline_mng->getOpSkeletonLhsPipeline().push_back(
202  new OpL2LhsPenalty(side_fe_ptr));
203 
204  // Push operators to the Pipeline for Boundary
205  pipeline_mng->getOpBoundaryLhsPipeline().push_back(
206  new OpL2LhsPenalty(side_fe_ptr));
207  pipeline_mng->getOpBoundaryRhsPipeline().push_back(
208  new OpL2BoundaryRhs(side_fe_ptr, u_exact));
209 
210  MoFEMFunctionReturn(0);
211 }

◆ boundaryCondition()

MoFEMErrorCode Poisson2DiscontGalerkin::boundaryCondition ( )
private

[Setup problem]

[Boundary condition]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 165 of file poisson_2d_dis_galerkin.cpp.

165  {
166  MoFEMFunctionBegin;
167  MoFEMFunctionReturn(0);
168 }

◆ checkResults()

MoFEMErrorCode Poisson2DiscontGalerkin::checkResults ( )
private

[Solve system]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 261 of file poisson_2d_dis_galerkin.cpp.

261  {
262  MoFEMFunctionBegin;
263 
264  auto pipeline_mng = mField.getInterface<PipelineManager>();
265  pipeline_mng->getDomainRhsFE().reset();
266  pipeline_mng->getDomainLhsFE().reset();
267  pipeline_mng->getSkeletonRhsFE().reset();
268  pipeline_mng->getSkeletonLhsFE().reset();
269  pipeline_mng->getBoundaryRhsFE().reset();
270  pipeline_mng->getBoundaryLhsFE().reset();
271 
272  auto rule = [](int, int, int p) -> int { return 2 * p; };
273  CHKERR pipeline_mng->setDomainRhsIntegrationRule(rule);
274  auto rule_2 = [this](int, int, int) { return 2 * oRder; };
275  CHKERR pipeline_mng->setSkeletonRhsIntegrationRule(rule_2);
276 
277  auto add_base_ops = [&](auto &pipeline) {
278  auto det_ptr = boost::make_shared<VectorDouble>();
279  auto jac_ptr = boost::make_shared<MatrixDouble>();
280  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
281  pipeline.push_back(new OpCalculateHOJacForFace(jac_ptr));
282  pipeline.push_back(new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
283  pipeline.push_back(new OpSetInvJacL2ForFace(inv_jac_ptr));
284  };
285 
286  auto u_vals_ptr = boost::make_shared<VectorDouble>();
287  auto u_grad_ptr = boost::make_shared<MatrixDouble>();
288 
289  add_base_ops(pipeline_mng->getOpDomainRhsPipeline());
290  pipeline_mng->getOpDomainRhsPipeline().push_back(
291  new OpCalculateScalarFieldValues(domainField, u_vals_ptr));
292  pipeline_mng->getOpDomainRhsPipeline().push_back(
293  new OpCalculateScalarFieldGradient<2>(domainField, u_grad_ptr));
294 
295  enum NORMS { L2 = 0, SEMI_NORM, H1, LAST_NORM };
296  std::array<int, LAST_NORM> error_indices;
297  for (int i = 0; i != LAST_NORM; ++i)
298  error_indices[i] = i;
299  auto l2_vec = createVectorMPI(
300  mField.get_comm(), (!mField.get_comm_rank()) ? LAST_NORM : 0, LAST_NORM);
301 
302  auto error_op = new DomainEleOp(domainField, DomainEleOp::OPROW);
303  error_op->doWorkRhsHook = [&](DataOperator *op_ptr, int side, EntityType type,
304  EntitiesFieldData::EntData &data) {
305  MoFEMFunctionBegin;
306  auto o = static_cast<DomainEleOp *>(op_ptr);
307 
308  FTensor::Index<'i', 2> i;
309 
310  if (const size_t nb_dofs = data.getIndices().size()) {
311 
312  const int nb_integration_pts = o->getGaussPts().size2();
313  auto t_w = o->getFTensor0IntegrationWeight();
314  auto t_val = getFTensor0FromVec(*u_vals_ptr);
315  auto t_grad = getFTensor1FromMat<2>(*u_grad_ptr);
316  auto t_coords = o->getFTensor1CoordsAtGaussPts();
317 
318  std::array<double, LAST_NORM> error;
319  std::fill(error.begin(), error.end(), 0);
320 
321  for (int gg = 0; gg != nb_integration_pts; ++gg) {
322 
323  const double alpha = t_w * o->getMeasure();
324  const double diff =
325  t_val - u_exact(t_coords(0), t_coords(1), t_coords(2));
326 
327  auto t_exact_grad = u_grad_exact(t_coords(0), t_coords(1), t_coords(2));
328 
329  const double diff_grad =
330  (t_grad(i) - t_exact_grad(i)) * (t_grad(i) - t_exact_grad(i));
331 
332  error[L2] += alpha * pow(diff, 2);
333  error[SEMI_NORM] += alpha * diff_grad;
334 
335  ++t_w;
336  ++t_val;
337  ++t_grad;
338  ++t_coords;
339  }
340 
341  error[H1] = error[L2] + error[SEMI_NORM];
342 
343  CHKERR VecSetValues(l2_vec, LAST_NORM, error_indices.data(), error.data(),
344  ADD_VALUES);
345  }
346 
347  MoFEMFunctionReturn(0);
348  };
349 
350  auto side_fe_ptr = boost::make_shared<FaceSideEle>(mField);
351  add_base_ops(side_fe_ptr->getOpPtrVector());
352  side_fe_ptr->getOpPtrVector().push_back(
353  new OpCalculateScalarFieldValues(domainField, u_vals_ptr));
354  std::array<VectorDouble, 2> side_vals;
355  std::array<double, 2> area_map;
356 
357  auto side_vals_op = new DomainEleOp(domainField, DomainEleOp::OPROW);
358  side_vals_op->doWorkRhsHook = [&](DataOperator *op_ptr, int side,
359  EntityType type,
360  EntitiesFieldData::EntData &data) {
361  MoFEMFunctionBegin;
362  auto o = static_cast<FaceSideOp *>(op_ptr);
363 
364  const auto nb_in_loop = o->getFEMethod()->nInTheLoop;
365  area_map[nb_in_loop] = o->getMeasure();
366  side_vals[nb_in_loop] = *u_vals_ptr;
367  if (!nb_in_loop) {
368  area_map[1] = 0;
369  side_vals[1].clear();
370  }
371 
372  MoFEMFunctionReturn(0);
373  };
374  side_fe_ptr->getOpPtrVector().push_back(side_vals_op);
375 
376  auto do_work_rhs_error = [&](DataOperator *op_ptr, int side, EntityType type,
377  EntitiesFieldData::EntData &data) {
378  MoFEMFunctionBegin;
379  auto o = static_cast<BoundaryEleOp *>(op_ptr);
380 
381  CHKERR o->loopSideFaces("dFE", *side_fe_ptr);
382  const auto in_the_loop = side_fe_ptr->nInTheLoop;
383 
384 #ifndef NDEBUG
385  const std::array<std::string, 2> ele_type_name = {"BOUNDARY", "SKELETON"};
386  MOFEM_LOG("SELF", Sev::noisy)
387  << "do_work_rhs_error in_the_loop " << ele_type_name[in_the_loop];
388 #endif
389 
390  const double s = o->getMeasure() / (area_map[0] + area_map[1]);
391  const double p = penalty * s;
392 
393  constexpr std::array<int, 2> sign_array{1, -1};
394 
395  std::array<double, LAST_NORM> error;
396  std::fill(error.begin(), error.end(), 0);
397 
398  const int nb_integration_pts = o->getGaussPts().size2();
399 
400  if (!in_the_loop) {
401  side_vals[1].resize(nb_integration_pts, false);
402  auto t_coords = o->getFTensor1CoordsAtGaussPts();
403  auto t_val_m = getFTensor0FromVec(side_vals[1]);
404  for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
405  t_val_m = u_exact(t_coords(0), t_coords(1), t_coords(2));
406  ++t_coords;
407  ++t_val_m;
408  }
409  }
410 
411  auto t_val_p = getFTensor0FromVec(side_vals[0]);
412  auto t_val_m = getFTensor0FromVec(side_vals[1]);
413  auto t_w = o->getFTensor0IntegrationWeight();
414 
415  for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
416 
417  const double alpha = o->getMeasure() * t_w;
418  const auto diff = t_val_p - t_val_m;
419  error[SEMI_NORM] += alpha * p * diff * diff;
420 
421  ++t_w;
422  ++t_val_p;
423  ++t_val_m;
424  }
425 
426 
427  error[H1] = error[SEMI_NORM];
428  CHKERR VecSetValues(l2_vec, LAST_NORM, error_indices.data(), error.data(),
429  ADD_VALUES);
430 
431  MoFEMFunctionReturn(0);
432  };
433 
434  auto skeleton_error_op = new BoundaryEleOp(NOSPACE, BoundaryEleOp::OPSPACE);
435  skeleton_error_op->doWorkRhsHook = do_work_rhs_error;
436  auto boundary_error_op = new BoundaryEleOp(NOSPACE, BoundaryEleOp::OPSPACE);
437  boundary_error_op->doWorkRhsHook = do_work_rhs_error;
438 
439  pipeline_mng->getOpDomainRhsPipeline().push_back(error_op);
440  pipeline_mng->getOpSkeletonLhsPipeline().push_back(skeleton_error_op);
441  pipeline_mng->getOpBoundaryRhsPipeline().push_back(boundary_error_op);
442 
443  CHKERR pipeline_mng->loopFiniteElements();
444 
445  CHKERR VecAssemblyBegin(l2_vec);
446  CHKERR VecAssemblyEnd(l2_vec);
447 
448  if (mField.get_comm_rank() == 0) {
449  const double *array;
450  CHKERR VecGetArrayRead(l2_vec, &array);
451  MOFEM_LOG_C("SELF", Sev::inform, "Error Norm L2 %6.4e",
452  std::sqrt(array[L2]));
453  MOFEM_LOG_C("SELF", Sev::inform, "Error Norm Energetic %6.4e",
454  std::sqrt(array[SEMI_NORM]));
455  MOFEM_LOG_C("SELF", Sev::inform, "Error Norm H1 %6.4e",
456  std::sqrt(array[H1]));
457 
458  if(is_test) {
459  constexpr double eps = 1e-12;
460  if (std::sqrt(array[H1]) > eps)
461  SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID, "Error is too big");
462  }
463 
464  CHKERR VecRestoreArrayRead(l2_vec, &array);
465  const MoFEM::Problem *problem_ptr;
466  CHKERR DMMoFEMGetProblemPtr(simpleInterface->getDM(), &problem_ptr);
467  MOFEM_LOG_C("SELF", Sev::inform, "Nb. DOFs %d",
468  problem_ptr->getNbDofsRow());
469  }
470 
471 
472 
473  MoFEMFunctionReturn(0);
474 }

◆ outputResults()

MoFEMErrorCode Poisson2DiscontGalerkin::outputResults ( )
private

[Output results]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 477 of file poisson_2d_dis_galerkin.cpp.

477  {
478  MoFEMFunctionBegin;
479 
480  auto pipeline_mng = mField.getInterface<PipelineManager>();
481  pipeline_mng->getDomainLhsFE().reset();
482  pipeline_mng->getSkeletonRhsFE().reset();
483  pipeline_mng->getSkeletonLhsFE().reset();
484  pipeline_mng->getBoundaryRhsFE().reset();
485  pipeline_mng->getBoundaryLhsFE().reset();
486 
487  auto post_proc_fe = boost::make_shared<PostProcEle>(mField);
488 
489  auto u_ptr = boost::make_shared<VectorDouble>();
490  post_proc_fe->getOpPtrVector().push_back(
491  new OpCalculateScalarFieldValues(domainField, u_ptr));
492 
493  using OpPPMap = OpPostProcMapInMoab<SPACE_DIM, SPACE_DIM>;
494 
495  post_proc_fe->getOpPtrVector().push_back(
496 
497  new OpPPMap(
498 
499  post_proc_fe->getPostProcMesh(), post_proc_fe->getMapGaussPts(),
500 
501  {{"U", u_ptr}},
502 
503  {},
504 
505  {},
506 
507  {})
508 
509  );
510 
511  pipeline_mng->getDomainRhsFE() = post_proc_fe;
512  CHKERR pipeline_mng->loopFiniteElements();
513  CHKERR post_proc_fe->writeFile("out_result.h5m");
514 
515  MoFEMFunctionReturn(0);
516 }

◆ readMesh()

MoFEMErrorCode Poisson2DiscontGalerkin::readMesh ( )
private

[Read mesh]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 104 of file poisson_2d_dis_galerkin.cpp.

104  {
105  MoFEMFunctionBegin;
106 
107  CHKERR mField.getInterface(simpleInterface);
108  CHKERR simpleInterface->getOptions();
109 
110  // Only L2 field is set in this example. Two lines bellow forces simple
111  // interface to creat lower dimension (edge) elements, despite that fact that
112  // there is no field spanning on such elements. We need them for DG method.
113  simpleInterface->getAddSkeletonFE() = true;
114  simpleInterface->getAddBoundaryFE() = true;
115 
116  CHKERR simpleInterface->loadFile();
117 
118  MoFEMFunctionReturn(0);
119 }

◆ runProgram()

MoFEMErrorCode Poisson2DiscontGalerkin::runProgram ( )

[Output results]

[Run program]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 520 of file poisson_2d_dis_galerkin.cpp.

520  {
521  MoFEMFunctionBegin;
522 
523  CHKERR readMesh();
524  CHKERR setupProblem();
525  CHKERR boundaryCondition();
526  CHKERR assembleSystem();
527  CHKERR setIntegrationRules();
528  CHKERR solveSystem();
529  CHKERR outputResults();
530  CHKERR checkResults();
531 
532  MoFEMFunctionReturn(0);
533 }

◆ setIntegrationRules()

MoFEMErrorCode Poisson2DiscontGalerkin::setIntegrationRules ( )
private

[Assemble system]

[Set integration rules]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 215 of file poisson_2d_dis_galerkin.cpp.

215  {
216  MoFEMFunctionBegin;
217 
218  auto rule_lhs = [](int, int, int p) -> int { return 2 * p; };
219  auto rule_rhs = [](int, int, int p) -> int { return 2 * p; };
220  auto rule_2 = [this](int, int, int) { return 2 * oRder; };
221 
222  auto pipeline_mng = mField.getInterface<PipelineManager>();
223  CHKERR pipeline_mng->setDomainLhsIntegrationRule(rule_lhs);
224  CHKERR pipeline_mng->setDomainRhsIntegrationRule(rule_rhs);
225 
226  CHKERR pipeline_mng->setSkeletonLhsIntegrationRule(rule_2);
227  CHKERR pipeline_mng->setSkeletonRhsIntegrationRule(rule_2);
228  CHKERR pipeline_mng->setBoundaryLhsIntegrationRule(rule_2);
229  CHKERR pipeline_mng->setBoundaryRhsIntegrationRule(rule_2);
230 
231  MoFEMFunctionReturn(0);
232 }

◆ setupProblem()

MoFEMErrorCode Poisson2DiscontGalerkin::setupProblem ( )
private

[Read mesh]

[Setup problem]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 123 of file poisson_2d_dis_galerkin.cpp.

123  {
124  MoFEMFunctionBegin;
125 
126  CHKERR PetscOptionsGetInt(PETSC_NULL, "", "-order", &oRder, PETSC_NULL);
127  CHKERR PetscOptionsGetScalar(PETSC_NULL, "", "-penalty", &penalty,
128  PETSC_NULL);
129  CHKERR PetscOptionsGetScalar(PETSC_NULL, "", "-phi", &phi, PETSC_NULL);
130  CHKERR PetscOptionsGetScalar(PETSC_NULL, "", "-nitsche", &nitsche,
131  PETSC_NULL);
132  PetscOptionsGetBool(PETSC_NULL, "", "-is_test", &is_test, PETSC_NULL);
133 
134  MOFEM_LOG("WORLD", Sev::inform) << "Set order: " << oRder;
135  MOFEM_LOG("WORLD", Sev::inform) << "Set penalty: " << penalty;
136  MOFEM_LOG("WORLD", Sev::inform) << "Set phi: " << phi;
137  MOFEM_LOG("WORLD", Sev::inform) << "Set nitche: " << nitsche;
138  MOFEM_LOG("WORLD", Sev::inform) << "Set test: " << (is_test == PETSC_TRUE);
139 
140  CHKERR simpleInterface->addDomainField(domainField, L2,
141  AINSWORTH_LOBATTO_BASE, 1);
142  CHKERR simpleInterface->setFieldOrder(domainField, oRder);
143  CHKERR simpleInterface->setUp();
144 
145  // This is only for debigging and experimentation, to see boundary and edge
146  // elements.
147  auto save_shared = [&](auto meshset, std::string prefix) {
148  MoFEMFunctionBegin;
149  auto file_name =
150  prefix + "_" +
151  boost::lexical_cast<std::string>(mField.get_comm_rank()) + ".vtk";
152  CHKERR mField.get_moab().write_file(file_name.c_str(), "VTK", "", &meshset,
153  1);
154  MoFEMFunctionReturn(0);
155  };
156 
157  // CHKERR save_shared(simpleInterface->getBoundaryMeshSet(), "bdy");
158  // CHKERR save_shared(simpleInterface->getSkeletonMeshSet(), "skel");
159 
160  MoFEMFunctionReturn(0);
161 }

◆ solveSystem()

MoFEMErrorCode Poisson2DiscontGalerkin::solveSystem ( )
private

[Set integration rules]

[Solve system]

Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 236 of file poisson_2d_dis_galerkin.cpp.

236  {
237  MoFEMFunctionBegin;
238 
239  auto pipeline_mng = mField.getInterface<PipelineManager>();
240 
241  auto ksp_solver = pipeline_mng->createKSP();
242  CHKERR KSPSetFromOptions(ksp_solver);
243  CHKERR KSPSetUp(ksp_solver);
244 
245  // Create RHS and solution vectors
246  auto dm = simpleInterface->getDM();
247  auto F = createDMVector(dm);
248  auto D = vectorDuplicate(F);
249 
250  CHKERR KSPSolve(ksp_solver, F, D);
251 
252  // Scatter result data on the mesh
253  CHKERR VecGhostUpdateBegin(D, INSERT_VALUES, SCATTER_FORWARD);
254  CHKERR VecGhostUpdateEnd(D, INSERT_VALUES, SCATTER_FORWARD);
255  CHKERR DMoFEMMeshToLocalVector(dm, D, INSERT_VALUES, SCATTER_REVERSE);
256 
257  MoFEMFunctionReturn(0);
258 }

Member Data Documentation

◆ domainField

std::string Poisson2DiscontGalerkin::domainField
private
Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 96 of file poisson_2d_dis_galerkin.cpp.

◆ mField

MoFEM::Interface& Poisson2DiscontGalerkin::mField
private
Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 92 of file poisson_2d_dis_galerkin.cpp.

◆ oRder

int Poisson2DiscontGalerkin::oRder
private
Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 97 of file poisson_2d_dis_galerkin.cpp.

◆ simpleInterface

Simple* Poisson2DiscontGalerkin::simpleInterface
private
Examples
poisson_2d_dis_galerkin.cpp.

Definition at line 93 of file poisson_2d_dis_galerkin.cpp.

The documentation for this struct was generated from the following file:
NOSPACE
@ NOSPACE
Definition: definitions.h:83
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.
Definition: UnknownInterface.hpp:93
MoFEM::EntitiesFieldData::EntData
Data on single entity (This is passed as argument to DataOperator::doWork)
Definition: EntitiesFieldData.hpp:127
Poisson2DiscontGalerkin::solveSystem
MoFEMErrorCode solveSystem()
[Set integration rules]
Definition: poisson_2d_dis_galerkin.cpp:236
Poisson2DiscontGalerkin::readMesh
MoFEMErrorCode readMesh()
[Read mesh]
Definition: poisson_2d_dis_galerkin.cpp:104
MoFEM::PipelineManager::getDomainRhsFE
boost::shared_ptr< FEMethod > & getDomainRhsFE()
Definition: PipelineManager.hpp:405
H1
@ H1
continuous field
Definition: definitions.h:85
Poisson2DiscontGalerkin::checkResults
MoFEMErrorCode checkResults()
[Solve system]
Definition: poisson_2d_dis_galerkin.cpp:261
MoFEM::DataOperator
base operator to do operations at Gauss Pt. level
Definition: DataOperators.hpp:24
u_grad_exact
auto u_grad_exact
Definition: poisson_2d_dis_galerkin.cpp:49
MoFEM::ForcesAndSourcesCore::UserDataOperator::OPSPACE
@ OPSPACE
operator do Work is execute on space data
Definition: ForcesAndSourcesCore.hpp:570
MoFEM::CoreInterface::get_comm
virtual MPI_Comm & get_comm() const =0
L2
@ L2
field with C-1 continuity
Definition: definitions.h:88
phi
static double phi
Definition: poisson_2d_dis_galerkin.cpp:29
MoFEM::CoreInterface::get_comm_rank
virtual int get_comm_rank() const =0
MoFEM::Simple::loadFile
MoFEMErrorCode loadFile(const std::string options, const std::string mesh_file_name, LoadFileFunc loadFunc=defaultLoadFileFunc)
Load mesh file.
Definition: Simple.cpp:194
MoFEM::PipelineManager
PipelineManager interface.
Definition: PipelineManager.hpp:24
MoFEM::DMoFEMMeshToLocalVector
PetscErrorCode DMoFEMMeshToLocalVector(DM dm, Vec l, InsertMode mode, ScatterMode scatter_mode)
set local (or ghosted) vector values on mesh for partition only
Definition: DMMoFEM.cpp:527
Poisson2DiscontGalerkinOperators::OpL2LhsPenalty
Operator the left hand side matrix.
Definition: PoissonDiscontinousGalerkin.hpp:115
OpDomainGradGrad
FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpGradGrad< BASE_DIM, FIELD_DIM, SPACE_DIM > OpDomainGradGrad
Definition: poisson_2d_dis_galerkin.cpp:36
DomainEleOp
DomainEle::UserDataOperator DomainEleOp
Definition: poisson_2d_dis_galerkin.cpp:17
MoFEM::VecSetValues
MoFEMErrorCode VecSetValues(Vec V, const EntitiesFieldData::EntData &data, const double *ptr, InsertMode iora)
Assemble PETSc vector.
Definition: EntitiesFieldData.hpp:1576
Poisson2DiscontGalerkinOperators::OpL2BoundaryRhs
Operator to evaluate Dirichlet boundary conditions using DG.
Definition: PoissonDiscontinousGalerkin.hpp:274
MoFEM::OpCalculateScalarFieldGradient
Get field gradients at integration pts for scalar filed rank 0, i.e. vector field.
Definition: UserDataOperators.hpp:1294
Poisson2DiscontGalerkin::setupProblem
MoFEMErrorCode setupProblem()
[Read mesh]
Definition: poisson_2d_dis_galerkin.cpp:123
MoFEM::OpSetInvJacL2ForFace
Definition: UserDataOperators.hpp:2929
MoFEM::Simple::getOptions
MoFEMErrorCode getOptions()
get options
Definition: Simple.cpp:180
u_exact
auto u_exact
Definition: poisson_2d_dis_galerkin.cpp:42
MoFEM::Simple::getAddSkeletonFE
bool & getAddSkeletonFE()
Get the addSkeletonFE.
Definition: Simple.hpp:425
Poisson2DiscontGalerkin::simpleInterface
Simple * simpleInterface
Definition: poisson_2d_dis_galerkin.cpp:93
MoFEM::Simple::getDM
MoFEMErrorCode getDM(DM *dm)
Get DM.
Definition: Simple.cpp:667
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:535
Poisson2DiscontGalerkin::domainField
std::string domainField
Definition: poisson_2d_dis_galerkin.cpp:96
MoFEM::createDMVector
auto createDMVector(DM dm)
Get smart vector from DM.
Definition: DMMoFEM.hpp:1018
Poisson2DiscontGalerkin::oRder
int oRder
Definition: poisson_2d_dis_galerkin.cpp:97
penalty
static double penalty
Definition: poisson_2d_dis_galerkin.cpp:28
MoFEM::CoreInterface::get_moab
virtual moab::Interface & get_moab()=0
FaceSideOp
MoFEM::FaceElementForcesAndSourcesCore::UserDataOperator
default operator for TRI element
Definition: FaceElementForcesAndSourcesCore.hpp:94
MOFEM_LOG_C
#define MOFEM_LOG_C(channel, severity, format,...)
Definition: LogManager.hpp:311
convert.type
type
Definition: convert.py:64
MoFEM::Problem::getNbDofsRow
DofIdx getNbDofsRow() const
Definition: ProblemsMultiIndices.hpp:376
OpPPMap
OpPostProcMapInMoab< SPACE_DIM, SPACE_DIM > OpPPMap
Definition: photon_diffusion.cpp:29
MoFEM::getFTensor0FromVec
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
Definition: Templates.hpp:135
MoFEM::OpCalculateScalarFieldValues
Get value at integration points for scalar field.
Definition: UserDataOperators.hpp:82
MoFEM::Simple::addDomainField
MoFEMErrorCode addDomainField(const std::string &name, const FieldSpace space, const FieldApproximationBase base, const FieldCoefficientsNumber nb_of_coefficients, const TagType tag_type=MB_TAG_SPARSE, const enum MoFEMTypes bh=MF_ZERO, int verb=-1)
Add field on domain.
Definition: Simple.cpp:264
MoFEM::PipelineManager::createKSP
SmartPetscObj< KSP > createKSP(SmartPetscObj< DM > dm=nullptr)
Create KSP (linear) solver.
Definition: PipelineManager.cpp:49
MoFEM::Simple::getAddBoundaryFE
bool & getAddBoundaryFE()
Get the addSkeletonFE.
Definition: Simple.hpp:436
Poisson2DiscontGalerkinOperators::OpCalculateSideData
Operator tp collect data from elements on the side of Edge/Face.
Definition: PoissonDiscontinousGalerkin.hpp:38
OpDomainSource
FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpSource< BASE_DIM, FIELD_DIM > OpDomainSource
Definition: poisson_2d_dis_galerkin.cpp:38
nitsche
static double nitsche
Definition: poisson_2d_dis_galerkin.cpp:31
AINSWORTH_LOBATTO_BASE
@ AINSWORTH_LOBATTO_BASE
Definition: definitions.h:62
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
FTensor::Index< 'i', 2 >
MoFEM::Simple::setFieldOrder
MoFEMErrorCode setFieldOrder(const std::string field_name, const int order, const Range *ents=NULL)
Set field order.
Definition: Simple.cpp:473
MoFEM::PipelineManager::getDomainLhsFE
boost::shared_ptr< FEMethod > & getDomainLhsFE()
Definition: PipelineManager.hpp:401
DomainEleOp
MOFEM_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition: LogManager.hpp:308
MoFEM::vectorDuplicate
SmartPetscObj< Vec > vectorDuplicate(Vec vec)
Create duplicate vector of smart vector.
Definition: PetscSmartObj.hpp:217
BoundaryEleOp
BoundaryEle::UserDataOperator BoundaryEleOp
Definition: poisson_2d_dis_galerkin.cpp:21
Poisson2DiscontGalerkin::boundaryCondition
MoFEMErrorCode boundaryCondition()
[Setup problem]
Definition: poisson_2d_dis_galerkin.cpp:165
Poisson2DiscontGalerkin::assembleSystem
MoFEMErrorCode assembleSystem()
[Boundary condition]
Definition: poisson_2d_dis_galerkin.cpp:173
eps
static const double eps
Definition: check_base_functions_derivatives_on_tet.cpp:11
Poisson2DiscontGalerkin::mField
MoFEM::Interface & mField
Definition: poisson_2d_dis_galerkin.cpp:92
MoFEM::DMMoFEMGetProblemPtr
PetscErrorCode DMMoFEMGetProblemPtr(DM dm, const MoFEM::Problem **problem_ptr)
Get pointer to problem data structure.
Definition: DMMoFEM.cpp:430
MoFEM::createVectorMPI
auto createVectorMPI(MPI_Comm comm, PetscInt n, PetscInt N)
Create MPI Vector.
Definition: PetscSmartObj.hpp:198
MoFEM::OpInvertMatrix
Definition: UserDataOperators.hpp:3254
MoFEM::OpCalculateHOJacForFace
OpCalculateHOJacForFaceImpl< 2 > OpCalculateHOJacForFace
Definition: HODataOperators.hpp:264
ReactionDiffusionEquation::D
const double D
diffusivity
Definition: reaction_diffusion.cpp:20
Poisson2DiscontGalerkin::outputResults
MoFEMErrorCode outputResults()
[Output results]
Definition: poisson_2d_dis_galerkin.cpp:477
MOFEM_ATOM_TEST_INVALID
@ MOFEM_ATOM_TEST_INVALID
Definition: definitions.h:40
MoFEM::PetscOptionsGetScalar
PetscErrorCode PetscOptionsGetScalar(PetscOptions *, const char pre[], const char name[], PetscScalar *dval, PetscBool *set)
Definition: DeprecatedPetsc.hpp:162
MoFEM::Problem
keeps basic data about problem
Definition: ProblemsMultiIndices.hpp:54
is_test
PetscBool is_test
Definition: poisson_2d_dis_galerkin.cpp:40
source
auto source
Definition: poisson_2d_dis_galerkin.cpp:61
MoFEM::Simple::setUp
MoFEMErrorCode setUp(const PetscBool is_partitioned=PETSC_TRUE)
Setup problem.
Definition: Simple.cpp:611
convert.int
int
Definition: convert.py:64
MoFEM::PetscOptionsGetInt
PetscErrorCode PetscOptionsGetInt(PetscOptions *, const char pre[], const char name[], PetscInt *ivalue, PetscBool *set)
Definition: DeprecatedPetsc.hpp:142
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:416
Poisson2DiscontGalerkin::setIntegrationRules
MoFEMErrorCode setIntegrationRules()
[Assemble system]
Definition: poisson_2d_dis_galerkin.cpp:215
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:346
F
@ F
Definition: free_surface.cpp:394
MoFEM::PetscOptionsGetBool
PetscErrorCode PetscOptionsGetBool(PetscOptions *, const char pre[], const char name[], PetscBool *bval, PetscBool *set)
Definition: DeprecatedPetsc.hpp:182
MoFEM::OpPostProcMapInMoab
Post post-proc data at points from hash maps.
Definition: PostProcBrokenMeshInMoabBase.hpp:698
Generated by Doxygen 1.8.17 and hosted at University of Glasgow