sparse matrix for big, variable sparse matrices. More...

#include <sparsematrix.h>

Classes
class	const_row_iterator

class	row_iterator

Public Types
enum	{ rows_sorted =true }

typedef AlgebraicConnection< TValueType >	connection

typedef ConstMatrixRow< this_type >	const_row_type

typedef MatrixRow< this_type >	row_type

typedef SparseMatrix< value_type >	this_type

typedef TValueType	value_type

Public Member Functions
template<typename M >
void	add (const M &mat)

void	add_matrix_row (size_t row, connection *c, size_t nr)

template<typename Vector_type >
bool	apply (Vector_type &res, const Vector_type &x) const
	calculate res = A x

template<typename vector_t >
void	apply_ignore_zero_rows (vector_t &dest, const number &beta1, const vector_t &w1) const
	calculated dest = beta1Aw1 . For empty rows, dest will not be changed

template<typename Vector_type >
bool	apply_transposed (Vector_type &res, const Vector_type &x) const
	calculate res = A.T x

template<typename vector_t >
void	apply_transposed_ignore_zero_rows (vector_t &dest, const number &beta1, const vector_t &w1) const
	calculated dest = beta1Aw1 . For empty cols of A (=empty rows of A^T), dest will not be changed

template<typename vector_t >
void	axpy (vector_t &dest, const number &alpha1, const vector_t &v1, const number &beta1, const vector_t &w1) const
	calculate dest = alpha1v1 + beta1A*w1 (A = this matrix)

template<typename vector_t >
void	axpy_transposed (vector_t &dest, const number &alpha1, const vector_t &v1, const number &beta1, const vector_t &w1) const
	calculate dest = alpha1v1 + beta1A^T*w1 (A = this matrix)

row_iterator	begin_row (size_t r)

const_row_iterator	begin_row (size_t r) const

void	check_rc (size_t r, size_t c) const

void	check_row_modifiable (size_t r) const

void	clear_and_free ()
	Clears the matrix vectors and frees their memory.

void	clear_retain_structure ()

int	col (size_t i) const

void	copy_crs (size_t &numRows, size_t &numCols, std::vector< value_type > &argValues, std::vector< int > &argRowStart, std::vector< int > &argColInd) const

void	defragment ()

void	defragment () const

row_iterator	end_row (size_t r)

const_row_iterator	end_row (size_t r) const

template<typename M >
void	get (M &mat) const
	get local matrix

row_iterator	get_connection (size_t r, size_t c)

const_row_iterator	get_connection (size_t r, size_t c) const

row_iterator	get_connection (size_t r, size_t c, bool &bFound)

const_row_iterator	get_connection (size_t r, size_t c, bool &bFound) const

void	get_crs (size_t &numRows, size_t &numCols, value_type &pValues, size_t pRowStart, size_t *pColInd, size_t &nnz) const

row_iterator	get_iterator_or_next (size_t r, size_t c)

row_type	get_row (size_t r)

const_row_type	get_row (size_t r) const

void	get_values (std::vector< value_type > &argValues) const

bool	has_connection (size_t r, size_t c) const

bool	is_isolated (size_t i) const
	check for isolated condition of an index

int	iters () const

template<typename vector_t >
void	mat_mult_add_row (size_t row, typename vector_t::value_type &dest, double alpha, const vector_t &v) const
	calculates dest += alpha * A[row, .] v;

template<typename Vector_type >
bool	matmul_minus (Vector_type &res, const Vector_type &x) const
	calculate res -= A x

size_t	num_cols () const
	returns the number of cols

size_t	num_connections (size_t i) const
	returns number of connections of row row.

size_t	num_rows () const
	returns number of rows

value_type &	operator() (size_t r, size_t c)

const value_type &	operator() (size_t r, size_t c) const

SparseMatrix< value_type > &	operator*= (double d)

SparseMatrix< value_type > &	operator= (const SparseMatrix< value_type > &B)

void	p () const

void	pr (size_t row) const

void	print (const char *const name=NULL) const

void	print_to_file (const char *filename) const

void	printrow (size_t row) const

void	printtype () const

void	resize_and_clear (size_t newRows, size_t newCols)
	resizes the SparseMatrix

void	resize_and_keep_values (size_t newRows, size_t newCols)

void	scale (double d)

template<typename M >
void	set (const M &mat)
	set local matrix

void	set (double a)
	set matrix to Id*a

void	set_as_copy_of (const SparseMatrix< value_type > &B, double scale=1.0)
	create/recreate this as a copy of SparseMatrix B

void	set_as_transpose_of (const SparseMatrix< value_type > &B, double scale=1.0)
	write in a empty SparseMatrix (this) the transpose SparseMatrix of B.

void	set_as_transpose_of2 (const SparseMatrix< value_type > &B, double scale=1.0)

void	set_matrix_row (size_t row, connection *c, size_t nr)

	SparseMatrix ()
	constructor for empty SparseMatrix

size_t	total_num_connections () const
	returns the total number of connections

virtual	~SparseMatrix ()
	destructor

Public Attributes
std::vector< int >	nrOfRowIterators

Protected Member Functions
void	check_fragmentation () const

void	copyToNewSize (size_t newSize)

void	copyToNewSize (size_t newSize, size_t maxCols)

int	get_index (int r, int c)

int	get_index_const (int r, int c) const

int	get_index_internal (size_t row, int col) const

int	get_nnz_max_cols (size_t maxCols)

Protected Attributes
bool	bNeedsValues

std::vector< int >	cols

size_t	fragmented

int	iIterators

int	m_numCols

int	maxValues

size_t	nnz

std::vector< int >	rowEnd

std::vector< int >	rowMax

std::vector< int >	rowStart

std::vector< value_type >	values

Private Member Functions
void	add_iterator (size_t row) const

void	assureValuesSize (size_t s)

void	check_row (size_t row, int i) const

size_t	get_nnz () const

void	remove_iterator (size_t row) const

	SparseMatrix (SparseMatrix &)
	disallow copy operator

Friends
std::ostream &	operator<< (std::ostream &out, const SparseMatrix &m)

Detailed Description

template<typename TValueType>
class ug::SparseMatrix< TValueType >

sparse matrix for big, variable sparse matrices.

SparseMatrix

matrix is stored independent row-wise When doing discretisation, use the add set and get methods for dealing with submatrices of A. For other things you can use the row iterators or operator()-methods.

See also: matrixrow, CreateAsMultiplyOf

Parameters

T	blocktype
T	blocktype

Member Typedef Documentation

◆ connection

template<typename TValueType >

typedef AlgebraicConnection<TValueType> ug::SparseMatrix< TValueType >::connection

◆ const_row_type

template<typename TValueType >

typedef ConstMatrixRow<this_type> ug::SparseMatrix< TValueType >::const_row_type

◆ row_type

template<typename TValueType >

typedef MatrixRow<this_type> ug::SparseMatrix< TValueType >::row_type

◆ this_type

template<typename TValueType >

typedef SparseMatrix<value_type> ug::SparseMatrix< TValueType >::this_type

◆ value_type

template<typename TValueType >

typedef TValueType ug::SparseMatrix< TValueType >::value_type

Member Enumeration Documentation

◆ anonymous enum

template<typename TValueType >

anonymous enum

Enumerator
rows_sorted

Constructor & Destructor Documentation

◆ SparseMatrix() [1/2]

template<typename T >

SparseMatrix< T >::SparseMatrix ( )

constructor for empty SparseMatrix

References PROFILE_SPMATRIX.

◆ ~SparseMatrix()

template<typename TValueType >

virtual ug::SparseMatrix< TValueType >::~SparseMatrix ( )

inlinevirtual

destructor

◆ SparseMatrix() [2/2]

template<typename TValueType >

ug::SparseMatrix< TValueType >::SparseMatrix ( SparseMatrix< TValueType > & )

private

disallow copy operator

Member Function Documentation

◆ add()

template<typename T >

template<typename M >

void SparseMatrix< T >::add ( const M & mat )

Add a local matrix

The local matrix type must declare the following members:

num_rows()
num_cols()
row_index(size_t i)
col_index(size_t j)
operator()(size_t i, size_t j) so that mat(i,j) will go to SparseMat(mat.row_index(i), mat.col_index(j))
Parameters

mat the whole local matrix type

◆ add_iterator()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::add_iterator ( size_t row ) const

inlineprivate

References ug::SparseMatrix< TValueType >::iIterators, and ug::SparseMatrix< TValueType >::nrOfRowIterators.

Referenced by ug::SparseMatrix< TValueType >::const_row_iterator::const_row_iterator(), ug::SparseMatrix< TValueType >::const_row_iterator::const_row_iterator(), ug::SparseMatrix< TValueType >::const_row_iterator::const_row_iterator(), ug::SparseMatrix< TValueType >::row_iterator::row_iterator(), ug::SparseMatrix< TValueType >::row_iterator::row_iterator(), and ug::SparseMatrix< TValueType >::row_iterator::row_iterator().

◆ add_matrix_row()

template<typename T >

void SparseMatrix< T >::add_matrix_row	(	size_t	row,
		connection *	c,
		size_t	nr
	)

adds the connections c to the matrixrow row. if c has a connection con with con.iIndex=i, and the matrix already has a connection (row, i), the function will set A(row,i) += con.dValue. otherwise the connection A(row, i) is created and set to con.dValue.

Parameters

row	row to add to
c	connections ("row") to be added the row.
nr	number of connections in array c.

Returns: true on success.

Note: you may use double connections in c.

References ug::AlgebraicConnection< T >::iIndex.

◆ apply()

template<typename TValueType >

template<typename Vector_type >

bool ug::SparseMatrix< TValueType >::apply	(	Vector_type &	res,
		const Vector_type &	x
	)		const

inline

calculate res = A x

References ug::SparseMatrix< TValueType >::axpy().

◆ apply_ignore_zero_rows()

template<typename T >

template<typename vector_t >

void SparseMatrix< T >::apply_ignore_zero_rows	(	vector_t &	dest,
		const number &	beta1,
		const vector_t &	w1
	)		const

calculated dest = beta1*A*w1 . For empty rows, dest will not be changed

References ug::MatMult(), ug::MatMultAdd(), and num_rows().

◆ apply_transposed()

template<typename TValueType >

template<typename Vector_type >

bool ug::SparseMatrix< TValueType >::apply_transposed	(	Vector_type &	res,
		const Vector_type &	x
	)		const

inline

calculate res = A.T x

References ug::SparseMatrix< TValueType >::axpy_transposed().

◆ apply_transposed_ignore_zero_rows()

template<typename T >

template<typename vector_t >

void SparseMatrix< T >::apply_transposed_ignore_zero_rows	(	vector_t &	dest,
		const number &	beta1,
		const vector_t &	w1
	)		const

calculated dest = beta1*A*w1 . For empty cols of A (=empty rows of A^T), dest will not be changed

References begin_row(), end_row(), ug::MatMultTransposedAdd(), and num_rows().

◆ assureValuesSize()

template<typename T >

void SparseMatrix< T >::assureValuesSize ( size_t s )

private

References s.

◆ axpy()

template<typename T >

template<typename vector_t >

void SparseMatrix< T >::axpy	(	vector_t &	dest,
		const number &	alpha1,
		const vector_t &	v1,
		const number &	beta1,
		const vector_t &	w1
	)		const

calculate dest = alpha1*v1 + beta1*A*w1 (A = this matrix)

References num_rows(), and PROFILE_SPMATRIX.

Referenced by ug::SparseMatrix< TValueType >::apply(), and ug::SparseMatrix< TValueType >::matmul_minus().

◆ axpy_transposed()

template<typename T >

template<typename vector_t >

void SparseMatrix< T >::axpy_transposed	(	vector_t &	dest,
		const number &	alpha1,
		const vector_t &	v1,
		const number &	beta1,
		const vector_t &	w1
	)		const

calculate dest = alpha1*v1 + beta1*A^T*w1 (A = this matrix)

References ug::MatMultTransposedAdd(), num_rows(), PROFILE_SPMATRIX, and ug::VecScaleAssign().

Referenced by ug::SparseMatrix< TValueType >::apply_transposed(), and ug::MatMultTransposedAdd().

◆ begin_row() [1/2]

template<typename TValueType >

row_iterator ug::SparseMatrix< TValueType >::begin_row ( size_t r )

inline

References ug::SparseMatrix< TValueType >::rowStart.

Referenced by boost::adjacent_vertices(), ug::SparseMatrix< TValueType >::set_as_copy_of(), ug::MapSparseMatrix< TValueType >::set_as_transpose_of(), and ug::SparseMatrix< TValueType >::set_as_transpose_of2().

◆ begin_row() [2/2]

template<typename TValueType >

const_row_iterator ug::SparseMatrix< TValueType >::begin_row ( size_t r ) const

inline

References ug::SparseMatrix< TValueType >::rowStart.

◆ check_fragmentation()

template<typename T >

void SparseMatrix< T >::check_fragmentation ( ) const

protected

References defragment().

◆ check_rc()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::check_rc	(	size_t	r,
		size_t	c
	)		const

inline

References ug::SparseMatrix< TValueType >::num_cols(), ug::SparseMatrix< TValueType >::num_rows(), and UG_ASSERT.

Referenced by ug::SparseMatrix< TValueType >::get_connection(), ug::SparseMatrix< TValueType >::get_connection(), ug::SparseMatrix< TValueType >::get_connection(), ug::SparseMatrix< TValueType >::get_iterator_or_next(), ug::SparseMatrix< TValueType >::has_connection(), ug::SparseMatrix< TValueType >::operator()(), and ug::SparseMatrix< TValueType >::operator()().

◆ check_row()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::check_row	(	size_t	row,
		int	i
	)		const

inlineprivate

References ug::SparseMatrix< TValueType >::rowEnd, ug::SparseMatrix< TValueType >::rowStart, and UG_ASSERT.

Referenced by ug::SparseMatrix< TValueType >::row_iterator::check(), and ug::SparseMatrix< TValueType >::const_row_iterator::check().

◆ check_row_modifiable()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::check_row_modifiable ( size_t r ) const

inline

References ug::SparseMatrix< TValueType >::nrOfRowIterators, and UG_ASSERT.

◆ clear_and_free()

template<typename T >

void SparseMatrix< T >::clear_and_free ( )

Clears the matrix vectors and frees their memory.

◆ clear_retain_structure()

template<typename T >

void SparseMatrix< T >::clear_retain_structure ( )

◆ col()

template<typename TValueType >

int ug::SparseMatrix< TValueType >::col ( size_t i ) const

inline

References ug::SparseMatrix< TValueType >::cols.

◆ copy_crs()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::copy_crs	(	size_t &	numRows,
		size_t &	numCols,
		std::vector< value_type > &	argValues,
		std::vector< int > &	argRowStart,
		std::vector< int > &	argColInd
	)		const

inline

copies the matrix to the standard CRS format

Parameters

numRows	(out) num rows of A
numCols	(out) num rows of A
argValues	(out) value_type vector with non-zero values
argRowStart	(out) row i is from argRowStart[i] to argRowStart[i+1]
argColInd	(out) argColInd[i] is colum index of nonzero i

References ug::SparseMatrix< TValueType >::cols, ug::SparseMatrix< TValueType >::defragment(), ug::SparseMatrix< TValueType >::num_cols(), ug::SparseMatrix< TValueType >::num_rows(), ug::SparseMatrix< TValueType >::rowStart, and ug::SparseMatrix< TValueType >::values.

◆ copyToNewSize() [1/2]

template<typename TValueType >

void ug::SparseMatrix< TValueType >::copyToNewSize ( size_t newSize )

inlineprotected

References ug::SparseMatrix< TValueType >::copyToNewSize(), and ug::SparseMatrix< TValueType >::num_cols().

Referenced by ug::SparseMatrix< TValueType >::copyToNewSize(), and ug::SparseMatrix< TValueType >::defragment().

◆ copyToNewSize() [2/2]

template<typename T >

void SparseMatrix< T >::copyToNewSize	(	size_t	newSize,
		size_t	maxCols
	)

protected

◆ defragment() [1/2]

template<typename TValueType >

void ug::SparseMatrix< TValueType >::defragment ( )

inline

References ug::SparseMatrix< TValueType >::copyToNewSize(), ug::SparseMatrix< TValueType >::nnz, ug::SparseMatrix< TValueType >::num_cols(), and ug::SparseMatrix< TValueType >::num_rows().

Referenced by ug::SparseMatrix< TValueType >::copy_crs(), ug::SparseMatrix< TValueType >::get_crs(), and ug::SparseMatrix< TValueType >::get_values().

◆ defragment() [2/2]

template<typename TValueType >

void ug::SparseMatrix< TValueType >::defragment ( ) const

inline

◆ end_row() [1/2]

template<typename TValueType >

row_iterator ug::SparseMatrix< TValueType >::end_row ( size_t r )

inline

References ug::SparseMatrix< TValueType >::rowEnd.

Referenced by boost::adjacent_vertices(), ug::SparseMatrix< TValueType >::get_connection(), ug::SparseMatrix< TValueType >::get_connection(), ug::SparseMatrix< TValueType >::get_iterator_or_next(), ug::SparseMatrix< TValueType >::set_as_copy_of(), ug::MapSparseMatrix< TValueType >::set_as_transpose_of(), and ug::SparseMatrix< TValueType >::set_as_transpose_of2().

◆ end_row() [2/2]

template<typename TValueType >

const_row_iterator ug::SparseMatrix< TValueType >::end_row ( size_t r ) const

inline

References ug::SparseMatrix< TValueType >::rowEnd.

◆ get()

template<typename T >

template<typename M >

void SparseMatrix< T >::get ( M & mat ) const

get local matrix

See also: add

◆ get_connection() [1/4]

template<typename TValueType >

row_iterator ug::SparseMatrix< TValueType >::get_connection	(	size_t	r,
		size_t	c
	)

inline

Parameters

r	index of the row
c	index of the column

Returns: a row_iterator to the connection A(r,c)

Remarks: creates connection if necessary.

References ug::SparseMatrix< TValueType >::bNeedsValues, ug::SparseMatrix< TValueType >::check_rc(), and ug::SparseMatrix< TValueType >::get_index().

◆ get_connection() [2/4]

template<typename TValueType >

const_row_iterator ug::SparseMatrix< TValueType >::get_connection	(	size_t	r,
		size_t	c
	)		const

inline

Parameters

r	index of the row
c	index of the column

Returns: a const_row_iterator to the connection A(r,c) if existing, otherwise end_row(row)

References ug::SparseMatrix< TValueType >::get_connection().

◆ get_connection() [3/4]

template<typename TValueType >

row_iterator ug::SparseMatrix< TValueType >::get_connection	(	size_t	r,
		size_t	c,
		bool &	bFound
	)

inline

Parameters

r	index of the row
c	index of the column

Returns: a row_iterator to the connection A(r,c) if existing, otherwise end_row(row)

References ug::SparseMatrix< TValueType >::check_rc(), ug::SparseMatrix< TValueType >::end_row(), and ug::SparseMatrix< TValueType >::get_index_const().

◆ get_connection() [4/4]

template<typename TValueType >

const_row_iterator ug::SparseMatrix< TValueType >::get_connection	(	size_t	r,
		size_t	c,
		bool &	bFound
	)		const

inline

Parameters

r	index of the row
c	index of the column

Returns: a const_row_iterator to the connection A(r,c) if existing, otherwise end_row(row)

References ug::SparseMatrix< TValueType >::check_rc(), ug::SparseMatrix< TValueType >::end_row(), and ug::SparseMatrix< TValueType >::get_index_const().

Referenced by ug::SparseMatrix< TValueType >::get_connection(), and ug::SparseMatrix< TValueType >::has_connection().

◆ get_crs()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::get_crs	(	size_t &	numRows,
		size_t &	numCols,
		value_type *&	pValues,
		size_t *	pRowStart,
		size_t *	pColInd,
		size_t &	nnz
	)		const

inline

returns pointers to CRS format. note that these are only valid as long as the matrix is not modified.

Parameters

numRows	(out) num rows of A
numCols	(out) num rows of A
pValues	(out) value_type vector with non-zero values
pRowStart	(out) row i is from pRowStart[i] to pRowStart[i+1]
pColInd	(out) pColInd[i] is colum index of nonzero i

References ug::SparseMatrix< TValueType >::defragment(), ug::SparseMatrix< TValueType >::nnz, ug::SparseMatrix< TValueType >::num_cols(), ug::SparseMatrix< TValueType >::num_rows(), ug::SparseMatrix< TValueType >::total_num_connections(), UG_THROW, and ug::SparseMatrix< TValueType >::values.

◆ get_index()

template<typename T >

int SparseMatrix< T >::get_index	(	int	r,
		int	c
	)

protected

Referenced by ug::SparseMatrix< TValueType >::get_connection(), and ug::SparseMatrix< TValueType >::operator()().

◆ get_index_const()

template<typename T >

int SparseMatrix< T >::get_index_const	(	int	r,
		int	c
	)		const

protected

Referenced by ug::SparseMatrix< TValueType >::get_connection(), ug::SparseMatrix< TValueType >::get_connection(), and ug::SparseMatrix< TValueType >::operator()().

◆ get_index_internal()

template<typename T >

int SparseMatrix< T >::get_index_internal	(	size_t	row,
		int	col
	)		const

protected

References UG_ASSERT.

Referenced by ug::SparseMatrix< TValueType >::get_iterator_or_next().

◆ get_iterator_or_next()

template<typename TValueType >

row_iterator ug::SparseMatrix< TValueType >::get_iterator_or_next	(	size_t	r,
		size_t	c
	)

inline

Parameters

r	index of the row
c	index of the column

Returns: a const_row_iterator to the connection A(r,c) if existing, otherwise end_row(row)

References ug::SparseMatrix< TValueType >::check_rc(), ug::SparseMatrix< TValueType >::end_row(), ug::SparseMatrix< TValueType >::get_index_internal(), ug::SparseMatrix< TValueType >::maxValues, ug::SparseMatrix< TValueType >::rowEnd, and ug::SparseMatrix< TValueType >::rowStart.

◆ get_nnz()

template<typename TValueType >

size_t ug::SparseMatrix< TValueType >::get_nnz ( ) const

inlineprivate

References ug::SparseMatrix< TValueType >::nnz.

◆ get_nnz_max_cols()

template<typename T >

int SparseMatrix< T >::get_nnz_max_cols ( size_t maxCols )

protected

References num_rows().

◆ get_row() [1/2]

template<typename TValueType >

row_type ug::SparseMatrix< TValueType >::get_row ( size_t r )

inline

◆ get_row() [2/2]

template<typename TValueType >

const_row_type ug::SparseMatrix< TValueType >::get_row ( size_t r ) const

inline

◆ get_values()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::get_values ( std::vector< value_type > & argValues ) const

inline

assigns a reference to the values vector to argument vector

Parameters

argValues vector to be assigned

References ug::SparseMatrix< TValueType >::defragment(), and ug::SparseMatrix< TValueType >::values.

◆ has_connection()

template<typename TValueType >

bool ug::SparseMatrix< TValueType >::has_connection	(	size_t	r,
		size_t	c
	)		const

inline

References ug::SparseMatrix< TValueType >::check_rc(), and ug::SparseMatrix< TValueType >::get_connection().

◆ is_isolated()

template<typename T >

bool SparseMatrix< T >::is_isolated ( size_t i ) const

inline

check for isolated condition of an index

Parameters

i

Returns: true if only A[i,i] != 0.0

References begin_row(), end_row(), num_rows(), and UG_ASSERT.

◆ iters()

template<typename TValueType >

int ug::SparseMatrix< TValueType >::iters ( ) const

inline

References ug::SparseMatrix< TValueType >::iIterators.

◆ mat_mult_add_row()

template<typename T >

template<typename vector_t >

void SparseMatrix< T >::mat_mult_add_row	(	size_t	row,
		typename vector_t::value_type &	dest,
		double	alpha,
		const vector_t &	v
	)		const

inline

calculates dest += alpha * A[row, .] v;

◆ matmul_minus()

template<typename TValueType >

template<typename Vector_type >

bool ug::SparseMatrix< TValueType >::matmul_minus	(	Vector_type &	res,
		const Vector_type &	x
	)		const

inline

calculate res -= A x

References ug::SparseMatrix< TValueType >::axpy().

◆ num_cols()

template<typename TValueType >

size_t ug::SparseMatrix< TValueType >::num_cols ( ) const

inline

returns the number of cols

References ug::SparseMatrix< TValueType >::m_numCols.

Referenced by ug::SparseMatrix< TValueType >::check_rc(), ug::SparseMatrix< TValueType >::copy_crs(), ug::SparseMatrix< TValueType >::copyToNewSize(), ug::SparseMatrix< TValueType >::defragment(), ug::SparseMatrix< TValueType >::get_crs(), boost::num_vertices(), ug::SparseMatrix< T >::set_as_copy_of(), ug::SparseMatrix< TValueType >::set_as_copy_of(), ug::MapSparseMatrix< TValueType >::set_as_transpose_of(), and ug::SparseMatrix< TValueType >::set_as_transpose_of2().

◆ num_connections()

template<typename TValueType >

size_t ug::SparseMatrix< TValueType >::num_connections ( size_t i ) const

inline

returns number of connections of row row.

References ug::SparseMatrix< TValueType >::rowEnd, and ug::SparseMatrix< TValueType >::rowStart.

◆ num_rows()

template<typename TValueType >

size_t ug::SparseMatrix< TValueType >::num_rows ( ) const

inline

returns number of rows

References ug::SparseMatrix< TValueType >::rowEnd.

Referenced by boost::adjacent_vertices(), ug::SparseMatrix< TValueType >::check_rc(), ug::SparseMatrix< TValueType >::copy_crs(), ug::SparseMatrix< TValueType >::defragment(), ug::SparseMatrix< TValueType >::get_crs(), boost::num_vertices(), boost::out_edges(), ug::SparseMatrix< T >::set_as_copy_of(), ug::SparseMatrix< TValueType >::set_as_copy_of(), ug::MapSparseMatrix< TValueType >::set_as_transpose_of(), ug::SparseMatrix< TValueType >::set_as_transpose_of2(), and boost::vertices().

◆ operator()() [1/2]

template<typename TValueType >

value_type & ug::SparseMatrix< TValueType >::operator()	(	size_t	r,
		size_t	c
	)

inline

operator() (size_t r, size_t c) const access or create connection (r, c)

Parameters

r	row
c	column

Note: (r,c) is added to sparsity pattern if not already there use operator()(r,c,bConnectionFound) to prevent

Returns: SparseMat(r, c)=0.0 if connection created, otherwise SparseMat(r, c)

References ug::SparseMatrix< TValueType >::check_rc(), ug::SparseMatrix< TValueType >::cols, ug::SparseMatrix< TValueType >::get_index(), ug::SparseMatrix< TValueType >::rowEnd, ug::SparseMatrix< TValueType >::rowStart, UG_ASSERT, and ug::SparseMatrix< TValueType >::values.

◆ operator()() [2/2]

template<typename TValueType >

const value_type & ug::SparseMatrix< TValueType >::operator()	(	size_t	r,
		size_t	c
	)		const

inline

operator() (size_t r, size_t c) const access connection (r, c)

Parameters

r	row
c	column

Note: if connection (r, c) is not there, returns 0.0

Returns: SparseMat(r, c)

References ug::SparseMatrix< TValueType >::check_rc(), ug::SparseMatrix< TValueType >::cols, ug::SparseMatrix< TValueType >::get_index_const(), ug::SparseMatrix< TValueType >::rowEnd, ug::SparseMatrix< TValueType >::rowStart, UG_ASSERT, and ug::SparseMatrix< TValueType >::values.

◆ operator*=()

template<typename TValueType >

SparseMatrix< value_type > & ug::SparseMatrix< TValueType >::operator*= ( double d )

inline

References ug::SparseMatrix< TValueType >::scale().

◆ operator=()

template<typename TValueType >

SparseMatrix< value_type > & ug::SparseMatrix< TValueType >::operator= ( const SparseMatrix< value_type > & B )

inline

References ug::SparseMatrix< TValueType >::set_as_copy_of().

◆ p()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::p ( ) const

inline

References ug::SparseMatrix< TValueType >::print().

◆ pr()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::pr ( size_t row ) const

inline

References ug::SparseMatrix< TValueType >::printrow().

◆ print_to_file()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::print_to_file ( const char * filename ) const

◆ remove_iterator()

template<typename TValueType >

void ug::SparseMatrix< TValueType >::remove_iterator ( size_t row ) const

inlineprivate

References ug::SparseMatrix< TValueType >::iIterators, ug::SparseMatrix< TValueType >::nrOfRowIterators, and UG_ASSERT.

Referenced by ug::SparseMatrix< TValueType >::const_row_iterator::~const_row_iterator(), and ug::SparseMatrix< TValueType >::row_iterator::~row_iterator().

◆ resize_and_clear()

template<typename T >

void SparseMatrix< T >::resize_and_clear	(	size_t	newRows,
		size_t	newCols
	)

resizes the SparseMatrix

Parameters

newRows	new nr of rows
newCols	new nr of cols

Returns

References PROFILE_SPMATRIX.

◆ resize_and_keep_values()

template<typename T >

void SparseMatrix< T >::resize_and_keep_values	(	size_t	newRows,
		size_t	newCols
	)

References PROFILE_SPMATRIX.

◆ scale()

template<typename T >

void SparseMatrix< T >::scale ( double d )

References begin_row(), end_row(), and num_rows().

Referenced by ug::SparseMatrix< TValueType >::operator*=().

◆ set() [1/2]

template<typename T >

template<typename M >

void SparseMatrix< T >::set ( const M & mat )

set local matrix

See also: add

◆ set() [2/2]

template<typename T >

void SparseMatrix< T >::set ( double a )

set matrix to Id*a

References begin_row(), end_row(), num_rows(), and PROFILE_SPMATRIX.

◆ set_as_copy_of()

template<typename TValueType >

void SparseMatrix< T >::set_as_copy_of	(	const SparseMatrix< value_type > &	B,
		double	scale = `1.0`
	)

create/recreate this as a copy of SparseMatrix B

Parameters

B	the matrix of which to create a copy of
scale	an optional scaling

Returns: true on success

References ug::SparseMatrix< TValueType >::begin_row(), ug::SparseMatrix< TValueType >::end_row(), ug::SparseMatrix< TValueType >::num_cols(), ug::SparseMatrix< TValueType >::num_rows(), and ug::SparseMatrix< TValueType >::row_iterator::value().

Referenced by ug::SparseMatrix< TValueType >::operator=().

◆ set_as_transpose_of()

template<typename T >

void SparseMatrix< T >::set_as_transpose_of	(	const SparseMatrix< value_type > &	B,
		double	scale = `1.0`
	)

write in a empty SparseMatrix (this) the transpose SparseMatrix of B.

Parameters

B	the matrix of which to create the transpose of
scale	an optional scaling

Returns: true on success

◆ set_as_transpose_of2()

template<typename T >

void SparseMatrix< T >::set_as_transpose_of2	(	const SparseMatrix< value_type > &	B,
		double	scale = `1.0`
	)

References ug::SparseMatrix< TValueType >::begin_row(), ug::SparseMatrix< TValueType >::bNeedsValues, ug::SparseMatrix< TValueType >::end_row(), ug::iszero(), ug::SparseMatrix< TValueType >::maxValues, ug::SparseMatrix< TValueType >::num_cols(), num_cols(), ug::SparseMatrix< TValueType >::num_rows(), num_rows(), and PROFILE_SPMATRIX.

◆ set_matrix_row()

template<typename T >

void SparseMatrix< T >::set_matrix_row	(	size_t	row,
		connection *	c,
		size_t	nr
	)

set a row of the matrix. all previous content in this row is destroyed (

See also: add_matrix_row).

Parameters

row	index of the row to set
c	pointer to a array of sorted connections of size nr
nr	number of connections in c

Remarks: will sort array c

References ug::AlgebraicConnection< T >::iIndex.

◆ total_num_connections()

template<typename TValueType >

size_t ug::SparseMatrix< TValueType >::total_num_connections ( ) const

inline

returns the total number of connections

References ug::SparseMatrix< TValueType >::nnz.

Referenced by ug::SparseMatrix< TValueType >::get_crs().

Friends And Related Symbol Documentation

◆ operator<<

template<typename TValueType >

std::ostream & operator<<	(	std::ostream &	out,
		const SparseMatrix< TValueType > &	m
	)

friend

Member Data Documentation

◆ bNeedsValues

template<typename TValueType >

bool ug::SparseMatrix< TValueType >::bNeedsValues

protected

Referenced by ug::SparseMatrix< TValueType >::get_connection(), and ug::SparseMatrix< TValueType >::set_as_transpose_of2().

◆ cols

template<typename TValueType >

std::vector<int> ug::SparseMatrix< TValueType >::cols

protected

Referenced by ug::SparseMatrix< TValueType >::col(), ug::SparseMatrix< TValueType >::copy_crs(), ug::SparseMatrix< TValueType >::row_iterator::index(), ug::SparseMatrix< TValueType >::const_row_iterator::index(), ug::SparseMatrix< TValueType >::operator()(), and ug::SparseMatrix< TValueType >::operator()().

◆ fragmented

template<typename TValueType >

size_t ug::SparseMatrix< TValueType >::fragmented

protected

◆ iIterators

template<typename TValueType >

int ug::SparseMatrix< TValueType >::iIterators

mutableprotected

Referenced by ug::SparseMatrix< TValueType >::add_iterator(), ug::SparseMatrix< TValueType >::iters(), and ug::SparseMatrix< TValueType >::remove_iterator().

◆ m_numCols

template<typename TValueType >

int ug::SparseMatrix< TValueType >::m_numCols

protected

Referenced by ug::SparseMatrix< TValueType >::num_cols().

◆ maxValues

template<typename TValueType >

int ug::SparseMatrix< TValueType >::maxValues

protected

Referenced by ug::SparseMatrix< TValueType >::get_iterator_or_next(), and ug::SparseMatrix< TValueType >::set_as_transpose_of2().

◆ nnz

template<typename TValueType >

size_t ug::SparseMatrix< TValueType >::nnz

protected

Referenced by ug::SparseMatrix< TValueType >::defragment(), ug::SparseMatrix< TValueType >::get_crs(), ug::SparseMatrix< TValueType >::get_nnz(), and ug::SparseMatrix< TValueType >::total_num_connections().

◆ nrOfRowIterators

template<typename TValueType >

std::vector<int> ug::SparseMatrix< TValueType >::nrOfRowIterators

mutable

Referenced by ug::SparseMatrix< TValueType >::add_iterator(), ug::SparseMatrix< TValueType >::check_row_modifiable(), and ug::SparseMatrix< TValueType >::remove_iterator().

◆ rowEnd

template<typename TValueType >

std::vector<int> ug::SparseMatrix< TValueType >::rowEnd

protected

Referenced by ug::SparseMatrix< TValueType >::check_row(), ug::SparseMatrix< TValueType >::end_row(), ug::SparseMatrix< TValueType >::end_row(), ug::SparseMatrix< TValueType >::get_iterator_or_next(), ug::SparseMatrix< TValueType >::num_connections(), ug::SparseMatrix< TValueType >::num_rows(), ug::SparseMatrix< TValueType >::operator()(), and ug::SparseMatrix< TValueType >::operator()().

◆ rowMax

template<typename TValueType >

std::vector<int> ug::SparseMatrix< TValueType >::rowMax

protected

◆ rowStart

template<typename TValueType >

std::vector<int> ug::SparseMatrix< TValueType >::rowStart

protected

Referenced by ug::SparseMatrix< TValueType >::begin_row(), ug::SparseMatrix< TValueType >::begin_row(), ug::SparseMatrix< TValueType >::check_row(), ug::SparseMatrix< TValueType >::copy_crs(), ug::SparseMatrix< TValueType >::get_iterator_or_next(), ug::SparseMatrix< TValueType >::num_connections(), ug::SparseMatrix< TValueType >::operator()(), and ug::SparseMatrix< TValueType >::operator()().

◆ values

template<typename TValueType >

std::vector<value_type> ug::SparseMatrix< TValueType >::values

protected

Referenced by ug::SparseMatrix< TValueType >::copy_crs(), ug::SparseMatrix< TValueType >::get_crs(), ug::SparseMatrix< TValueType >::get_values(), ug::SparseMatrix< TValueType >::operator()(), ug::SparseMatrix< TValueType >::operator()(), ug::SparseMatrix< TValueType >::row_iterator::value(), and ug::SparseMatrix< TValueType >::const_row_iterator::value().

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

ugbase/lib_algebra/cpu_algebra/sparsematrix.h
ugbase/lib_algebra/cpu_algebra/sparsematrix_impl.h
ugbase/lib_algebra/cpu_algebra/sparsematrix_print.h

Classes

Public Types

Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Friends

Detailed Description

Member Typedef Documentation

◆ connection

◆ const_row_type

◆ row_type

◆ this_type

◆ value_type

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ SparseMatrix() [1/2]

◆ ~SparseMatrix()

◆ SparseMatrix() [2/2]

Member Function Documentation

◆ add()

◆ add_iterator()

◆ add_matrix_row()

◆ apply()

◆ apply_ignore_zero_rows()

◆ apply_transposed()

◆ apply_transposed_ignore_zero_rows()

◆ assureValuesSize()

◆ axpy()

◆ axpy_transposed()

◆ begin_row() [1/2]

◆ begin_row() [2/2]

◆ check_fragmentation()

◆ check_rc()

◆ check_row()

◆ check_row_modifiable()

◆ clear_and_free()

◆ clear_retain_structure()

◆ col()

◆ copy_crs()

◆ copyToNewSize() [1/2]

◆ copyToNewSize() [2/2]

◆ defragment() [1/2]

◆ defragment() [2/2]

◆ end_row() [1/2]

◆ end_row() [2/2]

◆ get()

◆ get_connection() [1/4]

◆ get_connection() [2/4]

◆ get_connection() [3/4]

◆ get_connection() [4/4]

◆ get_crs()

◆ get_index()

◆ get_index_const()

◆ get_index_internal()

◆ get_iterator_or_next()

◆ get_nnz()

◆ get_nnz_max_cols()

◆ get_row() [1/2]

◆ get_row() [2/2]

◆ get_values()

◆ has_connection()

◆ is_isolated()

◆ iters()

◆ mat_mult_add_row()

◆ matmul_minus()

◆ num_cols()

◆ num_connections()

◆ num_rows()

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ operator*=()

◆ operator=()

◆ p()

◆ pr()

◆ print_to_file()

◆ remove_iterator()

◆ resize_and_clear()