#include <spmh.h>

Inheritance diagram for CompactSparseMatrixHandler_tpl< off >:

Collaboration diagram for CompactSparseMatrixHandler_tpl< off >:

Classes
class	const_iterator

Public Member Functions
	CompactSparseMatrixHandler_tpl (const integer &n, const integer &nn, std::vector< doublereal > &x, const std::vector< integer > &i, const std::vector< integer > &p)

virtual	~CompactSparseMatrixHandler_tpl (void)

CompactSparseMatrixHandler_tpl < off >::const_iterator	begin (void) const

const CompactSparseMatrixHandler_tpl < off >::const_iterator &	end (void) const

Public Member Functions inherited from CompactSparseMatrixHandler
	CompactSparseMatrixHandler (const integer &n, const integer &nn, std::vector< doublereal > &x, const std::vector< integer > &i, const std::vector< integer > &p)

virtual	~CompactSparseMatrixHandler ()

virtual CompactSparseMatrixHandler *	Copy (void) const =0

void	AddUnchecked (const CompactSparseMatrixHandler &m)

virtual const doublereal *	pdGetMat (void) const

void	Reset (void)

virtual integer	MakeCompressedColumnForm (doublereal const Ax, integer const Ai, integer *const Ap, int offset=0) const

virtual integer	MakeCompressedColumnForm (std::vector< doublereal > &Ax, std::vector< integer > &Ai, std::vector< integer > &Ap, int offset=0) const

virtual integer	MakeIndexForm (doublereal const Ax, integer const Arow, integer const Acol, integer const AcolSt, int offset=0) const

virtual integer	MakeIndexForm (std::vector< doublereal > &Ax, std::vector< integer > &Arow, std::vector< integer > &Acol, std::vector< integer > &AcolSt, int offset=0) const

Public Member Functions inherited from SparseMatrixHandler
const integer	Nz () const

	SparseMatrixHandler (const integer &n, const integer &nn=0)

virtual	~SparseMatrixHandler (void)

integer	iGetNumRows (void) const

integer	iGetNumCols (void) const

virtual VectorHandler &	GetCol (integer icol, VectorHandler &out) const =0

Public Member Functions inherited from MatrixHandler
virtual	~MatrixHandler (void)

virtual void	Resize (integer, integer)=0

virtual void	ResizeReset (integer, integer)

virtual doublereal *	pdGetMat (void)

virtual integer *	piGetRows (void) const

virtual integer *	piGetCols (void) const

virtual integer	PacMat (void)

virtual void	PutCoef (integer iRow, integer iCol, const doublereal &dCoef)

virtual void	IncCoef (integer iRow, integer iCol, const doublereal &dCoef)

virtual void	DecCoef (integer iRow, integer iCol, const doublereal &dCoef)

virtual const doublereal &	dGetCoef (integer iRow, integer iCol) const

virtual const doublereal &	operator() (integer iRow, integer iCol) const =0

virtual doublereal &	operator() (integer iRow, integer iCol)=0

virtual MatrixHandler &	operator= (const MatrixHandler &MH)

virtual MatrixHandler &	operator+= (const SubMatrixHandler &SubMH)

virtual MatrixHandler &	operator-= (const SubMatrixHandler &SubMH)

virtual MatrixHandler &	operator+= (const VariableSubMatrixHandler &SubMH)

virtual MatrixHandler &	operator-= (const VariableSubMatrixHandler &SubMH)

virtual MatrixHandler &	ScalarMul (const doublereal &d)

virtual MatrixHandler &	MatMatMul (MatrixHandler &out, const MatrixHandler &in) const

virtual MatrixHandler &	MatTMatMul (MatrixHandler &out, const MatrixHandler &in) const

virtual MatrixHandler &	MatMatIncMul (MatrixHandler &out, const MatrixHandler &in) const

virtual MatrixHandler &	MatTMatIncMul (MatrixHandler &out, const MatrixHandler &in) const

virtual MatrixHandler &	MatMatDecMul (MatrixHandler &out, const MatrixHandler &in) const

virtual MatrixHandler &	MatTMatDecMul (MatrixHandler &out, const MatrixHandler &in) const

virtual VectorHandler &	MatVecMul (VectorHandler &out, const VectorHandler &in) const

