FullMatrixHandler Class Reference

#include <fullmh.h>

Inheritance diagram for FullMatrixHandler:

Collaboration diagram for FullMatrixHandler:

Classes
class	const_iterator

Public Member Functions
	FullMatrixHandler (const FullMatrixHandler &)
FullMatrixHandler &	operator= (const FullMatrixHandler &)
FullMatrixHandler::const_iterator	begin (void) const
const FullMatrixHandler::const_iterator &	end (void) const
	FullMatrixHandler (doublereal *pd, doublereal **ppd, integer iSize, integer iNR, integer iNC, integer iMaxCols=0)
	FullMatrixHandler (integer iNR, integer iNC=0)
	FullMatrixHandler (void)
virtual	~FullMatrixHandler (void)
virtual void	Resize (integer iNewRows, integer iNewCols)
void	Detach (void)
void	Reset (void)
void	Attach (integer iNewRows, integer iNewCols, doublereal *pd, doublereal **ppd, integer iMSize=0, integer iMaxC=0)
const doublereal *	pdGetMat (void) const
doublereal *	pdGetMat (void)
virtual const doublereal *	pdGetVec (integer iCol) const
virtual doublereal *	pdGetVec (integer iCol)
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 integer	iGetNumRows (void) const
virtual integer	iGetNumCols (void) const
virtual doublereal &	operator() (integer iRow, integer iCol)
virtual const doublereal &	operator() (integer iRow, integer iCol) const
virtual MatrixHandler &	operator+= (const SubMatrixHandler &SubMH)
virtual MatrixHandler &	operator-= (const SubMatrixHandler &SubMH)
virtual MatrixHandler &	operator+= (const VariableSubMatrixHandler &SubMH)
virtual MatrixHandler &	operator-= (const VariableSubMatrixHandler &SubMH)
void	MatMul (const FullMatrixHandler &m1, const FullMatrixHandler &m2)
void	Add (integer iRow, integer iCol, const FullMatrixHandler &source)
void	Sub (integer iRow, integer iCol, const FullMatrixHandler &source)
void	Put (integer iRow, integer iCol, const FullMatrixHandler &source)
void	Add (integer iRow, integer iCol, const FullMatrixHandler &source, const doublereal dCoef)
void	Sub (integer iRow, integer iCol, const FullMatrixHandler &source, const doublereal dCoef)
void	Put (integer iRow, integer iCol, const FullMatrixHandler &source, const doublereal dCoef)
void	AddT (integer iRow, integer iCol, const FullMatrixHandler &source)
void	SubT (integer iRow, integer iCol, const FullMatrixHandler &source)
void	PutT (integer iRow, integer iCol, const FullMatrixHandler &source)
void	AddT (integer iRow, integer iCol, const FullMatrixHandler &source, const doublereal dCoef)
void	SubT (integer iRow, integer iCol, const FullMatrixHandler &source, const doublereal dCoef)
void	PutT (integer iRow, integer iCol, const FullMatrixHandler &source, const doublereal dCoef)
void	Add (integer iRow, integer iCol, const Vec3 &v)
void	Sub (integer iRow, integer iCol, const Vec3 &v)
void	Put (integer iRow, integer iCol, const Vec3 &v)
void	AddT (integer iRow, integer iCol, const Vec3 &v)
void	SubT (integer iRow, integer iCol, const Vec3 &v)
void	PutT (integer iRow, integer iCol, const Vec3 &v)
void	Add (integer iRow, integer iCol, const Mat3x3 &m)
void	AddT (integer iRow, integer iCol, const Mat3x3 &m)
void	Sub (integer iRow, integer iCol, const Mat3x3 &m)
void	SubT (integer iRow, integer iCol, const Mat3x3 &m)
void	Put (integer iRow, integer iCol, const Mat3x3 &m)
void	PutT (integer iRow, integer iCol, const Mat3x3 &m)
void	Add (integer iRow, integer iCol, const Mat3xN &m)
void	AddT (integer iRow, integer iCol, const Mat3xN &m)
void	Sub (integer iRow, integer iCol, const Mat3xN &m)
void	SubT (integer iRow, integer iCol, const Mat3xN &m)
void	Put (integer iRow, integer iCol, const Mat3xN &m)
void	PutT (integer iRow, integer iCol, const Mat3xN &m)
void	Add (integer iRow, integer iCol, const MatNx3 &m)
void	Sub (integer iRow, integer iCol, const MatNx3 &m)
void	Put (integer iRow, integer iCol, const MatNx3 &m)
void	CopyMatrixRow (integer dest_row, const FullMatrixHandler &source, integer source_row)
void	CopyMatrixBlock (integer dest_row, integer dest_col, const FullMatrixHandler &source, integer source_start_row, integer source_end_row, integer source_start_col, integer source_end_col)
Public Member Functions inherited from MatrixHandler
virtual	~MatrixHandler (void)
virtual void	ResizeReset (integer, integer)
virtual integer *	piGetRows (void) const
virtual integer *	piGetCols (void) const
virtual integer	PacMat (void)
virtual MatrixHandler &	operator= (const MatrixHandler &MH)
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
void	CreateColRow (integer iNR, integer iNC)
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
bool	bOwnsMemory
integer	iNumRows
integer	iNumCols
integer	iRawSize
integer	iMaxCols
doublereal *	pdRaw
doublereal *	pdRawm1
doublereal **	ppdCols
doublereal **	ppdColsm1

Private Attributes
const_iterator	m_end

Friends
class	FullSubMatrixHandler
class	SparseSubMatrixHandler
std::ostream &	operator<< (std::ostream &out, const FullMatrixHandler &m)
std::ostream &	Write (std::ostream &out, const FullMatrixHandler &m, const char *s, const char *s2)

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

Detailed Description

Definition at line 44 of file fullmh.h.

Constructor & Destructor Documentation

FullMatrixHandler::FullMatrixHandler

(

const FullMatrixHandler &

MH

)

Definition at line 129 of file fullmh.cc.

References iNumCols, iNumRows, pdRaw, and Resize().

: bOwnsMemory(true),
iNumRows(0), iNumCols(0), iRawSize(0), iMaxCols(0),
pdRaw(NULL), pdRawm1(NULL),
ppdCols(NULL), ppdColsm1(NULL), m_end(*this, true)
{
        if (MH.pdRaw != 0) {
                Resize(MH.iNumRows, MH.iNumCols);

                for (integer i = 0; i < iNumRows*iNumCols; i++) {
                        pdRaw[i] = MH.pdRaw[i];
                }
        }
}

Here is the call graph for this function:

FullMatrixHandler::FullMatrixHandler	(	doublereal *	pd,
		doublereal **	ppd,
		integer	iSize,
		integer	iNR,
		integer	iNC,
		integer	iMaxCols = 0
	)

Definition at line 75 of file fullmh.cc.

References ASSERT, CreateColRow(), iMaxCols, iNumCols, iNumRows, pdRawm1, and ppdColsm1.

: bOwnsMemory(false),
iNumRows(iNR), iNumCols(iNC), iRawSize(iSize), iMaxCols(iMaxC),
pdRaw(pd), pdRawm1(0),
ppdCols(ppd), ppdColsm1(0), m_end(*this, true)
{
        if (iMaxCols == 0) {
                iMaxCols = iNumCols;
        }

        ASSERT(pd);
        pdRawm1 = pd - 1;
        ASSERT(ppd);
        ppdColsm1 = ppd - 1;

        CreateColRow(iNumRows, iNumCols);

#ifdef DEBUG
        IsValid();
#endif /* DEBUG */
}

Here is the call graph for this function:

FullMatrixHandler::FullMatrixHandler	(	integer	iNR,
		integer	iNC = 0
	)

Definition at line 100 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, Reset(), and Resize().

: bOwnsMemory(true),
iNumRows(iNR), iNumCols(iNC ? iNC : iNumRows),
iRawSize(iNumRows*iNumCols), iMaxCols(iNumCols),
pdRaw(NULL), pdRawm1(NULL),
ppdCols(NULL), ppdColsm1(NULL), m_end(*this, true)
{
        ASSERT(iNumRows > 0);
        ASSERT(iNumCols > 0);
        Resize(iNumRows, iNumCols);
        Reset();

#ifdef DEBUG
        IsValid();
#endif /* DEBUG */
}

Here is the call graph for this function:

FullMatrixHandler::FullMatrixHandler

(

void

)

Definition at line 119 of file fullmh.cc.

References NO_OP.

: bOwnsMemory(true),
iNumRows(0), iNumCols(0), iRawSize(0), iMaxCols(0),
pdRaw(NULL), pdRawm1(NULL),
ppdCols(NULL), ppdColsm1(NULL), m_end(*this, true)
{
        NO_OP;
}

FullMatrixHandler::~FullMatrixHandler ( void  )

virtual

Definition at line 167 of file fullmh.cc.

References Detach().

{
        Detach();
}

Here is the call graph for this function:

Member Function Documentation

void FullMatrixHandler::Add	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source
	)

Definition at line 1039 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

Referenced by FullSubMatrixHandler::Add(), and ModuleNonsmoothNode::mbs_get_force_frictional().

{
        integer nr = source.iGetNumRows();
        integer nc = source.iGetNumCols();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ir, ic) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] += source.ppdColsm1[ic][ir];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::Add	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source,
		const doublereal	dCoef
	)

Definition at line 1150 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

{
        integer nr = source.iGetNumRows();
        integer nc = source.iGetNumCols();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ir, ic) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] += dCoef*source.ppdColsm1[ic][ir];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::Add	(	integer	iRow,
		integer	iCol,
		const Vec3 &	v
	)

Definition at line 1490 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, Vec3::pGetVec(), ppdColsm1, V1, V2, and V3.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols);
#endif /* DEBUG */

        /* Nota: assume che la matrice sia organizzata
         * "per colonne" (stile FORTRAN)
         */
        const doublereal* pdFrom = v.pGetVec();

        ppdColsm1[iCol][iRow] += pdFrom[V1];

        ppdColsm1[iCol][iRow + 1] += pdFrom[V2];

        ppdColsm1[iCol][iRow + 2] += pdFrom[V3];
}

Here is the call graph for this function:

void FullMatrixHandler::Add	(	integer	iRow,
		integer	iCol,
		const Mat3x3 &	m
	)

Definition at line 1849 of file fullmh.cc.

References ASSERT, c, iNumCols, iNumRows, M11, M12, M13, M21, M22, M23, M31, M32, M33, Mat3x3::pGetMat(), and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        /* Nota: assume che la matrice sia organizzata
         * "per colonne" (stile FORTRAN) e assume
         * che la matrice Mat3x3 sia organizzata anch'essa per colonne */
        const doublereal* pdFrom = m.pGetMat();

        ppdColsm1[iCol][iRow] += pdFrom[M11];
        ppdColsm1[iCol][iRow + 1] += pdFrom[M21];
        ppdColsm1[iCol][iRow + 2] += pdFrom[M31];

        ppdColsm1[iCol + 1][iRow] += pdFrom[M12];
        ppdColsm1[iCol + 1][iRow + 1] += pdFrom[M22];
        ppdColsm1[iCol + 1][iRow + 2] += pdFrom[M32];

        ppdColsm1[iCol + 2][iRow] += pdFrom[M13];
        ppdColsm1[iCol + 2][iRow + 1] += pdFrom[M23];
        ppdColsm1[iCol + 2][iRow + 2] += pdFrom[M33];

#if 0
        for (unsigned c = 0; c < 2; c++) {
                for (unsigned r = 0; r < 2; r++) {
                        ppdColsm1[iCol + c][iRow + r] += pdFrom[r];
                }
                pdFrom += 3;
        }
#endif
}

Here is the call graph for this function:

void FullMatrixHandler::Add	(	integer	iRow,
		integer	iCol,
		const Mat3xN &	m
	)

Definition at line 2080 of file fullmh.cc.

References ASSERT, c, Mat3xN::iGetNumCols(), iNumCols, iNumRows, and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - m.iGetNumCols() + 1);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = 3; r > 0; r--) {
                for (integer c = m.iGetNumCols(); c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] += m(r, c);
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::Add	(	integer	iRow,
		integer	iCol,
		const MatNx3 &	m
	)

Definition at line 2236 of file fullmh.cc.

References ASSERT, c, MatNx3::iGetNumRows(), iNumCols, iNumRows, and ppdColsm1.

{
#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - m.iGetNumRows() + 1);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = m.iGetNumRows(); r > 0; r--) {
                for (integer c = 3; c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] += m(r, c);
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::AddT	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source
	)

Definition at line 1264 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

Referenced by FullSubMatrixHandler::AddT().

{
        integer nr = source.iGetNumCols();
        integer nc = source.iGetNumRows();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ic, ir) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] += source.ppdColsm1[ir][ic];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::AddT	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source,
		const doublereal	dCoef
	)

Definition at line 1375 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

{
        integer nr = source.iGetNumCols();
        integer nc = source.iGetNumRows();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ic, ir) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] += dCoef*source.ppdColsm1[ir][ic];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::AddT	(	integer	iRow,
		integer	iCol,
		const Vec3 &	v
	)

Definition at line 1577 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, Vec3::pGetVec(), ppdColsm1, V1, V2, and V3.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        /* Nota: assume che la matrice sia organizzata
         * "per colonne" (stile FORTRAN)
         */
        const doublereal* pdFrom = v.pGetVec();

        ppdColsm1[iCol][iRow] += pdFrom[V1];

        ppdColsm1[iCol + 1][iRow] += pdFrom[V2];

        ppdColsm1[iCol + 2][iRow] += pdFrom[V3];
}

Here is the call graph for this function:

void FullMatrixHandler::AddT	(	integer	iRow,
		integer	iCol,
		const Mat3x3 &	m
	)

Definition at line 1892 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, M11, M12, M13, M21, M22, M23, M31, M32, M33, Mat3x3::pGetMat(), and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        /* Nota: assume che la matrice sia organizzata
         * "per colonne" (stile FORTRAN) e assume
         * che la matrice Mat3x3 sia organizzata anch'essa per colonne */
        const doublereal* pdFrom = m.pGetMat();

        ppdColsm1[iCol][iRow] += pdFrom[M11];
        ppdColsm1[iCol][iRow + 1] += pdFrom[M12];
        ppdColsm1[iCol][iRow + 2] += pdFrom[M13];

        ppdColsm1[iCol + 1][iRow] += pdFrom[M21];
        ppdColsm1[iCol + 1][iRow + 1] += pdFrom[M22];
        ppdColsm1[iCol + 1][iRow + 2] += pdFrom[M23];

        ppdColsm1[iCol + 2][iRow] += pdFrom[M31];
        ppdColsm1[iCol + 2][iRow + 1] += pdFrom[M32];
        ppdColsm1[iCol + 2][iRow + 2] += pdFrom[M33];
}

Here is the call graph for this function:

void FullMatrixHandler::AddT	(	integer	iRow,
		integer	iCol,
		const Mat3xN &	m
	)

Definition at line 2158 of file fullmh.cc.

References ASSERT, c, Mat3xN::iGetNumCols(), iNumCols, iNumRows, and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - m.iGetNumCols() + 1);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = m.iGetNumCols(); r > 0; r--) {
                for (integer c = 3; c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] += m(c, r);
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::Attach	(	integer	iNewRows,
		integer	iNewCols,
		doublereal *	pd,
		doublereal **	ppd,
		integer	iMSize = 0,
		integer	iMaxC = 0
	)

Definition at line 290 of file fullmh.cc.

References bOwnsMemory, Detach(), iMaxCols, iNumCols, iNumRows, iRawSize, MBDYN_EXCEPT_ARGS, pdRaw, pdRawm1, ppdCols, and ppdColsm1.

{
        if (bOwnsMemory || pdRaw != NULL) {
                Detach();
                bOwnsMemory = false;
        }

        iNumRows = iNewRows;
        iNumCols = iNewCols;
        if (iMSize > 0) {
                if (iMSize < iNumRows*iNumCols) {
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                iRawSize = iMSize;

        } else {
                iRawSize = iNumRows*iNumCols;
        }

        if (iMaxC > 0) {
                if (iMaxC < iNumCols) {
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                iMaxCols = iMaxC;

        } else {
                iMaxCols = iNumCols;
        }
        pdRaw = pd;
        pdRawm1 = pdRaw - 1;
        ppdCols = ppd;
        ppdColsm1 = ppdCols - 1;

#ifdef DEBUG
        IsValid();
#endif /* DEBUG */
}

Here is the call graph for this function:

FullMatrixHandler::const_iterator FullMatrixHandler::begin ( void  ) const

inline

Definition at line 102 of file fullmh.h.

Referenced by main().

                                                          {
                return const_iterator(*this);
        };

void FullMatrixHandler::CopyMatrixBlock	(	integer	dest_row,
		integer	dest_col,
		const FullMatrixHandler &	source,
		integer	source_start_row,
		integer	source_end_row,
		integer	source_start_col,
		integer	source_end_col
	)

Definition at line 2319 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), and ppdColsm1.

{
        ASSERT(source_start_row >= 1);
        ASSERT(source_end_row <= source.iGetNumRows());
        ASSERT(source_start_row <= source_end_row);
        ASSERT(dest_row >= 1);
        ASSERT(dest_row + (source_end_row - source_start_row) <= iGetNumRows());

        ASSERT(source_start_col >= 1);
        ASSERT(source_end_col <= source.iGetNumCols());
        ASSERT(source_start_col <= source_end_col);
        ASSERT(dest_col >= 1);
        ASSERT(dest_col + (source_end_col - source_start_col) <= iGetNumCols());

        for (integer ir = source_start_row; ir <= source_end_row; ir++) {
                integer row = dest_row + (ir - source_start_row);
                for (integer ic = source_start_col; ic <= source_end_col; ic++) {
                        integer col = dest_col + (ic - source_start_col);
                        // dest(row, col) = source(ir, ic);
                        ppdColsm1[col][row] = source.ppdColsm1[ic][ir];
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::CopyMatrixRow	(	integer	dest_row,
		const FullMatrixHandler &	source,
		integer	source_row
	)

Definition at line 2301 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), and ppdColsm1.

Referenced by Shell4EASANS::AssRes().

{
        integer nc = iGetNumCols();

        ASSERT(dest_row >= 1);
        ASSERT(dest_row <= iGetNumRows());
        ASSERT(source_row >= 1);
        ASSERT(source_row <= source.iGetNumRows());
        ASSERT(nc == source.iGetNumCols());

        for (integer ic = 1; ic <= nc; ic++) {
                // dest(dest_row, ic) = source(source_row, ic);
                ppdColsm1[ic][dest_row] = source.ppdColsm1[ic][source_row];
        }
}

Here is the call graph for this function:

void FullMatrixHandler::CreateColRow ( integer  iNR, integer  iNC  )

protected

Definition at line 58 of file fullmh.cc.

References ASSERT, c, pdRawm1, and ppdColsm1.

Referenced by FullMatrixHandler(), and Resize().

{
        ASSERT(iNC > 0);
        ASSERT(iNC > 0);
        ASSERT(pdRawm1);
        ASSERT(ppdColsm1);
        
        doublereal *pd = pdRawm1 + iNC*iNR;

        for (integer c = iNC; c > 0; c--) {
                pd -= iNR;
                ppdColsm1[c] = pd;
        }

        ASSERT(pd == pdRawm1);
}

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

inlinevirtual

Reimplemented from MatrixHandler.

Reimplemented in FullSubMatrixHandler.

Definition at line 202 of file fullmh.h.

References ASSERT, iNumCols, iNumRows, and ppdColsm1.

Referenced by FullSubMatrixHandler::DecCoef().

                                                                     {
#ifdef DEBUG
                IsValid();
                ASSERT(iRow > 0 && iRow <= iNumRows);
                ASSERT(iCol > 0 && iCol <= iNumCols);
#endif /* DEBUG */

                ppdColsm1[iCol][iRow] -= dCoef;
        };

void FullMatrixHandler::Detach

(

void

)

Definition at line 268 of file fullmh.cc.

References bOwnsMemory, iMaxCols, iNumCols, iNumRows, iRawSize, pdRaw, pdRawm1, ppdCols, ppdColsm1, and SAFEDELETEARR.

Referenced by Attach(), and ~FullMatrixHandler().

{
        if (bOwnsMemory) {
                if (pdRaw != NULL) {
                        SAFEDELETEARR(pdRaw);
                }
                if (ppdCols != NULL) {
                        SAFEDELETEARR(ppdCols);
                }
        }
        iRawSize = iMaxCols = iNumRows = iNumCols = 0;
        pdRaw = pdRawm1 = NULL;
        ppdCols = ppdColsm1 = NULL;

        /* so Resize() can be safely invoked */
        bOwnsMemory = true;
}

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

inlinevirtual

Reimplemented from MatrixHandler.

Reimplemented in FullSubMatrixHandler.

Definition at line 214 of file fullmh.h.

References ASSERT, iNumCols, iNumRows, and ppdColsm1.

Referenced by ModuleFMU::AssJac(), and FullSubMatrixHandler::dGetCoef().

                                                   {
#ifdef DEBUG
                IsValid();
                ASSERT(iRow > 0 && iRow <= iNumRows);
                ASSERT(iCol > 0 && iCol <= iNumCols);
#endif /* DEBUG */

                return ppdColsm1[iCol][iRow];
        };

const FullMatrixHandler::const_iterator& FullMatrixHandler::end ( void  ) const

inline

Definition at line 106 of file fullmh.h.

References m_end.

Referenced by main().

                                                               {
                return m_end;
        };

virtual integer FullMatrixHandler::iGetNumCols ( void  ) const

inlinevirtual

Implements MatrixHandler.

Reimplemented in FullSubMatrixHandler.

Definition at line 229 of file fullmh.h.

References iNumCols.

Referenced by Add(), AddT(), FullSubMatrixHandler::AddTo(), SparseSubMatrixHandler::AddTo(), FullSubMatrixHandler::AddToT(), SparseSubMatrixHandler::AddToT(), CopyMatrixBlock(), CopyMatrixRow(), MatMatMul_base(), MatMul(), MatTMatMul_base(), MatTVecMul_base(), MatVecMul_base(), DataManager::OutputEigFullMatrices(), Put(), PutT(), ReadMembraneConstLaw(), ReadShellConstLaw(), Sub(), FullSubMatrixHandler::SubFrom(), SparseSubMatrixHandler::SubFrom(), FullSubMatrixHandler::SubFromT(), SparseSubMatrixHandler::SubFromT(), and SubT().

                                                {
                return iNumCols;
        };

virtual integer FullMatrixHandler::iGetNumRows ( void  ) const

inlinevirtual

Implements MatrixHandler.

Reimplemented in FullSubMatrixHandler.

Definition at line 225 of file fullmh.h.

References iNumRows.

Referenced by Add(), AddT(), FullSubMatrixHandler::AddTo(), SparseSubMatrixHandler::AddTo(), FullSubMatrixHandler::AddToT(), SparseSubMatrixHandler::AddToT(), CopyMatrixBlock(), CopyMatrixRow(), MatMatMul_base(), MatMul(), MatTMatMul_base(), MatTVecMul_base(), MatVecMul_base(), DataManager::OutputEigFullMatrices(), Put(), PutT(), ReadMembraneConstLaw(), ReadShellConstLaw(), Sub(), FullSubMatrixHandler::SubFrom(), SparseSubMatrixHandler::SubFrom(), FullSubMatrixHandler::SubFromT(), SparseSubMatrixHandler::SubFromT(), and SubT().

                                                {
                return iNumRows;
        };

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

inlinevirtual

Reimplemented from MatrixHandler.

Reimplemented in FullSubMatrixHandler.

Definition at line 190 of file fullmh.h.

References ASSERT, iNumCols, iNumRows, and ppdColsm1.

Referenced by FullSubMatrixHandler::IncCoef().

                                                                     {
#ifdef DEBUG
                IsValid();
                ASSERT(iRow > 0 && iRow <= iNumRows);
                ASSERT(iCol > 0 && iCol <= iNumCols);
#endif /* DEBUG */

                ppdColsm1[iCol][iRow] += dCoef;
        };

MatrixHandler & FullMatrixHandler::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 601 of file fullmh.cc.

References c, MatrixHandler::DecCoef(), MatrixHandler::iGetNumCols(), iGetNumCols(), MatrixHandler::iGetNumRows(), iGetNumRows(), MatrixHandler::IncCoef(), MatrixHandler::MatMatMul_base(), MBDYN_EXCEPT_ARGS, pdRaw, ppdColsm1, and MatrixHandler::PutCoef().

{
        const FullMatrixHandler *pin = dynamic_cast<const FullMatrixHandler *>(&in);
        if (pin == 0) {
                // if input matrix is not FullMatrixHandler, use generic function
                return MatrixHandler::MatMatMul_base(op, out, in);
        }

        integer out_nc = out.iGetNumCols();
        integer out_nr = out.iGetNumRows();
        integer in_nr = in.iGetNumRows();

        if (out_nr != iGetNumRows()
                || out_nc != in.iGetNumCols()
                || in_nr != iGetNumCols())
        {
                const char *strop;

                if (op == &MatrixHandler::IncCoef) {
                        strop = "+=";
                } else if (op == &MatrixHandler::DecCoef) {
                        strop = "-=";
                } else if (op == &MatrixHandler::PutCoef) {
                        strop = "=";
                } else {
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                silent_cerr("FullMatrixHandler::MatMatMul_base: size mismatch "
                        "out(" << out_nr << ", " << out_nc << ") "
                        << strop << " this(" << iGetNumRows() << ", " << iGetNumCols() << ") "
                        "* in(" << in_nr << ", " << in.iGetNumCols() << ")"
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        FullMatrixHandler *pout = dynamic_cast<FullMatrixHandler *>(&out);
        if (pout == 0) {
                // if output matrix is not FullMatrixHandler,
                // optimize coefficient computation
                for (integer c = 1; c <= out_nc; c++) {
                        for (integer r = 1; r <= out_nr; r++) {
                                doublereal d = 0.;

                                for (integer k = 1; k <= in_nr; k++) {
                                        // this(r, k) * in(k, c)
                                        d += ppdColsm1[k][r]*pin->ppdColsm1[c][k];
                                }

                                (out.*op)(r, c, d);
                        }
                }

                return out;
        }

        // if all matrices are FullMatrixHandler,
        // either use LAPACK or optimize coefficient computation and store
#ifdef USE_LAPACK
        doublereal ALPHA, BETA;
        if (op == &MatrixHandler::IncCoef) {
                BETA = 1.;
                ALPHA = 1.;

        } else if (op == &MatrixHandler::DecCoef) {
                BETA = 1.;
                ALPHA = -1.;

        } else if (op == &MatrixHandler::PutCoef) {
                BETA = 0.;
                ALPHA = 1.;

        } else {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        const char *TRANSA = "N";
        const char *TRANSB = "N";
        __FC_DECL__(dgemm)(TRANSA, TRANSB, &out_nr, &out_nc, &in_nr,
                &ALPHA,
                this->pdRaw, &out_nr,
                pin->pdRaw, &in_nr,
                &BETA,
                pout->pdRaw, &out_nr);

#else // ! USE_LAPACK
        if (op == &MatrixHandler::IncCoef) {
                for (integer c = 1; c <= out_nc; c++) {
                        for (integer r = 1; r <= out_nr; r++) {
                                doublereal d = 0.;

                                for (integer k = 1; k <= in_nr; k++) {
                                        // this(r, k) * in(k, c)
                                        d += ppdColsm1[k][r]*pin->ppdColsm1[c][k];
                                }

                                // (out.*op)(r, c, d);
                                pout->ppdColsm1[c][r] += d;
                        }
                }

        } else if (op == &MatrixHandler::DecCoef) {
                for (integer c = 1; c <= out_nc; c++) {
                        for (integer r = 1; r <= out_nr; r++) {
                                doublereal d = 0.;

                                for (integer k = 1; k <= in_nr; k++) {
                                        // this(r, k) * in(k, c)
                                        d += ppdColsm1[k][r]*pin->ppdColsm1[c][k];
                                }

                                // (out.*op)(r, c, d);
                                pout->ppdColsm1[c][r] -= d;
                        }
                }

        } else if (op == &MatrixHandler::PutCoef) {
                for (integer c = 1; c <= out_nc; c++) {
                        for (integer r = 1; r <= out_nr; r++) {
                                doublereal d = 0.;

                                for (integer k = 1; k <= in_nr; k++) {
                                        // this(r, k) * in(k, c)
                                        d += ppdColsm1[k][r]*pin->ppdColsm1[c][k];
                                }

                                // (out.*op)(r, c, d);
                                pout->ppdColsm1[c][r] = d;
                        }
                }

        } else {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif // ! USE_LAPACK

        return out;
}

Here is the call graph for this function:

void FullMatrixHandler::MatMul	(	const FullMatrixHandler &	m1,
		const FullMatrixHandler &	m2
	)

Definition at line 433 of file fullmh.cc.

References iGetNumCols(), iGetNumRows(), iNumCols, iNumRows, MBDYN_EXCEPT_ARGS, and ppdColsm1.

Referenced by ModuleNonsmoothNode::mbs_get_force_frictional().

{
        if (m1.iGetNumRows() != iGetNumRows()
                || m2.iGetNumRows() != m1.iGetNumCols()
                || m2.iGetNumCols() != iGetNumCols())
        {
                silent_cerr("FullMatrixHandler::MatMul: size mismatch "
                        "this(" << iGetNumRows() << ", " << iGetNumCols() << ") "
                        "= m1(" << m1.iGetNumRows() << ", " << m1.iGetNumCols() << ") "
                        "* m2(" << m2.iGetNumRows() << ", " << m2.iGetNumCols() << ")"
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        integer nc = m1.iGetNumCols();

        for (integer iRow = 1; iRow <= iNumRows; iRow++) {
                for (integer iCol = 1; iCol <= iNumCols; iCol++) {
                        ppdColsm1[iCol][iRow] = 0.;
                        for (integer iK = 1; iK <= nc; iK++) {
                                // ppdColsm1[iCol][iRow] += m1(iRow, iK)*m2(iK, iCol);
                                ppdColsm1[iCol][iRow] += m1.ppdColsm1[iK][iRow]*m2.ppdColsm1[iCol][iK];
                        }
                }
        }
}

Here is the call graph for this function:

MatrixHandler & FullMatrixHandler::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 744 of file fullmh.cc.

References c, MatrixHandler::DecCoef(), MatrixHandler::iGetNumCols(), iGetNumCols(), MatrixHandler::iGetNumRows(), iGetNumRows(), MatrixHandler::IncCoef(), MatrixHandler::MatTMatMul_base(), MBDYN_EXCEPT_ARGS, pdRaw, ppdColsm1, and MatrixHandler::PutCoef().

{
        const FullMatrixHandler *pin = dynamic_cast<const FullMatrixHandler *>(&in);
        if (pin == 0) {
                // if input matrix is not FullMatrixHandler, use generic function
                return MatrixHandler::MatTMatMul_base(op, out, in);
        }

        integer out_nc = out.iGetNumCols();
        integer out_nr = out.iGetNumRows();
        integer in_nr = in.iGetNumRows();

        if (out_nr != iGetNumCols()
                || out_nc != in.iGetNumCols()
                || in_nr != iGetNumRows())
        {
                const char *strop;

                if (op == &MatrixHandler::IncCoef) {
                        strop = "+=";
                } else if (op == &MatrixHandler::DecCoef) {
                        strop = "-=";
                } else if (op == &MatrixHandler::PutCoef) {
                        strop = "=";
                } else {
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                silent_cerr("FullMatrixHandler::MatTMatMul_base: size mismatch "
                        "out(" << out_nr << ", " << out_nc << ") "
                        << strop << " this(" << iGetNumRows() << ", " << iGetNumCols() << ")^T "
                        "* in(" << in_nr << ", " << in.iGetNumCols() << ")"
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        FullMatrixHandler *pout = dynamic_cast<FullMatrixHandler *>(&out);
        if (pout == 0) {
                // if output matrix is not FullMatrixHandler,
                // optimize coefficient computation
                for (integer c = 1; c <= out_nc; c++) {
                        for (integer r = 1; r <= out_nr; r++) {
                                doublereal d = 0.;

                                for (integer k = 1; k <= in_nr; k++) {
                                        // this(k, r) * in(k, c)
                                        d += ppdColsm1[r][k]*pin->ppdColsm1[c][k];
                                }

                                (out.*op)(r, c, d);
                        }
                }

                return out;
        }

        // if all matrices are FullMatrixHandler,
        // either use LAPACK or optimize coefficient computation and store
#ifdef USE_LAPACK
        doublereal ALPHA, BETA;
        if (op == &MatrixHandler::IncCoef) {
                BETA = 1.;
                ALPHA = 1.;

        } else if (op == &MatrixHandler::DecCoef) {
                BETA = 1.;
                ALPHA = -1.;

        } else if (op == &MatrixHandler::PutCoef) {
                BETA = 0.;
                ALPHA = 1.;

        } else {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        const char *TRANSA = "T";
        const char *TRANSB = "N";
        __FC_DECL__(dgemm)(TRANSA, TRANSB, &out_nr, &out_nc, &in_nr,
                &ALPHA,
                this->pdRaw, &in_nr,
                pin->pdRaw, &in_nr,
                &BETA,
                pout->pdRaw, &out_nr);

#else // ! USE_LAPACK
        if (op == &MatrixHandler::IncCoef) {
                for (integer c = 1; c <= out_nc; c++) {
                        for (integer r = 1; r <= out_nr; r++) {
                                doublereal d = 0.;

                                for (integer k = 1; k <= in_nr; k++) {
                                        // this(k, r) * in(k, c)
                                        d += ppdColsm1[r][k]*pin->ppdColsm1[c][k];
                                }

                                // (out.*op)(r, c, d);
                                pout->ppdColsm1[c][r] += d;
                        }
                }

        } else if (op == &MatrixHandler::DecCoef) {
                for (integer c = 1; c <= out_nc; c++) {
                        for (integer r = 1; r <= out_nr; r++) {
                                doublereal d = 0.;

                                for (integer k = 1; k <= in_nr; k++) {
                                        // this(k, r) * in(k, c)
                                        d += ppdColsm1[r][k]*pin->ppdColsm1[c][k];
                                }

                                // (out.*op)(r, c, d);
                                pout->ppdColsm1[c][r] -= d;
                        }
                }

        } else if (op == &MatrixHandler::PutCoef) {
                for (integer c = 1; c <= out_nc; c++) {
                        for (integer r = 1; r <= out_nr; r++) {
                                doublereal d = 0.;

                                for (integer k = 1; k <= in_nr; k++) {
                                        // this(k, r) * in(k, c)
                                        d += ppdColsm1[r][k]*pin->ppdColsm1[c][k];
                                }

                                // (out.*op)(r, c, d);
                                pout->ppdColsm1[c][r] = d;
                        }
                }

        } else {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif // ! USE_LAPACK

        return out;
}

Here is the call graph for this function:

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

protectedvirtual

Reimplemented from MatrixHandler.

Definition at line 963 of file fullmh.cc.

References ASSERT, VectorHandler::DecCoef(), iGetNumCols(), iGetNumRows(), VectorHandler::iGetSize(), VectorHandler::IncCoef(), MBDYN_EXCEPT_ARGS, MyVectorHandler::pdVecm1, ppdColsm1, and VectorHandler::PutCoef().

{
        integer nr = iGetNumCols();
        integer nc = iGetNumRows();

        ASSERT(nc == in.iGetSize());
        ASSERT(nr == out.iGetSize());

        const MyVectorHandler *pin = dynamic_cast<const MyVectorHandler *>(&in);
        if (pin == 0) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ic, ir) * in(ic)
                                d += ppdColsm1[ir][ic]*in(ic);
                        }
                        (out.*op)(ir, d);
                }

                return out;
        }

        MyVectorHandler *pout = dynamic_cast<MyVectorHandler *>(&out);
        if (pout == 0) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ic, ir) * in(ic)
                                d += ppdColsm1[ir][ic]*pin->pdVecm1[ic];
                        }
                        (out.*op)(ir, d);
                }

                return out;
        }

        if (op == &VectorHandler::IncCoef) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ic, ir) * in(ic)
                                d += ppdColsm1[ir][ic]*pin->pdVecm1[ic];
                        }
                        pout->pdVecm1[ir] += d;
                }

        } else if (op == &VectorHandler::DecCoef) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ic, ir) * in(ic)
                                d += ppdColsm1[ir][ic]*pin->pdVecm1[ic];
                        }
                        pout->pdVecm1[ir] -= d;
                }

        } else if (op == &VectorHandler::PutCoef) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ic, ir) * in(ic)
                                d += ppdColsm1[ir][ic]*pin->pdVecm1[ic];
                        }
                        pout->pdVecm1[ir] = d;
                }

        } else {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        return out;
}

