Concrete DG diffusion operator acting on element-based unknowns. More...

#include <Op_Diff_DG_Elem.h>

Inheritance diagram for Op_Diff_DG_Elem:

Collaboration diagram for Op_Diff_DG_Elem:

Public Member Functions
virtual void	calculer_flux_bord (const DoubleTab &inco) const =delete
void	modifier_pour_Cl (Matrice_Morse &la_matrice, DoubleTab &secmem) const override
	DOES NOTHING - to override in derived classes.
void	completer () override
	Finalizes operator setup after all associations have been made.
void	dimensionner (Matrice_Morse &mat) const override
	Builds the sparsity pattern of the DG diffusion matrix in a Matrice_Morse.
int	has_interface_blocs () const override
void	dimensionner_blocs (matrices_t matrices, const tabs_t &semi_impl={}) const override
	Sizes the block matrices used by the interface_blocs assembly mechanism.
void	ajouter_blocs (matrices_t matrices, DoubleTab &secmem, const tabs_t &semi_impl={}) const override
	Assembles the SIP diffusion operator into the matrix and right-hand side.
void	dimensionner_termes_croises (Matrice_Morse &, const Probleme_base &autre_pb, int nl, int nc) const override
void	ajouter_termes_croises (const DoubleTab &inco, const Probleme_base &autre_pb, const DoubleTab &autre_inco, DoubleTab &resu) const override
void	contribuer_termes_croises (const DoubleTab &inco, const Probleme_base &autre_pb, const DoubleTab &autre_inco, Matrice_Morse &matrice) const override
void	contribuer_au_second_membre (DoubleTab &resu) const override
	Adds boundary condition contributions to the right-hand side.
Public Member Functions inherited from Op_Diff_DG_base
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.
double	calculer_dt_stab () const override
	Computes the maximum stable explicit time step for the diffusion operator.
void	associer_diffusivite (const Champ_base &diffu) override
const Champ_base &	diffusivite () const override
void	mettre_a_jour (double t) override
	DOES NOTHING - to override in derived classes.
void	update_nu () const
	Updates the cached effective diffusivity field nu_ by combining molecular and turbulent contributions.
double	nu (int i, int compo) const
DoubleTab &	calculer (const DoubleTab &, DoubleTab &) const override
	Computes the diffusion operator applied to inco and stores the result in resu.
int	impr (Sortie &os) const override
	DOES NOTHING - to override in derived classes.