virtual VectorHandler &	MatTVecMul (VectorHandler &out, const VectorHandler &in) const

virtual VectorHandler &	MatVecIncMul (VectorHandler &out, const VectorHandler &in) const

virtual VectorHandler &	MatTVecIncMul (VectorHandler &out, const VectorHandler &in) const

virtual VectorHandler &	MatVecDecMul (VectorHandler &out, const VectorHandler &in) const

virtual VectorHandler &	MatTVecDecMul (VectorHandler &out, const VectorHandler &in) const

virtual doublereal	ConditionNumber (enum Norm_t eNorm=NORM_1) const

virtual doublereal	Norm (enum Norm_t eNorm=NORM_1) const

Protected Member Functions
MatrixHandler &	MatMatMul_base (void(MatrixHandler::*op)(integer iRow, integer iCol, const doublereal &dCoef), MatrixHandler &out, const MatrixHandler &in) const

MatrixHandler &	MatTMatMul_base (void(MatrixHandler::*op)(integer iRow, integer iCol, const doublereal &dCoef), MatrixHandler &out, const MatrixHandler &in) const

virtual VectorHandler &	MatVecMul_base (void(VectorHandler::*op)(integer iRow, const doublereal &dCoef), VectorHandler &out, const VectorHandler &in) const

virtual VectorHandler &	MatTVecMul_base (void(VectorHandler::*op)(integer iRow, const doublereal &dCoef), VectorHandler &out, const VectorHandler &in) const

Protected Attributes
CompactSparseMatrixHandler_tpl < off >::const_iterator	m_end

Protected Attributes inherited from CompactSparseMatrixHandler
bool	bMatDuplicate

std::vector< doublereal > &	Ax

const std::vector< integer > &	Ai

const std::vector< integer > &	Ap

Protected Attributes inherited from SparseMatrixHandler
integer	NRows

integer	NCols

integer	NZ

Additional Inherited Members
Public Types inherited from MatrixHandler
enum	Norm_t { NORM_1, NORM_INF }

Detailed Description

template<int off>
class CompactSparseMatrixHandler_tpl< off >

Definition at line 210 of file spmh.h.

Constructor & Destructor Documentation

template<int off>

CompactSparseMatrixHandler_tpl< off >::CompactSparseMatrixHandler_tpl	(	const integer &	n,
		const integer &	nn,
		std::vector< doublereal > &	x,
		const std::vector< integer > &	i,
		const std::vector< integer > &	p
	)

Definition at line 148 of file spmh.cc.

References NO_OP.

 : CompactSparseMatrixHandler(n, nn, x, i, p),
 m_end(*this, true)
 {
         NO_OP;
 }

template<int off>

CompactSparseMatrixHandler_tpl< off >::~CompactSparseMatrixHandler_tpl ( void )

virtual

Definition at line 161 of file spmh.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

template<int off>

CompactSparseMatrixHandler_tpl<off>::const_iterator CompactSparseMatrixHandler_tpl< off >::begin ( void ) const

inline

Definition at line 262 of file spmh.h.

Referenced by main().

                                                                             {
                 return CompactSparseMatrixHandler_tpl<off>::const_iterator(*this);
         };

template<int off>

const CompactSparseMatrixHandler_tpl<off>::const_iterator& CompactSparseMatrixHandler_tpl< off >::end ( void ) const

inline

Definition at line 266 of file spmh.h.

References CompactSparseMatrixHandler_tpl< off >::m_end.

Referenced by main().

                                                                                  {
                 return m_end;
         };

template<int off>

MatrixHandler & CompactSparseMatrixHandler_tpl< off >::MatMatMul_base	(	void(MatrixHandler::*)(integer iRow, integer iCol, const doublereal &dCoef)	op,
		MatrixHandler &	out,
		const MatrixHandler &	in
	)		const

protectedvirtual

Reimplemented from MatrixHandler.

Definition at line 169 of file spmh.cc.