Here is the call graph for this function:

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

protectedvirtual

Reimplemented from MatrixHandler.

Definition at line 887 of file fullmh.cc.

References ASSERT, VectorHandler::DecCoef(), iGetNumCols(), iGetNumRows(), VectorHandler::iGetSize(), VectorHandler::IncCoef(), MBDYN_EXCEPT_ARGS, MyVectorHandler::pdVecm1, ppdColsm1, and VectorHandler::PutCoef().

{
        integer nr = iGetNumRows();
        integer nc = iGetNumCols();

        ASSERT(nc == in.iGetSize());
        ASSERT(nr == out.iGetSize());

        const MyVectorHandler *pin = dynamic_cast<const MyVectorHandler *>(&in);
        if (pin == 0) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ir, ic) * in(ic)
                                d += ppdColsm1[ic][ir]*in(ic);
                        }
                        (out.*op)(ir, d);
                }

                return out;
        }

        MyVectorHandler *pout = dynamic_cast<MyVectorHandler *>(&out);
        if (pout == 0) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ir, ic) * in(ic)
                                d += ppdColsm1[ic][ir]*pin->pdVecm1[ic];
                        }
                        (out.*op)(ir, d);
                }

                return out;
        }

        if (op == &VectorHandler::IncCoef) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ir, ic) * in(ic)
                                d += ppdColsm1[ic][ir]*pin->pdVecm1[ic];
                        }
                        pout->pdVecm1[ir] += d;
                }

        } else if (op == &VectorHandler::DecCoef) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ir, ic) * in(ic)
                                d += ppdColsm1[ic][ir]*pin->pdVecm1[ic];
                        }
                        pout->pdVecm1[ir] -= d;
                }

        } else if (op == &VectorHandler::PutCoef) {
                for (integer ir = 1; ir <= nr; ir++) {
                        doublereal d = 0.;
                        for (integer ic = 1; ic <= nc; ic++) {
                                // out(ir) ? this(ir, ic) * in(ic)
                                d += ppdColsm1[ic][ir]*pin->pdVecm1[ic];
                        }
                        pout->pdVecm1[ir] = d;
                }

        } else {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        return out;
}

Here is the call graph for this function:

virtual doublereal& FullMatrixHandler::operator() ( integer  iRow, integer  iCol  )

inlinevirtual

Implements MatrixHandler.

Reimplemented in VariableSubMatrixHandler, and FullSubMatrixHandler.

Definition at line 234 of file fullmh.h.

References ASSERT, iNumCols, iNumRows, and ppdColsm1.

Referenced by FullSubMatrixHandler::operator()().

                                                 {
#ifdef DEBUG
                IsValid();
                ASSERT(iRow > 0);
                ASSERT(iRow <= iNumRows);
                ASSERT(iCol > 0);
                ASSERT(iCol <= iNumCols);
#endif /* DEBUG */

                return ppdColsm1[iCol][iRow];
        };

virtual const doublereal& FullMatrixHandler::operator() ( integer  iRow, integer  iCol  ) const

inlinevirtual

Implements MatrixHandler.

Reimplemented in VariableSubMatrixHandler, and FullSubMatrixHandler.

Definition at line 247 of file fullmh.h.

References ASSERT, iNumCols, iNumRows, and ppdColsm1.

                                                       {
#ifdef DEBUG
                IsValid();
                ASSERT(iRow > 0);
                ASSERT(iRow <= iNumRows);
                ASSERT(iCol > 0);
                ASSERT(iCol <= iNumCols);
#endif /* DEBUG */

                return ppdColsm1[iCol][iRow];
        };

MatrixHandler & FullMatrixHandler::operator+= ( const SubMatrixHandler &  SubMH )

virtual

Reimplemented from MatrixHandler.

Definition at line 405 of file fullmh.cc.

References SubMatrixHandler::AddTo().

{
        return SubMH.AddTo(*this);
}

Here is the call graph for this function:

MatrixHandler & FullMatrixHandler::operator+= ( const VariableSubMatrixHandler &  SubMH )

virtual

Reimplemented from MatrixHandler.

Definition at line 419 of file fullmh.cc.

References VariableSubMatrixHandler::AddTo().

{
        return SubMH.AddTo(*this);
}

Here is the call graph for this function:

MatrixHandler & FullMatrixHandler::operator-= ( const SubMatrixHandler &  SubMH )

virtual

Reimplemented from MatrixHandler.

Definition at line 412 of file fullmh.cc.

References SubMatrixHandler::SubFrom().

{
        return SubMH.SubFrom(*this);
}

Here is the call graph for this function:

MatrixHandler & FullMatrixHandler::operator-= ( const VariableSubMatrixHandler &  SubMH )

virtual

Reimplemented from MatrixHandler.

Definition at line 426 of file fullmh.cc.

References VariableSubMatrixHandler::SubFrom().

{
        return SubMH.SubFrom(*this);
}

Here is the call graph for this function:

FullMatrixHandler & FullMatrixHandler::operator=

(

const FullMatrixHandler &

MH

)

Definition at line 146 of file fullmh.cc.

References iNumCols, iNumRows, MBDYN_EXCEPT_ARGS, pdRaw, and Resize().

{
        if (pdRaw == 0) {
                Resize(MH.iNumRows, MH.iNumCols);

        } else if (iNumRows != MH.iNumRows || iNumCols != MH.iNumCols) {
                silent_cerr("FullMatrixHandler::operator = (const FullMatrixHandler&): "
                        "incompatible size (" << iNumRows << ", " << iNumCols << ") "
                        " <> (" << MH.iNumRows << ", " << MH.iNumCols << ")"
                        << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        for (integer i = 0; i < iNumRows*iNumCols; i++) {
                pdRaw[i] = MH.pdRaw[i];
        }

        return *this;
}

Here is the call graph for this function:

const doublereal* FullMatrixHandler::pdGetMat ( void  ) const

inlinevirtual

Reimplemented from MatrixHandler.

Definition at line 144 of file fullmh.h.

References ASSERT, and pdRaw.

Referenced by ModuleNonsmoothNode::mbs_get_force_frictional(), and DataManager::OutputEigFullMatrices().

                                               {
                ASSERT(pdRaw != NULL);
                return pdRaw;
        };

doublereal* FullMatrixHandler::pdGetMat ( void  )

inlinevirtual

Reimplemented from MatrixHandler.

Definition at line 149 of file fullmh.h.

References ASSERT, and pdRaw.

                                   {
                ASSERT(pdRaw != NULL);
                return pdRaw;
        };

virtual const doublereal* FullMatrixHandler::pdGetVec ( integer  iCol ) const

inlinevirtual

Definition at line 154 of file fullmh.h.

References ASSERT, iNumCols, pdRaw, and ppdColsm1.

                                                                      {
#ifdef DEBUG
                ASSERT(pdRaw != NULL);
                ASSERT(ppdColsm1 != NULL);
                ASSERT(iCol > 0 && iCol <= iNumCols);
                ASSERT(ppdColsm1[iCol] != NULL);
#endif /* DEBUG */

                return &ppdColsm1[iCol][1];
        };

virtual doublereal* FullMatrixHandler::pdGetVec ( integer  iCol )

inlinevirtual

Definition at line 165 of file fullmh.h.

References ASSERT, iNumCols, pdRaw, and ppdColsm1.

                                                          {
#ifdef DEBUG
                ASSERT(pdRaw != NULL);
                ASSERT(ppdColsm1 != NULL);
                ASSERT(iCol > 0 && iCol <= iNumCols);
                ASSERT(ppdColsm1[iCol] != NULL);
#endif /* DEBUG */

                return &ppdColsm1[iCol][1];
        };

void FullMatrixHandler::Put	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source
	)

Definition at line 1113 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

Referenced by Shell4EASANS::AssRes(), ModuleNonsmoothNode::mbs_get_force_frictional(), and FullSubMatrixHandler::Put().

{
        integer nr = source.iGetNumRows();
        integer nc = source.iGetNumCols();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ir, ic) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] = source.ppdColsm1[ic][ir];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::Put	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source,
		const doublereal	dCoef
	)

