#include <mh.h>

Inheritance diagram for MatrixHandler:

Classes
class	ErrGeneric

class	ErrMatrixIsSingular

class	ErrRebuildMatrix

Public Types
enum	Norm_t { NORM_1, NORM_INF }

Public Member Functions
virtual	~MatrixHandler (void)

virtual void	Resize (integer, integer)=0

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 void	Reset (void)=0

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

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

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

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

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

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

virtual integer	iGetNumRows (void) const =0

virtual integer	iGetNumCols (void) const =0

virtual MatrixHandler &	operator= (const MatrixHandler &MH)

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

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

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

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

virtual MatrixHandler &	ScalarMul (const doublereal &d)

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

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

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

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

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

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

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

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

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

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

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

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

virtual doublereal	ConditionNumber (enum Norm_t eNorm=NORM_1) const

virtual doublereal	Norm (enum Norm_t eNorm=NORM_1) const

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

virtual 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

Detailed Description

Definition at line 60 of file mh.h.

Member Enumeration Documentation

enum MatrixHandler::Norm_t

Enumerator
NORM_1
NORM_INF

Definition at line 62 of file mh.h.

62 { NORM_1, NORM_INF };

MatrixHandler::NORM_INF

Definition: mh.h:62

MatrixHandler::NORM_1

Definition: mh.h:62

Constructor & Destructor Documentation

MatrixHandler::~MatrixHandler ( void )

virtual

Definition at line 49 of file mh.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

doublereal MatrixHandler::ConditionNumber ( enum Norm_t eNorm = NORM_1 ) const

virtual

Definition at line 451 of file mh.cc.

References ASSERT, MBDYN_EXCEPT_ARGS, NORM_1, and NORM_INF.

Referenced by LineSearchSolver::Jacobian(), and NewtonRaphsonSolver::Solve().

 {
 #if HAVE_CONDITION_NUMBER
         integer M;
         integer N;
         std::vector<doublereal> A;
 
         LapackCopyMatrix(*this, A, M, N);
 
         const doublereal ANORM = LapackMatrixNorm(A, M, N, eNorm);
 
 #ifdef DEBUG
         std::cerr << "ANORM=" << ANORM << std::endl;
         std::cerr << "A=" << std::endl;
         for (int i = 0; i < M; ++i) {
                 for (int j = 0; j < N; ++j) {
                         std::cerr << A[j * M + i] << '\t';
                 }
                 std::cerr << std::endl;
         }
 #endif // DEBUG
 
         integer INFO = 0;
         std::vector<integer> IPIV(std::min(M,N));
 
         __FC_DECL__(dgetrf)(&M, &N, &A[0], &M, &IPIV[0], &INFO);
 
         ASSERT(INFO == 0); // should not fail because the Jacobian has already been factorised before
 
         std::vector<doublereal> WORK(4*N);
         std::vector<integer> IWORK(N);
         doublereal RCOND = 0.;
         char norm;
 
         switch (eNorm) {
         case NORM_1:
                 norm = '1';
                 break;
         case NORM_INF:
                 norm = 'I';
                 break;
         default:
                 ASSERT(0);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         __FC_DECL__(dgecon)(&norm, &N, &A[0], &M, &ANORM, &RCOND, &WORK[0], &IWORK[0], &INFO);
 
         ASSERT(INFO == 0); // should not fail
 
 #ifdef DEBUG
         std::cerr << "RCOND=" << RCOND << std::endl;
 #endif // DEBUG
 
         return 1. / RCOND;
 #else // ! HAVE_CONDITION_NUMBER
         silent_cerr("Condition number is not available in this version of MBDyn" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 #endif // ! HAVE_CONDITION_NUMBER
 }

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

virtual

Reimplemented in SparseSubMatrixHandler, SchurMatrixHandlerUm, FullSubMatrixHandler, FullMatrixHandler, SpMapMatrixHandler, and SchurMatrixHandler.

Definition at line 379 of file mh.cc.

References operator()().

Referenced by SchurMatrixHandler::DecCoef(), SchurMatrixHandlerUm::DecCoef(), MatMatDecMul(), MatMatMul_base(), FullMatrixHandler::MatMatMul_base(), MatTMatDecMul(), MatTMatMul_base(), and FullMatrixHandler::MatTMatMul_base().

                                                                     {
         operator()(ix, iy) -= inc;
 }

Here is the call graph for this function:

const doublereal & MatrixHandler::dGetCoef	(	integer	iRow,
		integer	iCol
	)		const

virtual

Reimplemented in SparseSubMatrixHandler, SchurMatrixHandlerUm, FullSubMatrixHandler, FullMatrixHandler, SpMapMatrixHandler, and SchurMatrixHandler.

Definition at line 389 of file mh.cc.

References operator()().

                                                     {
         return operator()(ix, iy);
 }

Here is the call graph for this function:

virtual integer MatrixHandler::iGetNumCols ( void ) const

pure virtual

Implemented in SparseSubMatrixHandler, FullSubMatrixHandler, FullMatrixHandler, SparseMatrixHandler, NaiveMatrixHandler, and SchurMatrixHandler.

Referenced by FullSubMatrixHandler::AddTo(), SparseSubMatrixHandler::AddTo(), FullSubMatrixHandler::AddToT(), SparseSubMatrixHandler::AddToT(), DataManager::AssJac(), DirCColMatrixHandler< off >::FakeThirdOrderMulAndSumWithShift(), CColMatrixHandler< off >::FakeThirdOrderMulAndSumWithShift(), SpMapMatrixHandler::FakeThirdOrderMulAndSumWithShift(), StepNIntegrator::Jacobian(), NaiveMatrixHandler::MatMatMul_base(), MatMatMul_base(), CompactSparseMatrixHandler_tpl< off >::MatMatMul_base(), SpMapMatrixHandler::MatMatMul_base(), FullMatrixHandler::MatMatMul_base(), NaivePermMatrixHandler::MatMatMul_base(), NaiveMatrixHandler::MatTMatMul_base(), MatTMatMul_base(), CompactSparseMatrixHandler_tpl< off >::MatTMatMul_base(), SpMapMatrixHandler::MatTMatMul_base(), FullMatrixHandler::MatTMatMul_base(), NaivePermMatrixHandler::MatTMatMul_base(), MatTVecMul_base(), MatVecMul_base(), DirCColMatrixHandler< off >::MulAndSumWithShift(), CColMatrixHandler< off >::MulAndSumWithShift(), SpMapMatrixHandler::MulAndSumWithShift(), NaiveMatrixHandler::operator+=(), NaiveMatrixHandler::operator-=(), operator<<(), operator=(), DataManager::OutputEigenvectors(), MatrixScaleBase::PrepareCols(), ReportMatScale(), ScalarMul(), SetupSystem(), LineSearchSolver::Solve(), FullSubMatrixHandler::SubFrom(), SparseSubMatrixHandler::SubFrom(), FullSubMatrixHandler::SubFromT(), and SparseSubMatrixHandler::SubFromT().

virtual integer MatrixHandler::iGetNumRows ( void ) const

pure virtual

Implemented in SparseSubMatrixHandler, FullSubMatrixHandler, FullMatrixHandler, SparseMatrixHandler, NaiveMatrixHandler, and SchurMatrixHandler.

Referenced by FullSubMatrixHandler::AddTo(), SparseSubMatrixHandler::AddTo(), FullSubMatrixHandler::AddToT(), SparseSubMatrixHandler::AddToT(), DataManager::AssJac(), DirCColMatrixHandler< off >::FakeThirdOrderMulAndSumWithShift(), CColMatrixHandler< off >::FakeThirdOrderMulAndSumWithShift(), SpMapMatrixHandler::FakeThirdOrderMulAndSumWithShift(), StepNIntegrator::Jacobian(), NaiveMatrixHandler::MatMatMul_base(), MatMatMul_base(), CompactSparseMatrixHandler_tpl< off >::MatMatMul_base(), SpMapMatrixHandler::MatMatMul_base(), FullMatrixHandler::MatMatMul_base(), NaivePermMatrixHandler::MatMatMul_base(), NaiveMatrixHandler::MatTMatMul_base(), MatTMatMul_base(), CompactSparseMatrixHandler_tpl< off >::MatTMatMul_base(), SpMapMatrixHandler::MatTMatMul_base(), FullMatrixHandler::MatTMatMul_base(), NaivePermMatrixHandler::MatTMatMul_base(), MatTVecMul_base(), MatVecMul_base(), DirCColMatrixHandler< off >::MulAndSumWithShift(), CColMatrixHandler< off >::MulAndSumWithShift(), SpMapMatrixHandler::MulAndSumWithShift(), NaiveMatrixHandler::operator+=(), NaiveMatrixHandler::operator-=(), operator<<(), operator=(), DataManager::OutputEigenvectors(), MatrixScaleBase::PrepareRows(), ReportMatScale(), ScalarMul(), SetupSystem(), LineSearchSolver::Solve(), FullSubMatrixHandler::SubFrom(), SparseSubMatrixHandler::SubFrom(), FullSubMatrixHandler::SubFromT(), and SparseSubMatrixHandler::SubFromT().

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

virtual

Reimplemented in SparseSubMatrixHandler, SchurMatrixHandlerUm, FullSubMatrixHandler, FullMatrixHandler, SpMapMatrixHandler, and SchurMatrixHandler.

Definition at line 374 of file mh.cc.

References operator()().

Referenced by DirCColMatrixHandler< off >::FakeThirdOrderMulAndSumWithShift(), CColMatrixHandler< off >::FakeThirdOrderMulAndSumWithShift(), SpMapMatrixHandler::FakeThirdOrderMulAndSumWithShift(), SchurMatrixHandler::IncCoef(), SchurMatrixHandlerUm::IncCoef(), MatMatIncMul(), MatMatMul(), MatMatMul_base(), FullMatrixHandler::MatMatMul_base(), MatTMatIncMul(), MatTMatMul(), MatTMatMul_base(), FullMatrixHandler::MatTMatMul_base(), DirCColMatrixHandler< off >::MulAndSumWithShift(), CColMatrixHandler< off >::MulAndSumWithShift(), and SpMapMatrixHandler::MulAndSumWithShift().

                                                                     {
         operator()(ix, iy) += inc;
 }

Here is the call graph for this function:

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

virtual

Definition at line 262 of file mh.cc.

References DecCoef(), and MatMatMul_base().

 {
         return MatMatMul_base(&MatrixHandler::DecCoef, out, in);
 }

Here is the call graph for this function:

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

virtual

Definition at line 250 of file mh.cc.

References IncCoef(), and MatMatMul_base().

 {
         return MatMatMul_base(&MatrixHandler::IncCoef, out, in);
 }

Here is the call graph for this function:

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

virtual

Definition at line 232 of file mh.cc.

References IncCoef(), MatMatMul_base(), and Reset().

Referenced by Membrane4EAS::AssJac(), Shell4EAS::AssJac(), Shell4EASANS::AssJac(), Shell4EASANS::AssRes(), main(), Shell4EAS::Shell4EAS(), and Shell4EASANS::Shell4EASANS().

 {
         /* Put is implemented resetting out first, then passing IncCoef()
          * so that out-of-order assignments work */
         out.Reset();
         return MatMatMul_base(&MatrixHandler::IncCoef, out, in);
 }

Here is the call graph for this function:

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

protectedvirtual

Reimplemented in NaivePermMatrixHandler, FullMatrixHandler, SpMapMatrixHandler, CompactSparseMatrixHandler_tpl< off >, and NaiveMatrixHandler.

Definition at line 135 of file mh.cc.

References c, DecCoef(), iGetNumCols(), iGetNumRows(), IncCoef(), MBDYN_EXCEPT_ARGS, operator()(), and PutCoef().

Referenced by MatMatDecMul(), MatMatIncMul(), MatMatMul(), and FullMatrixHandler::MatMatMul_base().

 {
         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("MatrixHandler::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);
         }
 
         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++) {
                                 d += this->operator()(r, k)*in(k, c);
                         }
 
                         (out.*op)(r, c, d);
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

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

virtual

Definition at line 268 of file mh.cc.

References DecCoef(), and MatTMatMul_base().

 {
         return MatTMatMul_base(&MatrixHandler::DecCoef, out, in);
 }

Here is the call graph for this function:

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

virtual

Definition at line 256 of file mh.cc.

References IncCoef(), and MatTMatMul_base().

 {
         return MatTMatMul_base(&MatrixHandler::IncCoef, out, in);
 }

Here is the call graph for this function:

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

virtual

Definition at line 241 of file mh.cc.

References IncCoef(), MatTMatMul_base(), and Reset().

Referenced by main(), Shell4EAS::Shell4EAS(), and Shell4EASANS::Shell4EASANS().

 {
         /* Put is implemented resetting out first, then passing IncCoef()
          * so that out-of-order assignments work */
         out.Reset();
         return MatTMatMul_base(&MatrixHandler::IncCoef, out, in);
 }

Here is the call graph for this function:

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

protectedvirtual

Reimplemented in NaivePermMatrixHandler, FullMatrixHandler, SpMapMatrixHandler, CompactSparseMatrixHandler_tpl< off >, and NaiveMatrixHandler.

Definition at line 183 of file mh.cc.

References c, DecCoef(), iGetNumCols(), iGetNumRows(), IncCoef(), MBDYN_EXCEPT_ARGS, operator()(), and PutCoef().

Referenced by MatTMatDecMul(), MatTMatIncMul(), MatTMatMul(), and FullMatrixHandler::MatTMatMul_base().

 {
         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("MatrixHandler::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);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         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++) {
                                 d += this->operator()(k, r)*in(k, c);
                         }
 
                         (out.*op)(r, c, d);
                 }
         }
 
         return out;
 }

