Generic matrix. More...

#include <GeneralMatrix.h>

Inheritance diagram for GeneralMatrix:

Detailed Description

Generic matrix.

Definition at line 22 of file GeneralMatrix.h.

Public Member Functions
	GeneralMatrix ()=default
	Base Constructor.

virtual void	zero ()=0
	Zero the matrix elements.

virtual void	mult (const double b, double prod) const =0
	Multiply A*b and write result to prod.

virtual void	leftMult (const double const b, double const prod) const =0
	Multiply b*A and write result to prod.

virtual int	factor ()=0
	Factors the A matrix, overwriting A.

virtual int	factorQR ()
	Factors the A matrix using the QR algorithm, overwriting A.

virtual double	rcondQR ()
	Returns an estimate of the inverse of the condition number for the matrix.

virtual double	rcond (double a1norm)=0
	Returns an estimate of the inverse of the condition number for the matrix.

virtual void	useFactorAlgorithm (int fAlgorithm)
	Change the way the matrix is factored.

virtual int	factorAlgorithm () const =0
	Return the factor algorithm used.

virtual double	oneNorm () const =0
	Calculate the one norm of the matrix.

virtual size_t	nRows () const =0
	Return the number of rows in the matrix.

virtual void	clearFactorFlag ()
	clear the factored flag

virtual int	solve (double *b, size_t nrhs=1, size_t ldb=0)=0
	Solves the Ax = b system returning x in the b spot.

virtual bool	factored () const
	true if the current factorization is up to date with the matrix

virtual double *	ptrColumn (size_t j)=0
	Return a pointer to the top of column j, columns are assumed to be contiguous in memory.

virtual double &	operator() (size_t i, size_t j)=0
	Index into the (i,j) element.

virtual double	operator() (size_t i, size_t j) const =0
	Constant Index into the (i,j) element.

virtual double const	colPts ()=0
	Return a vector of const pointers to the columns.

virtual size_t	checkRows (double &valueSmall) const =0
	Check to see if we have any zero rows in the Jacobian.

virtual size_t	checkColumns (double &valueSmall) const =0
	Check to see if we have any zero columns in the Jacobian.

Protected Attributes
int	m_factored = false
	Indicates whether the matrix is factored.

Constructor & Destructor Documentation

◆ GeneralMatrix()

GeneralMatrix ( )

default

Base Constructor.

Member Function Documentation

◆ zero()

virtual void zero ( )

pure virtual

Zero the matrix elements.

Implemented in BandMatrix.

◆ mult()

virtual void mult	(	const double *	b,
		double *	prod
	)		const

pure virtual

Multiply A*b and write result to prod.

Parameters

b	Vector to do the rh multiplication
prod	OUTPUT vector to receive the result

Implemented in BandMatrix.

◆ leftMult()

virtual void leftMult	(	const double *const	b,
		double *const	prod
	)		const

pure virtual

Multiply b*A and write result to prod.

Parameters

b	Vector to do the lh multiplication
prod	OUTPUT vector to receive the result

Implemented in BandMatrix.

◆ factor()

virtual int factor ( )

pure virtual

Factors the A matrix, overwriting A.

We flip m_factored boolean to indicate that the matrix is now A-1.

Implemented in BandMatrix.

◆ factorQR()

virtual int factorQR ( )

inlinevirtual

Factors the A matrix using the QR algorithm, overwriting A.

we set m_factored to 2 to indicate the matrix is now QR factored

Returns: the info variable from LAPACK

Definition at line 59 of file GeneralMatrix.h.

◆ rcondQR()

virtual double rcondQR ( )

inlinevirtual

Returns an estimate of the inverse of the condition number for the matrix.

The matrix must have been previously factored using the QR algorithm

Returns: the inverse of the condition number

Definition at line 69 of file GeneralMatrix.h.

◆ rcond()

virtual double rcond ( double a1norm )

pure virtual

Returns an estimate of the inverse of the condition number for the matrix.

The matrix must have been previously factored using the LU algorithm

Parameters

a1norm Norm of the matrix

Returns: the inverse of the condition number

Implemented in BandMatrix.

◆ useFactorAlgorithm()

virtual void useFactorAlgorithm ( int fAlgorithm )

inlinevirtual

Change the way the matrix is factored.

Parameters

fAlgorithm integer 0 LU factorization 1 QR factorization

Definition at line 88 of file GeneralMatrix.h.

◆ factorAlgorithm()

virtual int factorAlgorithm ( ) const

pure virtual

Return the factor algorithm used.

Implemented in BandMatrix.

◆ oneNorm()

virtual double oneNorm ( ) const

pure virtual

Calculate the one norm of the matrix.

Implemented in BandMatrix.

◆ nRows()

virtual size_t nRows ( ) const

pure virtual

Return the number of rows in the matrix.

Implemented in BandMatrix.

◆ clearFactorFlag()

virtual void clearFactorFlag ( )

inlinevirtual

clear the factored flag

Definition at line 102 of file GeneralMatrix.h.

◆ solve()

virtual int solve	(	double *	b,
		size_t	nrhs = `1`,
		size_t	ldb = `0`
	)

pure virtual

Solves the Ax = b system returning x in the b spot.

Parameters

b	Vector for the RHS of the equation system
nrhs	Number of right-hand sides to solve, default 1
ldb	Leading dimension of the right-hand side array. Defaults to nRows()

Implemented in BandMatrix.

◆ factored()

virtual bool factored ( ) const

inlinevirtual

true if the current factorization is up to date with the matrix

Definition at line 116 of file GeneralMatrix.h.

◆ ptrColumn()

virtual double * ptrColumn ( size_t j )

pure virtual

Return a pointer to the top of column j, columns are assumed to be contiguous in memory.

Parameters

j	Value of the column

Returns: a pointer to the top of the column

Implemented in BandMatrix.

◆ operator()() [1/2]

virtual double & operator()	(	size_t	i,
		size_t	j
	)

pure virtual

Index into the (i,j) element.

Parameters

i	row
j	column

Returns: a changeable reference to the matrix entry

Implemented in BandMatrix.

◆ operator()() [2/2]

virtual double operator()	(	size_t	i,
		size_t	j
	)		const

pure virtual

Constant Index into the (i,j) element.

Parameters

i	row
j	column

Returns: an unchangeable reference to the matrix entry

Implemented in BandMatrix.

◆ colPts()

virtual double *const * colPts ( )

pure virtual

Return a vector of const pointers to the columns.

Note the value of the pointers are protected by their being const. However, the value of the matrix is open to being changed.

Returns: a vector of pointers to the top of the columns of the matrices.

Implemented in BandMatrix.

◆ checkRows()

virtual size_t checkRows ( double & valueSmall ) const

pure virtual

Check to see if we have any zero rows in the Jacobian.

This utility routine checks to see if any rows are zero. The smallest row is returned along with the largest coefficient in that row

Parameters

valueSmall OUTPUT value of the largest coefficient in the smallest row

Returns: index of the row that is most nearly zero

Implemented in BandMatrix.

◆ checkColumns()

virtual size_t checkColumns ( double & valueSmall ) const

pure virtual

Check to see if we have any zero columns in the Jacobian.

This utility routine checks to see if any columns are zero. The smallest column is returned along with the largest coefficient in that column

Parameters

valueSmall OUTPUT value of the largest coefficient in the smallest column

Returns: index of the column that is most nearly zero

Implemented in BandMatrix.

Member Data Documentation

◆ m_factored

int m_factored = false

protected

Indicates whether the matrix is factored.

0 for unfactored; Non-zero values indicate a particular factorization (LU=1, QR=2).

Definition at line 176 of file GeneralMatrix.h.

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

GeneralMatrix.h

Detailed Description

Public Member Functions

Protected Attributes

Constructor & Destructor Documentation

◆ GeneralMatrix()

Member Function Documentation

◆ zero()

◆ mult()

◆ leftMult()

◆ factor()

◆ factorQR()

◆ rcondQR()

◆ rcond()

◆ useFactorAlgorithm()

◆ factorAlgorithm()

◆ oneNorm()

◆ nRows()

◆ clearFactorFlag()

◆ solve()

◆ factored()

◆ ptrColumn()

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ colPts()

◆ checkRows()

◆ checkColumns()

Member Data Documentation

◆ m_factored