Definition at line 1226 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

{
        integer nr = source.iGetNumRows();
        integer nc = source.iGetNumCols();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ir, ic) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] = dCoef*source.ppdColsm1[ic][ir];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::Put	(	integer	iRow,
		integer	iCol,
		const Vec3 &	v
	)

Definition at line 1548 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, Vec3::pGetVec(), ppdColsm1, V1, V2, and V3.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols);
#endif /* DEBUG */

        /* Nota: assume che la matrice sia organizzata
         * "per colonne" (stile FORTRAN)
         */
        const doublereal* pdFrom = v.pGetVec();

        ppdColsm1[iCol][iRow] = pdFrom[V1];

        ppdColsm1[iCol][iRow + 1] = pdFrom[V2];

        ppdColsm1[iCol][iRow + 2] = pdFrom[V3];
}

Here is the call graph for this function:

void FullMatrixHandler::Put	(	integer	iRow,
		integer	iCol,
		const Mat3x3 &	m
	)

Definition at line 2005 of file fullmh.cc.

References ASSERT, c, iNumCols, iNumRows, M11, M12, M13, M21, M22, M23, M31, M32, M33, Mat3x3::pGetMat(), and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1
         */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        const doublereal* pdFrom = m.pGetMat();

        ppdColsm1[iCol][iRow] = pdFrom[M11];
        ppdColsm1[iCol][iRow + 1] = pdFrom[M21];
        ppdColsm1[iCol][iRow + 2] = pdFrom[M31];

        ppdColsm1[iCol + 1][iRow] = pdFrom[M12];
        ppdColsm1[iCol + 1][iRow + 1] = pdFrom[M22];
        ppdColsm1[iCol + 1][iRow + 2] = pdFrom[M32];

        ppdColsm1[iCol + 2][iRow] = pdFrom[M13];
        ppdColsm1[iCol + 2][iRow + 1] = pdFrom[M23];
        ppdColsm1[iCol + 2][iRow + 2] = pdFrom[M33];

#if 0
        for (unsigned c = 0; c < 2; c++) {
                for (unsigned r = 0; r < 2; r++) {
                        ppdColsm1[iCol + c][iRow + r] = pdFrom[r];
                }
                pdFrom += 3;
        }
#endif
}

Here is the call graph for this function:

void FullMatrixHandler::Put	(	integer	iRow,
		integer	iCol,
		const Mat3xN &	m
	)

Definition at line 2132 of file fullmh.cc.

References ASSERT, c, Mat3xN::iGetNumCols(), iNumCols, iNumRows, and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - m.iGetNumCols() + 1);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = 3; r > 0; r--) {
                for (integer c = m.iGetNumCols(); c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] = m(r, c);
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::Put	(	integer	iRow,
		integer	iCol,
		const MatNx3 &	m
	)

Definition at line 2280 of file fullmh.cc.

References ASSERT, c, MatNx3::iGetNumRows(), iNumCols, iNumRows, and ppdColsm1.

{
#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - m.iGetNumRows() + 1);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = m.iGetNumRows(); r > 0; r--) {
                for (integer c = 3; c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] = m(r, c);
                }
        }
}

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from MatrixHandler.

Reimplemented in FullSubMatrixHandler.

Definition at line 178 of file fullmh.h.

References ASSERT, iNumCols, iNumRows, and ppdColsm1.

Referenced by fmu2::GetDirectionalDerivatives(), FullSubMatrixHandler::PutCoef(), and ReadAerodynamicModal().

                                                                     {
#ifdef DEBUG
                IsValid();
                ASSERT(iRow > 0 && iRow <= iNumRows);
                ASSERT(iCol > 0 && iCol <= iNumCols);
#endif /* DEBUG */

                ppdColsm1[iCol][iRow] = dCoef;
        };

void FullMatrixHandler::PutT	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source
	)

Definition at line 1338 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

Referenced by Shell4EASANS::AssRes(), ModuleNonsmoothNode::mbs_get_force_frictional(), and FullSubMatrixHandler::PutT().

{
        integer nr = source.iGetNumCols();
        integer nc = source.iGetNumRows();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ic, ir) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] = source.ppdColsm1[ir][ic];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::PutT	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source,
		const doublereal	dCoef
	)

Definition at line 1451 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

{
        integer nr = source.iGetNumCols();
        integer nc = source.iGetNumRows();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ic, ir) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] = dCoef*source.ppdColsm1[ir][ic];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::PutT	(	integer	iRow,
		integer	iCol,
		const Vec3 &	v
	)

Definition at line 1635 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, Vec3::pGetVec(), ppdColsm1, V1, V2, and V3.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        /* Nota: assume che la matrice sia organizzata
         * "per colonne" (stile FORTRAN)
         */
        const doublereal* pdFrom = v.pGetVec();

        ppdColsm1[iCol][iRow] = pdFrom[V1];

        ppdColsm1[iCol + 1][iRow] = pdFrom[V2];

        ppdColsm1[iCol + 2][iRow] = pdFrom[V3];
}

Here is the call graph for this function:

void FullMatrixHandler::PutT	(	integer	iRow,
		integer	iCol,
		const Mat3x3 &	m
	)

Definition at line 2046 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, M11, M12, M13, M21, M22, M23, M31, M32, M33, Mat3x3::pGetMat(), and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1
         */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        const doublereal* pdFrom = m.pGetMat();

        ppdColsm1[iCol][iRow] = pdFrom[M11];
        ppdColsm1[iCol][iRow + 1] = pdFrom[M12];
        ppdColsm1[iCol][iRow + 2] = pdFrom[M13];

        ppdColsm1[iCol + 1][iRow] = pdFrom[M21];
        ppdColsm1[iCol + 1][iRow + 1] = pdFrom[M22];
        ppdColsm1[iCol + 1][iRow + 2] = pdFrom[M23];

        ppdColsm1[iCol + 2][iRow] = pdFrom[M31];
        ppdColsm1[iCol + 2][iRow + 1] = pdFrom[M32];
        ppdColsm1[iCol + 2][iRow + 2] = pdFrom[M33];
}

Here is the call graph for this function:

void FullMatrixHandler::PutT	(	integer	iRow,
		integer	iCol,
		const Mat3xN &	m
	)

Definition at line 2210 of file fullmh.cc.