Here is the call graph for this function:

VectorHandler & MatrixHandler::MatTVecDecMul	(	VectorHandler &	out,
		const VectorHandler &	in
	)		const

virtual

Definition at line 368 of file mh.cc.

References VectorHandler::DecCoef(), and MatTVecMul_base().

Referenced by LineSearchSolver::Solve().

 {
         return MatTVecMul_base(&VectorHandler::DecCoef, out, in);
 }

Here is the call graph for this function:

VectorHandler & MatrixHandler::MatTVecIncMul	(	VectorHandler &	out,
		const VectorHandler &	in
	)		const

virtual

Definition at line 356 of file mh.cc.

References VectorHandler::IncCoef(), and MatTVecMul_base().

 {
         return MatTVecMul_base(&VectorHandler::IncCoef, out, in);
 }

Here is the call graph for this function:

VectorHandler & MatrixHandler::MatTVecMul	(	VectorHandler &	out,
		const VectorHandler &	in
	)		const

virtual

Definition at line 341 of file mh.cc.

References VectorHandler::IncCoef(), MatTVecMul_base(), and VectorHandler::Reset().

Referenced by main(), ModalMappingExt::Recv(), StructMappingExtForce::RecvFromFileDes(), and StructMembraneMappingExtForce::RecvFromFileDes().

 {
         /* Put is implemented resetting out first, then passing IncCoef()
          * so that out-of-order assignments work */
         out.Reset();
         return MatTVecMul_base(&VectorHandler::IncCoef, out, in);
 }

