#include <schurmh.h>

Inheritance diagram for SchurMatrixHandlerUm:

Collaboration diagram for SchurMatrixHandlerUm:

Classes
class	ErrGeneric

Public Member Functions
	SchurMatrixHandlerUm (int LocSize, int IntSize, MatrixHandler pBM, integer pGlobToLoc)

	~SchurMatrixHandlerUm (void)

void	MatEFCReset (void)

void	Reset (void)

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

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

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

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

doublereal *	GetECol (const integer iCol) const

doublereal *	GetEColSol (const integer iCol) const

void	CompLocSchur (void)

VectorHandler &	CompNewf (VectorHandler &f, const VectorHandler &g) const

void	PrintMatrix (void)

Public Member Functions inherited from SchurMatrixHandler
	SchurMatrixHandler (int LocSize, int IntSize, MatrixHandler pBM, integer pGlobToLoc, doublereal *pdEv=0)

virtual	~SchurMatrixHandler (void)

virtual integer	iGetNumRows (void) const

virtual integer	iGetNumCols (void) const

MatrixHandler *	GetBMat (void)

void	SetBMat (MatrixHandler *pBM)

doublereal *	GetCMat (void)

virtual void	Resize (integer, integer)

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

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

virtual VectorHandler &	CompNewg (VectorHandler &g, const VectorHandler &f) 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

Private Attributes
doublereal *	pdEs

MyVectorHandler *	pEs

bool	Eflag

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

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

Protected Attributes inherited from SchurMatrixHandler
integer	LSize

integer	ISize

MatrixHandler *	pB

MyVectorHandler *	pE

doublereal *	pdE

SpMapMatrixHandler *	pF

MyVectorHandler *	pC

doublereal *	pdC

integer *	pGTL

bool	bExtpdE

Detailed Description

Definition at line 771 of file schurmh.h.

Constructor & Destructor Documentation

SchurMatrixHandlerUm::SchurMatrixHandlerUm	(	int	LocSize,
		int	IntSize,
		MatrixHandler *	pBM,
		integer *	pGlobToLoc
	)

Definition at line 211 of file schurmh.cc.

References MyVectorHandler::Attach(), SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, SchurMatrixHandler::pdE, pdEs, SchurMatrixHandler::pE, pEs, SAFENEWARR, and SAFENEWWITHCONSTRUCTOR.

 : SchurMatrixHandler(LocSize, IntSize, pBM, pGlobToLoc, 0),
 pdEs(0),
 pEs(0),
 Eflag(1)
 {
         SAFENEWARR(pdEs, doublereal, LSize*(ISize + 1));
         pdE = pdEs + LSize;
         pE->Attach(LSize*ISize, pdE, LSize*ISize);
         SAFENEWWITHCONSTRUCTOR(pEs,
                         MyVectorHandler,
                         MyVectorHandler(LSize*ISize, pdEs));
 }

Here is the call graph for this function:

SchurMatrixHandlerUm::~SchurMatrixHandlerUm ( void )

Definition at line 227 of file schurmh.cc.

References SchurMatrixHandler::pdE, pdEs, SchurMatrixHandler::pE, SAFEDELETE, and SAFEDELETEARR.

 {
         if (pE != 0) {
                 SAFEDELETE(pE);
         }
 
         if (pdEs != 0) {
                 SAFEDELETEARR(pdEs);
         }
 
         pdE = 0;
 }

Member Function Documentation

void SchurMatrixHandlerUm::CompLocSchur ( void )

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 1052 of file schurmh.h.

References Eflag, SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, SchurMatrixHandler::pdC, pdEs, and SchurMatrixHandler::pF.

 {
         Eflag = false;
         for (int j = 0; j < ISize; j++) {
                 int iColc = j * ISize;
                 int iCole = j * LSize;
 
                 for (int k = 0; k < LSize; k++) {
                         for (int i = 0; i < ISize; i++) {
                                 pdC[i + iColc] -=
                                         pF->operator()(i + 1, k + 1) * pdEs[k + iCole];
                         }
                 }
         }
 }