References ASSERT, c, Mat3xN::iGetNumCols(), iNumCols, iNumRows, and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - m.iGetNumCols() + 1);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = m.iGetNumCols(); r > 0; r--) {
                for (integer c = 3; c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] = m(c, r);
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::Reset ( void  )

virtual

Implements MatrixHandler.

Reimplemented in FullSubMatrixHandler.

Definition at line 44 of file fullmh.cc.

References c, iNumCols, iNumRows, and ppdColsm1.

Referenced by ModuleFMU::AssJac(), Shell4EASANS::AssRes(), FullMatrixHandler(), main(), ReadMembraneConstLaw(), ReadShellConstLaw(), and FullSubMatrixHandler::Reset().

{
#ifdef DEBUG
        IsValid();
#endif /* DEBUG */

        for (integer c = iNumCols; c > 0; c--) {
                for (integer r = iNumRows; r > 0; r--) {
                        ppdColsm1[c][r] = 0.;
                }
        }
}

void FullMatrixHandler::Resize ( integer  iNewRows, integer  iNewCols  )

virtual

Implements MatrixHandler.

Reimplemented in FullSubMatrixHandler.

Definition at line 174 of file fullmh.cc.

References bOwnsMemory, CreateColRow(), iMaxCols, iNumCols, iNumRows, iRawSize, MBDYN_EXCEPT_ARGS, pdRaw, pdRawm1, ppdCols, ppdColsm1, SAFEDELETEARR, and SAFENEWARR.

Referenced by ModuleFMU::AssJac(), FullMatrixHandler(), ModuleNonsmoothNode::ModuleNonsmoothNode(), operator=(), FullSubMatrixHandler::Resize(), and Gmres::Solve().

{
        if (iNewRows < 0 || iNewCols < 0) {
                silent_cerr("Negative size!" << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        integer iSize = iNewRows*iNewCols;
        if (bOwnsMemory) {
                if (pdRaw != NULL) {
                        if (iSize > iRawSize) {
                                /* crea */
                                doublereal* pd = NULL;
                                SAFENEWARR(pd, doublereal, iSize);
                                for (integer i = iRawSize; i-- > 0; ) {
                                        pd[i] = pdRaw[i];
                                }
                                iRawSize = iSize;
                                SAFEDELETEARR(pdRaw);
                                pdRaw = pd;
                                pdRawm1 = pd - 1;
                        }

                        if (iNewCols > iMaxCols) {
                                /* crea */
                                SAFEDELETEARR(ppdCols);
                                SAFENEWARR(ppdCols, doublereal*, iNewCols);
                                iMaxCols = iNewCols;
                                ppdColsm1 = ppdCols-1;
                                CreateColRow(iNewRows, iNewCols);
                        }

                        /* aggiorna iNumCols, iNumRows */
                        if (iNewRows != iNumRows || iNewCols != iNumCols) {
                                CreateColRow(iNewRows, iNewCols);
                        }

                        iNumRows = iNewRows;
                        iNumCols = iNewCols;

                } else {
                        /* crea tutto */
                        iRawSize = iSize;
                        iMaxCols = iNewCols;
                        iNumRows = iNewRows;
                        iNumCols = iNewCols;
                        SAFENEWARR(pdRaw, doublereal, iRawSize);
                        pdRawm1 = pdRaw-1;
                        SAFENEWARR(ppdCols, doublereal*, iNumCols);
                        ppdColsm1 = ppdCols-1;
                        CreateColRow(iNumRows, iNumCols);
                }

        } else {
                if (pdRaw != NULL) {
                        if (iSize > iRawSize || iNewCols > iMaxCols) {
                                if (iSize > iRawSize) {
                                        /* errore */
                                        silent_cerr("Can't resize to "
                                                << iSize
                                                << ": larger than max size "
                                                << iRawSize << std::endl);
                                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);

                                } else if (iNewCols > iMaxCols) {
                                        silent_cerr("Can't resize to "
                                                << iNewCols
                                                << " cols: larger than max "
                                                "num cols "
                                                << iMaxCols << std::endl);
                                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                }

                        } else {
                                /* aggiorna */
                                if (iNewRows != iNumRows ||
                                                iNewCols != iNumCols) {
                                        CreateColRow(iNewRows, iNewCols);
                                }
                                iNumRows = iNewRows;
                                iNumCols = iNewCols;
                        }

                } else {
                        /* errore */
                        silent_cerr("internal error!" << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::Sub	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source
	)

Definition at line 1076 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

Referenced by ModuleNonsmoothNode::mbs_get_force_frictional(), and FullSubMatrixHandler::Sub().

{
        integer nr = source.iGetNumRows();
        integer nc = source.iGetNumCols();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ir, ic) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] -= source.ppdColsm1[ic][ir];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::Sub	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source,
		const doublereal	dCoef
	)

Definition at line 1188 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

{
        integer nr = source.iGetNumRows();
        integer nc = source.iGetNumCols();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ir, ic) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] -= dCoef*source.ppdColsm1[ic][ir];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::Sub	(	integer	iRow,
		integer	iCol,
		const Vec3 &	v
	)

Definition at line 1519 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, Vec3::pGetVec(), ppdColsm1, V1, V2, and V3.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols);
#endif /* DEBUG */

        /* Nota: assume che la matrice sia organizzata
         * "per colonne" (stile FORTRAN)
         */
        const doublereal* pdFrom = v.pGetVec();

        ppdColsm1[iCol][iRow] -= pdFrom[V1];

        ppdColsm1[iCol][iRow + 1] -= pdFrom[V2];

        ppdColsm1[iCol][iRow + 2] -= pdFrom[V3];
}

Here is the call graph for this function:

void FullMatrixHandler::Sub	(	integer	iRow,
		integer	iCol,
		const Mat3x3 &	m
	)

Definition at line 1929 of file fullmh.cc.

References ASSERT, c, iNumCols, iNumRows, M11, M12, M13, M21, M22, M23, M31, M32, M33, Mat3x3::pGetMat(), and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1
         */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        const doublereal* pdFrom = m.pGetMat();

        ppdColsm1[iCol][iRow] -= pdFrom[M11];
        ppdColsm1[iCol][iRow + 1] -= pdFrom[M21];
        ppdColsm1[iCol][iRow + 2] -= pdFrom[M31];

        ppdColsm1[iCol + 1][iRow] -= pdFrom[M12];
        ppdColsm1[iCol + 1][iRow + 1] -= pdFrom[M22];
        ppdColsm1[iCol + 1][iRow + 2] -= pdFrom[M32];

        ppdColsm1[iCol + 2][iRow] -= pdFrom[M13];
        ppdColsm1[iCol + 2][iRow + 1] -= pdFrom[M23];
        ppdColsm1[iCol + 2][iRow + 2] -= pdFrom[M33];

#if 0
        for (unsigned c = 0; c < 2; c++) {
                for (unsigned r = 0; r < 2; r++) {
                        ppdColsm1[iCol + c][iRow + r] -= pdFrom[r];
                }
                pdFrom += 3;
        }
#endif
}

Here is the call graph for this function:

void FullMatrixHandler::Sub	(	integer	iRow,
		integer	iCol,
		const Mat3xN &	m
	)

Definition at line 2106 of file fullmh.cc.

References ASSERT, c, Mat3xN::iGetNumCols(), iNumCols, iNumRows, and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - m.iGetNumCols() + 1);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = 3; r > 0; r--) {
                for (integer c = m.iGetNumCols(); c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] -= m(r, c);
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::Sub	(	integer	iRow,
		integer	iCol,
		const MatNx3 &	m
	)

Definition at line 2258 of file fullmh.cc.

References ASSERT, c, MatNx3::iGetNumRows(), iNumCols, iNumRows, and ppdColsm1.

{
#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - m.iGetNumRows() + 1);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = m.iGetNumRows(); r > 0; r--) {
                for (integer c = 3; c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] -= m(r, c);
                }
        }
}

Here is the call graph for this function:

void FullMatrixHandler::SubT	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source
	)

Definition at line 1301 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

Referenced by FullSubMatrixHandler::SubT().

{
        integer nr = source.iGetNumCols();
        integer nc = source.iGetNumRows();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ic, ir) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] -= source.ppdColsm1[ir][ic];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::SubT	(	integer	iRow,
		integer	iCol,
		const FullMatrixHandler &	source,
		const doublereal	dCoef
	)

Definition at line 1413 of file fullmh.cc.

References ASSERT, iGetNumCols(), iGetNumRows(), MBDYN_EXCEPT_ARGS, and ppdColsm1.

{
        integer nr = source.iGetNumCols();
        integer nc = source.iGetNumRows();
        iRow--;
        iCol--;

#ifdef DEBUG
        IsValid();

        ASSERT(iRow >= 0);
        ASSERT(iCol >= 0);
        ASSERT(iRow + nr <= iGetNumRows());
        ASSERT(iCol + nc <= iGetNumCols());
#endif // DEBUG
#if 0
        if (iRow < 0 || iCol < 0
                || iRow + nr > iGetNumRows()
                || iCol + nc > iGetNumCols())
        {
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif

        for (integer ir = 1; ir <= nr; ir++) {
                for (integer ic = 1; ic <= nc; ic++) {
                        // IncCoef(ir + iRow, ic + iCol, source(ic, ir) * dCoef);
                        ppdColsm1[ic + iCol][ir + iRow] -= dCoef*source.ppdColsm1[ir][ic];
                }
        }

        return;
}

Here is the call graph for this function:

void FullMatrixHandler::SubT	(	integer	iRow,
		integer	iCol,
		const Vec3 &	v
	)

Definition at line 1606 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, Vec3::pGetVec(), ppdColsm1, V1, V2, and V3.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        /* Nota: assume che la matrice sia organizzata
         * "per colonne" (stile FORTRAN)
         */
        const doublereal* pdFrom = v.pGetVec();

        ppdColsm1[iCol][iRow] -= pdFrom[V1];

        ppdColsm1[iCol + 1][iRow] -= pdFrom[V2];

        ppdColsm1[iCol + 2][iRow] -= pdFrom[V3];
}

Here is the call graph for this function:

void FullMatrixHandler::SubT	(	integer	iRow,
		integer	iCol,
		const Mat3x3 &	m
	)

Definition at line 1970 of file fullmh.cc.

