approx_sphere.cpp File Reference

#include <MoFEM.hpp>

Go to the source code of this file.

Classes
struct	ElementsAndOps< DIM >
struct	ElementsAndOps< 2 >
struct	OpRhs
struct	OpLhs
struct	OpError
struct	ApproxSphere

Typedefs
using	DomainEle = ElementsAndOps< FM_DIM >::DomainEle
using	DomainEleOp = DomainEle::UserDataOperator
using	EntData = EntitiesFieldData::EntData
using	AssemblyDomainEleOp = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::OpBase

Functions
int	main (int argc, char *argv[])
	[Postprocess results]

Variables
static char	help [] = "...\n\n"
constexpr int	FM_DIM = 2
constexpr double	a = 1
constexpr double	a2 = a * a
constexpr double	a4 = a2 * a2
constexpr double	A = 6371220
FTensor::Index< 'i', 3 >	i
FTensor::Index< 'j', 3 >	j
FTensor::Index< 'k', 3 >	k
auto	res_J
auto	res_J_dx
auto	lhs_J_dx2

Typedef Documentation

AssemblyDomainEleOp

using AssemblyDomainEleOp = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::OpBase

Definition at line 27 of file approx_sphere.cpp.

DomainEle

using DomainEle = ElementsAndOps<FM_DIM>::DomainEle

Definition at line 23 of file approx_sphere.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 24 of file approx_sphere.cpp.

EntData

using EntData = EntitiesFieldData::EntData

Definition at line 25 of file approx_sphere.cpp.

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

[Postprocess results]

[Register MoFEM discrete manager in PETSc]

[Register MoFEM discrete manager in PETSc

[Create MoAB]

< mesh database

< mesh database interface

[Create MoAB]

[Create MoFEM]

< finite element database

< finite element database insterface

[Create MoFEM]

[ApproxSphere]

[ApproxSphere]

Definition at line 473 of file approx_sphere.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  auto core_log = logging::core::get();
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "EXAMPLE"));
  LogManager::setLog("EXAMPLE");
  MOFEM_LOG_TAG("EXAMPLE", "example");
 
  try {
 
    //! [Register MoFEM discrete manager in PETSc]
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
    //! [Register MoFEM discrete manager in PETSc
 
    //! [Create MoAB]
    moab::Core mb_instance;              ///< mesh database
    moab::Interface &moab = mb_instance; ///< mesh database interface
    //! [Create MoAB]
 
    //! [Create MoFEM]
    MoFEM::Core core(moab);           ///< finite element database
    MoFEM::Interface &m_field = core; ///< finite element database insterface
    //! [Create MoFEM]
 
    //! [ApproxSphere]
    ApproxSphere ex(m_field);
    CHKERR ex.runProblem();
    //! [ApproxSphere]
  }
  CATCH_ERRORS;
 
  MoFEM::Core::Finalize();
}

Variable Documentation

a

constexpr double a = 1

constexpr

Examples

mofem/atom_tests/matrix_function.cpp, mofem/atom_tests/prism_polynomial_approximation.cpp, mofem/atom_tests/quad_polynomial_approximation.cpp, mofem/tutorials/adv-0/plastic.cpp, mofem/tutorials/adv-0/src/PlasticOpsGeneric.hpp, mofem/tutorials/adv-5/src/SeepageOps.hpp, mofem/tutorials/cor-0to1/src/UnsaturatedFlow.hpp, mofem/tutorials/cor-2to5/analytical_poisson_field_split.cpp, mofem/tutorials/cor-2to5/src/PoissonOperators.hpp, mofem/tutorials/cor-6/simple_elasticity.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/scl-1/src/poisson_2d_homogeneous.hpp, mofem/tutorials/scl-6/heat_equation.cpp, mofem/tutorials/scl-9/heat_method.cpp, mofem/tutorials/vec-2/src/HenckyOps.hpp, mofem/tutorials/vec-4/approx_sphere.cpp, mofem/tutorials/vec-4/shallow_wave.cpp, mofem/tutorials/vec-5/free_surface.cpp, mofem/tutorials/vec-6/plate.cpp, mofem/tutorials/vec-7/adjoint.cpp, mofem/users_modules/adolc-plasticity/adolc_plasticity.cpp, mofem/users_modules/basic_finite_elements/src/HookeElement.hpp, mofem/users_modules/basic_finite_elements/src/impl/HookeElement.cpp, mofem/users_modules/bone_remodelling/src/impl/Remodeling.cpp, mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp, and mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp.

Definition at line 30 of file approx_sphere.cpp.

A

constexpr double A = 6371220

constexpr

Examples

mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 34 of file approx_sphere.cpp.

a2

constexpr double a2 = a * a

constexpr

Definition at line 31 of file approx_sphere.cpp.

a4

constexpr double a4 = a2 * a2

constexpr

Definition at line 32 of file approx_sphere.cpp.

FM_DIM

constexpr int FM_DIM = 2

constexpr

Examples

mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 15 of file approx_sphere.cpp.

help

char help[] = "...\n\n"

static

Definition at line 11 of file approx_sphere.cpp.

i

FTensor::Index<'i', 3> i

Examples

mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 36 of file approx_sphere.cpp.

j

FTensor::Index<'j', 3> j

Examples

mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 37 of file approx_sphere.cpp.

k

FTensor::Index<'k', 3> k

Examples

mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 38 of file approx_sphere.cpp.

lhs_J_dx2

auto lhs_J_dx2

Initial value:

= [](const double x, const double y, const double z) {
  const double res = res_J(x, y, z);
  return FTensor::Tensor2<double, 3, 3>{
 
      (res * 2 + (4 * x * x)),
      (4 * y * x),
      (4 * z * x),
 
      (4 * x * y),
      (2 * res + (4 * y * y)),
      (4 * z * y),
 
      (4 * x * z),
      (4 * y * z),
      (2 * res + (4 * z * z))};
}

Examples

mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 51 of file approx_sphere.cpp.

                                                                    {
  const double res = res_J(x, y, z);
  return FTensor::Tensor2<double, 3, 3>{
 
      (res * 2 + (4 * x * x)),
      (4 * y * x),
      (4 * z * x),
 
      (4 * x * y),
      (2 * res + (4 * y * y)),
      (4 * z * y),
 
      (4 * x * z),
      (4 * y * z),
      (2 * res + (4 * z * z))};
};

res_J

auto res_J

Initial value:

= [](const double x, const double y, const double z) {
  const double res = (x * x + y * y + z * z - a2);
  return res;
}

Examples

mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 40 of file approx_sphere.cpp.

                                                                {
  const double res = (x * x + y * y + z * z - a2);
  return res;
};

res_J_dx

auto res_J_dx

Initial value:

= [](const double x, const double y, const double z) {
  const double res = res_J(x, y, z);
  return FTensor::Tensor1<double, 3>{res * (2 * x), res * (2 * y),
                                     res * (2 * z)};
}

Examples

mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 45 of file approx_sphere.cpp.

                                                                   {
  const double res = res_J(x, y, z);
  return FTensor::Tensor1<double, 3>{res * (2 * x), res * (2 * y),
                                     res * (2 * z)};
};