#include <Matrice_Dense.h>

Inheritance diagram for Matrice_Dense:

Collaboration diagram for Matrice_Dense:

Public Member Functions
	Matrice_Dense ()
	Matrice_Dense (const int nb_lines, const int nb_cols)
void	dimensionner (const int nb_lines, const int nb_cols)
void	read_from_file (const Nom &filename)
void	convert_to_morse_matrix (Matrice_Morse &morse_matrix) const
const double &	operator() (const int line, const int col) const
double &	operator() (const int line, const int col)
void	build_matrix_from_coefficients_line_by_line (const DoubleVect &coefficients)
void	build_matrix_from_coefficients_column_by_column (const DoubleVect &coefficients)
bool	is_the_same (const Matrice_Dense &other_matrix, const double tol=1e-14) const
void	build_the_transposed (Matrice_Dense &transposed) const
void	set_coefficient (const int i, const int j, const double value)
Sortie &	imprimer_formatte (Sortie &s) const override
int	ordre () const override
	If square matrix, returns number of lines, otherwise 0.
int	nb_lignes () const override
	Return local number of lines (=size on the current proc).
int	nb_colonnes () const override
	Return local number of columns (=size on the current proc).
DoubleVect &	ajouter_multvect_ (const DoubleVect &x, DoubleVect &r) const override
	Matrix-vector multiply-accumulate (saxpy) operation.
DoubleVect &	ajouter_multvectT_ (const DoubleVect &x, DoubleVect &r) const override
	Matrix-vector multiply-accumulate (saxpy) operation using the transposed matrix.
DoubleTab &	ajouter_multTab_ (const DoubleTab &x, DoubleTab &r) const override
void	scale (const double x) override
void	clean () override
void	get_stencil (Stencil &stencil) const override
void	inverse ()
void	solve (const ArrOfDouble &b, ArrOfDouble &x)
	Solves the linear system A*x = b using LU factorization.
void	multiplyToRight (const Matrice_Dense &B, Matrice_Dense &RES) const
DoubleTab &	coeffs ()
Public Member Functions inherited from Matrice_Base
virtual DoubleVect &	ajouter_multvect (const DoubleVect &x, DoubleVect &r) const
	Matrix-vector multiply-accumulate operation (saxpy).
virtual DoubleVect &	ajouter_multvectT (const DoubleVect &x, DoubleVect &r) const
	Matrix-vector multiply-accumulate operation (saxpy).
virtual DoubleTab &	ajouter_multTab (const DoubleTab &, DoubleTab &r) const
	NOT IMPLEMENTED Matrix-matrix multiply-accumulate operation (saxpy).
virtual DoubleVect &	multvect_ (const DoubleVect &, DoubleVect &) const
virtual DoubleVect &	multvect (const DoubleVect &, DoubleVect &) const
	Multiplication of a vector by the matrix.
virtual DoubleVect &	multvectT_ (const DoubleVect &, DoubleVect &) const
virtual DoubleVect &	multvectT (const DoubleVect &, DoubleVect &) const
	Multiplication of a vector by the transposed matrix.
virtual DoubleTab &	multTab (const DoubleTab &, DoubleTab &r) const
	NOT IMPLEMENTED Multiplication of a matrix represented by an array by the matrix.
virtual Sortie &	imprimer (Sortie &) const
virtual void	get_symmetric_stencil (Stencil &stencil) const
virtual void	get_stencil_and_coefficients (Stencil &stencil, StencilCoeffs &coefficients) const
virtual void	get_stencil_and_coeff_ptrs (Stencil &stencil, std::vector< const double * > &coeff_ptr) const
virtual void	get_symmetric_stencil_and_coefficients (Stencil &stencil, StencilCoeffs &coefficients) const
int	get_stencil_size () const
virtual void	build_stencil ()
void	set_stencil (const Stencil &stencil)
bool	is_stencil_up_to_date () 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 ()

Friends
Matrice_Dense	operator+ (const Matrice_Dense &A, const Matrice_Dense &B)
Matrice_Dense	operator* (const double &a, const Matrice_Dense &B)

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
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 inherited from Matrice_Base
bool	is_stencil_up_to_date_ = false
Stencil	stencil_

Detailed Description

Definition at line 24 of file Matrice_Dense.h.

Constructor & Destructor Documentation

◆ Matrice_Dense() [1/2]

Matrice_Dense::Matrice_Dense ( )

Definition at line 64 of file Matrice_Dense.cpp.

◆ Matrice_Dense() [2/2]

Matrice_Dense::Matrice_Dense	(	const int	nb_lines,
		const int	nb_cols )

Definition at line 69 of file Matrice_Dense.cpp.

Member Function Documentation

◆ ajouter_multTab_()

DoubleTab & Matrice_Dense::ajouter_multTab_	(	const DoubleTab &	x,
		DoubleTab &	r ) const

overridevirtual

Implements Matrice_Base.

Definition at line 374 of file Matrice_Dense.cpp.

◆ ajouter_multvect_()

DoubleVect & Matrice_Dense::ajouter_multvect_	(	const DoubleVect &	x,
		DoubleVect &	resu ) const

overridevirtual

Matrix-vector multiply-accumulate (saxpy) operation.

Operation: resu = resu + Matrix_ * x

Implements Matrice_Base.

Definition at line 215 of file Matrice_Dense.cpp.

◆ ajouter_multvectT_()

DoubleVect & Matrice_Dense::ajouter_multvectT_	(	const DoubleVect &	x,
		DoubleVect &	resu ) const

overridevirtual

Matrix-vector multiply-accumulate (saxpy) operation using the transposed matrix.

Operation: resu = resu + A^{T}*x

Implements Matrice_Base.

Definition at line 299 of file Matrice_Dense.cpp.

◆ build_matrix_from_coefficients_column_by_column()

void Matrice_Dense::build_matrix_from_coefficients_column_by_column ( const DoubleVect & coefficients )

Definition at line 132 of file Matrice_Dense.cpp.

◆ build_matrix_from_coefficients_line_by_line()

void Matrice_Dense::build_matrix_from_coefficients_line_by_line ( const DoubleVect & coefficients )

Definition at line 115 of file Matrice_Dense.cpp.

◆ build_the_transposed()

void Matrice_Dense::build_the_transposed ( Matrice_Dense & transposed ) const

Definition at line 248 of file Matrice_Dense.cpp.

◆ clean()

void Matrice_Dense::clean ( )

overridevirtual

Reimplemented from Matrice_Base.

Definition at line 339 of file Matrice_Dense.cpp.

◆ coeffs()

DoubleTab & Matrice_Dense::coeffs ( )

inline

Definition at line 64 of file Matrice_Dense.h.

◆ convert_to_morse_matrix()

void Matrice_Dense::convert_to_morse_matrix ( Matrice_Morse & morse_matrix ) const

Definition at line 150 of file Matrice_Dense.cpp.

◆ dimensionner()

void Matrice_Dense::dimensionner	(	const int	nb_lines,
		const int	nb_cols )

Definition at line 181 of file Matrice_Dense.cpp.

◆ get_stencil()

void Matrice_Dense::get_stencil ( Stencil & stencil ) const

overridevirtual

Reimplemented from Matrice_Base.

Definition at line 352 of file Matrice_Dense.cpp.

◆ imprimer_formatte()

Sortie & Matrice_Dense::imprimer_formatte ( Sortie & s ) const

overridevirtual

Reimplemented from Matrice_Base.

Definition at line 46 of file Matrice_Dense.cpp.

◆ inverse()

void Matrice_Dense::inverse ( )

Definition at line 384 of file Matrice_Dense.cpp.

◆ is_the_same()

bool Matrice_Dense::is_the_same	(	const Matrice_Dense &	other_matrix,
		const double	tol = 1e-14 ) const

