Implementation of sparse matrix stored in compressed column storage. More...

#include <dyncomprow.h>

Inheritance diagram for oofem::DynCompRow:

Collaboration diagram for oofem::DynCompRow:

Public Member Functions
	DynCompRow (int n)
	Constructor. More...

	DynCompRow ()
	Constructor. More...

	DynCompRow (const DynCompRow &S)
	Copy constructor. More...

DynCompRow &	operator= (const DynCompRow &C)
	Assignment operator. More...

virtual	~DynCompRow ()
	Destructor. More...

SparseMtrx *	GiveCopy () const
	Returns a newly allocated copy of receiver. More...

virtual void	times (const FloatArray &x, FloatArray &answer) const
	Evaluates $y = A \cdot x$ . More...

virtual void	timesT (const FloatArray &x, FloatArray &answer) const
	Evaluates $y = A^{\mathrm{T}} \cdot x$ . More...

virtual void	times (double x)
	Multiplies receiver by scalar value. More...

virtual int	buildInternalStructure (EngngModel *, int, const UnknownNumberingScheme &)
	Builds internal structure of receiver. More...

virtual int	assemble (const IntArray &loc, const FloatMatrix &mat)
	Assembles sparse matrix from contribution of local elements. More...

virtual int	assemble (const IntArray &rloc, const IntArray &cloc, const FloatMatrix &mat)
	Assembles sparse matrix from contribution of local elements. More...

virtual bool	canBeFactorized () const
	Determines, whether receiver can be factorized. More...

virtual void	zero ()
	Zeroes the receiver. More...

virtual const char *	giveClassName () const

virtual SparseMtrxType	giveType () const
	Sparse matrix type identification. More...

virtual bool	isAsymmetric () const
	Returns true if asymmetric. More...

virtual void	printStatistics () const
	Prints the receiver statistics (one-line) to stdout. More...

virtual double &	at (int i, int j)
	Returns coefficient at position (i,j). More...

virtual double	at (int i, int j) const
	Returns coefficient at position (i,j). More...

void	ILUPYourself (int part_fill=5, double drop_tol=1.e-8)
	Performs LU factorization on yourself; modifies receiver This routine computes the L and U factors of the ILU(p). More...

void	ILUPsolve (const FloatArray &x, FloatArray &y) const

void	ILUPtrans_solve (const FloatArray &x, FloatArray &y) const

const IntArray *	col_ind (int i) const
	Returns col index for i-th row. More...

const FloatArray *	row (int i) const
	Returns row values. More...

Public Member Functions inherited from oofem::SparseMtrx
	SparseMtrx (int n, int m)
	Constructor, creates (n,m) sparse matrix. More...

	SparseMtrx ()
	Constructor. More...

virtual	~SparseMtrx ()
	Destructor. More...

SparseMtrxVersionType	giveVersion ()
	Return receiver version. More...

void	checkBounds (int i, int j) const
	Checks size of receiver towards requested bounds. More...

int	giveNumberOfRows () const
	Returns number of rows of receiver. More...

int	giveNumberOfColumns () const
	Returns number of columns of receiver. More...

bool	isSquare () const
	Returns nonzero if receiver is square matrix. More...

bool	isNotEmpty () const
	Tests for empty matrix. More...

virtual void	times (const FloatMatrix &B, FloatMatrix &answer) const
	Evaluates $C = A^{\mathrm{T}} \cdot B$ . More...

virtual void	timesT (const FloatMatrix &B, FloatMatrix &answer) const
	Evaluates $C = A^{\mathrm{T}} \cdot B$ . More...

virtual void	add (double x, SparseMtrx &m)
	Adds x * m. More...

virtual void	addDiagonal (double x, FloatArray &m)
	Adds x * m (treats m as a diagonal matrix, stored as an array) More...

virtual int	buildInternalStructure (EngngModel *eModel, int n, int m, const IntArray &I, const IntArray &J)
	Builds internal structure of receiver based on I and J. More...

virtual int	buildInternalStructure (EngngModel *eModel, int di, const UnknownNumberingScheme &r_s, const UnknownNumberingScheme &c_s)
	Build internal structure of receiver. More...

virtual int	assembleBegin ()
	Starts assembling the elements. More...

virtual int	assembleEnd ()
	Returns when assemble is completed. More...

virtual SparseMtrx *	factorized ()
	Returns the receiver factorized. More...

virtual FloatArray *	backSubstitutionWith (FloatArray &y) const
	Computes the solution of linear system $A\cdot x = y$ where A is receiver. More...

virtual double	computeNorm () const
	Returns the norm of receiver. More...

virtual SparseMtrx *	giveSubMatrix (const IntArray &rows, const IntArray &cols)

virtual bool	isAllocatedAt (int i, int j) const
	Checks whether memory is allocated at position (i,j). More...

virtual void	toFloatMatrix (FloatMatrix &answer) const
	Converts receiving sparse matrix to a dense float matrix. More...

virtual void	printYourself () const
	Prints receiver to stdout. Works only for relatively small matrices. More...

virtual void	writeToFile (const char *fname) const
	Helpful for debugging, writes the matrix to given file. More...

std::string	errorInfo (const char *func) const
	Error printing helper. More...

FloatArray	operator* (const FloatArray &x) const
	IML compatibility, $A \cdot x$ . More...

FloatArray	trans_mult (const FloatArray &x) const
	IML compatibility, $A^{\mathrm{T}} \cdot x$ . More...

