CColMatrixHandler< off > Class Template Reference

#include <ccmh.h>

Inheritance diagram for CColMatrixHandler< off >:

Collaboration diagram for CColMatrixHandler< off >:

Public Member Functions
	CColMatrixHandler (std::vector< doublereal > &x, const std::vector< integer > &i, const std::vector< integer > &p)
virtual	~CColMatrixHandler ()
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 ir, integer ic)
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

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

Definition at line 46 of file ccmh.h.

Constructor & Destructor Documentation

template<int off>

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

Definition at line 41 of file ccmh.cc.

References NO_OP.

: CompactSparseMatrixHandler_tpl<off>(p.size() - 1, p.size() - 1, x, i, p)
{
        NO_OP;
}

template<int off>

CColMatrixHandler< off >::~CColMatrixHandler ( )

virtual

Definition at line 50 of file ccmh.cc.

References NO_OP.

{
        NO_OP;
}

Member Function Documentation

template<int off>

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

virtual

Implements CompactSparseMatrixHandler.

Definition at line 59 of file ccmh.cc.

References CompactSparseMatrixHandler::bMatDuplicate.

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

        return p;
}

template<int off>

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

Definition at line 132 of file ccmh.cc.

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

{
        silent_cerr("CColMatrixHandler<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 CColMatrixHandler<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];
                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 & CColMatrixHandler< off >::GetCol ( integer  icol, VectorHandler &  out  ) const

virtual

Implements SparseMatrixHandler.

Definition at line 82 of file ccmh.cc.

References SparseMatrixHandler::iGetNumCols(), and MBDYN_EXCEPT_ARGS.

{
        // NOTE: out must be zeroed by caller
        
        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 & CColMatrixHandler< off >::MulAndSumWithShift	(	MatrixHandler &	out,
		doublereal	s = 1.,
		integer	drow = 0,
		integer	dcol = 0
	)		const

Definition at line 104 of file ccmh.cc.

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

{
        silent_cerr("CColMatrixHandler<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 CColMatrixHandler<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& CColMatrixHandler< off >::operator() ( integer  i_row, integer  i_col  )

inlinevirtual

Implements MatrixHandler.

Definition at line 66 of file ccmh.h.

References ASSERTMSGBREAK, SparseMatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), and MBDYN_EXCEPT_ARGS.

                                                              {
                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");
                i_row--;
                integer row_begin = CompactSparseMatrixHandler_tpl<off>::Ap[i_col - 1] - off;
                integer row_end = CompactSparseMatrixHandler_tpl<off>::Ap[i_col] - off - 1;
                integer idx;
                integer row;

                if (row_begin == (CompactSparseMatrixHandler_tpl<off>::Ap[i_col] - off)
                                || (CompactSparseMatrixHandler_tpl<off>::Ai[row_begin] - off) > i_row
                                || (CompactSparseMatrixHandler_tpl<off>::Ai[row_end] - off) < i_row)
                {
                        /* matrix must be rebuilt */
                        throw MatrixHandler::ErrRebuildMatrix(MBDYN_EXCEPT_ARGS);
                }

                while (row_end >= row_begin) {
                        idx = (row_begin + row_end)/2;
                        row = CompactSparseMatrixHandler_tpl<off>::Ai[idx] - off;
                        if (i_row < row) {
                                row_end = idx - 1;
                        } else if (i_row > row) {
                                row_begin = idx + 1;
                        } else {
                                return CompactSparseMatrixHandler_tpl<off>::Ax[idx];
                        }
                }

                /* matrix must be rebuilt */
                throw MatrixHandler::ErrRebuildMatrix(MBDYN_EXCEPT_ARGS);
        };

Here is the call graph for this function:

template<int off>

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

inlinevirtual

Implements MatrixHandler.

Definition at line 103 of file ccmh.h.

References ASSERTMSGBREAK, SparseMatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), 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");
                i_row--;
                integer row_begin = CompactSparseMatrixHandler_tpl<off>::Ap[i_col - 1] - off;
                integer row_end = CompactSparseMatrixHandler_tpl<off>::Ap[i_col] - off - 1;
                integer idx;
                integer row;

                if (row_begin == CompactSparseMatrixHandler_tpl<off>::Ap[i_col] - off
                                || CompactSparseMatrixHandler_tpl<off>::Ai[row_begin] - off > i_row
                                || CompactSparseMatrixHandler_tpl<off>::Ai[row_end] - off < i_row)
                {
                        return ::Zero1;
                }

                while (row_end >= row_begin) {
                        idx = (row_begin + row_end)/2;
                        row = CompactSparseMatrixHandler_tpl<off>::Ai[idx] - off;
                        if (i_row < row) {
                                row_end = idx - 1;
                        } else if (i_row > row) {
                                row_begin = idx + 1;
                        } else {
                                return CompactSparseMatrixHandler_tpl<off>::Ax[idx];
                        }
                }

                return ::Zero1;
        };

Here is the call graph for this function:

template<int off>

void CColMatrixHandler< off >::Resize ( integer  ir, integer  ic  )

virtual

Implements MatrixHandler.

Definition at line 73 of file ccmh.cc.

References MBDYN_EXCEPT_ARGS.

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

---

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

• ccmh.h
• ccmh.cc