Tetra Namespace Reference

Classes
class	prm
class	Tet

Functions
Eigen::Matrix< double, NPI, 1 >	calc_alpha_eff (const double dt, const double alpha, Eigen::Ref< Eigen::Matrix< double, NPI, 1 >> uHeff)

Variables
const int	N = 4
const int	NPI = 5
const double	epsilon
constexpr double	A = 1. / 4.
constexpr double	B = 1. / 6.
constexpr double	C = 1. / 2.
constexpr double	D = -2. / 15.
constexpr double	E = 3. / 40.
constexpr double	u [NPI] = {A, B, B, B, C}
constexpr double	v [NPI] = {A, B, B, C, B}
constexpr double	w [NPI] = {A, B, C, B, B}
constexpr double	pds [NPI] = {D, E, E, E, E}
constexpr double	a [N][NPI]

Detailed Description

to grab altogether some constants for struct Tet

Function Documentation

calc_alpha_eff()

Eigen::Matrix< double, NPI, 1 > Tetra::calc_alpha_eff	(	const double	dt,
		const double	alpha,
		Eigen::Ref< Eigen::Matrix< double, NPI, 1 >>	uHeff
	)

to perform some second order corrections, an effective \( \alpha \) is computed here with a piecewise formula

Variable Documentation

A

constexpr double Tetra::A = 1. / 4.

constexpr

constant to build hat functions

a

constexpr double Tetra::a[N][NPI]

constexpr

Initial value:

= {{1. - u[0] - v[0] - w[0], 1. - u[1] - v[1] - w[1],
                               1. - u[2] - v[2] - w[2], 1. - u[3] - v[3] - w[3],
                               1. - u[4] - v[4] - w[4]},
                              {u[0], u[1], u[2], u[3], u[4]},
                              {v[0], v[1], v[2], v[3], v[4]},
                              {w[0], w[1], w[2], w[3], w[4]}}

constant matrix \( j:0..4 \) a[0][j] = 1.-u[j]-v[j]-w[j]; a[1][j] = u[j]; a[2][j] = v[j]; a[3][j] = w[j];

B

constexpr double Tetra::B = 1. / 6.

constexpr

constant to build hat functions

C

constexpr double Tetra::C = 1. / 2.

constexpr

constant to build hat functions

D

constexpr double Tetra::D = -2. / 15.

constexpr

constant to build hat functions

E

constexpr double Tetra::E = 3. / 40.

constexpr

constant to build hat functions

epsilon

const double Tetra::epsilon

Initial value:

=
        EPSILON

this constant is defined from a macro in config.h.in, it is used to check the validity of the tetrahedreon, a degeneracy test

N

const int Tetra::N = 4

number of sommits

NPI

const int Tetra::NPI = 5

number of weights

pds

constexpr double Tetra::pds[NPI] = {D, E, E, E, E}

constexpr

some constant weights to build hat functions

u

constexpr double Tetra::u[NPI] = {A, B, B, B, C}

constexpr

some constants to build hat functions

v

constexpr double Tetra::v[NPI] = {A, B, B, C, B}

constexpr

some constants to build hat functions

w

constexpr double Tetra::w[NPI] = {A, B, C, B, B}

constexpr

some constants to build hat functions