Protected Member Functions
double	operator() (int i, int j) const
	implements 0-based access More...

double &	operator() (int i, int j)
	implements 0-based access More...

int	giveColIndx (int row, int col) const
	returns the column index of given column at given row, else returns zero. More...

int	insertColInRow (int row, int col)
	insert column entry into row, preserving order of column indexes, returns the index of new row. More...

void	checkSizeTowards (IntArray &)

void	checkSizeTowards (const IntArray &rloc, const IntArray &cloc)

void	growTo (int)

void	qsortRow (IntArray &ind, IntArray &ir, FloatArray &val, int l, int r)

int	qsortRowPartition (IntArray &ind, IntArray &ir, FloatArray &val, int l, int r)

Protected Attributes
FloatArray **	rows_

IntArray **	colind_

IntArray	diag_rowptr_

int	base_

Protected Attributes inherited from oofem::SparseMtrx
int	nRows
	Number of rows. More...

int	nColumns
	Number of columns. More...

SparseMtrxVersionType	version
	Allows to track if receiver changes. More...

Additional Inherited Members
Public Types inherited from oofem::SparseMtrx
typedef long	SparseMtrxVersionType

Detailed Description

Implementation of sparse matrix stored in compressed column storage.

Designed to allow simple dynamic runtime grow of receiver.

Definition at line 57 of file dyncomprow.h.

Constructor & Destructor Documentation

oofem::DynCompRow::DynCompRow ( int n )

Constructor.

Before any operation an internal profile must be built.

See also: builInternalStructure

Definition at line 63 of file dyncomprow.C.

References colind_, oofem::SparseMtrx::nColumns, oofem::SparseMtrx::nRows, and rows_.

oofem::DynCompRow::DynCompRow ( void )

Constructor.

Before any operation an internal profile must be built.

See also: builInternalStructure

Definition at line 56 of file dyncomprow.C.

References colind_, and rows_.

Referenced by GiveCopy().

oofem::DynCompRow::DynCompRow ( const DynCompRow & S )

Copy constructor.

Definition at line 76 of file dyncomprow.C.

References colind_, oofem::SparseMtrx::nColumns, oofem::SparseMtrx::nRows, rows_, and oofem::SparseMtrx::version.

oofem::DynCompRow::~DynCompRow ( )

virtual

Destructor.

Definition at line 104 of file dyncomprow.C.

References colind_, oofem::SparseMtrx::nRows, and rows_.

Member Function Documentation

int oofem::DynCompRow::assemble	(	const IntArray &	loc,
		const FloatMatrix &	mat
	)

virtual

Assembles sparse matrix from contribution of local elements.

This method for each element adds its contribution to itself. Mapping between local element contribution and its global position is given by local code numbers of element.

Parameters

loc	Location array. The values corresponding to zero loc array value are not assembled.
mat	Contribution to be assembled using loc array.

Returns: Zero iff successful.

Implements oofem::SparseMtrx.

Definition at line 321 of file dyncomprow.C.

References at(), oofem::IntArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::giveNumberOfRows(), oofem::IntArray::giveSize(), OOFEM_ERROR, and oofem::SparseMtrx::version.

int oofem::DynCompRow::assemble	(	const IntArray &	rloc,
		const IntArray &	cloc,
		const FloatMatrix &	mat
	)

virtual

Assembles sparse matrix from contribution of local elements.

This method for each element adds its contribution to itself. Mapping between local element contribution and its global position is given by row and column local code numbers.

Parameters

rloc	Row location array. The values corresponding to zero loc array value are not assembled.
cloc	Column location array. The values corresponding to zero loc array value are not assembled.
mat	Contribution to be assembled using rloc and cloc arrays. The rloc position determines the row, the cloc position determines the corresponding column.

Returns: Zero iff successful.

Implements oofem::SparseMtrx.

Definition at line 352 of file dyncomprow.C.

References oofem::FloatArray::at(), oofem::IntArray::giveSize(), insertColInRow(), rows_, and oofem::SparseMtrx::version.

double & oofem::DynCompRow::at	(	int	i,
		int	j
	)

virtual

Returns coefficient at position (i,j).

Implements oofem::SparseMtrx.

Definition at line 406 of file dyncomprow.C.

References oofem::FloatArray::at(), giveColIndx(), OOFEM_ERROR, rows_, and oofem::SparseMtrx::version.

Referenced by assemble(), giveColIndx(), ILUPsolve(), ILUPtrans_solve(), ILUPYourself(), insertColInRow(), times(), and timesT().

double oofem::DynCompRow::at	(	int	i,
		int	j
	)		const

virtual

Returns coefficient at position (i,j).

Implements oofem::SparseMtrx.

Definition at line 421 of file dyncomprow.C.

References oofem::FloatArray::at(), giveColIndx(), oofem::SparseMtrx::nColumns, oofem::SparseMtrx::nRows, OOFEM_ERROR, and rows_.

int oofem::DynCompRow::buildInternalStructure	(	EngngModel *	eModel,
		int	di,
		const UnknownNumberingScheme &	s
	)

virtual

Builds internal structure of receiver.

This method determines the internal profile of sparse matrix, allocates necessary space for storing nonzero coefficients and initializes receiver. In general, the profile of sparse matrix is determined using one (or more) loop over local code numbers of elements. This method must be called before any operation, like assembly, zeroing, or multiplication.

Parameters