#include <spmapmh.h>

Inheritance diagram for SpMapMatrixHandler:

Collaboration diagram for SpMapMatrixHandler:

Classes
class	const_iterator

Public Member Functions
SpMapMatrixHandler::const_iterator	begin (void) const

const SpMapMatrixHandler::const_iterator &	end (void) const

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

virtual	~SpMapMatrixHandler (void)

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

void	IncCoef (integer ix, integer iy, const doublereal &inc)

void	DecCoef (integer ix, integer iy, const doublereal &inc)

void	PutCoef (integer ix, integer iy, const doublereal &val)

const doublereal &	dGetCoef (integer ix, integer iy) const

const doublereal &	operator() (integer ix, integer iy) const

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

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

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

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

void	Reset (void)

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 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 const doublereal *	pdGetMat (void) const

virtual doublereal *	pdGetMat (void)

virtual integer *	piGetRows (void) const

virtual integer *	piGetCols (void) const

virtual integer	PacMat (void)

virtual MatrixHandler &	operator= (const MatrixHandler &MH)

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

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

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

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

virtual MatrixHandler &	ScalarMul (const doublereal &d)

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

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

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

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

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

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

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

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

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

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

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

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

virtual doublereal	ConditionNumber (enum Norm_t eNorm=NORM_1) const

virtual doublereal	Norm (enum Norm_t eNorm=NORM_1) const

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

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

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

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

Private Types
typedef std::map< integer, doublereal >	row_cont_type

Private Member Functions
	SpMapMatrixHandler (const SpMapMatrixHandler &)

Private Attributes
std::vector< row_cont_type >	col_indices

const_iterator	m_end

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

Protected Attributes inherited from SparseMatrixHandler
integer	NRows

integer	NCols

integer	NZ

Detailed Description

Definition at line 76 of file spmapmh.h.

Member Typedef Documentation

typedef std::map<integer, doublereal> SpMapMatrixHandler::row_cont_type

private

Definition at line 78 of file spmapmh.h.

Constructor & Destructor Documentation

SpMapMatrixHandler::SpMapMatrixHandler ( const SpMapMatrixHandler & )

private

SpMapMatrixHandler::SpMapMatrixHandler	(	const integer &	n = `0`,
		const integer &	nn = `0`
	)

Definition at line 68 of file spmapmh.cc.

References col_indices, and SparseMatrixHandler::NCols.

 : SparseMatrixHandler(n, nn), m_end(*this, true)
 {
         col_indices.resize(NCols);
 }

SpMapMatrixHandler::~SpMapMatrixHandler ( void )

virtual

Definition at line 74 of file spmapmh.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

SpMapMatrixHandler::const_iterator SpMapMatrixHandler::begin ( void ) const

inline

Definition at line 116 of file spmapmh.h.

Referenced by FakeThirdOrderMulAndSumWithShift(), GetCol(), main(), MakeCompressedColumnForm(), MakeIndexForm(), MatMatMul_base(), MatTMatMul_base(), MatTVecMul_base(), MatVecMul_base(), MulAndSumWithShift(), DataManager::OutputEigSparseMatrices(), and Reset().

                                                            {
                 return const_iterator(*this);
         };

void SpMapMatrixHandler::DecCoef	(	integer	ix,
		integer	iy,
		const doublereal &	inc
	)

inlinevirtual

Reimplemented from MatrixHandler.

Definition at line 170 of file spmapmh.h.

References ASSERTMSGBREAK, SparseMatrixHandler::NCols, SparseMatrixHandler::NRows, and operator()().

Referenced by SchurMatrixHandler::DecCoef(), and SchurMatrixHandlerUm::DecCoef().

                                                                     {
                 ASSERTMSGBREAK(ix > 0 && ix <= NRows,
                                 "Error in SpMapMatrixHandler::DecCoef(), "
                                 "row index out of range");
                 ASSERTMSGBREAK(iy > 0 && iy <= NCols,
                                 "Error in SpMapMatrixHandler::DecCoef(), "
                                 "col index out of range");
 #ifdef MBDYN_X_KEEP_SPARSITY
                 /* try to keep sparsity */
                 if (inc != 0.) {
 #endif /* MBDYN_X_KEEP_SPARSITY */
                         operator()(ix,iy) -= inc;
 #ifdef MBDYN_X_KEEP_SPARSITY
                 }
 #endif /* MBDYN_X_KEEP_SPARSITY */
         };

Here is the call graph for this function:

