MatrixScale< T > Class Template Reference

#include <dgeequ.h>

Inheritance diagram for MatrixScale< T >:

Collaboration diagram for MatrixScale< T >:

Public Member Functions
	MatrixScale (const SolutionManager::ScaleOpt &scale)
virtual	~MatrixScale ()
T &	ScaleMatrix (T &mh) const
bool	ComputeScaleFactors (const T &mh)
Public Member Functions inherited from MatrixScaleBase
	MatrixScaleBase (const SolutionManager::ScaleOpt &scale)
virtual	~MatrixScaleBase ()
VectorHandler &	ScaleRightHandSide (VectorHandler &bVH) const
VectorHandler &	ScaleSolution (VectorHandler &xVH) const
std::ostream &	Report (std::ostream &os) const
const std::vector< doublereal > &	GetRowScale () const
const std::vector< doublereal > &	GetColScale () const
bool	bGetInitialized () const

Static Public Member Functions
static MatrixScale< T > *	Allocate (const SolutionManager::ScaleOpt &scale)

Protected Member Functions
virtual bool	ComputeScaleFactors (const T &mh, std::vector< doublereal > &rowScale, std::vector< doublereal > &colScale)=0
Protected Member Functions inherited from MatrixScaleBase
MatrixHandler::Norm_t	GetCondNumNorm () const
virtual std::ostream &	vReport (std::ostream &os) const =0
void	Prepare (const MatrixHandler &mh, integer &nrows, integer &ncols)
void	PrepareRows (const MatrixHandler &mh, integer &nrows)
void	PrepareCols (const MatrixHandler &mh, integer &ncols)
bool	bReport () const

Additional Inherited Members
Protected Attributes inherited from MatrixScaleBase
std::vector< doublereal >	rowScale
std::vector< doublereal >	colScale
doublereal	dCondBefore
doublereal	dCondAfter
const unsigned	uFlags
bool	bOK

Detailed Description

template<typename T>
class MatrixScale< T >

Definition at line 79 of file dgeequ.h.

Constructor & Destructor Documentation

template<typename T >

MatrixScale< T >::MatrixScale ( const SolutionManager::ScaleOpt &  scale )

inlineexplicit

Definition at line 261 of file dgeequ.h.

        :MatrixScaleBase(scale)
{

}

template<typename T >

MatrixScale< T >::~MatrixScale ( )

virtual

Definition at line 268 of file dgeequ.h.

{

}

Member Function Documentation

template<typename T >

MatrixScale< T > * MatrixScale< T >::Allocate ( const SolutionManager::ScaleOpt &  scale )

static

Definition at line 313 of file dgeequ.h.

References SolutionManager::ScaleOpt::algorithm, ASSERT, MBDYN_EXCEPT_ARGS, SAFENEWWITHCONSTRUCTOR, SolutionManager::SCALEA_COL_MAX, SolutionManager::SCALEA_COL_SUM, SolutionManager::SCALEA_ITERATIVE, SolutionManager::SCALEA_LAPACK, SolutionManager::SCALEA_ROW_MAX, SolutionManager::SCALEA_ROW_MAX_COL_MAX, and SolutionManager::SCALEA_ROW_SUM.

Referenced by NaiveSparseSolutionManager::GetMatrixScale().

{
        MatrixScale<T>* pMatScale = 0;

        switch (scale.algorithm) {
        case SolutionManager::SCALEA_ROW_MAX:
                SAFENEWWITHCONSTRUCTOR(pMatScale,
                                                           RowMaxMatrixScale<T>,
                                                           RowMaxMatrixScale<T>(scale));
                break;

        case SolutionManager::SCALEA_ROW_SUM:
                SAFENEWWITHCONSTRUCTOR(pMatScale,
                                                           RowSumMatrixScale<T>,
                                                           RowSumMatrixScale<T>(scale));
                break;

        case SolutionManager::SCALEA_COL_MAX:
                        SAFENEWWITHCONSTRUCTOR(pMatScale,
                                                                   ColMaxMatrixScale<T>,
                                                                   ColMaxMatrixScale<T>(scale));
                        break;

        case SolutionManager::SCALEA_COL_SUM:
                SAFENEWWITHCONSTRUCTOR(pMatScale,
                                                           ColSumMatrixScale<T>,
                                                           ColSumMatrixScale<T>(scale));
                break;

        case SolutionManager::SCALEA_LAPACK:
                SAFENEWWITHCONSTRUCTOR(pMatScale,
                                                           LapackMatrixScale<T>,
                                                           LapackMatrixScale<T>(scale));
                break;

        case SolutionManager::SCALEA_ITERATIVE:
                SAFENEWWITHCONSTRUCTOR(pMatScale,
                                                           IterativeMatrixScale<T>,
                                                           IterativeMatrixScale<T>(scale));
                break;

        case SolutionManager::SCALEA_ROW_MAX_COL_MAX:
                SAFENEWWITHCONSTRUCTOR(pMatScale,
                                                           RowMaxColMaxMatrixScale<T>,
                                                           RowMaxColMaxMatrixScale<T>(scale));
                break;
        default:
                ASSERT(0);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        return pMatScale;
}

template<typename T >

bool MatrixScale< T >::ComputeScaleFactors ( const T &  mh )

inline

Definition at line 305 of file dgeequ.h.

Referenced by ScaleMatrix(), and NaiveSparseSolutionManager::ScaleMatrixAndRightHandSide().

{
        bOK = ComputeScaleFactors(mh, rowScale, colScale);

        return bOK;
}

template<typename T>

virtual bool MatrixScale< T >::ComputeScaleFactors ( const T &  mh, std::vector< doublereal > &  rowScale, std::vector< doublereal > &  colScale  )

protectedpure virtual

Implemented in RowMaxColMaxMatrixScale< T >, IterativeMatrixScale< T >, LapackMatrixScale< T >, ColMaxMatrixScale< T >, ColSumMatrixScale< T >, RowMaxMatrixScale< T >, and RowSumMatrixScale< T >.

template<typename T >

T & MatrixScale< T >::ScaleMatrix ( T &  mh ) const

inline

Definition at line 274 of file dgeequ.h.

References SolutionManager::SCALEF_COND_NUM.

Referenced by ScaleMatrix(), and NaiveSparseSolutionManager::ScaleMatrixAndRightHandSide().

{
        if (uFlags & SolutionManager::SCALEF_COND_NUM) {
                dCondBefore = mh.ConditionNumber(GetCondNumNorm());
        }

        const bool bScaleRows = !rowScale.empty();
        const bool bScaleCols = !colScale.empty();

        for (typename T::const_iterator i = mh.begin(); i != mh.end(); ++i) {
                doublereal dCoef = i->dCoef;

                if (bScaleRows) {
                        dCoef *= rowScale[i->iRow];
                }

                if (bScaleCols) {
                        dCoef *= colScale[i->iCol];
                }

                mh(i->iRow + 1, i->iCol + 1) = dCoef;
        }

        if (uFlags & SolutionManager::SCALEF_COND_NUM) {
                dCondAfter = mh.ConditionNumber(GetCondNumNorm());
        }

        return mh;
}

---

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

• dgeequ.h