#include <Assembleur_P_VDF.h>

Inheritance diagram for Assembleur_P_VDF:

Collaboration diagram for Assembleur_P_VDF:

Public Member Functions
void	associer_domaine_dis_base (const Domaine_dis_base &) override
void	associer_domaine_cl_dis_base (const Domaine_Cl_dis_base &) override
const Domaine_dis_base &	domaine_dis_base () const override
const Domaine_Cl_dis_base &	domaine_Cl_dis_base () const override
int	assembler (Matrice &) override
	Assembles the pressure matrix M such that MP = div(porosity grad(P)).
int	assembler_mat (Matrice &, const DoubleVect &, int incr_pression, int resoudre_en_u) override
int	assembler_rho_variable (Matrice &, const Champ_Don_base &rho) override
	Assembles the pressure matrix M such that MP = div(porosity/rho grad(P)).
int	assembler_QC (const DoubleTab &, Matrice &) override
	Assembles the pressure matrix for a quasi-compressible fluid.
int	modifier_secmem (DoubleTab &) override
	Modifies the right-hand side to apply boundary conditions.
int	modifier_solution (DoubleTab &) override
void	completer (const Equation_base &) override
void	dimensionner_continuite (matrices_t matrices, int aux_only=0) const override
void	assembler_continuite (matrices_t matrices, DoubleTab &secmem, int aux_only=0) const override
DoubleTab	norme_continuite () const override
Public Member Functions inherited from Assembleur_base
int	set_resoudre_increment_pression (int flag)
	Sets the value of the resoudre_increment_pression_ flag.
int	get_resoudre_increment_pression () const
	Returns the value of the resoudre_increment_pression_ flag (0 or 1) Returns -1 if the flag has not been initialized.
int	set_resoudre_en_u (int flag)
	Sets the value of the resoudre_en_u__ flag.
int	get_resoudre_en_u () const
	Returns the value of the resoudre_en_u_ flag (0 or 1) Returns -1 if the flag has not been initialized.
virtual void	modifier_secmem_pour_incr_p (const DoubleTab &press, const double fac, DoubleTab &incr) const
virtual void	corriger_vitesses (const DoubleTab &dP, DoubleTab &dv) 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 ()

Protected Member Functions
int	construire (Matrice &la_matrice)
	Determines the nonzero entries of the matrix and prepares the storage.
int	remplir (Matrice &la_matrice, const DoubleVect &volumes_entrelaces, const Champ_Don_base *rho_ptr)
	Computes the coefficients of the pressure matrix with a rho field.
void	modifier_secmem_pression_imposee (const Neumann_sortie_libre &cond_lim, const Front_VF &frontiere_vf, DoubleTab &secmem)
	Modifies the right-hand side of the pressure solver for a "Neumann_sortie_libre" condition.
void	modifier_secmem_vitesse_imposee (const Entree_fluide_vitesse_imposee &cond_lim, const Front_VF &frontiere_vf, DoubleTab &secmem)
	Modifies the right-hand side of the pressure system for an imposed velocity boundary condition.
int	liste_faces_periodiques (ArrOfInt &faces)
	Fills the array faces with the list of indices of the periodic faces in the face_voisins array.
	OBS_PTR (Domaine_VDF) le_dom_VDF
	OBS_PTR (Domaine_Cl_VDF) le_dom_Cl_VDF
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 Attributes
ArrOfDouble	les_coeff_pression
int	has_P_ref

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)
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

Detailed Description

Definition at line 28 of file Assembleur_P_VDF.h.

Member Function Documentation

◆ assembler()

int Assembleur_P_VDF::assembler ( Matrice & matrice )

overridevirtual

Assembles the pressure matrix M such that M*P = div(porosity * grad(P)).

and computes the coefficients for modifier_secmem.

Implements Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 669 of file Assembleur_P_VDF.cpp.

◆ assembler_continuite()

void Assembleur_P_VDF::assembler_continuite	(	matrices_t	matrices,
		DoubleTab &	secmem,
		int	aux_only = 0 ) const

overridevirtual

Reimplemented from Assembleur_base.

Definition at line 779 of file Assembleur_P_VDF.cpp.

◆ assembler_mat()

int Assembleur_P_VDF::assembler_mat	(	Matrice &	matrice,
		const DoubleVect &	volumes_entrelaces,
		int	incr_pression,
		int	resoudre_en_u )

overridevirtual

Reimplemented from Assembleur_base.

Definition at line 648 of file Assembleur_P_VDF.cpp.

◆ assembler_QC()

int Assembleur_P_VDF::assembler_QC	(	const DoubleTab &	tab_rho,
		Matrice &	matrice )

overridevirtual

Assembles the pressure matrix for a quasi-compressible fluid.

The matrix M is such that M*P = div( porosity * grad(P) ). The resoudre_increment_pression flag is set to zero if not yet assigned.

Parameters

tab_rho Density array.

Returns: Always returns 1.

Reimplemented from Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 724 of file Assembleur_P_VDF.cpp.

◆ assembler_rho_variable()

int Assembleur_P_VDF::assembler_rho_variable	(	Matrice &	matrice,
		const Champ_Don_base &	rho )

overridevirtual

Assembles the pressure matrix M such that M*P = div(porosity/rho * grad(P)).

and computes the coefficients for modifier_secmem.

Parameters

matrice	The matrix to assemble. Constraint: either the matrix has not yet been typed (in which case it is "constructed"), or it is the same as from the previous call.
rho	Density field.

Reimplemented from Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 691 of file Assembleur_P_VDF.cpp.

◆ associer_domaine_cl_dis_base()

void Assembleur_P_VDF::associer_domaine_cl_dis_base ( const Domaine_Cl_dis_base & le_dom_Cl_dis )

overridevirtual

Implements Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 818 of file Assembleur_P_VDF.cpp.

◆ associer_domaine_dis_base()

void Assembleur_P_VDF::associer_domaine_dis_base ( const Domaine_dis_base & le_dom_dis )

overridevirtual

Implements Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 813 of file Assembleur_P_VDF.cpp.

◆ completer()

void Assembleur_P_VDF::completer ( const Equation_base & Eqn )

overridevirtual

Implements Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 823 of file Assembleur_P_VDF.cpp.

◆ construire()

int Assembleur_P_VDF::construire ( Matrice & la_matrice )

protected

Determines the nonzero entries of the matrix and prepares the storage.

Sparse matrix of size nb_elements (rows) * nb_elem_tot (columns) Stored as a block matrix composed of two Morse matrices:

Square symmetric matrix nb_elements * nb_elements (contains the terms M(i,j) where i and j are indices of real elements)
Rectangular matrix nb_elements * (nb_elem_tot - nb_elem) (contains the terms M(i,j) where i is real and j is virtual)

Definition at line 105 of file Assembleur_P_VDF.cpp.

◆ dimensionner_continuite()

void Assembleur_P_VDF::dimensionner_continuite	(	matrices_t	matrices,
		int	aux_only = 0 ) const

overridevirtual

Reimplemented from Assembleur_base.

Definition at line 768 of file Assembleur_P_VDF.cpp.

◆ domaine_Cl_dis_base()

const Domaine_Cl_dis_base & Assembleur_P_VDF::domaine_Cl_dis_base ( ) const

overridevirtual

Implements Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 808 of file Assembleur_P_VDF.cpp.

◆ domaine_dis_base()

const Domaine_dis_base & Assembleur_P_VDF::domaine_dis_base ( ) const

overridevirtual

Implements Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 803 of file Assembleur_P_VDF.cpp.

◆ liste_faces_periodiques()

int Assembleur_P_VDF::liste_faces_periodiques ( ArrOfInt & faces )

protected

Fills the array faces with the list of indices of the periodic faces in the face_voisins array.

Each periodic face appears twice in face_voisins (each face corresponds to the opposite face). Only the one of the two with the smaller index in the list of faces of each periodic boundary is stored in the array faces. Return value: number of periodic faces (equal to the size of the faces array).

Definition at line 56 of file Assembleur_P_VDF.cpp.

◆ modifier_secmem()

int Assembleur_P_VDF::modifier_secmem ( DoubleTab & secmem )

overridevirtual

Modifies the right-hand side to apply boundary conditions.

The supported conditions are: Neumann_sortie_libre, Entree_fluide_vitesse_imposee, Dirichlet_paroi_defilante (nothing to do), Dirichlet_paroi_fixe (nothing to do), Symetrie (nothing to do)

Implements Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 476 of file Assembleur_P_VDF.cpp.

◆ modifier_secmem_pression_imposee()

void Assembleur_P_VDF::modifier_secmem_pression_imposee	(	const Neumann_sortie_libre &	cond_lim,
		const Front_VF &	frontiere_vf,
		DoubleTab &	secmem )

protected

Modifies the right-hand side of the pressure solver for a "Neumann_sortie_libre" condition.

Computation in "pressure increment" mode: add the pressure increment, i.e. zero (unsteady boundary condition not supported) Computation in "pressure" mode: Add the term Pimpose * surface / volume_entrelace to the right-hand side in the discretization of the pressure at the boundary (between an element elem0 and a fictitious exterior element with imposed pressure): grad P = (P(elem0) - Pimpose) * surface / volume_entrelace

Definition at line 539 of file Assembleur_P_VDF.cpp.

◆ modifier_secmem_vitesse_imposee()

void Assembleur_P_VDF::modifier_secmem_vitesse_imposee	(	const Entree_fluide_vitesse_imposee &	cond_lim,
		const Front_VF &	frontiere_vf,
		DoubleTab &	secmem )

protected

Modifies the right-hand side of the pressure system for an imposed velocity boundary condition.

If solving in pressure increment mode, ... otherwise nothing to do.

Definition at line 583 of file Assembleur_P_VDF.cpp.

◆ modifier_solution()

int Assembleur_P_VDF::modifier_solution ( DoubleTab & pression )

overridevirtual

Implements Assembleur_base.

Reimplemented in AssembleurPVDF_PF.

Definition at line 629 of file Assembleur_P_VDF.cpp.

◆ norme_continuite()

DoubleTab Assembleur_P_VDF::norme_continuite ( ) const

overridevirtual

Reimplemented from Assembleur_base.

Definition at line 795 of file Assembleur_P_VDF.cpp.

◆ OBS_PTR() [1/2]

Assembleur_P_VDF::OBS_PTR ( Domaine_Cl_VDF )

protected

◆ OBS_PTR() [2/2]

Assembleur_P_VDF::OBS_PTR ( Domaine_VDF )

protected

◆ remplir()

int Assembleur_P_VDF::remplir	(	Matrice &	la_matrice,
		const DoubleVect &	volumes_entrelaces,
		const Champ_Don_base *	rho_ptr )

protected

Computes the coefficients of the pressure matrix with a rho field.

If rho_ptr == 0, compute the matrix -div( porosity * grad P ), otherwise compute -div( porosity/rho grad P ) and *rho_ptr must be a Champ_Fonc_Face_VDF.

Definition at line 279 of file Assembleur_P_VDF.cpp.

Member Data Documentation

◆ has_P_ref

int Assembleur_P_VDF::has_P_ref

protected

Definition at line 60 of file Assembleur_P_VDF.h.

◆ les_coeff_pression

ArrOfDouble Assembleur_P_VDF::les_coeff_pression

protected

Definition at line 57 of file Assembleur_P_VDF.h.

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

/home/docs/trust-code/src/VDF/Solveurs/Assembleur_P_VDF.h
/home/docs/trust-code/src/VDF/Solveurs/Assembleur_P_VDF.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ assembler()

◆ assembler_continuite()

◆ assembler_mat()

◆ assembler_QC()

◆ assembler_rho_variable()

◆ associer_domaine_cl_dis_base()

◆ associer_domaine_dis_base()

◆ completer()

◆ construire()

◆ dimensionner_continuite()

◆ domaine_Cl_dis_base()

◆ domaine_dis_base()

◆ liste_faces_periodiques()

◆ modifier_secmem()

◆ modifier_secmem_pression_imposee()

◆ modifier_secmem_vitesse_imposee()

◆ modifier_solution()

◆ norme_continuite()

◆ OBS_PTR() [1/2]

◆ OBS_PTR() [2/2]

◆ remplir()

Member Data Documentation

◆ has_P_ref

◆ les_coeff_pression