References ASSERT, iNumCols, iNumRows, M11, M12, M13, M21, M22, M23, M31, M32, M33, Mat3x3::pGetMat(), and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1
         */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - 2);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        const doublereal* pdFrom = m.pGetMat();

        ppdColsm1[iCol][iRow] -= pdFrom[M11];
        ppdColsm1[iCol][iRow + 1] -= pdFrom[M12];
        ppdColsm1[iCol][iRow + 2] -= pdFrom[M13];

        ppdColsm1[iCol + 1][iRow] -= pdFrom[M21];
        ppdColsm1[iCol + 1][iRow + 1] -= pdFrom[M22];
        ppdColsm1[iCol + 1][iRow + 2] -= pdFrom[M23];

        ppdColsm1[iCol + 2][iRow] -= pdFrom[M31];
        ppdColsm1[iCol + 2][iRow + 1] -= pdFrom[M32];
        ppdColsm1[iCol + 2][iRow + 2] -= pdFrom[M33];
}

Here is the call graph for this function:

void FullMatrixHandler::SubT	(	integer	iRow,
		integer	iCol,
		const Mat3xN &	m
	)

Definition at line 2184 of file fullmh.cc.

References ASSERT, c, Mat3xN::iGetNumCols(), iNumCols, iNumRows, and ppdColsm1.

{
        /* iRow e iCol sono gli indici effettivi di riga e colonna
         * es. per il primo coefficiente:
         *     iRow = 1, iCol = 1 */

#ifdef DEBUG
        IsValid();

        ASSERT(iRow > 0);
        ASSERT(iRow <= iNumRows - m.iGetNumCols() + 1);
        ASSERT(iCol > 0);
        ASSERT(iCol <= iNumCols - 2);
#endif /* DEBUG */

        --iRow;
        --iCol;
        for (int r = m.iGetNumCols(); r > 0; r--) {
                for (integer c = 3; c > 0; c--) {
                        ppdColsm1[iCol + c][iRow + r] -= m(c, r);
                }
        }
}

Here is the call graph for this function:

Friends And Related Function Documentation

friend class FullSubMatrixHandler

friend

Definition at line 50 of file fullmh.h.

std::ostream& operator<< ( std::ostream &  out, const FullMatrixHandler &  m  )

friend

Definition at line 352 of file fullmh.cc.

{
#ifdef HAVE_FMTFLAGS_IN_IOS
        std::ios::fmtflags oldbits = out.setf(std::ios::scientific);
#else /* !HAVE_FMTFLAGS_IN_IOS */
        long oldbits = out.setf(ios::scientific);
#endif /* !HAVE_FMTFLAGS_IN_IOS */

        out << "<FullMatrixHandler> n. rows: " << m.iNumRows
                << "; n. cols: " << m.iNumCols << std::endl;
        for (int i = 1; i <= m.iNumRows; i++) {
                out << "Row " << std::setw(8) << i;
                for (int j = 1; j <= m.iNumCols; j++) {
                        out << std::setw(10) << std::setprecision(2)
                                << m.ppdColsm1[j][i];
                }
                out << std::endl;
        }

        out.flags(oldbits);
        return out;
}

friend class SparseSubMatrixHandler

friend

Definition at line 51 of file fullmh.h.

std::ostream& Write ( std::ostream &  out, const FullMatrixHandler &  m, const char *  s, const char *  s2  )

friend

Definition at line 376 of file fullmh.cc.

{
#ifdef HAVE_FMTFLAGS_IN_IOS
        std::ios::fmtflags oldbits = out.setf(std::ios::scientific);
#else /* !HAVE_FMTFLAGS_IN_IOS */
        long oldbits = out.setf(ios::scientific);
#endif /* !HAVE_FMTFLAGS_IN_IOS */

        if (s2 == NULL) {
                s2 = s;
        }

        for (int i = 1; i <= m.iNumRows; i++) {
                for (int j = 1; j <= m.iNumCols; j++) {
                        out << std::setw(20) << std::setprecision(12)
                                << m.ppdColsm1[j][i] << s;
                }
                out << s2;
        }

        out.flags(oldbits);
        return out;
}

Member Data Documentation

bool FullMatrixHandler::bOwnsMemory

protected

Definition at line 61 of file fullmh.h.

Referenced by Attach(), FullSubMatrixHandler::Attach(), Detach(), FullSubMatrixHandler::FullSubMatrixHandler(), Resize(), and FullSubMatrixHandler::~FullSubMatrixHandler().

integer FullMatrixHandler::iMaxCols

protected

Definition at line 67 of file fullmh.h.

Referenced by Attach(), Detach(), FullMatrixHandler(), Resize(), and FullSubMatrixHandler::Resize().

integer FullMatrixHandler::iNumCols

protected

Definition at line 64 of file fullmh.h.

Referenced by Add(), AddT(), FullSubMatrixHandler::AddTo(), FullSubMatrixHandler::AddToT(), Attach(), FullSubMatrixHandler::Attach(), DecCoef(), Detach(), dGetCoef(), FullMatrixHandler(), FullSubMatrixHandler::iGetColIndex(), iGetNumCols(), FullSubMatrixHandler::iGetNumCols(), IncCoef(), MatMul(), operator()(), operator<<(), operator=(), pdGetVec(), Put(), PutCoef(), FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutCross(), FullSubMatrixHandler::PutDiag(), PutT(), FullMatrixHandler::const_iterator::reset(), Reset(), FullSubMatrixHandler::Reset(), Resize(), FullSubMatrixHandler::Resize(), Sub(), FullSubMatrixHandler::SubFrom(), FullSubMatrixHandler::SubFromT(), SubT(), and Write().

integer FullMatrixHandler::iNumRows

protected

Definition at line 63 of file fullmh.h.

Referenced by Add(), AddT(), FullSubMatrixHandler::AddTo(), FullSubMatrixHandler::AddToT(), Attach(), FullSubMatrixHandler::Attach(), DecCoef(), Detach(), dGetCoef(), FullMatrixHandler(), FullSubMatrixHandler::FullSubMatrixHandler(), iGetNumRows(), FullSubMatrixHandler::iGetNumRows(), FullSubMatrixHandler::iGetRowIndex(), IncCoef(), MatMul(), operator()(), operator<<(), operator=(), Put(), PutCoef(), FullSubMatrixHandler::PutCross(), FullSubMatrixHandler::PutDiag(), FullSubMatrixHandler::PutRowIndex(), PutT(), FullMatrixHandler::const_iterator::reset(), Reset(), FullSubMatrixHandler::Reset(), Resize(), FullSubMatrixHandler::Resize(), Sub(), FullSubMatrixHandler::SubFrom(), FullSubMatrixHandler::SubFromT(), SubT(), and Write().

integer FullMatrixHandler::iRawSize

protected

Definition at line 66 of file fullmh.h.

Referenced by Attach(), Detach(), Resize(), and FullSubMatrixHandler::Resize().

const_iterator FullMatrixHandler::m_end

private

Definition at line 99 of file fullmh.h.

Referenced by end().

doublereal* FullMatrixHandler::pdRaw

protected

Definition at line 69 of file fullmh.h.

Referenced by Attach(), Detach(), FullMatrixHandler(), MatMatMul_base(), MatTMatMul_base(), operator=(), pdGetMat(), pdGetVec(), FullMatrixHandler::const_iterator::reset(), and Resize().

doublereal* FullMatrixHandler::pdRawm1

protected

Definition at line 70 of file fullmh.h.

Referenced by Attach(), CreateColRow(), Detach(), FullMatrixHandler(), and Resize().

doublereal** FullMatrixHandler::ppdCols

protected

Definition at line 71 of file fullmh.h.

Referenced by Attach(), Detach(), NaiveMatrixHandler::operator+=(), NaiveMatrixHandler::operator-=(), and Resize().

doublereal** FullMatrixHandler::ppdColsm1

protected

Definition at line 72 of file fullmh.h.

Referenced by Add(), AddT(), FullSubMatrixHandler::AddTo(), SparseSubMatrixHandler::AddTo(), FullSubMatrixHandler::AddToT(), SparseSubMatrixHandler::AddToT(), Attach(), CopyMatrixBlock(), CopyMatrixRow(), CreateColRow(), DecCoef(), FullSubMatrixHandler::DecCoef(), Detach(), dGetCoef(), FullSubMatrixHandler::dGetCoef(), FullMatrixHandler(), IncCoef(), FullSubMatrixHandler::IncCoef(), MatMatMul_base(), MatMul(), MatTMatMul_base(), MatTVecMul_base(), MatVecMul_base(), operator()(), operator<<(), pdGetVec(), Put(), PutCoef(), FullSubMatrixHandler::PutCoef(), FullSubMatrixHandler::PutCross(), FullSubMatrixHandler::PutDiag(), PutT(), Reset(), Resize(), Sub(), FullSubMatrixHandler::SubFrom(), SparseSubMatrixHandler::SubFrom(), FullSubMatrixHandler::SubFromT(), SparseSubMatrixHandler::SubFromT(), SubT(), and Write().

---

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

• fullmh.h
• fullmh.cc