#include <naivewrap.h>

Inheritance diagram for NaiveSparsePermSolutionManager< T >:

Collaboration diagram for NaiveSparsePermSolutionManager< T >:

Public Member Functions
	NaiveSparsePermSolutionManager (const integer Dim, const doublereal dMP=1.e-9, const ScaleOpt &scale=ScaleOpt())

virtual	~NaiveSparsePermSolutionManager (void)

virtual void	Solve (void)

virtual void	MatrInitialize (void)

Public Member Functions inherited from NaiveSparseSolutionManager
	NaiveSparseSolutionManager (const integer Dim, const doublereal dMP=1.e-9, const ScaleOpt &scale=ScaleOpt())

virtual	~NaiveSparseSolutionManager (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	SolveT (void)

doublereal *	pdSetResVec (doublereal *pd)

doublereal *	pdSetSolVec (doublereal *pd)

bool	bGetConditionNumber (doublereal &dCond) const

Protected Types
enum	{ PERM_NO, PERM_INTERMEDIATE, PERM_READY }

Protected Member Functions
virtual void	MatrReset (void)

Protected Member Functions inherited from NaiveSparseSolutionManager
template<class MH >
void	ScaleMatrixAndRightHandSide (MH &mh)

template<typename MH >
MatrixScale< MH > &	GetMatrixScale ()

void	ScaleSolution (void)

Protected Attributes
enum NaiveSparsePermSolutionManager:: { ... }	ePermState

std::vector< integer >	perm

std::vector< integer >	invperm

Protected Attributes inherited from NaiveSparseSolutionManager
NaiveMatrixHandler *	A

MyVectorHandler	VH

ScaleOpt	scale

MatrixScaleBase *	pMatScale

Protected Attributes inherited from SolutionManager
LinearSolver *	pLS

Private Member Functions
void	ComputePermutation (void)

void	BackPerm (void)

template<>
void	ComputePermutation (void)

Private Attributes
const doublereal	dMinPiv

MyVectorHandler	TmpH

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

template<class T>
class NaiveSparsePermSolutionManager< T >

Definition at line 124 of file naivewrap.h.

Member Enumeration Documentation

template<class T>

anonymous enum

protected

Enumerator
PERM_NO
PERM_INTERMEDIATE
PERM_READY

Definition at line 133 of file naivewrap.h.

              {
                 PERM_NO,
                 PERM_INTERMEDIATE,
                 PERM_READY
         } ePermState;

Constructor & Destructor Documentation

template<class T >

NaiveSparsePermSolutionManager< T >::NaiveSparsePermSolutionManager	(	const integer	Dim,
		const doublereal	dMP = `1.e-9`,
		const ScaleOpt &	scale = `ScaleOpt()`
	)

Definition at line 281 of file naivewrap.cc.

References NaiveSparseSolutionManager::A, NaiveSparsePermSolutionManager< T >::invperm, NaiveSparsePermSolutionManager< T >::MatrInitialize(), NaivePermMatrixHandler::NaivePermMatrixHandler(), NaiveSparsePermSolutionManager< T >::perm, SolutionManager::pLS, SAFEDELETE, and SAFENEWWITHCONSTRUCTOR.

 : NaiveSparseSolutionManager(Dim, dMP, scale),
 dMinPiv(dMP < 0 ? 0 : dMP),
 TmpH(Dim),
 ePermState(PERM_NO)
 {
         perm.resize(Dim, 0);
         invperm.resize(Dim, 0);
 
         // replace matrix handler
         SAFEDELETE(A);
         A = 0;
         SAFENEWWITHCONSTRUCTOR(A, NaivePermMatrixHandler,
                         NaivePermMatrixHandler(Dim, perm, invperm));
 
         dynamic_cast<NaiveSolver *>(pLS)->SetMat(A);
 
         MatrInitialize();
 }

Here is the call graph for this function:

template<class T >

NaiveSparsePermSolutionManager< T >::~NaiveSparsePermSolutionManager ( void )

virtual

Definition at line 305 of file naivewrap.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

template<class T >

void NaiveSparsePermSolutionManager< T >::BackPerm ( void )

private

Definition at line 323 of file naivewrap.cc.

References ASSERT, and NaivePermMatrixHandler::invperm.

 {
         /* NOTE: use whatever is stored in pLS - someone could
          * trick us into using its memory */
         doublereal *pd = pLS->pdGetResVec();
 
         ASSERT(pd != TmpH.pdGetVec());
         
         for (integer i = 0; i < A->iGetNumCols(); i++) {
                 pd[invperm[i]] = TmpH(i + 1);
         }
 }

template<class T>

void NaiveSparsePermSolutionManager< T >::ComputePermutation ( void )

private

template<>

void NaiveSparsePermSolutionManager< Colamd_ordering >::ComputePermutation ( void )

private

Definition at line 390 of file naivewrap.cc.

References COLAMD_KNOBS, COLAMD_STATS, NaivePermMatrixHandler::invperm, mbdyn_colamd(), mbdyn_colamd_recommended(), mbdyn_colamd_set_defaults(), MBDYN_EXCEPT_ARGS, and NaivePermMatrixHandler::perm.

 {
         std::vector<integer> Ai;
         A->MakeCCStructure(Ai, invperm);
         doublereal knobs[COLAMD_KNOBS];
         integer stats[COLAMD_STATS];
         integer Alen = mbdyn_colamd_recommended(Ai.size(), A->iGetNumRows(),
                         A->iGetNumCols());
         Ai.resize(Alen);
         mbdyn_colamd_set_defaults(knobs);
         if (!mbdyn_colamd(A->iGetNumRows(), A->iGetNumCols(), Alen,
                 &Ai[0], &invperm[0], knobs, stats))
         {
                 silent_cerr("colamd permutation failed" << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
         for (integer i = 0; i < A->iGetNumRows(); i++) {
                 perm[invperm[i]] = i;
         }
         ePermState = PERM_INTERMEDIATE;
 }

Here is the call graph for this function:

template<class T >

void NaiveSparsePermSolutionManager< T >::MatrInitialize ( void )

virtual

Reimplemented from SolutionManager.

Definition at line 371 of file naivewrap.cc.

References NaivePermMatrixHandler::invperm, and NaivePermMatrixHandler::perm.

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

 {
         ePermState = PERM_NO;
         for (integer i = 0; i < A->iGetNumRows(); i++) {
                 perm[i] = i;
                 invperm[i] = i;
         }
 
         MatrReset();
 }

template<class T >

void NaiveSparsePermSolutionManager< T >::MatrReset ( void )

protectedvirtual

Reimplemented from NaiveSparseSolutionManager.

Definition at line 312 of file naivewrap.cc.

References NaiveSparseSolutionManager::MatrReset().

 {
         if (ePermState == PERM_INTERMEDIATE) {
                 ePermState = PERM_READY;
         }
 
         NaiveSparseSolutionManager::MatrReset();
 }

Here is the call graph for this function:

template<class T >

void NaiveSparsePermSolutionManager< T >::Solve ( void )

virtual

Reimplemented from NaiveSparseSolutionManager.

Definition at line 340 of file naivewrap.cc.

References ASSERT.

 {
         doublereal *pd = 0;
 
         // Will throw std::bad_cast if the data type does not match
         ScaleMatrixAndRightHandSide(dynamic_cast<NaivePermMatrixHandler&>(*A));
 
         if (ePermState == PERM_NO) {
                 ComputePermutation();
 
         } else if (ePermState == PERM_READY) {
                 /* We need to use local storage to allow BackPerm();
                  * save and restore original pointer */
                 pd = pLS->pdSetSolVec(TmpH.pdGetVec());
         }
 
         pLS->Solve();
 
         if (ePermState == PERM_READY) {
                 BackPerm();
 
                 ASSERT(pd != 0);
                 pLS->pdSetSolVec(pd);
         }
 
         ScaleSolution();
 }