#include <dirccmh.h>

Inheritance diagram for DirCColMatrixHandler< off >:

Collaboration diagram for DirCColMatrixHandler< off >:

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

virtual	~DirCColMatrixHandler ()

CompactSparseMatrixHandler *	Copy (void) const

doublereal &	operator() (integer i_row, integer i_col)

const doublereal &	operator() (integer i_row, integer i_col) const

void	Resize (integer n, integer nn)

VectorHandler &	GetCol (integer icol, VectorHandler &out) const

MatrixHandler &	MulAndSumWithShift (MatrixHandler &out, doublereal s=1., integer drow=0, integer dcol=0) const

MatrixHandler &	FakeThirdOrderMulAndSumWithShift (MatrixHandler &out, std::vector< bool > b, doublereal s=1., integer drow=0, integer dcol=0) const

Public Member Functions inherited from 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)

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 ()

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

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

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 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

Private Member Functions
	DirCColMatrixHandler (const DirCColMatrixHandler &)

Private Attributes
std::vector< integer * >	pindices

std::vector< integer >	indices

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

Protected Member Functions inherited from CompactSparseMatrixHandler_tpl< off >
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 inherited from CompactSparseMatrixHandler_tpl< off >
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

Detailed Description

template<int off>
class DirCColMatrixHandler< off >

Definition at line 46 of file dirccmh.h.

Constructor & Destructor Documentation

template<int off>

DirCColMatrixHandler< off >::DirCColMatrixHandler ( const DirCColMatrixHandler< off > & )

private

template<int off>

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

Definition at line 41 of file dirccmh.cc.

References SparseMatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), DirCColMatrixHandler< off >::indices, and DirCColMatrixHandler< off >::pindices.

 : CompactSparseMatrixHandler_tpl<off>(p.size() - 1, p.size() - 1, x, i, p),
   pindices(SparseMatrixHandler::iGetNumCols() + 1),
 indices(SparseMatrixHandler::iGetNumRows()*SparseMatrixHandler::iGetNumCols(), -1)
 {
         for (integer col = 1; col <= SparseMatrixHandler::iGetNumCols(); col++) {
                 pindices[col] = &indices[(col - 1)*SparseMatrixHandler::iGetNumRows()] - 1;
 
                 integer row_begin = p[col - 1] - off, row_end = p[col] - off;
 
                 for (integer r = row_begin; r < row_end; r++) {
                         pindices[col][i[r] - off + 1] = r;
                 }
         }
 }

Here is the call graph for this function:

template<int off>

DirCColMatrixHandler< off >::~DirCColMatrixHandler ( )

virtual

Definition at line 60 of file dirccmh.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

template<int off>

CompactSparseMatrixHandler * DirCColMatrixHandler< off >::Copy ( void ) const

virtual

Implements CompactSparseMatrixHandler.

Definition at line 69 of file dirccmh.cc.

References CompactSparseMatrixHandler::bMatDuplicate.

 {
         std::vector<doublereal> *pax =
                 new std::vector<doublereal>(CompactSparseMatrixHandler_tpl<off>::Ax);
         DirCColMatrixHandler<off> *p =
                 new DirCColMatrixHandler<off>(*pax, CompactSparseMatrixHandler_tpl<off>::Ai,
                         CompactSparseMatrixHandler_tpl<off>::Ap);
         p->bMatDuplicate = true;
 
         return p;
 }

template<int off>

MatrixHandler & DirCColMatrixHandler< off >::FakeThirdOrderMulAndSumWithShift	(	MatrixHandler &	out,
		std::vector< bool >	b,
		doublereal	s = `1.`,
		integer	drow = `0`,
		integer	dcol = `0`
	)		const

Definition at line 145 of file dirccmh.cc.

References SparseMatrixHandler::iGetNumCols(), MatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), MatrixHandler::iGetNumRows(), MatrixHandler::IncCoef(), and MBDYN_EXCEPT_ARGS.

 {
         silent_cerr("DirCColMatrixHandler<off>::FakeThirdOrderMulAndSumWithShift "
                         "called" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);            
         if ((out.iGetNumCols() < SparseMatrixHandler::iGetNumCols() + dcol)
                         || (out.iGetNumRows() < SparseMatrixHandler::iGetNumRows() + drow))
         {
                 silent_cerr("Assertion fault "
                                 "in DirCColMatrixHandler<off>::MulAndSumWithShift"
                                 << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
         drow = drow + 1;
         for (integer col = 0; col < SparseMatrixHandler::iGetNumCols(); col++) {
                 integer idx = CompactSparseMatrixHandler_tpl<off>::Ap[col] - off;
                 integer idxe = CompactSparseMatrixHandler_tpl<off>::Ap[col + 1] - off;
                 integer newcol = col + dcol + 1;
                 for (; idx < idxe; idx++) {
                         if (b[CompactSparseMatrixHandler_tpl<off>::Ai[idx] - off]) {
                                 out.IncCoef(CompactSparseMatrixHandler_tpl<off>::Ai[idx] - off + drow,
                                         newcol, CompactSparseMatrixHandler_tpl<off>::Ax[idx]*s);
                         }
                 }
         }
         return out;     
 }

Here is the call graph for this function:

template<int off>

VectorHandler & DirCColMatrixHandler< off >::GetCol	(	integer	icol,
		VectorHandler &	out
	)		const

virtual

Implements SparseMatrixHandler.

Definition at line 92 of file dirccmh.cc.

References SparseMatrixHandler::iGetNumCols(), and MBDYN_EXCEPT_ARGS.

 {
         /*
          * Note: we assume out has been reset
          */
         
         if (icol > SparseMatrixHandler::iGetNumCols()) {
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         integer idx = CompactSparseMatrixHandler_tpl<off>::Ap[icol - 1] - off;
         integer idxe = CompactSparseMatrixHandler_tpl<off>::Ap[icol] - off;
 
         for ( ; idx < idxe; idx++) {
                 out(CompactSparseMatrixHandler_tpl<off>::Ai[idx] - off + 1) =
                         CompactSparseMatrixHandler_tpl<off>::Ax[idx];
         }
 
         return out;
 }

Here is the call graph for this function:

template<int off>

MatrixHandler & DirCColMatrixHandler< off >::MulAndSumWithShift	(	MatrixHandler &	out,
		doublereal	s = `1.`,
		integer	drow = `0`,
		integer	dcol = `0`
	)		const

Definition at line 116 of file dirccmh.cc.

References SparseMatrixHandler::iGetNumCols(), MatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), MatrixHandler::iGetNumRows(), MatrixHandler::IncCoef(), and MBDYN_EXCEPT_ARGS.

 {
         silent_cerr("DirCColMatrixHandler<off>::MulAndSumWithShift called"
                         << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);            
         if ((out.iGetNumCols() < SparseMatrixHandler::iGetNumCols() + dcol)
                 || (out.iGetNumRows() < SparseMatrixHandler::iGetNumRows() + drow))
         {
                 silent_cerr("Assertion fault "
                                 "in DirCColMatrixHandler<off>::MulAndSumWithShift"
                                 << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
         drow = drow + 1;
         for (integer col = 0; col < SparseMatrixHandler::iGetNumCols(); col++) {
                 integer idx = CompactSparseMatrixHandler_tpl<off>::Ap[col] - off;
                 integer idxe = CompactSparseMatrixHandler_tpl<off>::Ap[col+1] - off;
                 integer newcol = col + dcol + 1;
                 for (; idx < idxe; idx++) {
                         out.IncCoef(CompactSparseMatrixHandler_tpl<off>::Ai[idx] - off + drow,
                                 newcol, CompactSparseMatrixHandler_tpl<off>::Ax[idx]*s);
                 }
         }
         return out;
 }

Here is the call graph for this function:

template<int off>

doublereal& DirCColMatrixHandler< off >::operator()	(	integer	i_row,
		integer	i_col
	)

inlinevirtual

Implements MatrixHandler.

Definition at line 71 of file dirccmh.h.

References ASSERTMSGBREAK, SparseMatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), MBDYN_EXCEPT_ARGS, and DirCColMatrixHandler< off >::pindices.

                                                               {
                 ASSERTMSGBREAK(i_row > 0 && i_row <= SparseMatrixHandler::iGetNumRows(),
                                 "Error in CColMatrixHandler::operator(), "
                                 "row index out of range");
                 ASSERTMSGBREAK(i_col > 0 && i_col <= SparseMatrixHandler::iGetNumCols(),
                                 "Error in CColMatrixHandler::operator(), "
                                 "col index out of range");
 
                 integer idx = pindices[i_col][i_row];
                 if (idx == -1) {
                         /* matrix must be rebuilt */
                         throw MatrixHandler::ErrRebuildMatrix(MBDYN_EXCEPT_ARGS);
                 }
 
                 return CompactSparseMatrixHandler_tpl<off>::Ax[idx];
         };

Here is the call graph for this function:

template<int off>

const doublereal& DirCColMatrixHandler< off >::operator()	(	integer	i_row,
		integer	i_col
	)		const

inlinevirtual

Implements MatrixHandler.

Definition at line 88 of file dirccmh.h.

References ASSERTMSGBREAK, SparseMatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), DirCColMatrixHandler< off >::pindices, and Zero1.

                                                                            {
                 ASSERTMSGBREAK(i_row > 0 && i_row <= SparseMatrixHandler::iGetNumRows(),
                                 "Error in CColMatrixHandler::operator(), "
                                 "row index out of range");
                 ASSERTMSGBREAK(i_col > 0 && i_col <= SparseMatrixHandler::iGetNumCols(),
                                 "Error in CColMatrixHandler::operator(), "
                                 "col index out of range");
 
                 integer idx = pindices[i_col][i_row];
                 if (idx == -1) {
                         /* matrix must be rebuilt */
                         return ::Zero1;
 
                 }
 
                 return CompactSparseMatrixHandler_tpl<off>::Ax[idx];
         };

Here is the call graph for this function:

template<int off>

void DirCColMatrixHandler< off >::Resize	(	integer	n,
		integer	nn
	)

virtual

Implements MatrixHandler.

Definition at line 83 of file dirccmh.cc.

References MBDYN_EXCEPT_ARGS.

 {
         silent_cerr("DirCColMatrixHandler<off>::Resize called" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

Member Data Documentation

template<int off>

std::vector<integer> DirCColMatrixHandler< off >::indices

private

Definition at line 54 of file dirccmh.h.

Referenced by DirCColMatrixHandler< off >::DirCColMatrixHandler().

template<int off>

std::vector<integer *> DirCColMatrixHandler< off >::pindices

private

Definition at line 53 of file dirccmh.h.

Referenced by DirCColMatrixHandler< off >::DirCColMatrixHandler(), and DirCColMatrixHandler< off >::operator()().

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

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

template<int off> class DirCColMatrixHandler< off >

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

template<int off>
class DirCColMatrixHandler< off >