NaiveSparseSolutionManager Class Reference

#include <naivewrap.h>

Inheritance diagram for NaiveSparseSolutionManager:

Collaboration diagram for NaiveSparseSolutionManager:

Public Member Functions
	NaiveSparseSolutionManager (const integer Dim, const doublereal dMP=1.e-9, const ScaleOpt &scale=ScaleOpt())
virtual	~NaiveSparseSolutionManager (void)
virtual void	MatrReset (void)
virtual void	Solve (void)
virtual MatrixHandler *	pMatHdl (void) const
virtual MyVectorHandler *	pResHdl (void) const
virtual MyVectorHandler *	pSolHdl (void) const
Public Member Functions inherited from SolutionManager
	SolutionManager (void)
virtual	~SolutionManager (void)
virtual void	LinkToSolution (VectorHandler &XCurr, VectorHandler &XPrimeCurr)
virtual void	MatrInitialize (void)
virtual void	SolveT (void)
doublereal *	pdSetResVec (doublereal *pd)
doublereal *	pdSetSolVec (doublereal *pd)
bool	bGetConditionNumber (doublereal &dCond) const

Protected Member Functions
template<class MH >
void	ScaleMatrixAndRightHandSide (MH &mh)
template<typename MH >
MatrixScale< MH > &	GetMatrixScale ()
void	ScaleSolution (void)

Protected Attributes
NaiveMatrixHandler *	A
MyVectorHandler	VH
ScaleOpt	scale
MatrixScaleBase *	pMatScale
Protected Attributes inherited from SolutionManager
LinearSolver *	pLS

Additional Inherited Members
Public Types inherited from SolutionManager
enum	ScaleWhen { SCALEW_NEVER = 0, SCALEW_ONCE, SCALEW_ALWAYS }
enum	ScaleAlgorithm {   SCALEA_NONE, SCALEA_UNDEF, SCALEA_ROW_MAX, SCALEA_ROW_SUM,   SCALEA_COL_MAX, SCALEA_COL_SUM, SCALEA_LAPACK, SCALEA_ITERATIVE,   SCALEA_ROW_MAX_COL_MAX }
enum	ScaleFlags {   SCALEF_DEFAULT = 0x0u, SCALEF_WARN = 0x1u, SCALEF_VERBOSE = 0x2u, SCALEF_COND_NUM_1 = 0x4u,   SCALEF_COND_NUM_INF = 0x8u, SCALEF_COND_NUM = SCALEF_COND_NUM_1 | SCALEF_COND_NUM_INF }

Detailed Description

Definition at line 75 of file naivewrap.h.

Constructor & Destructor Documentation

NaiveSparseSolutionManager::NaiveSparseSolutionManager	(	const integer	Dim,
		const doublereal	dMP = 1.e-9,
		const ScaleOpt &	scale = ScaleOpt()
	)

Definition at line 151 of file naivewrap.cc.

References A, SolutionManager::ScaleOpt::algorithm, NaiveMatrixHandler::NaiveMatrixHandler(), NaiveMatrixHandler::NaiveSolver, MyVectorHandler::pdGetVec(), LinearSolver::pdSetResVec(), LinearSolver::pdSetSolVec(), SolutionManager::pLS, SAFENEWWITHCONSTRUCTOR, scale, SolutionManager::SCALEA_NONE, SolutionManager::SCALEA_UNDEF, SolutionManager::SCALEW_NEVER, LinearSolver::SetSolutionManager(), VH, and SolutionManager::ScaleOpt::when.

: A(0),
VH(Dim),
scale(s),
pMatScale(0)
{
        switch (scale.algorithm) {
        case SCALEA_NONE:
        case SCALEA_UNDEF:
                scale.when = SCALEW_NEVER; // Default scaling is not supported for this solver
                break;

        default:
                // Allocate MatrixScale<T> on demand
                ;
        }

        SAFENEWWITHCONSTRUCTOR(A, NaiveMatrixHandler, NaiveMatrixHandler(Dim));
        SAFENEWWITHCONSTRUCTOR(pLS, NaiveSolver, NaiveSolver(Dim, dMP, A));

        pLS->pdSetResVec(VH.pdGetVec());
        pLS->pdSetSolVec(VH.pdGetVec());

        pLS->SetSolutionManager(this);
}

Here is the call graph for this function:

NaiveSparseSolutionManager::~NaiveSparseSolutionManager ( void  )

virtual

Definition at line 178 of file naivewrap.cc.

References A, pMatScale, and SAFEDELETE.

{
        if (A != 0) {
                SAFEDELETE(A);
                A = 0;
        }

        if (pMatScale != 0) {
                SAFEDELETE(pMatScale);
        }
}

Member Function Documentation

template<typename MH >

MatrixScale< MH > & NaiveSparseSolutionManager::GetMatrixScale ( )

protected

Definition at line 223 of file naivewrap.cc.

References MatrixScale< T >::Allocate(), pMatScale, and scale.

{
        if (pMatScale == 0) {
                pMatScale = MatrixScale<MH>::Allocate(scale);
        }

        // Will throw std::bad_cast if the type does not match
        return dynamic_cast<MatrixScale<MH>&>(*pMatScale);
}

Here is the call graph for this function:

void NaiveSparseSolutionManager::MatrReset ( void  )

virtual

Implements SolutionManager.

Reimplemented in NaiveSparsePermSolutionManager< T >.

Definition at line 191 of file naivewrap.cc.

References SolutionManager::pLS, and LinearSolver::Reset().

Referenced by NaiveSparsePermSolutionManager< T >::MatrReset().

{
        pLS->Reset();
}

Here is the call graph for this function:

MatrixHandler * NaiveSparseSolutionManager::pMatHdl ( void  ) const

virtual

Implements SolutionManager.

Definition at line 256 of file naivewrap.cc.

References A.

{
        return A;
}

MyVectorHandler * NaiveSparseSolutionManager::pResHdl ( void  ) const

virtual

Implements SolutionManager.

Definition at line 263 of file naivewrap.cc.

References VH.

{
        return &VH;
}

MyVectorHandler * NaiveSparseSolutionManager::pSolHdl ( void  ) const

virtual

Implements SolutionManager.

Definition at line 270 of file naivewrap.cc.

References VH.

{
        return &VH;
}

template<class MH >

void NaiveSparseSolutionManager::ScaleMatrixAndRightHandSide ( MH &  mh )

protected

Definition at line 199 of file naivewrap.cc.

References MatrixScaleBase::bGetInitialized(), LinearSolver::bReset(), MatrixScale< T >::ComputeScaleFactors(), SolutionManager::pLS, MatrixScaleBase::Report(), scale, MatrixScale< T >::ScaleMatrix(), MatrixScaleBase::ScaleRightHandSide(), SolutionManager::SCALEW_ALWAYS, SolutionManager::SCALEW_NEVER, VH, and SolutionManager::ScaleOpt::when.

Referenced by Solve().

{
        if (scale.when != SCALEW_NEVER) {
                MatrixScale<MH>& rMatScale = GetMatrixScale<MH>();

                if (pLS->bReset()) {
                        if (!rMatScale.bGetInitialized()
                                || scale.when == SolutionManager::SCALEW_ALWAYS) {
                                // (re)compute
                                rMatScale.ComputeScaleFactors(mh);
                        }
                        // in any case scale matrix and right-hand-side
                        rMatScale.ScaleMatrix(mh);

                        if (silent_err) {
                                rMatScale.Report(std::cerr);
                        }
                }

                rMatScale.ScaleRightHandSide(VH);
        }
}

Here is the call graph for this function:

void NaiveSparseSolutionManager::ScaleSolution ( void  )

protected

Definition at line 234 of file naivewrap.cc.

References ASSERT, pMatScale, scale, MatrixScaleBase::ScaleSolution(), SolutionManager::SCALEW_NEVER, VH, and SolutionManager::ScaleOpt::when.

Referenced by Solve().

{
        if (scale.when != SCALEW_NEVER) {
                ASSERT(pMatScale != 0);
                // scale solution
                pMatScale->ScaleSolution(VH);
        }
}

Here is the call graph for this function:

void NaiveSparseSolutionManager::Solve ( void  )

virtual

Implements SolutionManager.

Reimplemented in NaiveSparsePermSolutionManager< T >.

Definition at line 245 of file naivewrap.cc.

References A, SolutionManager::pLS, ScaleMatrixAndRightHandSide(), ScaleSolution(), and LinearSolver::Solve().

{
        ScaleMatrixAndRightHandSide(*A);

        pLS->Solve();

        ScaleSolution();
}

Here is the call graph for this function:

Member Data Documentation

NaiveMatrixHandler* NaiveSparseSolutionManager::A

mutableprotected

Definition at line 77 of file naivewrap.h.

Referenced by NaiveSparsePermSolutionManager< T >::NaiveSparsePermSolutionManager(), NaiveSparseSolutionManager(), pMatHdl(), Solve(), and ~NaiveSparseSolutionManager().

MatrixScaleBase* NaiveSparseSolutionManager::pMatScale

protected

Definition at line 81 of file naivewrap.h.

Referenced by GetMatrixScale(), ScaleSolution(), and ~NaiveSparseSolutionManager().

ScaleOpt NaiveSparseSolutionManager::scale

protected

Definition at line 80 of file naivewrap.h.

Referenced by GetMatrixScale(), NaiveSparseSolutionManager(), ScaleMatrixAndRightHandSide(), and ScaleSolution().

MyVectorHandler NaiveSparseSolutionManager::VH

mutableprotected

Definition at line 78 of file naivewrap.h.

Referenced by NaiveSparseSolutionManager(), pResHdl(), pSolHdl(), ScaleMatrixAndRightHandSide(), and ScaleSolution().

---

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

• naivewrap.h
• naivewrap.cc