Here is the call graph for this function:

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

protectedvirtual

Reimplemented in NaivePermMatrixHandler, SpMapMatrixHandler, FullMatrixHandler, CompactSparseMatrixHandler_tpl< off >, and NaiveMatrixHandler.

Definition at line 304 of file mh.cc.

References c, iGetNumCols(), iGetNumRows(), VectorHandler::iGetSize(), MBDYN_EXCEPT_ARGS, and operator()().

Referenced by MatTVecDecMul(), MatTVecIncMul(), and MatTVecMul().

 {
         integer nr = iGetNumRows();
         integer nc = iGetNumCols();
 
         if (out.iGetSize() != nc || in.iGetSize() != nr) {
                 silent_cerr("MatrixHandler::MatVecMul_base(): size mismatch "
                         "out(" << out.iGetSize() << ", 1) "
                         "= this(" << nr << ", " << nc << ")^T "
                         "* in(" << in.iGetSize() << ", 1)" << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         for (integer r = 1; r <= nc; r++) {
                 doublereal d = 0.;
 
                 for (integer c = 1; c <= nr; c++) {
                         d += this->operator()(c, r)*in(c);
                 }
                 (out.*op)(r, d);
         }
 
         return out;
 }

Here is the call graph for this function:

VectorHandler & MatrixHandler::MatVecDecMul	(	VectorHandler &	out,
		const VectorHandler &	in
	)		const

virtual

Definition at line 362 of file mh.cc.

References VectorHandler::DecCoef(), and MatVecMul_base().

Referenced by SchurMatrixHandler::CompNewg().

 {
         return MatVecMul_base(&VectorHandler::DecCoef, out, in);
 }

Here is the call graph for this function:

VectorHandler & MatrixHandler::MatVecIncMul	(	VectorHandler &	out,
		const VectorHandler &	in
	)		const

virtual

Definition at line 350 of file mh.cc.

References VectorHandler::IncCoef(), and MatVecMul_base().

Referenced by AerodynamicModal::AssVec(), and main().

 {
         return MatVecMul_base(&VectorHandler::IncCoef, out, in);
 }

Here is the call graph for this function:

VectorHandler & MatrixHandler::MatVecMul	(	VectorHandler &	out,
		const VectorHandler &	in
	)		const

virtual

Definition at line 332 of file mh.cc.

References VectorHandler::IncCoef(), MatVecMul_base(), and VectorHandler::Reset().

Referenced by main(), ModuleNonsmoothNode::mbs_get_force_frictional(), ModalMappingExt::Send(), StructMappingExtForce::SendToFileDes(), and StructMembraneMappingExtForce::SendToFileDes().

 {
         /* Put is implemented resetting out first, then passing IncCoef()
          * so that out-of-order assignments work */
         out.Reset();
         return MatVecMul_base(&VectorHandler::IncCoef, out, in);
 }

Here is the call graph for this function:

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

protectedvirtual

Reimplemented in NaivePermMatrixHandler, SpMapMatrixHandler, FullMatrixHandler, CompactSparseMatrixHandler_tpl< off >, and NaiveMatrixHandler.

Definition at line 275 of file mh.cc.

References c, iGetNumCols(), iGetNumRows(), VectorHandler::iGetSize(), MBDYN_EXCEPT_ARGS, and operator()().

Referenced by MatVecDecMul(), MatVecIncMul(), and MatVecMul().

 {
         integer nr = iGetNumRows();
         integer nc = iGetNumCols();
 
         if (out.iGetSize() != nr || in.iGetSize() != nc) {
                 silent_cerr("MatrixHandler::MatVecMul_base(): size mismatch "
                         "out(" << out.iGetSize() << ", 1) "
                         "= this(" << nr << ", " << nc << ") "
                         "* in(" << in.iGetSize() << ", 1)" << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         for (integer r = 1; r <= nr; r++) {
                 doublereal d = 0.;
 
                 for (integer c = 1; c <= nc; c++) {
                         d += this->operator()(r, c)*in(c);
                 }
                 (out.*op)(r, d);
         }
 
         return out;
 
 }

Here is the call graph for this function:

doublereal MatrixHandler::Norm ( enum Norm_t eNorm = NORM_1 ) const

virtual

Definition at line 512 of file mh.cc.

 {
         integer M;
         integer N;
         std::vector<doublereal> A;
 
         LapackCopyMatrix(*this, A, M, N);
 
         return LapackMatrixNorm(A, M, N, eNorm);
 }

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

pure virtual

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

pure virtual

Implemented in VariableSubMatrixHandler, SparseSubMatrixHandler, FullSubMatrixHandler, NaivePermMatrixHandler, FullMatrixHandler, NaiveMatrixHandler, SpMapMatrixHandler, SchurMatrixHandler, DirCColMatrixHandler< off >, and CColMatrixHandler< off >.

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

virtual

Reimplemented in FullMatrixHandler, and NaiveMatrixHandler.

Definition at line 92 of file mh.cc.

References SubMatrixHandler::AddTo().

 {
         return SubMH.AddTo(*this);
 }

Here is the call graph for this function:

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

virtual

Reimplemented in FullMatrixHandler, and NaiveMatrixHandler.

Definition at line 106 of file mh.cc.

References VariableSubMatrixHandler::AddTo().

 {
         return SubMH.AddTo(*this);
 }

Here is the call graph for this function:

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

virtual

Reimplemented in FullMatrixHandler, and NaiveMatrixHandler.

Definition at line 99 of file mh.cc.

References SubMatrixHandler::SubFrom().

 {
         return SubMH.SubFrom(*this);
 }

Here is the call graph for this function:

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

virtual

Reimplemented in FullMatrixHandler, and NaiveMatrixHandler.

Definition at line 113 of file mh.cc.

References VariableSubMatrixHandler::SubFrom().

 {
         return SubMH.SubFrom(*this);
 }

Here is the call graph for this function:

MatrixHandler & MatrixHandler::operator= ( const MatrixHandler & MH )

virtual

Definition at line 72 of file mh.cc.

References iGetNumCols(), iGetNumRows(), operator()(), and Resize().

 {
         integer nr = MH.iGetNumRows();
         integer nc = MH.iGetNumCols();
 
         Resize(nr, nc);
 
         for (integer i = 1; i <= nr; i++) { 
                 for (integer ii = 1; ii <= nc; ii++) { 
                         this->operator()(i, ii) = MH(i, ii);
                 }
         }
 
         return *this;
 }

Here is the call graph for this function:

integer MatrixHandler::PacMat ( void )

virtual

Definition at line 65 of file mh.cc.

 {
         return 0L;
 }

const doublereal * MatrixHandler::pdGetMat ( void ) const

inlinevirtual

Reimplemented in CompactSparseMatrixHandler, and FullMatrixHandler.

Definition at line 209 of file mh.h.

 {
         return NULL;
 }

doublereal * MatrixHandler::pdGetMat ( void )

inlinevirtual

Reimplemented in FullMatrixHandler.

Definition at line 215 of file mh.h.

 {
         return NULL;
 }

integer * MatrixHandler::piGetCols ( void ) const

inlinevirtual

Definition at line 229 of file mh.h.

 {
         return NULL;
 }

integer * MatrixHandler::piGetRows ( void ) const

inlinevirtual

Definition at line 222 of file mh.h.

 {
         return NULL;
 }

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

virtual

Reimplemented in SparseSubMatrixHandler, SchurMatrixHandlerUm, FullSubMatrixHandler, SpMapMatrixHandler, FullMatrixHandler, and SchurMatrixHandler.

Definition at line 384 of file mh.cc.

References operator()().

Referenced by NRTrim::DoTrim(), DataManager::InitialJointAssembly(), MatMatMul_base(), FullMatrixHandler::MatMatMul_base(), MatTMatMul_base(), FullMatrixHandler::MatTMatMul_base(), SchurMatrixHandler::PutCoef(), and SchurMatrixHandlerUm::PutCoef().

                                                                     {
         operator()(ix, iy) = val;
 }

Here is the call graph for this function:

virtual void MatrixHandler::Reset ( void )

pure virtual

Implemented in SparseSubMatrixHandler, SchurMatrixHandlerUm, SpMapMatrixHandler, FullSubMatrixHandler, CompactSparseMatrixHandler, FullMatrixHandler, NaiveMatrixHandler, and SchurMatrixHandler.

Referenced by DataManager::AssConstrJac(), DataManager::AssJac(), Solver::Eig(), ThirdOrderIntegrator::Jacobian(), main(), MatMatMul(), MatTMatMul(), SchurMatrixHandler::Reset(), SchurMatrixHandlerUm::Reset(), and ResizeReset().

virtual void MatrixHandler::Resize	(	integer	,
		integer
	)

pure virtual

Implemented in SparseSubMatrixHandler, SpMapMatrixHandler, FullSubMatrixHandler, CColMatrixHandler< off >, FullMatrixHandler, NaiveMatrixHandler, DirCColMatrixHandler< off >, SubMatrixHandler, and SchurMatrixHandler.

Referenced by operator=(), and ResizeReset().

void MatrixHandler::ResizeReset	(	integer	iNewRow,
		integer	iNewCol
	)

virtual

Reimplemented in SparseSubMatrixHandler, FullSubMatrixHandler, and SubMatrixHandler.

Definition at line 56 of file mh.cc.

References Reset(), and Resize().

 {
         Resize(iNewRow, iNewCol);
         Reset();
 }

Here is the call graph for this function:

MatrixHandler & MatrixHandler::ScalarMul ( const doublereal & d )

virtual

Definition at line 119 of file mh.cc.

References iGetNumCols(), iGetNumRows(), and operator()().

 {
         integer nr = iGetNumRows();
         integer nc = iGetNumCols();
 
         for (integer i = 1; i <= nr; i++) {
                 for (integer j = 1; j <= nc; j++) {
                         this->operator()(i, j) *= d;
                 }
         }
 
         return *this;
 }

Here is the call graph for this function:

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

mh.h
mh.cc

Classes

Public Types

Public Member Functions

Protected Member Functions

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation