#include <Op_Div_DG.h>

Inheritance diagram for Op_Div_DG:

Collaboration diagram for Op_Div_DG:

Public Member Functions
void	associer (const Domaine_dis_base &, const Domaine_Cl_dis_base &, const Champ_Inc_base &) override
void	completer () override
	Associates the operator with the domaine_dis, the domaine_Cl_dis, and the unknown of its equation.
DoubleTab &	calculer (const DoubleTab &, DoubleTab &) const override
int	impr (Sortie &os) const override
	DOES NOTHING - to override in derived classes.
void	volumique (DoubleTab &) const override
	Converts the divergence field from an integrated to a volumetric (per-unit-volume) form.
int	has_interface_blocs () const override
void	dimensionner_blocs (matrices_t matrices, const tabs_t &semi_impl={}) const override
	Sizes the velocity-pressure and pressure-pressure matrix blocks.
void	ajouter_blocs_ext (const DoubleTab &vit, matrices_t matrices, DoubleTab &secmem, const tabs_t &semi_impl={ }) const override
	Assembles the DG divergence operator and, if needed, the pressure stabilization term.
Public Member Functions inherited from Operateur_Div_base
DoubleVect &	multvect (const DoubleTab &, DoubleTab &) const
void	ajouter_blocs (matrices_t matrices, DoubleTab &secmem, const tabs_t &semi_impl={}) const override
DoubleTab &	ajouter (const DoubleTab &vit, DoubleTab &div) const override
Public Member Functions inherited from Operateur_base
virtual void	associer_champ (const Champ_Inc_base &, const std::string &nom_ch)
virtual void	associer_domaine_cl_dis (const Domaine_Cl_dis_base &)
virtual void	dimensionner (Matrice_Morse &) const
	DOES NOTHING - to override in derived classes.
virtual void	dimensionner_bloc_vitesse (Matrice_Morse &matrice) const
virtual void	modifier_pour_Cl (Matrice_Morse &, DoubleTab &) const
	DOES NOTHING - to override in derived classes.
virtual void	contribuer_a_avec (const DoubleTab &, Matrice_Morse &) const
	DOES NOTHING - to override in derived classes.
virtual void	contribuer_bloc_vitesse (const DoubleTab &, Matrice_Morse &) const
virtual void	contribuer_au_second_membre (DoubleTab &) const
	DOES NOTHING - to override in derived classes.
void	tester_contribuer_a_avec (const DoubleTab &, const Matrice_Morse &)
virtual void	dimensionner_termes_croises (Matrice_Morse &, const Probleme_base &autre_pb, int nl, int nc) const
virtual void	ajouter_termes_croises (const DoubleTab &inco, const Probleme_base &autre_pb, const DoubleTab &autre_inco, DoubleTab &resu) const
virtual void	contribuer_termes_croises (const DoubleTab &inco, const Probleme_base &autre_pb, const DoubleTab &autre_inco, Matrice_Morse &matrice) const
virtual double	calculer_dt_stab () const
	Computes dt_stab.
virtual void	calculer_dt_local (DoubleTab &) const
virtual void	mettre_a_jour (double temps)
	DOES NOTHING - to override in derived classes.
virtual void	abortTimeStep ()
virtual void	resetTime (double time)
void	associer_eqn (const Equation_base &)
	Associates an equation with the operator.
int	get_decal_temps () const
int	set_decal_temps (int)
int	get_nb_ss_pas_de_temps () const
int	set_nb_ss_pas_de_temps (int)
const Matrice &	get_matrice () const
Matrice &	set_matrice ()
const SolveurSys &	get_solveur () const
SolveurSys &	set_solveur ()
Entree &	lire_solveur (Entree &)
virtual int	systeme_invariant () const
virtual void	ajouter_contribution_explicite_au_second_membre (const Champ_Inc_base &inconnue, DoubleTab &derivee) const
const Champ_Inc_base &	mon_inconnue () const
bool	has_champ_inco () const
const std::string &	nom_inconnue () const
void	ouvrir_fichier (SFichier &os, const Nom &, const int flag=1) const
	Opening/creation of a file for printing an operator. To override in derived classes.
void	ouvrir_fichier_partage (EcrFicPartage &, const Nom &, const int flag=1) const
	Opening/creation of a shared file for printing an operator. To override in derived classes.
void	set_fichier (const Nom &)
const Nom	fichier () const
void	set_description (const Nom &nom)
const Nom	description () const
DoubleTab &	flux_bords ()
DoubleTab &	flux_bords () const
void	creer_champ (const Motcle &motlu) override
const Champ_base &	get_champ (const Motcle &nom) const override
bool	has_champ (const Motcle &nom, OBS_PTR(Champ_base) &ref_champ) const override
bool	has_champ (const Motcle &nom) const override
void	get_noms_champs_postraitables (Noms &nom, Option opt=NONE) const override
void	calculer_pour_post (Champ_base &espace_stockage, const Nom &option, int comp) const override
Motcle	get_localisation_pour_post (const Nom &option) const override
virtual void	ajouter_flux (const DoubleTab &inconnue, DoubleTab &contribution) const
virtual void	calculer_flux (const DoubleTab &inconnue, DoubleTab &flux) const
virtual void	preparer_calcul ()
bool	has_impr_file () const
Public Member Functions inherited from Objet_U
	~Objet_U () override
	Destructor. Removes the object from the list of objects registered in "memory".
int	numero () const
	Returns the index of the object in Memoire::data.
virtual int	duplique () const =0
virtual Sortie &	printOn (Sortie &) const
	Writes the object to an output stream. Virtual method to override.
virtual Entree &	readOn (Entree &)
	Reads an Objet_U from an input stream. Virtual method to override.
virtual unsigned	taille_memoire () const =0
virtual int	est_egal_a (const Objet_U &) const
	Returns 1 if x and *this are the same instance (same memory address).
virtual const Nom &	le_nom () const
	Returns the name of the Objet_U. Virtual method to override: returns "neant" in this implementation.
virtual void	nommer (const Nom &)
	Assigns a name to the Objet_U. Virtual method to override.
virtual int	reprendre (Entree &)
	Restores an Objet_U from an input stream. Virtual method to override.
virtual int	sauvegarder (Sortie &) const
	Saves an Objet_U to an output stream. Virtual method to override.
int	get_object_id () const
	Returns the unique identifier of the object (object_id_).
virtual const Type_info *	get_info () const
	Returns type information for the Objet_U.
const Nom &	que_suis_je () const
	Returns the string identifying the class.
const char *	le_type () const
	Returns the type name of the Objet_U.
virtual int	change_num (const int *const)
	Changes the internal number of the Objet_U.
virtual int	lire_motcle_non_standard (const Motcle &motlu, Entree &is)
	Reads non-simple-type parameters of an Objet_U from an input stream.
virtual int	associer_ (Objet_U &)
	Associates the Objet_U with another Objet_U. Virtual method to override.
const Interprete &	interprete () const
Interprete &	interprete ()
Public Member Functions inherited from Process
virtual	~Process ()
Public Member Functions inherited from MorEqn
void	associer_eqn (const Equation_base &)
	Associates an equation with the object.
const Equation_base &	equation () const
	Returns the reference to the equation pointed to by MorEqn::mon_equation.
Equation_base &	equation ()
	Returns the reference to the equation pointed to by MorEqn::mon_equation.
int	mon_equation_non_nul () const
virtual void	check_multiphase_compatibility () const
Public Member Functions inherited from Champs_compris_interface
virtual	~Champs_compris_interface ()
virtual bool	has_champ (const Motcle &nom, OBS_PTR(Champ_base)&ref_champ) const =0

Protected Member Functions
	OBS_PTR (Domaine_DG) le_dom_DG
	OBS_PTR (Domaine_Cl_DG) le_dcl_DG
	OBS_PTR (Op_Diff_DG_base) op_diff_
Protected Member Functions inherited from Operateur_base
	OBS_PTR (Champ_Inc_base) le_champ_inco
Protected Member Functions inherited from Objet_U
	Objet_U ()
	Default constructor: assigns a unique identifier to the object (object_id_) and registers the object in "memory" by giving it a num_obj number.
	Objet_U (const Objet_U &)
	Copy constructor.
const Objet_U &	operator= (const Objet_U &)
	Assignment operator: does nothing (the number and identifier are preserved).
virtual void	set_param (Param &) const
Protected Member Functions inherited from MorEqn
	OBS_PTR (Equation_base) mon_equation
virtual	~MorEqn ()

Additional Inherited Members
Static Public Member Functions inherited from Objet_U
static const Nom &	nom_du_cas ()
	Returns a constant reference to the case name. This method is static.
static Nom &	get_set_nom_du_cas ()
	Returns a non-constant reference to the case name (to allow modification). This method is static.
static const Type_info *	info ()
	Returns type information for the Objet_U (static version).
static const Objet_U &	self_cast (const Objet_U &)
	Method added for casting in Python.
static Objet_U &	self_cast (Objet_U &)
Static Public Member Functions inherited from Process
static int	me ()
	Returns the rank of the local processor in the current communication group. See Comm_Group::rank() and PE_Groups::current_group().
static int	nproc ()
	Returns the number of processors in the current group. See Comm_Group::nproc() and PE_Groups::current_group().
static bool	is_parallel ()
static void	exit (int exit_code=-1)
	Exit routine for TRUST within a Kokkos region.
static double	mp_sum (double)
	Computes the sum of x over all processors in the current group.
static float	mp_sum (float)
static trustIdType	mp_sum (trustIdType)
	Computes the sum of x over all processors in the current group.
static double	mp_max (double)
static double	mp_min (double)
static int	mp_max (int)
	Returns the maximum value of x across all processors in the current group.
static int	mp_min (int)
	Returns the minimum value of x across all processors in the current group.
static double	mp_sum_as_double (int v)
static trustIdType	mppartial_sum (trustIdType i)
	Computes the partial sum of x over processors 0 to me()-1 (returns 0 on processor 0).
template<typename T>
static void	mp_sum_for_each (T &arg1, T &arg2)
	C++14 compatible mp_sum_for_each: combine multiple mp_sum calls into one collective operation Usage: mp_sum_for_each(a, b); mp_sum_for_each(a, b, c); mp_sum_for_each(a, b, c, d); mp_sum_for_each(a, b, c, d, e); All arguments must be of the same type (double or int) and are modified in place. Supports 2-5 parameters.
template<typename T>
static void	mp_sum_for_each (T &arg1, T &arg2, T &arg3)
template<typename T>
static void	mp_sum_for_each (T &arg1, T &arg2, T &arg3, T &arg4)
template<typename T>
static void	mp_sum_for_each (T &arg1, T &arg2, T &arg3, T &arg4, T &arg5)
template<typename T>
static void	mp_max_for_each (T &arg1, T &arg2)
	C++14 compatible mp_max_for_each: combine multiple mp_max calls into one collective operation.
template<typename T>
static void	mp_max_for_each (T &arg1, T &arg2, T &arg3)
template<typename T>
static void	mp_max_for_each (T &arg1, T &arg2, T &arg3, T &arg4)
template<typename T>
static void	mp_max_for_each (T &arg1, T &arg2, T &arg3, T &arg4, T &arg5)
template<typename T>
static void	mp_min_for_each (T &arg1, T &arg2)
	C++14 compatible mp_min_for_each: combine multiple mp_min calls into one collective operation.
template<typename T>
static void	mp_min_for_each (T &arg1, T &arg2, T &arg3)
template<typename T>
static void	mp_min_for_each (T &arg1, T &arg2, T &arg3, T &arg4)
template<typename T>
static void	mp_min_for_each (T &arg1, T &arg2, T &arg3, T &arg4, T &arg5)
template<typename _TYPE_>
static void	mp_sum_for_each_item (TRUSTArray< _TYPE_ > &x, int n=-1)
template<typename _TYPE_>
static void	mp_max_for_each_item (TRUSTArray< _TYPE_ > &x, int n=-1)
template<typename _TYPE_>
static void	mp_min_for_each_item (TRUSTArray< _TYPE_ > &x, int n=-1)
static bool	mp_and (bool)
	Computes the logical AND of b across all processors in the current group.
static bool	mp_or (bool)
static int	check_int_overflow (trustIdType)
static int	je_suis_maitre ()
	Returns 1 if on the master processor of the current group (i.e. me() == 0), 0 otherwise.
static KOKKOS_INLINE_FUNCTION void	Kokkos_exit (const char *)
	Exit routine for TRUST within a Kokkos region.
static int	node_master ()
	Returns 1 if on the NUMA node master processor, 0 otherwise.
static void	exit (const Nom &message, int exit_code=-1)
static bool	is_sequential ()
static void	barrier ()
	Synchronizes all processors in the current group (waits until all processors have reached the barrier).
static void	abort ()
	Abort routine for TRUST on a fatal error.
static Sortie &	Journal (int message_level=0)
	Returns a static Sortie object used as an event journal.
static double	ram_processeur ()
static void	imprimer_ram_totale (int all_process=0)
static bool	force_single_file (const int ranks, const Nom &filename)
Public Attributes inherited from Operateur_base
int	col_width_
Static Public Attributes inherited from Objet_U
static double	precision_geom = 1e-10
static constexpr bool	HAS_POINTER = false
static int	dimension =0
static int	format_precision_geom =11
static int	axi =0
static int	bidim_axi =0
static int	DEACTIVATE_SIGINT_CATCH =0
static Type_info	info_obj
static bool	disable_TU =false
	Flag to disable or not the writing of the .TU files.
static bool	stat_per_proc_perf_log =false
	Flag to enable the writing of the statistics detailed per processor in _csv.TU file.
Static Public Attributes inherited from Process
static int	exception_sur_exit =0
static int	multiple_files =5120
Protected Attributes inherited from Operateur_Div_base
SFichier	Flux_div
Protected Attributes inherited from Operateur_base
int	decal_temps
int	nb_ss_pas_de_temps
SolveurSys	solveur
Matrice	matrice_
Nom	out_
Nom	description_
DoubleTab	flux_bords_
Champs_compris	champs_compris_
std::string	nom_inco_

Detailed Description

class Op_Div_DG

This class represents the divergence operator with DG discretization. It computes the divergence of a velocity field (P1NC-type).

See also: Operateur_Div_base

Definition at line 34 of file Op_Div_DG.h.

Member Function Documentation

◆ ajouter_blocs_ext()

void Op_Div_DG::ajouter_blocs_ext	(	const DoubleTab &	vit,
		matrices_t	matrices,
		DoubleTab &	secmem,
		const tabs_t &	semi_impl = { } ) const

overridevirtual

Assembles the DG divergence operator and, if needed, the pressure stabilization term.

The assembly proceeds in two independent parts:

1. Divergence part (always assembled)

Uses the velocity basis gradient and the pressure basis:

Volume term: for each element, integral of q_h * div(u_h) = integral of grad(q_h) . u_h (integration by parts) i.e., integral of grad(phi_p_j) . phi_v_i over the element.
Internal face jump term: for each internal face shared by elem0 and elem1, integral of [u_h . n]_f * {{q_h}} = 0.5 * integral of (phi_v_i0 - phi_v_i1) . n * (phi_p_j0 + phi_p_j1) The four (elem0/elem1) x (elem0/elem1) combinations are assembled with the appropriate sign.
Boundary face term: for Dirichlet boundaries, the boundary velocity contributes to the face normal flux. The commented-out block handles time-varying Dirichlet data on the RHS (not yet active).

2. Pressure stabilization part (only when order_v == order_p and matp is allocated)

Adds a Dohrmann-Bochev-type face stabilization to recover inf-sup stability for equal-order velocity-pressure pairs: nu_f * integral of [p_h]_f * [q_h]_f where nu_f is the harmonic mean of the diffusivities of the two adjacent elements, and the jump [.]_f is assembled as the four cross-element combinations.

Parameters

vit	Velocity field values.
matrices	Map of matrix name → Matrice_Morse pointer to accumulate into.
secmem	Right-hand side (continuity residual) to accumulate into.
semi_impl	Map of semi-implicit field values.

Reimplemented from Operateur_Div_base.

Definition at line 256 of file Op_Div_DG.cpp.

◆ associer()

void Op_Div_DG::associer	(	const Domaine_dis_base &	domaine_dis,
		const Domaine_Cl_dis_base &	domaine_Cl_dis,
		const Champ_Inc_base &	)

overridevirtual

Implements Operateur_base.

Definition at line 40 of file Op_Div_DG.cpp.

◆ calculer()

DoubleTab & Op_Div_DG::calculer	(	const DoubleTab &	vit,
		DoubleTab &	div ) const

overridevirtual

Reimplemented from Operateur_base.

Definition at line 575 of file Op_Div_DG.cpp.

◆ completer()

void Op_Div_DG::completer ( )

overridevirtual

Associates the operator with the domaine_dis, the domaine_Cl_dis, and the unknown of its equation.

Exceptions

no	associated equation

Reimplemented from Operateur_base.

Definition at line 46 of file Op_Div_DG.cpp.

◆ dimensionner_blocs()

void Op_Div_DG::dimensionner_blocs	(	matrices_t	matrices,
		const tabs_t &	semi_impl = {} ) const

overridevirtual

Sizes the velocity-pressure and pressure-pressure matrix blocks.

Builds two sparsity patterns depending on what is needed:

Velocity-pressure block (matv, always built when the "vitesse" matrix is present): Each pressure DOF of element T is coupled to all velocity DOFs of T and its face-neighbours, as given by the pre-computed stencil. The block has size_p rows and size_v columns.
Pressure-pressure block (matp, only built when order_v == order_p): Required for the equal-order pressure stabilization term. Each pressure DOF of element T is coupled to all pressure DOFs of T and its face-neighbours. When order_v != order_p, a minimal dummy pattern (one non-zero) is allocated to satisfy the matrix infrastructure without adding real entries.

Parameters

matrices	Map of matrix name → Matrice_Morse pointer to be sized.
semi_impl	Map of semi-implicit field names; if "vitesse" is present, returns immediately.

Reimplemented from Operateur_base.

Definition at line 71 of file Op_Div_DG.cpp.

◆ has_interface_blocs()

int Op_Div_DG::has_interface_blocs ( ) const

inlineoverridevirtual

Reimplemented from Operateur_base.

Definition at line 45 of file Op_Div_DG.h.

◆ impr()

int Op_Div_DG::impr ( Sortie & os ) const

overridevirtual

DOES NOTHING - to override in derived classes.

Prints the operator to an output stream.

Parameters

(Sortie& os)

Returns: (int) always returns 1

Reimplemented from Operateur_base.

Definition at line 581 of file Op_Div_DG.cpp.

◆ OBS_PTR() [1/3]

Op_Div_DG::OBS_PTR ( Domaine_Cl_DG )

protected

◆ OBS_PTR() [2/3]

Op_Div_DG::OBS_PTR ( Domaine_DG )

protected

◆ OBS_PTR() [3/3]

Op_Div_DG::OBS_PTR ( Op_Diff_DG_base )

protected

◆ volumique()

void Op_Div_DG::volumique ( DoubleTab & div ) const

overridevirtual

Converts the divergence field from an integrated to a volumetric (per-unit-volume) form.

Divides each element's divergence value by the element volume.

Parameters

div	The divergence field to scale in-place.

Implements Operateur_Div_base.

Definition at line 593 of file Op_Div_DG.cpp.

The documentation for this class was generated from the following files:

/home/docs/trust-code/src/DG/Operateurs/Op_Divers/Op_Div_DG.h
/home/docs/trust-code/src/DG/Operateurs/Op_Divers/Op_Div_DG.cpp

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ ajouter_blocs_ext()

◆ associer()

◆ calculer()

◆ completer()

◆ dimensionner_blocs()

◆ has_interface_blocs()

◆ impr()

◆ OBS_PTR() [1/3]

◆ OBS_PTR() [2/3]

◆ OBS_PTR() [3/3]

◆ volumique()