References ASSERT, MatrixHandler::iGetNumCols(), MatrixHandler::iGetNumRows(), SparseMatrixHandler::NCols, and SparseMatrixHandler::NRows.

 {
         ASSERT(in.iGetNumRows() == NCols);
         ASSERT(out.iGetNumRows() == NRows);
         ASSERT(in.iGetNumCols() == out.iGetNumCols());
 
         // NOTE: out must be zeroed by caller
 
         integer ncols_in = in.iGetNumCols();
 
         for (integer col_idx = 1; col_idx <= NCols; col_idx++) {
                 integer re = Ap[col_idx] - off;
                 integer ri = Ap[col_idx - 1] - off;
                 for ( ; ri < re; ri++) {
                         integer row_idx = Ai[ri] - off + 1;
                         const doublereal& d = Ax[ri];
                         for (integer col_in = 1; col_in <= ncols_in; col_in++) {
                                 (out.*op)(row_idx, col_in, d*in(col_idx, col_in));
                         }
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

template<int off>

MatrixHandler & CompactSparseMatrixHandler_tpl< off >::MatTMatMul_base	(	void(MatrixHandler::*)(integer iRow, integer iCol, const doublereal &dCoef)	op,
		MatrixHandler &	out,
		const MatrixHandler &	in
	)		const

protectedvirtual

Reimplemented from MatrixHandler.

Definition at line 199 of file spmh.cc.

References ASSERT, MatrixHandler::iGetNumCols(), MatrixHandler::iGetNumRows(), SparseMatrixHandler::NCols, and SparseMatrixHandler::NRows.

 {
         ASSERT(in.iGetNumRows() == NRows);
         ASSERT(out.iGetNumRows() == NCols);
         ASSERT(in.iGetNumCols() == out.iGetNumCols());
         // NOTE: out must be zeroed by caller
 
         integer ncols_in = in.iGetNumCols();
 
         for (integer col_idx = 1; col_idx <= NCols; col_idx++) {
                 integer re = Ap[col_idx] - off;
                 integer ri = Ap[col_idx - 1] - off;
                 for ( ; ri < re; ri++) {
                         integer row_idx = Ai[ri] - off + 1;
                         const doublereal& d = Ax[ri];
                         for (integer col_in = 1; col_in <= ncols_in; col_in++) {
                                 (out.*op)(col_idx, col_in, d*in(row_idx, col_in));
                         }
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

template<int off>

VectorHandler & CompactSparseMatrixHandler_tpl< off >::MatTVecMul_base	(	void(VectorHandler::*)(integer iRow, const doublereal &dCoef)	op,
		VectorHandler &	out,
		const VectorHandler &	in
	)		const

protectedvirtual

Reimplemented from MatrixHandler.

Definition at line 253 of file spmh.cc.

References ASSERT, VectorHandler::iGetSize(), SparseMatrixHandler::NCols, and SparseMatrixHandler::NRows.

 {
         ASSERT(in.iGetSize() == NRows);
         ASSERT(out.iGetSize() == NCols);
 
         // NOTE: out must be zeroed by caller
 
         for (integer col_idx = 1; col_idx <= NCols; col_idx++) {
                 integer re = Ap[col_idx] - off;
                 integer ri = Ap[col_idx - 1] - off;
                 for ( ; ri < re; ri++) {
                         integer row_idx = Ai[ri] - off + 1;
                         const doublereal& d = Ax[ri];
                         (out.*op)(col_idx, d*in(row_idx));
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

template<int off>

VectorHandler & CompactSparseMatrixHandler_tpl< off >::MatVecMul_base	(	void(VectorHandler::*)(integer iRow, const doublereal &dCoef)	op,
		VectorHandler &	out,
		const VectorHandler &	in
	)		const

protectedvirtual

Reimplemented from MatrixHandler.

Definition at line 228 of file spmh.cc.

References ASSERT, VectorHandler::iGetSize(), SparseMatrixHandler::NCols, and SparseMatrixHandler::NRows.

 {
         ASSERT(in.iGetSize() == NCols);
         ASSERT(out.iGetSize() == NRows);
 
         // NOTE: out must be zeroed by caller
 
         for (integer col_idx = 1; col_idx <= NCols; col_idx++) {
                 integer re = Ap[col_idx] - off;
                 integer ri = Ap[col_idx - 1] - off;
                 for ( ; ri < re; ri++) {
                         integer row_idx = Ai[ri] - off + 1;
                         const doublereal& d = Ax[ri];
                         (out.*op)(row_idx, d*in(col_idx));
                 }
         }
 
         return out;
 }