Public Member Functions inherited from Operateur_Diff_base
virtual void	associer_diffusivite_pour_pas_de_temps (const Champ_base &)
	Associates the true diffusivity in m^2/s (in QC for example, the operator is applied to rho*u, and the dynamic viscosity is then associated to.
virtual void	associer_diffusivite_volumique (const Champ_base &)
virtual void	calculer_borne_locale (DoubleVect &, double, double) const
virtual void	init_op_ext () const
virtual bool	is_turb () const
virtual const Correlation_base *	correlation_viscosite_turbulente () const
virtual void	calculer_von_mises (const DoubleTab &deplacement, DoubleTab &deformation, DoubleTab &contraintes, DoubleTab &von_mises) const
Public Member Functions inherited from Operateur_base
virtual DoubleTab &	ajouter (const DoubleTab &, DoubleTab &) const
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	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
void	tester_contribuer_a_avec (const DoubleTab &, const Matrice_Morse &)
virtual void	calculer_dt_local (DoubleTab &) const
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
Public Member Functions inherited from Support_Champ_Masse_Volumique
	Support_Champ_Masse_Volumique ()
	Constructor of the class.
virtual	~Support_Champ_Masse_Volumique ()
	Virtual destructor (to avoid warnings).
virtual void	associer_champ_masse_volumique (const Champ_base &)
	Method to be called during problem preparation to ask the object to take into account the density field passed as parameter.
virtual void	dissocier_champ_masse_volumique ()
	Cancels the reference to the density field.
virtual int	has_champ_masse_volumique () const
	Returns 1 if the density field has been associated, 0 otherwise.
virtual const Champ_base &	get_champ_masse_volumique () const
	Returns the density field.
const int &	support_ok () const
Public Member Functions inherited from Op_Diff_Turbulent_base
virtual	~Op_Diff_Turbulent_base ()
void	associer_diffusivite_turbulente (const Champ_Fonc_base &)
const Champ_Fonc_base &	diffusivite_turbulente () const
bool	has_diffusivite_turbulente () const

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_Diff_base
std::vector< const Operateur_Diff_base * >	op_ext
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 Member Functions inherited from Op_Diff_DG_base
	OBS_PTR (Domaine_DG) le_dom_dg_
	OBS_PTR (Domaine_Cl_DG) la_zcl_dg_
	OBS_PTR (Champ_base) diffusivite_
Protected Member Functions inherited from Operateur_Diff_base
virtual const Champ_base &	diffusivite_pour_pas_de_temps () const
	Returns the field corresponding to the true diffusivity of the medium used for the time step computation.
	OBS_PTR (Champ_base) diffusivite_pour_pas_de_temps_
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 ()
Protected Member Functions inherited from Support_Champ_Masse_Volumique
virtual void	declare_support_masse_volumique (int ok)
	The constructor of a derived class that uses the density field must call this function with the value 1.
Protected Attributes inherited from Op_Diff_DG_base
SFichier	Flux
SFichier	Flux_moment
SFichier	Flux_sum
int	nu_a_jour_ = 0
DoubleTab	nu_
bool	is_var_
bool	is_aniso_
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

Concrete DG diffusion operator acting on element-based unknowns.

This class implements the Symmetric Interior Penalty (SIP) Galerkin discretization of the diffusion operator for scalar and vector fields discretized on DG elements. It assembles both the global stiffness matrix and the right-hand side contribution through the interface_blocs mechanism (has_interface_blocs() returns 1).

The SIP formulation introduces three families of face integrals at every internal face shared by elements T0 and T1:

Consistency term: -0.5 * integral of { nu * grad(phi_i) } . n * [phi_j]
Symmetry term: -0.5 * integral of { nu * grad(phi_j) } . n * [phi_i] (transpose of consistency)
Penalty term: gamma * (eta_F / h_T) * integral of [phi_i] * [phi_j]

where gamma = 2*nu_0*nu_1/(nu_0+nu_1) is the harmonic mean of the diffusivities on each side, eta_F is the element-local penalty coefficient, and h_T is the minimum characteristic size of the two adjacent elements.

At boundary faces the same three terms are applied to a single element, with the boundary condition value contributing to the right-hand side through contribuer_au_second_membre(). Supported boundary conditions are:

Neumann / Neumann_paroi: flux is added directly to the RHS.
Dirichlet / Dirichlet_homogene: enforced weakly via SIP penalty + symmetry terms.

The stencil (sparsity pattern) of the assembled matrix couples every element to all its face-neighbours, as computed by dimensionner() and dimensionner_blocs().

Cross-problem coupling (dimensionner_termes_croises, ajouter_termes_croises, contribuer_termes_croises) is declared but not yet implemented (throws at runtime).

Note: The ordering of degrees of freedom inside an element follows: phi0.ex, phi1.ex, ..., phi_{nb_bfunc-1}.ex, phi0.ey, phi1.ey, ..., phi_{nb_bfunc-1}.ey, ... for vector fields, and phi0, phi1, ..., phi_{nb_bfunc-1} for scalar fields.

Definition at line 57 of file Op_Diff_DG_Elem.h.

Member Function Documentation

◆ ajouter_blocs()

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

overridevirtual

Assembles the SIP diffusion operator into the matrix and right-hand side.

This is the core assembly routine. It proceeds in three stages:

1. Volume integrals (stiffness term) For each element, the term integral of nu * grad(phi_i) . grad(phi_j) is integrated using the element quadrature rule and accumulated into the diagonal block of the matrix and into secmem.

2. Internal face integrals (SIP terms) For each internal face shared by elem0 and elem1:

Penalty term: gamma * (eta_F/h_T) * integral of phi_i * phi_j, added to the diagonal blocks (same-element test and trial) and subtracted from the off-diagonal blocks (cross-element), enforcing the jump penalization symmetrically.
Consistency + symmetry terms: 0.5 * integral of { nu * grad(phi_i) } . n * phi_j, assembled into all four (elem0/elem1) x (elem0/elem1) block combinations with the appropriate sign to yield a symmetric bilinear form.

3. Boundary face integrals (Dirichlet enforcement) For boundary faces flagged as Dirichlet (fcl flag > 5), the same penalty and consistency terms are applied to the single adjacent element, imposing the boundary condition weakly in the SIP sense. The Dirichlet value contribution to secmem is handled separately by contribuer_au_second_membre().

Both isotropic and anisotropic diffusivities are supported: for anisotropic cases, nu_F is computed as the normal projection of the diffusivity tensor onto the face.

The DOF ordering within an element is: phi0.ex, phi1.ex, ..., phi0.ey, phi1.ey, ... for vector fields (nb_bfunc DOFs per spatial direction).

Parameters

matrices	Map of matrix name → Matrice_Morse pointer to accumulate into.
secmem	Right-hand side array to accumulate into.
semi_impl	Map of semi-implicit field values (matrix assembly is skipped if the unknown is present here).

Reimplemented from Operateur_base.

Definition at line 236 of file Op_Diff_DG_Elem.cpp.

◆ ajouter_termes_croises()

void Op_Diff_DG_Elem::ajouter_termes_croises	(	const DoubleTab &	inco,
		const Probleme_base &	autre_pb,
		const DoubleTab &	autre_inco,
		DoubleTab &	resu ) const

overridevirtual

Reimplemented from Operateur_base.

Definition at line 559 of file Op_Diff_DG_Elem.cpp.

◆ calculer_flux_bord()

virtual void Op_Diff_DG_Elem::calculer_flux_bord ( const DoubleTab & inco ) const

virtualdelete

◆ completer()

void Op_Diff_DG_Elem::completer ( )

overridevirtual

Finalizes operator setup after all associations have been made.

Calls the parent completer(), then:

Casts the unknown field to Champ_Elem_DG and checks that the ghost cell layer is at least 1 element thick (required for face-neighbour communication in DG).
Initializes the face boundary flux array flux_bords_ with the correct number of components (2 for Transport_K_Epsilon, otherwise the line size of the unknown).
If the operator is a turbulent diffusion operator (name starts with "Op_Dift"), retrieves the turbulent conductivity from the turbulence model and registers it as the turbulent diffusivity via associer_diffusivite_turbulente().

Reimplemented from Op_Diff_DG_base.

Definition at line 58 of file Op_Diff_DG_Elem.cpp.

◆ contribuer_au_second_membre()

void Op_Diff_DG_Elem::contribuer_au_second_membre ( DoubleTab & resu ) const

overridevirtual

Adds boundary condition contributions to the right-hand side.

Loops over all boundary faces and, depending on the boundary condition type, adds the appropriate weak enforcement term to resu:

Neumann / Neumann_paroi: the imposed flux g_N is integrated against each basis function on the boundary face: resu(elem, i) += integral of phi_i * g_N Point-wise flux values are retrieved either from a Champ_front_var_instationnaire (using valeur_au_temps_et_au_point) or directly from Neumann::flux_impose().
Dirichlet (weak SIP enforcement): the boundary value g_D enters two terms:
- Symmetry: -integral of { nu * grad(phi_i) } . n * g_D
- Penalty: (nu_F * eta_F / h_T) * integral of phi_i * g_D The Dirichlet value is retrieved either point-wise (Champ_front_var_instationnaire) or face-wise (Dirichlet::val_imp_au_temps()).
Dirichlet_homogene: no contribution (homogeneous condition, nothing to add).

Both isotropic and anisotropic diffusivities are handled for the Dirichlet terms. The use of Champ_front_softanalytique is explicitly forbidden and triggers an error.

Parameters

resu	The right-hand side array to accumulate boundary contributions into.

Reimplemented from Operateur_base.

Definition at line 596 of file Op_Diff_DG_Elem.cpp.

◆ contribuer_termes_croises()

void Op_Diff_DG_Elem::contribuer_termes_croises	(	const DoubleTab &	inco,
		const Probleme_base &	autre_pb,
		const DoubleTab &	autre_inco,
		Matrice_Morse &	matrice ) const

overridevirtual

Reimplemented from Operateur_base.

Definition at line 565 of file Op_Diff_DG_Elem.cpp.

◆ dimensionner()

void Op_Diff_DG_Elem::dimensionner ( Matrice_Morse & la_matrice ) const

overridevirtual

Builds the sparsity pattern of the DG diffusion matrix in a Matrice_Morse.

The matrix couples each element to all its face-neighbours as given by the pre-computed sorted stencil (Domaine_DG::get_stencil_sorted()). The global index space is built from the BasisFunction index map indices_glob_elem, so that each element contributes a block of nb_basis_func * dim rows and the column range of a row spans nb_basis_func columns per neighbour element (including itself).

The method fills tab1 (row pointers) in a first pass, then tab2 (column indices) in a second pass, and finally marks the stencil as sorted.

Parameters

la_matrice The Matrice_Morse whose sparsity pattern is to be set.

Reimplemented from Operateur_base.

Definition at line 92 of file Op_Diff_DG_Elem.cpp.

◆ dimensionner_blocs()

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

overridevirtual

Sizes the block matrices used by the interface_blocs assembly mechanism.

If the unknown is treated semi-implicitly (its name appears in semi_impl), no matrix needs to be dimensioned and the method returns immediately. Otherwise, for each external operator registered in op_ext (used for monolithic thermal coupling), the corresponding sub-matrix is sized by calling dimensionner(). Cross-problem coupling (i > 0) is not yet implemented and throws at runtime.

Parameters

matrices	Map of matrix name → Matrice_Morse pointer to be sized.
semi_impl	Map of semi-implicit field names to their current values.

Reimplemented from Operateur_base.

Definition at line 174 of file Op_Diff_DG_Elem.cpp.

◆ dimensionner_termes_croises()

void Op_Diff_DG_Elem::dimensionner_termes_croises	(	Matrice_Morse &	matrice,
		const Probleme_base &	autre_pb,
		int	nl,
		int	nc ) const

overridevirtual

Reimplemented from Operateur_base.

Definition at line 553 of file Op_Diff_DG_Elem.cpp.

◆ has_interface_blocs()

int Op_Diff_DG_Elem::has_interface_blocs ( ) const

inlineoverridevirtual

Reimplemented from Operateur_base.

Definition at line 69 of file Op_Diff_DG_Elem.h.

◆ modifier_pour_Cl()

void Op_Diff_DG_Elem::modifier_pour_Cl	(	Matrice_Morse &	,
		DoubleTab &	) const

inlineoverridevirtual

DOES NOTHING - to override in derived classes.

Parameters

(Matrice_Morse&)	a matrix in Morse format
(DoubleTab&)	a double-valued array

Exceptions

method to override

Reimplemented from Operateur_base.

Definition at line 64 of file Op_Diff_DG_Elem.h.

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

/home/docs/trust-code/src/DG/Operateurs/Op_Diff_Dift/Op_Diff_DG_Elem.h
/home/docs/trust-code/src/DG/Operateurs/Op_Diff_Dift/Op_Diff_DG_Elem.cpp

Public Member Functions

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ ajouter_blocs()

◆ ajouter_termes_croises()

◆ calculer_flux_bord()

◆ completer()

◆ contribuer_au_second_membre()

◆ contribuer_termes_croises()

◆ dimensionner()

◆ dimensionner_blocs()

◆ dimensionner_termes_croises()

◆ has_interface_blocs()

◆ modifier_pour_Cl()