#include <stdlib.h>
#include <cmath>
#include <MoFEM.hpp>
#include <boost/math/quadrature/gauss_kronrod.hpp>

Classes
struct	ElementsAndOps< DIM >

struct	OpError

struct	PhotonDiffusion

Macros
#define	BOOST_MATH_GAUSS_NO_COMPUTE_ON_DEMAND

Typedefs
using	DomainEle = VolumeElementForcesAndSourcesCore
	[Define dimension]

using	DomainEleOp = DomainEle::UserDataOperator

using	PostProcEle = PostProcBrokenMeshInMoab< DomainEle >

using	OpDomainMass = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMass< 1, 1 >

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

using	OpDomainTimesScalarField = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpBaseTimesScalar< 1 >

using	OpDomainGradTimesVec = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpGradTimesTensor< 1, 1, SPACE_DIM >

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

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

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

constexpr int	SPACE_DIM = 3
	[Define dimension]

const double	n = 1.44
	refractive index of diffusive medium

const double	c = 30.
	speed of light (cm/ns)

const double	v = c / n
	phase velocity of light in medium (cm/ns)

double	mu_a
	absorption coefficient (cm^-1)

double	mu_sp
	scattering coefficient (cm^-1)

double	D

double	slab_thickness

double	beam_radius

double	beam_centre_x

double	beam_centre_y

double	flux_magnitude = 1e3
	impulse magnitude

double	initial_time

char	out_file_name [255] = "init_file.dat"

int	numHoLevels = 1

PetscBool	output_volume = PETSC_FALSE

Macro Definition Documentation

◆ BOOST_MATH_GAUSS_NO_COMPUTE_ON_DEMAND

#define BOOST_MATH_GAUSS_NO_COMPUTE_ON_DEMAND

Definition at line 11 of file initial_diffusion.cpp.

Typedef Documentation

◆ DomainEle

using DomainEle = VolumeElementForcesAndSourcesCore

[Define dimension]

Definition at line 23 of file initial_diffusion.cpp.

◆ DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 24 of file initial_diffusion.cpp.

◆ OpDomainGradGrad

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

Definition at line 29 of file initial_diffusion.cpp.

◆ OpDomainGradTimesVec

using OpDomainGradTimesVec = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<GAUSS>::OpGradTimesTensor<1, 1, SPACE_DIM>

Examples: mofem/tutorials/scl-8/radiation.cpp.

Definition at line 33 of file initial_diffusion.cpp.

◆ OpDomainMass

using OpDomainMass = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::BiLinearForm<GAUSS>::OpMass<1, 1>

Definition at line 27 of file initial_diffusion.cpp.

◆ OpDomainSource

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

Definition at line 35 of file initial_diffusion.cpp.

◆ OpDomainTimesScalarField

using OpDomainTimesScalarField = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<GAUSS>::OpBaseTimesScalar<1>

Examples: mofem/atom_tests/child_and_parent.cpp.

Definition at line 31 of file initial_diffusion.cpp.

◆ PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 25 of file initial_diffusion.cpp.

Function Documentation

◆ main()

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

Definition at line 399 of file initial_diffusion.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  // Add logging channel for example
  auto core_log = logging::core::get();
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "INITIAL"));
  LogManager::setLog("INITIAL");
  MOFEM_LOG_TAG("INITIAL", "initial_diffusion")
 
  // 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
    PhotonDiffusion heat_problem(m_field);
    CHKERR heat_problem.runProgram();
  }
  CATCH_ERRORS;
 
  // Finish work: cleaning memory, getting statistics, etc.
  MoFEM::Core::Finalize();
 
  return 0;
}