v0.14.0
approx_sphere.cpp
Go to the documentation of this file.
1 /**
2  * \file approx_sphere.cpp
3  * \example approx_sphere.cpp
4  *
5  */
6 
7 #include <MoFEM.hpp>
8 
9 using namespace MoFEM;
10 
11 static char help[] = "...\n\n";
12 
13 #include <BasicFiniteElements.hpp>
14 
15 constexpr int FM_DIM = 2;
16 
17 template <int DIM> struct ElementsAndOps {};
18 
19 template <> struct ElementsAndOps<2> {
20  using DomainEle = MoFEM::FaceElementForcesAndSourcesCore;
21 };
22 
23 using DomainEle = ElementsAndOps<FM_DIM>::DomainEle;
24 using DomainEleOp = DomainEle::UserDataOperator;
25 using EntData = EntitiesFieldData::EntData;
26 
27 using AssemblyDomainEleOp =
28  FormsIntegrators<DomainEleOp>::Assembly<PETSC>::OpBase;
29 
30 constexpr double a = 1;
31 constexpr double a2 = a * a;
32 constexpr double a4 = a2 * a2;
33 
34 constexpr double A = 6371220;
35 
36 FTensor::Index<'i', 3> i;
37 FTensor::Index<'j', 3> j;
38 FTensor::Index<'k', 3> k;
39 
40 auto res_J = [](const double x, const double y, const double z) {
41  const double res = (x * x + y * y + z * z - a2);
42  return res;
43 };
44 
45 auto res_J_dx = [](const double x, const double y, const double z) {
46  const double res = res_J(x, y, z);
47  return FTensor::Tensor1<double, 3>{res * (2 * x), res * (2 * y),
48  res * (2 * z)};
49 };
50 
51 auto lhs_J_dx2 = [](const double x, const double y, const double z) {
52  const double res = res_J(x, y, z);
53  return FTensor::Tensor2<double, 3, 3>{
54 
55  (res * 2 + (4 * x * x)),
56  (4 * y * x),
57  (4 * z * x),
58 
59  (4 * x * y),
60  (2 * res + (4 * y * y)),
61  (4 * z * y),
62 
63  (4 * x * z),
64  (4 * y * z),
65  (2 * res + (4 * z * z))};
66 };
67 
68 struct OpRhs : public AssemblyDomainEleOp {
69 
70  OpRhs(const std::string field_name, boost::shared_ptr<MatrixDouble> x_ptr,
71  boost::shared_ptr<MatrixDouble> dot_x_ptr)
72  : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
73  xPtr(x_ptr), xDotPtr(dot_x_ptr) {}
74 
75  MoFEMErrorCode iNtegrate(EntitiesFieldData::EntData &row_data) {
76  MoFEMFunctionBegin;
77 
78  auto t_w = getFTensor0IntegrationWeight();
79  auto t_row_base = row_data.getFTensor0N();
80 
81  auto t_x0 = getFTensor1CoordsAtGaussPts();
82  auto t_x = getFTensor1FromMat<3>(*xPtr);
83  auto t_dot_x = getFTensor1FromMat<3>(*xDotPtr);
84  auto t_normal = getFTensor1NormalsAtGaussPts();
85 
86  constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
87 
88  for (int gg = 0; gg != nbIntegrationPts; gg++) {
89 
90  constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
91  FTensor::Tensor1<double, 3> t_n{t_x0(0), t_x0(1), t_x0(2)};
92  t_n.normalize();
93  FTensor::Tensor2<double, 3, 3> t_P, t_Q;
94  t_P(i, j) = t_n(i) * t_n(j);
95  t_Q(i, j) = t_kd(i, j) - t_P(i, j);
96 
97  auto t_J_res = res_J_dx(t_x(0), t_x(1), t_x(2));
98 
99  const double alpha = t_w;
100  auto t_nf = getFTensor1FromArray<3, 3>(locF);
101  double l = std::sqrt(t_normal(i) * t_normal(i));
102 
103  FTensor::Tensor1<double, 3> t_res;
104  t_res(i) =
105  alpha * l * ((t_P(i, k) * t_J_res(k) + t_Q(i, k) * t_dot_x(k)));
106 
107  int rr = 0;
108  for (; rr != nbRows / 3; ++rr) {
109 
110  t_nf(j) += t_row_base * t_res(j);
111 
112  ++t_row_base;
113  ++t_nf;
114  }
115  for (; rr < nbRowBaseFunctions; ++rr) {
116  ++t_row_base;
117  }
118 
119  ++t_w;
120  ++t_x;
121  ++t_dot_x;
122  ++t_x0;
123  ++t_normal;
124  }
125 
126  MoFEMFunctionReturn(0);
127  }
128 
129 private:
130  boost::shared_ptr<MatrixDouble> xPtr;
131  boost::shared_ptr<MatrixDouble> xDotPtr;
132 };
133 
134 struct OpLhs : public AssemblyDomainEleOp {
135 
136  OpLhs(const std::string field_name, boost::shared_ptr<MatrixDouble> x_ptr,
137  boost::shared_ptr<MatrixDouble> dot_x_ptr)
138  : AssemblyDomainEleOp(field_name, field_name,
139  AssemblyDomainEleOp::OPROWCOL),
140  xPtr(x_ptr), xDotPtr(dot_x_ptr) {
141  this->sYmm = false;
142  }
143 
144  MoFEMErrorCode iNtegrate(EntitiesFieldData::EntData &row_data,
145  EntitiesFieldData::EntData &col_data) {
146  MoFEMFunctionBegin;
147 
148  auto t_w = getFTensor0IntegrationWeight();
149  auto t_row_base = row_data.getFTensor0N();
150 
151  auto t_x0 = getFTensor1CoordsAtGaussPts();
152  auto t_x = getFTensor1FromMat<3>(*xPtr);
153  auto t_dot_x = getFTensor1FromMat<3>(*xDotPtr);
154  auto t_normal = getFTensor1NormalsAtGaussPts();
155 
156  auto get_t_mat = [&](const int rr) {
157  return FTensor::Tensor2<FTensor::PackPtr<double *, 3>, 3, 3>{
158  &locMat(rr + 0, 0), &locMat(rr + 0, 1), &locMat(rr + 0, 2),
159 
160  &locMat(rr + 1, 0), &locMat(rr + 1, 1), &locMat(rr + 1, 2),
161 
162  &locMat(rr + 2, 0), &locMat(rr + 2, 1), &locMat(rr + 2, 2)};
163  };
164 
165  constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
166  const double ts_a = getTSa();
167 
168  for (int gg = 0; gg != nbIntegrationPts; gg++) {
169 
170  constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
171  FTensor::Tensor1<double, 3> t_n{t_x0(0), t_x0(1), t_x0(2)};
172  t_n.normalize();
173  FTensor::Tensor2<double, 3, 3> t_P, t_Q;
174  t_P(i, j) = t_n(i) * t_n(j);
175  t_Q(i, j) = t_kd(i, j) - t_P(i, j);
176 
177  auto t_J_lhs = lhs_J_dx2(t_x(0), t_x(1), t_x(2));
178  double l = std::sqrt(t_normal(i) * t_normal(i));
179 
180  const double alpha = t_w;
181  FTensor::Tensor2<double, 3, 3> t_lhs;
182  t_lhs(i, j) =
183  (alpha * l) * (t_P(i, k) * t_J_lhs(k, j) + t_Q(i, j) * ts_a);
184 
185  int rr = 0;
186  for (; rr != nbRows / 3; rr++) {
187 
188  auto t_col_base = col_data.getFTensor0N(gg, 0);
189  auto t_mat = get_t_mat(3 * rr);
190 
191  for (int cc = 0; cc != nbCols / 3; cc++) {
192 
193  const double rc = t_row_base * t_col_base;
194  t_mat(i, j) += rc * t_lhs(i, j);
195 
196  ++t_col_base;
197  ++t_mat;
198  }
199  ++t_row_base;
200  }
201 
202  for (; rr < nbRowBaseFunctions; ++rr)
203  ++t_row_base;
204 
205  ++t_w;
206  ++t_x;
207  ++t_x0;
208  ++t_normal;
209  }
210 
211  MoFEMFunctionReturn(0);
212  }
213 
214 private:
215  boost::shared_ptr<MatrixDouble> xPtr;
216  boost::shared_ptr<MatrixDouble> xDotPtr;
217 };
218 
219 struct OpError : public DomainEleOp {
220 
221  OpError(const std::string field_name, boost::shared_ptr<MatrixDouble> x_ptr)
222  : DomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
223  xPtr(x_ptr) {
224 
225  std::fill(&doEntities[MBEDGE], &doEntities[MBMAXTYPE], false);
226  }
227 
228  MoFEMErrorCode doWork(int side, EntityType type, EntData &data) {
229 
230  MoFEMFunctionBegin;
231 
232  auto t_w = getFTensor0IntegrationWeight();
233  auto t_x = getFTensor1FromMat<3>(*xPtr);
234  auto t_normal = getFTensor1NormalsAtGaussPts();
235  auto nb_integration_pts = getGaussPts().size2();
236 
237  double error = 0;
238 
239  for (int gg = 0; gg != nb_integration_pts; gg++) {
240 
241  double l = std::sqrt(t_normal(i) * t_normal(i));
242  error += t_w * l * std::abs((t_x(i) * t_x(i) - A * A));
243 
244  ++t_w;
245  ++t_x;
246  ++t_normal;
247  }
248 
249  CHKERR VecSetValue(errorVec, 0, error, ADD_VALUES);
250 
251  MoFEMFunctionReturn(0);
252  }
253 
254  static SmartPetscObj<Vec> errorVec;
255 
256 private:
257  boost::shared_ptr<MatrixDouble> xPtr;
258 };
259 
260 SmartPetscObj<Vec> OpError::errorVec;
261 
262 struct ApproxSphere {
263 
264  ApproxSphere(MoFEM::Interface &m_field) : mField(m_field) {}
265 
266  MoFEMErrorCode runProblem();
267 
268 private:
269  MoFEM::Interface &mField;
270 
271  MoFEMErrorCode getOptions();
272  MoFEMErrorCode readMesh();
273  MoFEMErrorCode setupProblem();
274  MoFEMErrorCode setOPs();
275  MoFEMErrorCode solveSystem();
276  MoFEMErrorCode outputResults();
277 };
278 
279 //! [Run programme]
280 MoFEMErrorCode ApproxSphere::runProblem() {
281  MoFEMFunctionBegin;
282  CHKERR getOptions();
283  CHKERR readMesh();
284  CHKERR setupProblem();
285  CHKERR setOPs();
286  CHKERR solveSystem();
287  CHKERR outputResults();
288  MoFEMFunctionReturn(0);
289 }
290 //! [Run programme]
291 
292 MoFEMErrorCode ApproxSphere::getOptions() {
293  MoFEMFunctionBeginHot;
294  MoFEMFunctionReturnHot(0);
295 }
296 
297 //! [Read mesh]
298 MoFEMErrorCode ApproxSphere::readMesh() {
299  MoFEMFunctionBegin;
300  auto simple = mField.getInterface<Simple>();
301  CHKERR simple->getOptions();
302  CHKERR simple->loadFile();
303  MoFEMFunctionReturn(0);
304 }
305 //! [Read mesh]
306 
307 //! [Set up problem]
308 MoFEMErrorCode ApproxSphere::setupProblem() {
309  MoFEMFunctionBegin;
310 
311  auto simple = mField.getInterface<Simple>();
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 
320  MoFEMFunctionReturn(0);
321 }
322 //! [Set up problem]
323 
324 //! [Push operators to pipeline]
325 MoFEMErrorCode ApproxSphere::setOPs() {
326  MoFEMFunctionBegin;
327  PipelineManager *pipeline_mng = mField.getInterface<PipelineManager>();
328 
329  auto integration_rule = [](int, int, int approx_order) {
330  return 3 * approx_order;
331  };
332  CHKERR pipeline_mng->setDomainRhsIntegrationRule(integration_rule);
333  CHKERR pipeline_mng->setDomainLhsIntegrationRule(integration_rule);
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 
356  MoFEMFunctionReturn(0);
357 }
358 //! [Push operators to pipeline]
359 
360 //! [Solve]
361 MoFEMErrorCode ApproxSphere::solveSystem() {
362  MoFEMFunctionBegin;
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();
372  MoFEM::SmartPetscObj<TS> ts;
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 = createDMVector(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 
390  MoFEMFunctionReturn(0);
391 }
392 
393 //! [Solve]
394 MoFEMErrorCode ApproxSphere::outputResults() {
395  MoFEMFunctionBegin;
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 
402  auto simple = mField.getInterface<Simple>();
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 
411  using OpPPMap = OpPostProcMapInMoab<3, 3>;
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 = [&]() {
440  MoFEMFunctionBegin;
441  MOFEM_LOG("EXAMPLE", Sev::inform) << "Create vec ";
442  OpError::errorVec = createVectorMPI(
443  mField.get_comm(), (!mField.get_comm_rank()) ? 1 : 0, 1);
444  VecZeroEntries(OpError::errorVec);
445  MoFEMFunctionReturn(0);
446  };
447 
448  error_fe->postProcessHook = [&]() {
449  MoFEMFunctionBegin;
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();
457  MoFEMFunctionReturn(0);
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");
469  MoFEMFunctionReturn(0);
470 }
471 //! [Postprocess results]
472 
473 int main(int argc, char *argv[]) {
474 
475  // Initialisation of MoFEM/PETSc and MOAB data structures
476  const char param_file[] = "param_file.petsc";
477  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
478 
479  auto core_log = logging::core::get();
480  core_log->add_sink(
481  LogManager::createSink(LogManager::getStrmWorld(), "EXAMPLE"));
482  LogManager::setLog("EXAMPLE");
483  MOFEM_LOG_TAG("EXAMPLE", "example");
484 
485  try {
486 
487  //! [Register MoFEM discrete manager in PETSc]
488  DMType dm_name = "DMMOFEM";
489  CHKERR DMRegister_MoFEM(dm_name);
490  //! [Register MoFEM discrete manager in PETSc
491 
492  //! [Create MoAB]
493  moab::Core mb_instance; ///< mesh database
494  moab::Interface &moab = mb_instance; ///< mesh database interface
495  //! [Create MoAB]
496 
497  //! [Create MoFEM]
498  MoFEM::Core core(moab); ///< finite element database
499  MoFEM::Interface &m_field = core; ///< finite element database insterface
500  //! [Create MoFEM]
501 
502  //! [ApproxSphere]
503  ApproxSphere ex(m_field);
504  CHKERR ex.runProblem();
505  //! [ApproxSphere]
506  }
507  CATCH_ERRORS;
508 
509  MoFEM::Core::Finalize();
510 }
help
static char help[]
Definition: approx_sphere.cpp:11
MoFEMFunctionReturnHot
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:447
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
MoFEM::CoreTmp< 0 >
Core (interface) class.
Definition: Core.hpp:82
H1
@ H1
continuous field
Definition: definitions.h:85
OpError
Definition: initial_diffusion.cpp:61
FTensor::Tensor1< double, 3 >
ApproxSphere::getOptions
MoFEMErrorCode getOptions()
[Run programme]
Definition: approx_sphere.cpp:292
ApproxSphere::ApproxSphere
ApproxSphere(MoFEM::Interface &m_field)
Definition: approx_sphere.cpp:264
OpError::OpError
OpError(const std::string field_name, boost::shared_ptr< MatrixDouble > x_ptr)
Definition: approx_sphere.cpp:221
OpLhs::xDotPtr
boost::shared_ptr< MatrixDouble > xDotPtr
Definition: approx_sphere.cpp:216
ApproxSphere::setupProblem
MoFEMErrorCode setupProblem()
[Read mesh]
Definition: approx_sphere.cpp:308
ApproxSphere::readMesh
MoFEMErrorCode readMesh()
[Read mesh]
Definition: approx_sphere.cpp:298
MoFEM::Exceptions::MoFEMErrorCode
PetscErrorCode MoFEMErrorCode
MoFEM/PETSc error code.
Definition: Exceptions.hpp:56
MoFEM::OpCalculateVectorFieldValues
Get values at integration pts for tensor filed rank 1, i.e. vector field.
Definition: UserDataOperators.hpp:466
MoFEM::PipelineManager::getOpDomainRhsPipeline
boost::ptr_deque< UserDataOperator > & getOpDomainRhsPipeline()
Get the Op Domain Rhs Pipeline object.
Definition: PipelineManager.hpp:773
FTensor::Kronecker_Delta
Kronecker Delta class.
Definition: Kronecker_Delta.hpp:15
OpError::errorVec
static SmartPetscObj< Vec > errorVec
Definition: approx_sphere.cpp:254
MoFEM::PETSC
@ PETSC
Definition: FormsIntegrators.hpp:104
MoFEM::PipelineManager
PipelineManager interface.
Definition: PipelineManager.hpp:24
OpRhs::OpRhs
OpRhs(const std::string field_name, boost::shared_ptr< MatrixDouble > x_ptr, boost::shared_ptr< MatrixDouble > dot_x_ptr)
Definition: approx_sphere.cpp:70
MoFEM.hpp
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
BasicFiniteElements.hpp
MoFEM::Projection10NodeCoordsOnField
Projection of edge entities with one mid-node on hierarchical basis.
Definition: Projection10NodeCoordsOnField.hpp:24
MoFEM::CoreTmp< 0 >::Finalize
static MoFEMErrorCode Finalize()
Checks for options to be called at the conclusion of the program.
Definition: Core.cpp:112
ApproxSphere::runProblem
MoFEMErrorCode runProblem()
[Run programme]
Definition: approx_sphere.cpp:280
OpBase
OpBaseImpl< PETSC, EdgeEleOp > OpBase
Definition: radiation.cpp:29
MoFEM::LogManager::createSink
static boost::shared_ptr< SinkType > createSink(boost::shared_ptr< std::ostream > stream_ptr, std::string comm_filter)
Create a sink object.
Definition: LogManager.cpp:298
OpRhs::iNtegrate
MoFEMErrorCode iNtegrate(EntitiesFieldData::EntData &row_data)
Class dedicated to integrate operator.
Definition: approx_sphere.cpp:75
MoFEM::Simple
Simple interface for fast problem set-up.
Definition: Simple.hpp:27
OpRhs::xPtr
boost::shared_ptr< MatrixDouble > xPtr
Definition: approx_sphere.cpp:130
FTensor::Tensor2< double, 3, 3 >
order
constexpr int order
Definition: dg_projection.cpp:18
MoFEM::DeprecatedCoreInterface
Deprecated interface functions.
Definition: DeprecatedCoreInterface.hpp:16
OpLhs::OpLhs
OpLhs(const std::string field_name, boost::shared_ptr< MatrixDouble > x_ptr, boost::shared_ptr< MatrixDouble > dot_x_ptr)
Definition: approx_sphere.cpp:136
MoFEM::Interface
DeprecatedCoreInterface Interface
Definition: Interface.hpp:1975
approx_order
static constexpr int approx_order
Definition: prism_polynomial_approximation.cpp:14
MoFEM::OpCalculateVectorFieldValuesFromPetscVecImpl
Approximate field values for given petsc vector.
Definition: UserDataOperators.hpp:595
MoFEM::EntitiesFieldData::EntData::getFTensor0N
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
Definition: EntitiesFieldData.hpp:1489
res_J_dx
auto res_J_dx
Definition: approx_sphere.cpp:45
j
FTensor::Index< 'j', 3 > j
Definition: approx_sphere.cpp:37
MoFEM::OpBaseImpl
Definition: FormsIntegrators.hpp:170
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:535
FM_DIM
constexpr int FM_DIM
Definition: approx_sphere.cpp:15
ApproxSphere::setOPs
MoFEMErrorCode setOPs()
[Set up problem]
Definition: approx_sphere.cpp:325
MoFEM::createDMVector
auto createDMVector(DM dm)
Get smart vector from DM.
Definition: DMMoFEM.hpp:1018
EntData
EntitiesFieldData::EntData EntData
Definition: approx_sphere.cpp:25
ApproxSphere::outputResults
MoFEMErrorCode outputResults()
[Solve]
Definition: approx_sphere.cpp:394
OpError::xPtr
boost::shared_ptr< MatrixDouble > xPtr
Definition: approx_sphere.cpp:257
main
int main(int argc, char *argv[])
[Postprocess results]
Definition: approx_sphere.cpp:473
MoFEM
implementation of Data Operators for Forces and Sources
Definition: Common.hpp:10
OpLhs
Definition: approx_sphere.cpp:134
a
constexpr double a
Definition: approx_sphere.cpp:30
i
FTensor::Index< 'i', 3 > i
Definition: approx_sphere.cpp:36
simple
void simple(double P1[], double P2[], double P3[], double c[], const int N)
Definition: acoustic.cpp:69
convert.type
type
Definition: convert.py:64
MoFEM::FormsIntegrators::Assembly
Assembly methods.
Definition: FormsIntegrators.hpp:302
OpPPMap
OpPostProcMapInMoab< SPACE_DIM, SPACE_DIM > OpPPMap
Definition: photon_diffusion.cpp:29
MoFEM::LogManager::getStrmWorld
static boost::shared_ptr< std::ostream > getStrmWorld()
Get the strm world object.
Definition: LogManager.cpp:344
MoFEM::DMRegister_MoFEM
PetscErrorCode DMRegister_MoFEM(const char sname[])
Register MoFEM problem.
Definition: DMMoFEM.cpp:47
MoFEM::FaceElementForcesAndSourcesCore
Face finite element.
Definition: FaceElementForcesAndSourcesCore.hpp:23
MoFEM::OpGetHONormalsOnFace
Calculate normals at Gauss points of triangle element.
Definition: HODataOperators.hpp:281
MoFEM::PipelineManager::setDomainRhsIntegrationRule
MoFEMErrorCode setDomainRhsIntegrationRule(RuleHookFun rule)
Definition: PipelineManager.hpp:530
MOFEM_LOG_TAG
#define MOFEM_LOG_TAG(channel, tag)
Tag channel.
Definition: LogManager.hpp:339
MoFEM::PipelineManager::getOpDomainLhsPipeline
boost::ptr_deque< UserDataOperator > & getOpDomainLhsPipeline()
Get the Op Domain Lhs Pipeline object.
Definition: PipelineManager.hpp:749
t_kd
constexpr auto t_kd
Definition: free_surface.cpp:137
lhs_J_dx2
auto lhs_J_dx2
Definition: approx_sphere.cpp:51
ApproxSphere::solveSystem
MoFEMErrorCode solveSystem()
[Push operators to pipeline]
Definition: approx_sphere.cpp:361
A
constexpr double A
Definition: approx_sphere.cpp:34
field_name
constexpr auto field_name
Definition: poisson_2d_homogeneous.cpp:13
FTensor::Index< 'i', 3 >
FTensor::Tensor1::normalize
Tensor1< T, Tensor_Dim > normalize()
Definition: Tensor1_value.hpp:77
OpLhs::iNtegrate
MoFEMErrorCode iNtegrate(EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
Integrate grad-grad operator.
Definition: approx_sphere.cpp:144
OpRhs
Definition: approx_sphere.cpp:68
MoFEM::PipelineManager::setDomainLhsIntegrationRule
MoFEMErrorCode setDomainLhsIntegrationRule(RuleHookFun rule)
Definition: PipelineManager.hpp:503
integration_rule
auto integration_rule
Definition: free_surface.cpp:185
ElementsAndOps
Definition: child_and_parent.cpp:18
DomainEleOp
ApproxSphere::mField
MoFEM::Interface & mField
Definition: approx_sphere.cpp:269
MoFEM::CoreTmp< 0 >::Initialize
static MoFEMErrorCode Initialize(int *argc, char ***args, const char file[], const char help[])
Initializes the MoFEM database PETSc, MOAB and MPI.
Definition: Core.cpp:72
MOFEM_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition: LogManager.hpp:308
CATCH_ERRORS
#define CATCH_ERRORS
Catch errors.
Definition: definitions.h:372
MoFEM::Core
CoreTmp< 0 > Core
Definition: Core.hpp:1094
OpError::doWork
MoFEMErrorCode doWork(int side, EntityType type, EntData &data)
Definition: approx_sphere.cpp:228
UserDataOperator
ForcesAndSourcesCore::UserDataOperator UserDataOperator
Definition: HookeElement.hpp:75
k
FTensor::Index< 'k', 3 > k
Definition: approx_sphere.cpp:38
AINSWORTH_LEGENDRE_BASE
@ AINSWORTH_LEGENDRE_BASE
Ainsworth Cole (Legendre) approx. base .
Definition: definitions.h:60
MoFEM::createVectorMPI
auto createVectorMPI(MPI_Comm comm, PetscInt n, PetscInt N)
Create MPI Vector.
Definition: PetscSmartObj.hpp:198
a2
constexpr double a2
Definition: approx_sphere.cpp:31
OpRhs::xDotPtr
boost::shared_ptr< MatrixDouble > xDotPtr
Definition: approx_sphere.cpp:131
ApproxSphere
Definition: approx_sphere.cpp:262
a4
constexpr double a4
Definition: approx_sphere.cpp:32
MoFEMFunctionBeginHot
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:440
DomainEle
ElementsAndOps< SPACE_DIM >::DomainEle DomainEle
Definition: child_and_parent.cpp:34
MoFEM::LogManager::setLog
static LoggerType & setLog(const std::string channel)
Set ans resset chanel logger.
Definition: LogManager.cpp:389
MoFEM::SmartPetscObj< Vec >
res_J
auto res_J
Definition: approx_sphere.cpp:40
MoFEM::DMoFEMLoopFiniteElements
PetscErrorCode DMoFEMLoopFiniteElements(DM dm, const char fe_name[], MoFEM::FEMethod *method, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
Executes FEMethod for finite elements in DM.
Definition: DMMoFEM.cpp:590
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
MoFEM::FieldBlas
Basic algebra on fields.
Definition: FieldBlas.hpp:21
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:346
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
OpLhs::xPtr
boost::shared_ptr< MatrixDouble > xPtr
Definition: approx_sphere.cpp:215
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