poisson_2d_dis_galerkin.cpp File Reference

#include <MoFEM.hpp>
#include <PoissonDiscontinousGalerkin.hpp>

Go to the source code of this file.

Classes
struct	Poisson2DiscontGalerkin

Typedefs
using	EntData = EntitiesFieldData::EntData
using	DomainEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::DomainEle
using	DomainEleOp = DomainEle::UserDataOperator
using	BoundaryEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::BoundaryEle
using	BoundaryEleOp = BoundaryEle::UserDataOperator
using	FaceSideEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::FaceSideEle
using	FaceSideOp = FaceSideEle::UserDataOperator
using	PostProcEle = PostProcBrokenMeshInMoab< DomainEle >
using	OpDomainGradGrad = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpGradGrad< BASE_DIM, FIELD_DIM, SPACE_DIM >
using	OpDomainSource = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpSource< BASE_DIM, FIELD_DIM >

Functions
int	main (int argc, char *argv[])
	[Run program]

Variables
constexpr int	BASE_DIM = 1
constexpr int	FIELD_DIM = 1
constexpr int	SPACE_DIM = 2
static double	penalty = 1e6
static double	phi
static double	nitsche = 1
PetscBool	is_test = PETSC_FALSE
auto	u_exact
auto	u_grad_exact
auto	source
static char	help [] = "...\n\n"

Typedef Documentation

BoundaryEle

using BoundaryEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::BoundaryEle

Definition at line 20 of file poisson_2d_dis_galerkin.cpp.

BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 22 of file poisson_2d_dis_galerkin.cpp.

DomainEle

using DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle

Definition at line 17 of file poisson_2d_dis_galerkin.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 18 of file poisson_2d_dis_galerkin.cpp.

EntData

using EntData = EntitiesFieldData::EntData

Definition at line 15 of file poisson_2d_dis_galerkin.cpp.

FaceSideEle

using FaceSideEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::FaceSideEle

Definition at line 23 of file poisson_2d_dis_galerkin.cpp.

FaceSideOp

using FaceSideOp = FaceSideEle::UserDataOperator

Examples

mofem/tutorials/scl-11/src/PoissonDiscontinousGalerkin.hpp.

Definition at line 25 of file poisson_2d_dis_galerkin.cpp.

OpDomainGradGrad

using OpDomainGradGrad = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::BiLinearForm<GAUSS>::OpGradGrad<BASE_DIM, FIELD_DIM, SPACE_DIM>

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 36 of file poisson_2d_dis_galerkin.cpp.

OpDomainSource

using OpDomainSource = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<GAUSS>::OpSource<BASE_DIM, FIELD_DIM>

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 38 of file poisson_2d_dis_galerkin.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 27 of file poisson_2d_dis_galerkin.cpp.

Function Documentation

main()

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

[Run program]

[Main]

Definition at line 536 of file poisson_2d_dis_galerkin.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  // Error handling
  try {
    // Register MoFEM discrete manager in PETSc
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
 
    // Create MOAB instance
    moab::Core mb_instance;              // mesh database
    moab::Interface &moab = mb_instance; // mesh database interface
 
    // Create MoFEM instance
    MoFEM::Core core(moab);           // finite element database
    MoFEM::Interface &m_field = core; // finite element interface
 
    // Run the main analysis
    Poisson2DiscontGalerkin poisson_problem(m_field);
    CHKERR poisson_problem.runProgram();
  }
  CATCH_ERRORS;
 
  // Finish work: cleaning memory, getting statistics, etc.
  MoFEM::Core::Finalize();
 
  return 0;
}

Variable Documentation

BASE_DIM

constexpr int BASE_DIM = 1

constexpr

Definition at line 11 of file poisson_2d_dis_galerkin.cpp.

FIELD_DIM

constexpr int FIELD_DIM = 1

constexpr

Definition at line 12 of file poisson_2d_dis_galerkin.cpp.

help

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

static

Definition at line 72 of file poisson_2d_dis_galerkin.cpp.

is_test

PetscBool is_test = PETSC_FALSE

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 41 of file poisson_2d_dis_galerkin.cpp.

nitsche

double nitsche = 1

static

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-11/src/PoissonDiscontinousGalerkin.hpp, and mofem/tutorials/vec-6/plate.cpp.

Definition at line 32 of file poisson_2d_dis_galerkin.cpp.

penalty

double penalty = 1e6

static

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-11/src/PoissonDiscontinousGalerkin.hpp, mofem/tutorials/vec-4/shallow_wave.cpp, and mofem/tutorials/vec-6/plate.cpp.

Definition at line 29 of file poisson_2d_dis_galerkin.cpp.

phi

double phi

static

Initial value:

=
    -1

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-11/src/PoissonDiscontinousGalerkin.hpp, mofem/tutorials/vec-4/shallow_wave.cpp, mofem/tutorials/vec-6/plate.cpp, mofem/users_modules/adolc-plasticity/src/ADOLCPlasticityMaterialModels.hpp, and mofem/users_modules/basic_finite_elements/src/impl/NavierStokesElement.cpp.

Definition at line 30 of file poisson_2d_dis_galerkin.cpp.

source

auto source

Initial value:

= [](const double x, const double y, const double) {
  if (is_test)
    return -(2 * x * x + 2 * y * y);
  else
    return 8 * M_PI * M_PI * cos(2 * x * M_PI) * cos(2 * y * M_PI);
}

Examples

mofem/atom_tests/test_broken_space.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, and mofem/tutorials/vec-6/plate.cpp.

Definition at line 62 of file poisson_2d_dis_galerkin.cpp.

                                                               {
  if (is_test)
    return -(2 * x * x + 2 * y * y);
  else
    return 8 * M_PI * M_PI * cos(2 * x * M_PI) * cos(2 * y * M_PI);
};

SPACE_DIM

constexpr int SPACE_DIM = 2

constexpr

Definition at line 13 of file poisson_2d_dis_galerkin.cpp.

u_exact

auto u_exact

Initial value:

= [](const double x, const double y, const double) {
  if (is_test)
    return x * x * y * y;
  else
    return cos(2 * x * M_PI) * cos(2 * y * M_PI);
}

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 43 of file poisson_2d_dis_galerkin.cpp.

                                                                {
  if (is_test)
    return x * x * y * y;
  else
    return cos(2 * x * M_PI) * cos(2 * y * M_PI);
};

u_grad_exact

auto u_grad_exact

Initial value:

= [](const double x, const double y, const double) {
  if (is_test)
    return FTensor::Tensor1<double, 2>{2 * x * y * y, 2 * x * x * y};
  else
    return FTensor::Tensor1<double, 2>{
 
        -2 * M_PI * cos(2 * M_PI * y) * sin(2 * M_PI * x),
        -2 * M_PI * cos(2 * M_PI * x) * sin(2 * M_PI * y)
 
    };
}

Examples

mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 50 of file poisson_2d_dis_galerkin.cpp.

                                                                     {
  if (is_test)
    return FTensor::Tensor1<double, 2>{2 * x * y * y, 2 * x * x * y};
  else
    return FTensor::Tensor1<double, 2>{
 
        -2 * M_PI * cos(2 * M_PI * y) * sin(2 * M_PI * x),
        -2 * M_PI * cos(2 * M_PI * x) * sin(2 * M_PI * y)
 
    };
};