VectorHandler & SchurMatrixHandlerUm::CompNewf	(	VectorHandler &	f,
		const VectorHandler &	g
	)		const

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 1069 of file schurmh.h.

References ASSERT, VectorHandler::DecCoef(), VectorHandler::iGetSize(), SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, and pdEs.

 {
 #ifdef DEBUG
         ASSERT(f.iGetSize() == LSize);
         ASSERT(g.iGetSize() == ISize);
 #endif /* DEBUG */
 
         for (int j = 0; j < ISize; j++) {
                 int iColx = j * LSize;
 
                 for (int i = 0; i < LSize; i++) {
                         if (pdEs[i + iColx] != 0) {
                                 f.DecCoef(i + 1, pdEs[i + iColx]*g(j + 1));
                         }
                 }
         }
         return f;
 }

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 948 of file schurmh.h.

References MatrixHandler::DecCoef(), SpMapMatrixHandler::DecCoef(), MyVectorHandler::DecCoef(), Eflag, SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, MBDYN_EXCEPT_ARGS, SchurMatrixHandler::pB, SchurMatrixHandler::pC, SchurMatrixHandler::pE, pEs, SchurMatrixHandler::pF, and SchurMatrixHandler::pGTL.

 {
 #ifdef DEBUG
         IsValid();
 #endif /* DEBUG */
 
         if (pGTL[iRow] > 0) {
                 if (pGTL[iCol] > 0) {
                         pB->DecCoef(pGTL[iRow], pGTL[iCol], dCoef);
                         return;
 
                 } else if (pGTL[iCol] < 0) {
                         if (Eflag) {
                                 pE->DecCoef(pGTL[iRow] - (pGTL[iCol] + 1)*LSize, dCoef);
 
                         } else {
                                 pEs->DecCoef(pGTL[iRow] - (pGTL[iCol] + 1)*LSize, dCoef);
                         }
                         return;
                 }
 
         } else if (pGTL[iRow] < 0) {
                 if (iCol > 0) {
                         pF->DecCoef(-pGTL[iRow], pGTL[iCol], dCoef);
                         return;
 
                 } else if (iCol < 0) {
                         pC->DecCoef(-pGTL[iRow] - (pGTL[iCol] + 1)*ISize, dCoef);
                         return;
                 }
         }
 
 #ifdef USE_MPI
         silent_cerr("SchurMatrixHandlerUm::DecCoef() "
                 "Process(" << MBDynComm.Get_rank() << "): "
                 "trying to operate on nonlocal indices "
                 << iRow << "," << iCol << std::endl);
 #else /* ! USE_MPI */
         silent_cerr("SchurMatrixHandlerUm::DecCoef() "
                 "trying to operate on nonlocal indices "
                 << iRow << "," << iCol << std::endl);
 #endif /* ! USE_MPI */
 
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 996 of file schurmh.h.

References Eflag, SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, SchurMatrixHandler::pB, SchurMatrixHandler::pC, SchurMatrixHandler::pE, pEs, SchurMatrixHandler::pF, SchurMatrixHandler::pGTL, and Zero1.

 {
 #ifdef DEBUG
         IsValid();
 #endif /* DEBUG */
 
         if (pGTL[iRow] > 0) {
                 if (pGTL[iCol] > 0) {
                         return pB->operator()(pGTL[iRow], pGTL[iCol]);
 
                 } else if (pGTL[iCol] < 0) {
                         if (Eflag) {
                                 return pE->operator()(pGTL[iRow] - (pGTL[iCol] + 1)*LSize);
 
                         } else {
                                 return pEs->operator()(pGTL[iRow] - (pGTL[iCol] + 1)*LSize);
                         }
                 }
 
         } else if (pGTL[iRow] < 0) {
                 if (pGTL[iCol] > 0) {
                         return pF->operator()(-pGTL[iRow], pGTL[iCol]);
 
                 } else if (pGTL[iCol] < 0) {
                         return pC->operator()(-pGTL[iRow] - (pGTL[iCol] + 1)*ISize);
                 }
         }
 
 #ifdef USE_MPI
         silent_cerr("SchurMatrixHandlerUm::dGetCoef() "
                 "Process(" << MBDynComm.Get_rank() << "): "
                 "trying to operate on nonlocal indices "
                 << iRow << "," << iCol << std::endl);
 #else /* ! USE_MPI */
         silent_cerr("SchurMatrixHandlerUm::dGetCoef() "
                 "trying to operate on nonlocal indices "
                 << iRow << "," << iCol << std::endl);
 #endif /* ! USE_MPI */
 
         return ::Zero1;
 }

doublereal * SchurMatrixHandlerUm::GetECol ( const integer iCol ) const

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 1039 of file schurmh.h.

References SchurMatrixHandler::LSize, and SchurMatrixHandler::pdE.

 {
         return &pdE[iCol*LSize];
 }

doublereal * SchurMatrixHandlerUm::GetEColSol ( const integer iCol ) const

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 1045 of file schurmh.h.

References SchurMatrixHandler::LSize, and pdEs.

 {
         return &pdEs[iCol*LSize];
 }

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

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 900 of file schurmh.h.

References Eflag, MatrixHandler::IncCoef(), SpMapMatrixHandler::IncCoef(), MyVectorHandler::IncCoef(), SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, MBDYN_EXCEPT_ARGS, SchurMatrixHandler::pB, SchurMatrixHandler::pC, SchurMatrixHandler::pE, pEs, SchurMatrixHandler::pF, and SchurMatrixHandler::pGTL.

 {
 #ifdef DEBUG
         IsValid();
 #endif /* DEBUG */
 
         if (pGTL[iRow] > 0) {
                 if (pGTL[iCol] > 0) {
                         pB->IncCoef(pGTL[iRow], pGTL[iCol], dCoef);
                         return;
 
                 } else if (pGTL[iCol] < 0) {
                         if (Eflag) {
                                 pE->IncCoef(pGTL[iRow] - (pGTL[iCol] + 1)*LSize, dCoef);
 
                         } else {
                                 pEs->IncCoef(pGTL[iRow] - (pGTL[iCol] + 1)*LSize, dCoef);
                         }
                         return;
                 }
 
         } else if (pGTL[iRow] < 0) {
                 if (pGTL[iCol] > 0) {
                         pF->IncCoef(-pGTL[iRow], pGTL[iCol], dCoef);
                         return;
 
                 } else if (pGTL[iCol] < 0) {
                         pC->IncCoef(-pGTL[iRow] - (pGTL[iCol] + 1)*ISize, dCoef);
                         return;
                 }
         }
 
 #ifdef USE_MPI
         silent_cerr("SchurMatrixHandlerUm::IncCoef() "
                 "Process(" << MBDynComm.Get_rank() << "): "
                 "trying to operate on nonlocal indices "
                 << iRow << "," << iCol << std::endl);
 #else /* ! USE_MPI */
         silent_cerr("SchurMatrixHandlerUm::IncCoef() "
                 "trying to operate on nonlocal indices "
                 << iRow << "," << iCol << std::endl);
 #endif /* ! USE_MPI */
 
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

Here is the call graph for this function:

void SchurMatrixHandlerUm::MatEFCReset ( void )

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 826 of file schurmh.h.

References Eflag, SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, SchurMatrixHandler::pC, pdEs, SchurMatrixHandler::pF, MyVectorHandler::Reset(), and SpMapMatrixHandler::Reset().

Referenced by Reset().

 {
 #ifdef DEBUG
         IsValid();
 #endif /* DEBUG */
 
         Eflag = true;
         for (int i = 0; i < LSize*(ISize + 1); i++) {
                 pdEs[i] = 0.;
         }
         pF->Reset();
         pC->Reset();
 }

Here is the call graph for this function:

void SchurMatrixHandlerUm::PrintMatrix ( void )

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 1089 of file schurmh.h.

References Eflag, SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, SchurMatrixHandler::pB, SchurMatrixHandler::pdC, SchurMatrixHandler::pdE, pdEs, and SchurMatrixHandler::pF.

 {
         silent_cout("Schur Matrix " << std::endl);
 
         for (int i = 0; i < LSize; i++) {
                 for (int j = 0; j < LSize; j++) {
                         silent_cout(pB->operator()(i + 1, j + 1) << " ");
                 }
 
                 for (int j = 0; j < ISize; j++) {
                         if (Eflag) {
                                 silent_cout(pdE[i + j*LSize] << " ");
                         } else {
                                 silent_cout(pdEs[i + j*LSize] << " ");
                         }
                 }
                 silent_cout(std::endl);
         }
 
         for (int i = 0; i < ISize; i++) {
                 for (int j = 0;j < LSize; j++) {
                         silent_cout(pF->operator()(i + 1, j + 1) << " ");
                 }
 
                 for (int j = 0; j < ISize; j++) {
                         silent_cout(pdC[i + j*ISize] << " ");
                 }
                 silent_cout(std::endl);
         }
 }

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

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 852 of file schurmh.h.

References Eflag, SchurMatrixHandler::ISize, SchurMatrixHandler::LSize, MBDYN_EXCEPT_ARGS, SchurMatrixHandler::pB, SchurMatrixHandler::pC, SchurMatrixHandler::pE, pEs, SchurMatrixHandler::pF, SchurMatrixHandler::pGTL, MatrixHandler::PutCoef(), MyVectorHandler::PutCoef(), and SpMapMatrixHandler::PutCoef().

 {
 #ifdef DEBUG
         IsValid();
 #endif /* DEBUG */
 
         if (pGTL[iRow] > 0) {
                 if (pGTL[iCol] > 0) {
                         pB->PutCoef(pGTL[iRow], pGTL[iCol], dCoef);
                         return;
 
                 } else if (pGTL[iCol] < 0 ) {
                         if (Eflag) {
                                 pE->PutCoef(pGTL[iRow] - (pGTL[iCol] + 1)*LSize, dCoef);
 
                         } else {
                                 pEs->PutCoef(pGTL[iRow] - (pGTL[iCol] + 1)*LSize, dCoef);
                         }
                         return;
                 }
 
         } else if (pGTL[iRow] < 0 ) {
                 if (pGTL[iCol] > 0) {
                         pF->PutCoef(-pGTL[iRow], pGTL[iCol], dCoef);
                         return;
 
                 } else if (pGTL[iCol] < 0) {
                         pC->PutCoef(-pGTL[iRow] - (pGTL[iCol] + 1)*ISize, dCoef);
                         return;
                 }
         }
 
 #ifdef USE_MPI
         silent_cerr("SchurMatrixHandlerUm::PutCoef() "
                 "Process(" << MBDynComm.Get_rank() << "): "
                 "trying to operate on nonlocal indices "
                 << iRow << "," << iCol << std::endl);
 #else /* ! USE_MPI */
         silent_cerr("SchurMatrixHandlerUm::PutCoef() "
                 "trying to operate on nonlocal indices "
                 << iRow << "," << iCol << std::endl);
 #endif /* ! USE_MPI */
 
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

Here is the call graph for this function:

void SchurMatrixHandlerUm::Reset ( void )

inlinevirtual

Reimplemented from SchurMatrixHandler.

Definition at line 840 of file schurmh.h.

References Eflag, MatEFCReset(), SchurMatrixHandler::pB, and MatrixHandler::Reset().

 {
 #ifdef DEBUG
         IsValid();
 #endif /* DEBUG */
 
         Eflag = true;
         pB->Reset();
         MatEFCReset();
 }

Here is the call graph for this function:

Member Data Documentation

bool SchurMatrixHandlerUm::Eflag

private

Definition at line 782 of file schurmh.h.

Referenced by CompLocSchur(), DecCoef(), dGetCoef(), IncCoef(), MatEFCReset(), PrintMatrix(), PutCoef(), and Reset().

doublereal* SchurMatrixHandlerUm::pdEs

private

Definition at line 780 of file schurmh.h.

Referenced by CompLocSchur(), CompNewf(), GetEColSol(), MatEFCReset(), PrintMatrix(), SchurMatrixHandlerUm(), and ~SchurMatrixHandlerUm().

MyVectorHandler* SchurMatrixHandlerUm::pEs

private

Definition at line 781 of file schurmh.h.

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

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

Classes

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation