The compressed row storage sparse matrix data structure. More...

#include <CRS.hpp>

Inheritance diagram for CRS< T >:

Collaboration diagram for CRS< T >:

Public Member Functions
	CRS ()
	Base constructor. More...

	CRS (std::string file, T zero=0)
	Base constructor. More...

	CRS (const ULI number_of_nonzeros, const ULI number_of_rows, const ULI number_of_cols, T zero)
	Base constructor which only initialises the internal arrays. More...

	CRS (CRS< T > &toCopy)
	Copy constructor. More...

	CRS (std::vector< Triplet< T > > input, ULI m, ULI n, T zero)
	Constructor which transforms a collection of input triplets to CRS format. More...

virtual void	load (std::vector< Triplet< T > > &input, ULI m, ULI n, T zero)

T &	random_access (ULI i, ULI j)
	Method which provides random matrix access to the stored sparse matrix. More...

virtual void	getFirstIndexPair (ULI &row, ULI &col)
	Returns the first nonzero index, per reference. More...

virtual void	zxa (const T __restrict__ x, T __restrict__ z)
	In-place z=xA function.

virtual void	zax (const T __restrict__ x, T __restrict__ z)
	In-place z=Ax function. More...

template<size_t k>
void	ZaX (const T __restrict__ const __restrict__ const X, T __restrict__ const __restrict__ const Z)

template<size_t k>
void	ZXa (const T __restrict__ const __restrict__ const X, T __restrict__ const __restrict__ const Z)

ULI *	rowJump ()
	Returns pointer to the row_start vector. More...

ULI *	columnIndices ()
	Returns pointer to the column index vector. More...

T *	values ()
	Returns pointer to the matrix nonzeros vector. More...

virtual	~CRS ()
	Base deconstructor. More...

virtual size_t	bytesUsed ()

Public Member Functions inherited from SparseMatrix< T, ULI >
	SparseMatrix ()
	Base constructor. More...

	SparseMatrix (const ULInzs, const ULInr, const ULInc, const T zero)
	Base constructor. More...

virtual	~SparseMatrix ()
	Base deconstructor. More...

void	loadFromFile (const std::string file, const T zero=0)
	Function which loads a matrix from a matrix market file. More...

virtual unsigned long int	m ()
	Queries the number of rows this matrix contains. More...

virtual unsigned long int	n ()
	Queries the number of columns this matrix contains. More...

virtual unsigned long int	nzs ()
	Queries the number of nonzeroes stored in this matrix. More...

virtual T *	mv (const T *x)
	Calculates and returns z=Ax. More...

Public Member Functions inherited from Matrix< T >
	Matrix ()
	Base constructor. More...

virtual	~Matrix ()
	Base deconstructor. More...

virtual void	zax (const T __restrict__ x, T __restrict__ z, const size_t k, const clockid_t clock_id=0, double *elapsed_time=NULL)
	Wrapper function to call the zax kernel multiple times successively, while timing the duration of the operation. More...

template<size_t k>
void	ZaX (const T __restrict__ const __restrict__ const X, T __restrict__ const __restrict__ const Z)
	In-place Z=AX function, where A is m x n, Z = m x k, and X is n x k. More...

template<size_t k>
void	ZXa (const T __restrict__ const __restrict__ const X, T __restrict__ const __restrict__ const Z)
	In-place Z=XA function, where A is m x n, Z = k x n, and X is k x m. More...

virtual void	zxa (const T __restrict__ x, T __restrict__ z, const unsigned long int repeat, const clockid_t clock_id=0, double *elapsed_time=NULL)
	Wrapper function to call the zxa kernel multiple times successively, while timing the operation duration. More...

Protected Member Functions
bool	find (const ULI col_index, const ULI search_start, const ULI search_end, ULI &ret)
	Helper function which finds a value with a given column index on a given subrange of indices. More...

Static Protected Member Functions
static int	compareTriplets (const void left, const void right)
	Sorts 1D columnwise.

Protected Attributes
ULI *	row_start
	Array keeping track of individual row starting indices. More...

T *	ds
	Array containing the actual nnz non-zeros. More...

ULI *	col_ind
	Array containing the column indeces corresponding to the elements in ds. More...

Protected Attributes inherited from SparseMatrix< T, ULI >
ULI	nor
	Number of rows. More...

ULI	noc
	Number of columns.

ULI	nnz
	Number of non-zeros. More...

Additional Inherited Members
Public Attributes inherited from SparseMatrix< T, ULI >
T	zero_element
	The element considered to be zero. More...

Detailed Description

template<typename T>
class CRS< T >

The compressed row storage sparse matrix data structure.

Constructor & Destructor Documentation

template<typename T>

CRS< T >::CRS ( )

inline

Base constructor.

template<typename T>

CRS< T >::CRS	(	std::string	file,
		T	zero = `0`
	)

inline

Base constructor.

Will read in from Matrix Market file.

See also: SparseMatrix::SparseMatrix( file, zero )

References SparseMatrix< T, ULI >::loadFromFile().

template<typename T>

CRS< T >::CRS	(	const ULI	number_of_nonzeros,
		const ULI	number_of_rows,
		const ULI	number_of_cols,
		T	zero
	)

inline

Base constructor which only initialises the internal arrays.

Note that to gain a valid CRS structure, these arrays have to be filled by some external mechanism (i.e., after calling this constructor, the internal arrays contain garbage, resuling in invalid datastructure).

Parameters

number_of_nonzeros	The number of non-zeros to be stored.
number_of_rows	The number of rows to be stored.
number_of_cols	The number of columns of the matrix.
zero	The element considered to be zero.

References CRS< T >::col_ind, CRS< T >::ds, SparseMatrix< T, ULI >::nnz, SparseMatrix< T, ULI >::noc, SparseMatrix< T, ULI >::nor, CRS< T >::row_start, and SparseMatrix< T, ULI >::zero_element.

template<typename T>

CRS< T >::CRS ( CRS< T > & toCopy )

inline

Copy constructor.

Parameters

toCopy reference to the CRS datastructure to copy.

References CRS< T >::col_ind, CRS< T >::ds, SparseMatrix< T, ULI >::nnz, SparseMatrix< T, ULI >::nor, CRS< T >::row_start, and SparseMatrix< T, ULI >::zero_element.

template<typename T>

CRS< T >::CRS	(	std::vector< Triplet< T > >	input,
		ULI	m,
		ULI	n,
		T	zero
	)

inline

Constructor which transforms a collection of input triplets to CRS format.

The input collection is considered to have at most one triplet with unique pairs of indeces. Unspecified behaviour occurs when this assumption is not met.

Parameters

input	The input collection.
m	Total number of rows.
n	Total number of columns.
zero	The element considered to be zero.

References CRS< T >::load().

template<typename T>

virtual CRS< T >::~CRS ( )

inlinevirtual

Base deconstructor.

References CRS< T >::col_ind, CRS< T >::ds, and CRS< T >::row_start.

Member Function Documentation

template<typename T>

virtual size_t CRS< T >::bytesUsed ( )

inlinevirtual

See also: Matrix::bytesUsed()

Implements Matrix< T >.

References SparseMatrix< T, ULI >::nnz, and SparseMatrix< T, ULI >::nor.

template<typename T>

ULI* CRS< T >::columnIndices ( )

inline

Returns pointer to the column index vector.

References CRS< T >::col_ind.

template<typename T>

bool CRS< T >::find	(	const ULI	col_index,
		const ULI	search_start,
		const ULI	search_end,
		ULI &	ret
	)

inlineprotected

Helper function which finds a value with a given column index on a given subrange of indices.

Parameters

col_index	The given column index.
search_start	The start index of the subrange (inclusive).
search_end	The end index of the subrange (exlusive).
ret	Reference to the variable where the return index is stored.