#include <MoFEM.hpp>

Classes
struct	LevelSet

struct	LevelSet::SideData
	data structure carrying information on skeleton on both sides. More...

struct	LevelSet::WrapperClass
	Wrapper executing stages while mesh refinement. More...

struct	LevelSet::WrapperClassInitalSolution
	Used to execute inital mesh approximation while mesh refinement. More...

struct	LevelSet::WrapperClassErrorProjection
	Use peculated errors on all levels while mesh projection. More...

struct	LevelSet::OpRhsDomain

struct	LevelSet::OpLhsDomain

struct	LevelSet::OpRhsSkeleton

struct	LevelSet::OpLhsSkeleton

Typedefs
using	EntData = EntitiesFieldData::EntData

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

using	DomainParentEle = PipelineManager::ElementsAndOpsByDim< FE_DIM >::DomainParentEle

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

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

using	DomainEleOp = DomainEle::UserDataOperator

using	BoundaryEleOp = BoundaryEle::UserDataOperator

using	FaceSideEleOp = FaceSideEle::UserDataOperator

using	PostProcEle = PostProcBrokenMeshInMoab< DomainEle >

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

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

constexpr int	FE_DIM = EXECUTABLE_DIMENSION
	[Define dimension]

constexpr int	SPACE_DIM = FE_DIM

constexpr int	DIM1 = 1

constexpr int	DIM2 = 1

FTensor::Index< 'I', DIM1 >	I

FTensor::Index< 'J', DIM1 >	J

constexpr AssemblyType	A = AssemblyType::PETSC

constexpr IntegrationType	G

constexpr FieldSpace	potential_velocity_space = FE_DIM == 2 ? H1 : HCURL

constexpr size_t	potential_velocity_field_dim = FE_DIM == 2 ? 1 : 3

constexpr bool	debug = true

constexpr int	nb_levels = 3

constexpr int	start_bit

constexpr int	current_bit
	dofs bit used to do calculations

constexpr int	skeleton_bit = 2 * start_bit + 2
	skeleton elements bit

constexpr int	aggregate_bit
	all bits for advection problem

constexpr int	projection_bit

constexpr int	aggregate_projection_bit
	all bits for projection problem

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

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

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

LevelSet *	level_set_raw_ptr = nullptr

MoFEM::TsCtx *	ts_ctx

Typedef Documentation

◆ BoundaryEle

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

Definition at line 40 of file level_set.cpp.

◆ BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 44 of file level_set.cpp.

◆ DomainEle

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

Definition at line 37 of file level_set.cpp.

◆ DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Examples: mofem/tutorials/adv-3/level_set.cpp.

Definition at line 43 of file level_set.cpp.

◆ DomainParentEle

using DomainParentEle = PipelineManager::ElementsAndOpsByDim<FE_DIM>::DomainParentEle

Definition at line 38 of file level_set.cpp.

◆ EntData

using EntData = EntitiesFieldData::EntData

Definition at line 36 of file level_set.cpp.

◆ FaceSideEle

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

Definition at line 41 of file level_set.cpp.

◆ FaceSideEleOp

using FaceSideEleOp = FaceSideEle::UserDataOperator

Definition at line 45 of file level_set.cpp.

◆ PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 47 of file level_set.cpp.

Function Documentation

◆ main()

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

Definition at line 464 of file level_set.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  try {
 
    // Create MoAB database
    moab::Core moab_core;
    moab::Interface &moab = moab_core;
 
    // Create MoFEM database and link it to MoAB
    MoFEM::Core mofem_core(moab);
    MoFEM::Interface &m_field = mofem_core;
 
    // Register DM Manager
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
 
    // Add logging channel for example
    auto core_log = logging::core::get();
    core_log->add_sink(
        LogManager::createSink(LogManager::getStrmWorld(), "LevelSet"));
    LogManager::setLog("LevelSet");
    MOFEM_LOG_TAG("LevelSet", "LevelSet");
 
    LevelSet level_set(m_field);
    CHKERR level_set.runProblem();
  }
  CATCH_ERRORS;
 
  // finish work cleaning memory, getting statistics, etc.
  MoFEM::Core::Finalize();
 
  return 0;
}