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	DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle
using	BoundaryEle
using	FaceSideEle
using	FaceSideOp = FaceSideEle::UserDataOperator
using	PostProcEle = PostProcBrokenMeshInMoab<DomainEle>
using	OpDomainGradGrad
using	OpDomainSource

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

Initial value:

 
    PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::BoundaryEle

Definition at line 20 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.

FaceSideEle

using FaceSideEle

Initial value:

 
    PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::FaceSideEle

Definition at line 23 of file poisson_2d_dis_galerkin.cpp.

FaceSideOp

using FaceSideOp = FaceSideEle::UserDataOperator

Examples

PoissonDiscontinousGalerkin.hpp.

Definition at line 25 of file poisson_2d_dis_galerkin.cpp.

OpDomainGradGrad

using OpDomainGradGrad

Initial value:

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

Examples

poisson_2d_dis_galerkin.cpp.

Definition at line 36 of file poisson_2d_dis_galerkin.cpp.

OpDomainSource

using OpDomainSource

Initial value:

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

Examples

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 538 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

int BASE_DIM = 1

constexpr

Definition at line 11 of file poisson_2d_dis_galerkin.cpp.

FIELD_DIM

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

poisson_2d_dis_galerkin.cpp.

Definition at line 41 of file poisson_2d_dis_galerkin.cpp.

nitsche

double nitsche = 1

static

Examples

PoissonDiscontinousGalerkin.hpp, plate.cpp, and poisson_2d_dis_galerkin.cpp.

Definition at line 32 of file poisson_2d_dis_galerkin.cpp.

penalty

double penalty = 1e6

static

Examples

PoissonDiscontinousGalerkin.hpp, plate.cpp, poisson_2d_dis_galerkin.cpp, and shallow_wave.cpp.

Definition at line 29 of file poisson_2d_dis_galerkin.cpp.

phi

double phi

static

Initial value:

=
    -1

Examples

ADOLCPlasticityMaterialModels.hpp, EshelbianPlasticity.cpp, NavierStokesElement.cpp, PoissonDiscontinousGalerkin.hpp, plate.cpp, poisson_2d_dis_galerkin.cpp, and shallow_wave.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

mixed_poisson.cpp, plate.cpp, poisson_2d_dis_galerkin.cpp, and test_broken_space.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

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

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

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)
 
    };
};