Definition at line 266 of file Matrice_Dense.cpp.

◆ multiplyToRight()

void Matrice_Dense::multiplyToRight	(	const Matrice_Dense &	B,
		Matrice_Dense &	RES ) const

Definition at line 453 of file Matrice_Dense.cpp.

◆ nb_colonnes()

int Matrice_Dense::nb_colonnes ( ) const

overridevirtual

Return local number of columns (=size on the current proc).

Implements Matrice_Base.

Definition at line 191 of file Matrice_Dense.cpp.

◆ nb_lignes()

int Matrice_Dense::nb_lignes ( ) const

overridevirtual

Return local number of lines (=size on the current proc).

Implements Matrice_Base.

Definition at line 186 of file Matrice_Dense.cpp.

◆ operator()() [1/2]

double & Matrice_Dense::operator()	(	const int	line,
		const int	col )

inline

Definition at line 83 of file Matrice_Dense.h.

◆ operator()() [2/2]

const double & Matrice_Dense::operator()	(	const int	line,
		const int	col ) const

inline

Definition at line 76 of file Matrice_Dense.h.

◆ ordre()

int Matrice_Dense::ordre ( ) const

overridevirtual

If square matrix, returns number of lines, otherwise 0.

Implements Matrice_Base.

Definition at line 196 of file Matrice_Dense.cpp.

◆ read_from_file()

void Matrice_Dense::read_from_file ( const Nom & filename )

Definition at line 82 of file Matrice_Dense.cpp.

◆ scale()

void Matrice_Dense::scale ( const double x )

overridevirtual

Implements Matrice_Base.

Definition at line 325 of file Matrice_Dense.cpp.

◆ set_coefficient()

void Matrice_Dense::set_coefficient	(	const int	i,
		const int	j,
		const double	value )

Definition at line 241 of file Matrice_Dense.cpp.

◆ solve()

void Matrice_Dense::solve	(	const ArrOfDouble &	b,
		ArrOfDouble &	x )

Solves the linear system A*x = b using LU factorization.

This method solves a system of linear equations using LAPACK routines:

DGETRF for LU factorization
DGETRS for solving the system using the factored matrix

Parameters

[in]	b	Right-hand side vector of the system
[out]	x	Solution vector of the system

Precondition: The matrix must be square; The size of vectors b and x must match the matrix dimensions

Exceptions

Process::exit if LU factorization or solve step fails

Definition at line 426 of file Matrice_Dense.cpp.

◆ operator*

Matrice_Dense operator*	(	const double &	a,
		const Matrice_Dense &	B )

friend

Definition at line 491 of file Matrice_Dense.cpp.

◆ operator+

Matrice_Dense operator+	(	const Matrice_Dense &	A,
		const Matrice_Dense &	B )

friend

Definition at line 476 of file Matrice_Dense.cpp.

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

/home/docs/trust-code/src/Kernel/Math/Matrices/Matrice_Dense.h
/home/docs/trust-code/src/Kernel/Math/Matrices/Matrice_Dense.cpp

Public Member Functions

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Matrice_Dense() [1/2]

◆ Matrice_Dense() [2/2]

Member Function Documentation

◆ ajouter_multTab_()

◆ ajouter_multvect_()

◆ ajouter_multvectT_()

◆ build_matrix_from_coefficients_column_by_column()

◆ build_matrix_from_coefficients_line_by_line()

◆ build_the_transposed()

◆ clean()

◆ coeffs()

◆ convert_to_morse_matrix()

◆ dimensionner()

◆ get_stencil()

◆ imprimer_formatte()

◆ inverse()

◆ is_the_same()

◆ multiplyToRight()

◆ nb_colonnes()

◆ nb_lignes()

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ ordre()

◆ read_from_file()

◆ scale()

◆ set_coefficient()

◆ solve()

Friends And Related Symbol Documentation

◆ operator*

◆ operator+