DG gradient operator acting on a DG pressure field to produce a velocity-space residual. More...

#include <Op_Grad_DG.h>

Inheritance diagram for Op_Grad_DG:

Collaboration diagram for Op_Grad_DG:

Public Member Functions
void	associer (const Domaine_dis_base &, const Domaine_Cl_dis_base &, const Champ_Inc_base &) override
	Associates the operator with a DG domain and its boundary conditions.
int	has_interface_blocs () const override
void	dimensionner_blocs (matrices_t matrices, const tabs_t &semi_impl={}) const override
	Sizes the velocity-pressure matrix block for the gradient operator.
void	ajouter_blocs (matrices_t matrices, DoubleTab &secmem, const tabs_t &semi_impl={ }) const override
	Assembles the DG gradient operator into the matrix and right-hand side.
DoubleTab &	calculer (const DoubleTab &, DoubleTab &) const override
int	impr (Sortie &os) const override
	DOES NOTHING - to override in derived classes.
Public Member Functions inherited from Operateur_Grad_base
void	dimensionner (Matrice_Morse &) const override
	DOES NOTHING - to override in derived classes.
void	contribuer_a_avec (const DoubleTab &, Matrice_Morse &) const override
	DOES NOTHING - to override in derived classes.
DoubleTab &	ajouter (const DoubleTab &inco, DoubleTab &secmem) const override
virtual DoubleVect &	multvect (const DoubleTab &, DoubleTab &) const
	Computation without boundary conditions ?
virtual void	calculer_flux_bords () const
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_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_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	completer ()
	Associates the operator with the domaine_dis, the domaine_Cl_dis, and the unknown of its equation.
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
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_Grad_base
SFichier	Flux_grad
SFichier	Flux_grad_moment
SFichier	Flux_grad_sum
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

DG gradient operator acting on a DG pressure field to produce a velocity-space residual.

This class implements the DG discretization of the pressure gradient operator used in the incompressible Navier-Stokes momentum equation. It assembles the coupling between pressure DOFs and velocity DOFs through the interface_blocs mechanism.

The weak formulation integrates -grad(p_h) against velocity test functions v_h:

Volume term: integral of grad(p_h) . v_h (element-wise, by parts)
Internal face term: integral of {{p_h}} * [v_h . n] (average-jump coupling)

This formulation is the transpose of the divergence operator assembled by Op_Div_DG in the continuous sense. However, Op_Grad_DG and Op_Div_DG use independent implementations: Op_Grad_DG integrates grad(phi_p) . phi_v on element volumes and {{phi_p}} * [phi_v . n] on faces, while Op_Div_DG integrates phi_p * div(phi_v) and [phi_v . n] * {{phi_p}}. Both produce the same matrix up to a sign when the two basis sets are compatible, but they are kept separate to allow independent tuning.

See also: Op_Div_DG, Operateur_Grad_base

Definition at line 45 of file Op_Grad_DG.h.

Member Function Documentation

◆ ajouter_blocs()

void Op_Grad_DG::ajouter_blocs	(	matrices_t	matrices,
		DoubleTab &	secmem,
		const tabs_t &	semi_impl = { } ) const

overridevirtual

Assembles the DG gradient operator into the matrix and right-hand side.

The assembly has two stages:

1. Volume term — for each element: integral of grad(phi_p_j) . phi_v_i over the element, accumulated as (*mat)(v_dof, p_dof) += coeff and secmem(elem, v_dof) -= coeff * p(elem, p_dof).

2. Internal face term — for each internal face shared by elem0 and elem1: The average-jump coupling integral of {{phi_p}} * [phi_v . n] is expanded into four pointwise products at face quadrature points:

coeff00: -0.5 * phi_p0 * n_d * phi_v0 (elem0 pressure, elem0 velocity)
coeff01: +0.5 * phi_p1 * n_d * phi_v0 (elem1 pressure, elem0 velocity)
coeff10: -0.5 * phi_p0 * n_d * phi_v1 (elem0 pressure, elem1 velocity)
coeff11: +0.5 * phi_p1 * n_d * phi_v1 (elem1 pressure, elem1 velocity)

The sign convention follows the outward normal of elem0: the jump [v.n] on the face is (v0 - v1).n/|f|, hence the minus sign on elem0-side contributions.

Parameters

matrices	Map of matrix name → Matrice_Morse pointer to accumulate into.
secmem	Right-hand side (momentum residual) to accumulate into.
semi_impl	Map of semi-implicit field values; if "pression" is present, its values are used instead of the current pressure field.

Reimplemented from Operateur_base.

Definition at line 166 of file Op_Grad_DG.cpp.

◆ associer()

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

overridevirtual

Associates the operator with a DG domain and its boundary conditions.

Parameters

domaine_dis	The discretized domain, expected to be a Domaine_DG.
domaine_Cl_dis	The boundary condition container, expected to be a Domaine_Cl_DG.

Implements Operateur_base.

Definition at line 41 of file Op_Grad_DG.cpp.

◆ calculer()

DoubleTab & Op_Grad_DG::calculer	(	const DoubleTab &	pressure,
		DoubleTab &	grad ) const

overridevirtual

Reimplemented from Operateur_base.

Definition at line 285 of file Op_Grad_DG.cpp.

◆ dimensionner_blocs()

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

overridevirtual

Sizes the velocity-pressure matrix block for the gradient operator.

Builds the sparsity pattern of the rectangular matrix mapping pressure DOFs to velocity DOFs. Each velocity DOF of element T is coupled to all pressure DOFs of T and its face-neighbours, as given by the pre-computed stencil. The resulting matrix has size_v rows (velocity global DOFs) and size_p columns (pressure global DOFs).

Returns immediately if "pression" is absent from matrices or is treated semi-implicitly.

Parameters

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

Reimplemented from Operateur_base.

Definition at line 60 of file Op_Grad_DG.cpp.

◆ has_interface_blocs()

int Op_Grad_DG::has_interface_blocs ( ) const

inlineoverridevirtual

Reimplemented from Operateur_base.

Definition at line 52 of file Op_Grad_DG.h.

◆ impr()

int Op_Grad_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 292 of file Op_Grad_DG.cpp.

◆ OBS_PTR() [1/2]

Op_Grad_DG::OBS_PTR ( Domaine_Cl_DG )

protected

◆ OBS_PTR() [2/2]

Op_Grad_DG::OBS_PTR ( Domaine_DG )

protected

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

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ ajouter_blocs()

◆ associer()

◆ calculer()

◆ dimensionner_blocs()

◆ has_interface_blocs()

◆ impr()

◆ OBS_PTR() [1/2]

◆ OBS_PTR() [2/2]