const doublereal& SpMapMatrixHandler::dGetCoef	(	integer	ix,
		integer	iy
	)		const

inlinevirtual

Reimplemented from MatrixHandler.

Definition at line 211 of file spmapmh.h.

References ASSERTMSGBREAK, col_indices, SparseMatrixHandler::NCols, SparseMatrixHandler::NRows, and Zero1.

                                                                  {
                 ASSERTMSGBREAK(ix > 0 && ix <= NRows,
                                 "Error in SpMapMatrixHandler::dGetCoef(), "
                                 "row index out of range");
                 ASSERTMSGBREAK(iy > 0 && iy <= NCols,
                                 "Error in SpMapMatrixHandler::dGetCoef(), "
                                 "col index out of range");
                 row_cont_type::iterator i;
                 row_cont_type& row = col_indices[iy - 1];
                 i = row.find(ix - 1);
                 if (i == row.end()) {
                         return ::Zero1;
 
                 } else {
                         return i->second;
                 }
         };

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

inline

Definition at line 120 of file spmapmh.h.

References m_end.

Referenced by main(), and DataManager::OutputEigSparseMatrices().

                                                                 {
                 return m_end;
         };

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

Definition at line 296 of file spmapmh.cc.

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

Referenced by ThirdOrderIntegrator::Jacobian().

 {
         if ((out.iGetNumCols() < iGetNumCols() + dcol)
                         || (out.iGetNumRows() < iGetNumRows() + drow))
         {
                 silent_cerr("Assertion fault "
                         "in SpMapMatrixHandler::MulAndSumWithShift"
                         << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         drow = drow + 1;
 
         for (integer col = 0; col < NCols; col++) {
                 row_cont_type::const_iterator ri, re;
                 re = col_indices[col].end();
                 integer newcol = col + dcol + 1;
                 for (ri = col_indices[col].begin(); ri != re; ++ri) {
                         if (b[ri->first]) {
                                 out.IncCoef(ri->first + drow,
                                                 newcol, ri->second*s);
                         }
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

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

virtual

Implements SparseMatrixHandler.

Definition at line 194 of file spmapmh.cc.

References begin(), col_indices, SparseMatrixHandler::iGetNumCols(), and MBDYN_EXCEPT_ARGS.

 {
         if (icol > iGetNumCols()) {
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         row_cont_type::const_iterator ri, re;
         re = col_indices[icol].end();
         for (ri = col_indices[icol].begin(); ri != re; ++ri) {
                 out(ri->first + 1) = ri->second;
         }
         return out;
 }

Here is the call graph for this function:

void SpMapMatrixHandler::IncCoef	(	integer	ix,
		integer	iy,
		const doublereal &	inc
	)

inlinevirtual

Reimplemented from MatrixHandler.

Definition at line 153 of file spmapmh.h.

References ASSERTMSGBREAK, SparseMatrixHandler::NCols, SparseMatrixHandler::NRows, and operator()().

Referenced by SchurMatrixHandler::IncCoef(), and SchurMatrixHandlerUm::IncCoef().

                                                                     {
                 ASSERTMSGBREAK(ix > 0 && ix <= NRows,
                                 "Error in SpMapMatrixHandler::IncCoef(), "
                                 "row index out of range");
                 ASSERTMSGBREAK(iy > 0 && iy <= NCols,
                                 "Error in SpMapMatrixHandler::IncCoef(), "
                                 "col index out of range");
 #ifdef MBDYN_X_KEEP_SPARSITY
                 /* try to keep sparsity */
                 if (inc != 0.) {
 #endif /* MBDYN_X_KEEP_SPARSITY */
                         operator()(ix,iy) += inc;
 #ifdef MBDYN_X_KEEP_SPARSITY
                 }
 #endif /* MBDYN_X_KEEP_SPARSITY */
         };

Here is the call graph for this function:

integer SpMapMatrixHandler::MakeCompressedColumnForm	(	doublereal *const	Ax,
		integer *const	Ai,
		integer *const	Ap,
		int	offset = `0`
	)		const

virtual

Implements SparseMatrixHandler.

Definition at line 80 of file spmapmh.cc.

References ASSERTMSGBREAK, begin(), col_indices, SparseMatrixHandler::NCols, SparseMatrixHandler::NZ, and SparseMatrixHandler::Nz().

Referenced by main(), and MakeCompressedColumnForm().

 {
         integer x_ptr = 0;
 
         row_cont_type::iterator ri;
         row_cont_type::const_iterator re;
 
         for (integer col = 0; col < NCols; col++) {
                 Ap[col] = x_ptr + offset;
                 re = col_indices[col].end();
                 for (ri = col_indices[col].begin(); ri != re; ++ri) {
                         Ax[x_ptr] = ri->second;
                         Ai[x_ptr] = ri->first + offset;
                         x_ptr++;
                 }
         }
 
         ASSERTMSGBREAK(x_ptr == NZ, "Error in "
                         "SpMapMatrixHandler::MakeCompressedColumnForm");
 
         Ap[NCols] = x_ptr + offset;
 
         return Nz();
 }

Here is the call graph for this function:

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

virtual

Implements SparseMatrixHandler.

Definition at line 107 of file spmapmh.cc.

References SparseMatrixHandler::iGetNumCols(), MakeCompressedColumnForm(), and SparseMatrixHandler::Nz().

 {
         Ax.resize(Nz());
         Ai.resize(Nz());
         Ap.resize(iGetNumCols() + 1);
 
         return MakeCompressedColumnForm(&Ax[0], &Ai[0], &Ap[0], offset);
 }

Here is the call graph for this function:

integer SpMapMatrixHandler::MakeIndexForm	(	doublereal *const	Ax,
		integer *const	Arow,
		integer *const	Acol,
		integer *const	AcolSt,
		int	offset = `0`
	)		const

virtual

Implements SparseMatrixHandler.

Definition at line 119 of file spmapmh.cc.

References ASSERTMSGBREAK, begin(), col_indices, SparseMatrixHandler::NCols, SparseMatrixHandler::NZ, and SparseMatrixHandler::Nz().

Referenced by MakeIndexForm().

 {
         integer x_ptr = 0;
 
         row_cont_type::iterator ri;
         row_cont_type::const_iterator re;
 
         for (integer col = 0; col < NCols; col++) {
                 Ap[col] = x_ptr + offset;
                 re = col_indices[col].end();
                 for (ri = col_indices[col].begin(); ri != re; ++ri) {
                         Ax[x_ptr] = ri->second;
                         Arow[x_ptr] = ri->first + offset;
                         Acol[x_ptr] = col + offset;
                         x_ptr++;
                 }
         }
 
         ASSERTMSGBREAK(x_ptr == NZ, "Error in "
                         "SpMapMatrixHandler::MakeIndexForm");
 
         Ap[NCols] = x_ptr + offset;
 
         return Nz();
 }

Here is the call graph for this function:

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

virtual

Implements SparseMatrixHandler.

Definition at line 148 of file spmapmh.cc.

References SparseMatrixHandler::iGetNumCols(), MakeIndexForm(), and SparseMatrixHandler::Nz().

 {
         Ax.resize(Nz());
         Arow.resize(Nz());
         Acol.resize(Nz());
         Ap.resize(iGetNumCols() + 1);
 
         return MakeIndexForm(&Ax[0], &Arow[0], &Acol[0], &Ap[0], offset);
 }

Here is the call graph for this function:

MatrixHandler & SpMapMatrixHandler::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 210 of file spmapmh.cc.

References begin(), col_indices, SparseMatrixHandler::iGetNumCols(), MatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), MatrixHandler::iGetNumRows(), MBDYN_EXCEPT_ARGS, and SparseMatrixHandler::NCols.

 {
         if ((in.iGetNumRows() != iGetNumCols())
                         || (in.iGetNumCols() != out.iGetNumCols())
                         || (out.iGetNumRows() != iGetNumRows()))
         {
                 silent_cerr("Assertion fault "
                         "in SpMapMatrixHandler::MatMatIncMul" << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         integer ncols_in = in.iGetNumCols();
         for (integer row_in = 0; row_in < NCols; row_in++) {
                 row_cont_type::const_iterator ri, re;
                 re = col_indices[row_in].end();
                 for (ri = col_indices[row_in].begin(); ri != re; ++ri) {
                         for (integer col_in = 1; col_in <= ncols_in;  col_in++) {
                                 (out.*op)(ri->first + 1, col_in,
                                                 ri->second*in(row_in + 1, col_in));
                         }
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

MatrixHandler & SpMapMatrixHandler::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 239 of file spmapmh.cc.

References begin(), col_indices, SparseMatrixHandler::iGetNumCols(), MatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), MatrixHandler::iGetNumRows(), MBDYN_EXCEPT_ARGS, and SparseMatrixHandler::NCols.

 {
         if ((in.iGetNumRows() != iGetNumCols())
                         || (in.iGetNumCols() != out.iGetNumCols())
                         || (out.iGetNumRows() != iGetNumRows()))
         {
                 silent_cerr("Assertion fault "
                         "in SpMapMatrixHandler::MatTMatMul_base" << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         integer ncols_in = in.iGetNumCols();
         for (integer row_out = 0; row_out < NCols; row_out++) {
                 row_cont_type::const_iterator ri, re;
                 re = col_indices[row_out].end();
                 for (ri = col_indices[row_out].begin(); ri != re; ++ri) {
                         for (integer col_in = 1; col_in <= ncols_in; col_in++) {
                                 (out.*op)(row_out + 1, col_in,
                                                 ri->second*in(ri->first + 1, col_in));
                         }
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

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

protectedvirtual

Reimplemented from MatrixHandler.

Definition at line 351 of file spmapmh.cc.

References begin(), col_indices, SparseMatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), VectorHandler::iGetSize(), MBDYN_EXCEPT_ARGS, and SparseMatrixHandler::NCols.

 {
         if (out.iGetSize() != iGetNumCols()
                         || in.iGetSize() != iGetNumRows())
         {
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         row_cont_type::const_iterator ri, re;
 
         for (integer col = 0; col < NCols; col++) {
                 doublereal d = 0.;
                 re = col_indices[col].end();
                 for (ri = col_indices[col].begin(); ri != re; ++ri) {
                         d += ri->second*in(ri->first + 1);
                 }
                 (out.*op)(col+1, d);
         }
 
         return out;
 }

Here is the call graph for this function:

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

protectedvirtual

Reimplemented from MatrixHandler.

Definition at line 327 of file spmapmh.cc.

References begin(), col_indices, SparseMatrixHandler::iGetNumCols(), SparseMatrixHandler::iGetNumRows(), VectorHandler::iGetSize(), MBDYN_EXCEPT_ARGS, and SparseMatrixHandler::NCols.

 {
         if (in.iGetSize() != iGetNumCols()
                         || out.iGetSize() != iGetNumRows())
         {
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         row_cont_type::const_iterator ri, re;
 
         for (integer col = 0; col < NCols; col++) {
                 re = col_indices[col].end();
                 for (ri = col_indices[col].begin(); ri != re; ++ri) {
                         doublereal d = ri->second*in(col + 1);
                         (out.*op)(ri->first + 1, d);
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

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

Definition at line 269 of file spmapmh.cc.

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

Referenced by ThirdOrderIntegrator::Jacobian().

 {
         if ((out.iGetNumCols() < iGetNumCols() + dcol)
                         || (out.iGetNumRows() < iGetNumRows() + drow))
         {
                 silent_cerr("Assertion fault "
                         "in SpMapMatrixHandler::MulAndSumWithShift"
                         << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         drow = drow + 1;
 
         for (integer col = 0; col < NCols; col++) {
                 row_cont_type::const_iterator ri, re;
                 re = col_indices[col].end();
                 integer newcol = col + dcol + 1;
                 for (ri = col_indices[col].begin(); ri != re; ++ri) {
                         out.IncCoef(ri->first + drow, newcol, ri->second*s);
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

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

inlinevirtual

Implements MatrixHandler.

Definition at line 129 of file spmapmh.h.

References ASSERTMSGBREAK, col_indices, SparseMatrixHandler::NCols, SparseMatrixHandler::NRows, and SparseMatrixHandler::NZ.

Referenced by DecCoef(), IncCoef(), and PutCoef().

                                                               {
                 ASSERTMSGBREAK(i_row > 0 && i_row <= NRows,
                                 "Error in SpMapMatrixHandler::operator(), "
                                 "row index out of range");
                 ASSERTMSGBREAK(i_col > 0 && i_col <= NCols, 
                                 "Error in SpMapMatrixHandler::operator(), "
                                 "col index out of range");
                 i_row--;
                 i_col--;
                 row_cont_type::iterator i;
                 row_cont_type& row = col_indices[i_col];
                 i = row.find(i_row);
                 if (i == row.end()) {
                         NZ++;
 
                         row[i_row] = 0.;
 
                         return row[i_row];
 
                 } else {
                         return i->second;
                 }
         };

const doublereal& SpMapMatrixHandler::operator()	(	integer	ix,
		integer	iy
	)		const

inlinevirtual

Implements MatrixHandler.

Definition at line 229 of file spmapmh.h.

References ASSERTMSGBREAK, col_indices, SparseMatrixHandler::NCols, SparseMatrixHandler::NRows, and Zero1.

                                                                      {
                 ASSERTMSGBREAK(ix > 0 && ix <= NRows,
                                 "Error in SpMapMatrixHandler::operator(), "
                                 "row index out of range");
                 ASSERTMSGBREAK(iy > 0 && iy <= NCols,
                                 "Error in SpMapMatrixHandler::operator(), "
                                 "col index out of range");
                 row_cont_type::iterator i;
                 row_cont_type& row = col_indices[iy - 1];
                 i = row.find(ix - 1);
                 if (i == row.end()) {
                         return ::Zero1;
 
                 } else {
                         return i->second;
                 }
         };

void SpMapMatrixHandler::PutCoef	(	integer	ix,
		integer	iy,
		const doublereal &	val
	)

inlinevirtual

Reimplemented from MatrixHandler.

Definition at line 187 of file spmapmh.h.

References ASSERTMSGBREAK, col_indices, SparseMatrixHandler::NCols, SparseMatrixHandler::NRows, and operator()().

Referenced by SchurMatrixHandler::PutCoef(), SchurMatrixHandlerUm::PutCoef(), and ReadAerodynamicModal().

                                                                     {
                 ASSERTMSGBREAK(ix - 1 < NRows,
                                 "Error in SpMapMatrixHandler::PutCoef(), "
                                 "row index out of range");
                 ASSERTMSGBREAK(iy - 1 < NCols,
                                 "Error in SpMapMatrixHandler::PutCoef(), "
                                 "col index out of range");
 #ifdef MBDYN_X_KEEP_SPARSITY
                 /* try to keep sparsity */
                 if (val != 0.) {
 #endif /* MBDYN_X_KEEP_SPARSITY */
                         operator()(ix,iy) = val;
 #ifdef MBDYN_X_KEEP_SPARSITY
                 } else {
                         row_cont_type::iterator i;
                         row_cont_type& row = col_indices[iy-1];
                         i = row.find(ix - 1);
                         if (i != row.end()) {
                                 i->second = val;
                         }
                 }
 #endif /* MBDYN_X_KEEP_SPARSITY */
         };

Here is the call graph for this function:

void SpMapMatrixHandler::Reset ( void )

virtual

Implements MatrixHandler.

Definition at line 161 of file spmapmh.cc.

References begin(), col_indices, and SparseMatrixHandler::NCols.

Referenced by main(), SchurMatrixHandler::MatEFCReset(), SchurMatrixHandlerUm::MatEFCReset(), and ReadAerodynamicModal().

 {
         row_cont_type::const_iterator re;
         row_cont_type::iterator ri;
         for (integer col = 0; col < NCols; col++) {
                 re = col_indices[col].end();
                 for (ri = col_indices[col].begin(); ri != re; ++ri) {
                         ri->second = 0.;
                 }
         }
 }

Here is the call graph for this function:

void SpMapMatrixHandler::Resize	(	integer	ir,
		integer	ic
	)

virtual

Implements MatrixHandler.

Definition at line 174 of file spmapmh.cc.

References col_indices, m_end, SparseMatrixHandler::NCols, SparseMatrixHandler::NRows, SparseMatrixHandler::NZ, and SpMapMatrixHandler::const_iterator::reset().

Referenced by ThirdOrderIntegrator::Advance().

 {
         if (nn == 0) {
                 nn = n;
         }
 
         for (integer col = 0; col < NCols; col++) {
                 col_indices[col].clear();
         }
 
         col_indices.resize(nn);
         NRows = n;
         NCols = nn;
         NZ = 0;
 
         m_end.reset();
 }

Here is the call graph for this function:

Member Data Documentation

std::vector<row_cont_type> SpMapMatrixHandler::col_indices

mutableprivate

Definition at line 79 of file spmapmh.h.

Referenced by dGetCoef(), FakeThirdOrderMulAndSumWithShift(), GetCol(), MakeCompressedColumnForm(), MakeIndexForm(), MatMatMul_base(), MatTMatMul_base(), MatTVecMul_base(), MatVecMul_base(), MulAndSumWithShift(), operator()(), PutCoef(), SpMapMatrixHandler::const_iterator::reset(), Reset(), Resize(), and SpMapMatrixHandler().

const_iterator SpMapMatrixHandler::m_end

private

Definition at line 113 of file spmapmh.h.

Referenced by end(), and Resize().

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

Classes

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation