#include <stepsol.h>

Inheritance diagram for StepNIntegrator:

Collaboration diagram for StepNIntegrator:

Public Member Functions
	StepNIntegrator (const integer MaxIt, const doublereal dT, const doublereal dSolutionTol, const integer stp, const bool bmod_res_test)

virtual	~StepNIntegrator (void)

virtual void	Residual (VectorHandler *pRes) const

virtual void	Jacobian (MatrixHandler *pJac) const

virtual void	Update (const VectorHandler *pSol) const

virtual doublereal	TestScale (const NonlinearSolverTest *pTest, doublereal &dAlgebraicEqu) const

Public Member Functions inherited from ImplicitStepIntegrator
	ImplicitStepIntegrator (const integer MaxIt, const doublereal dT, const doublereal dSolutionTol, const integer stp, const integer sts, const bool bmod_res_test)

virtual	~ImplicitStepIntegrator (void)

virtual void	EvalProd (doublereal Tau, const VectorHandler &f0, const VectorHandler &w, VectorHandler &z) const

Public Member Functions inherited from StepIntegrator
	StepIntegrator (const integer MaxIt, const doublereal dT, const doublereal dSolutionTol, const integer stp, const integer sts)

virtual	~StepIntegrator (void)

void	SetDataManager (DataManager *pDatMan)

virtual integer	GetIntegratorNumPreviousStates (void) const

virtual integer	GetIntegratorNumUnknownStates (void) const

virtual integer	GetIntegratorMaxIters (void) const

virtual doublereal	GetIntegratorDTol (void) const

virtual doublereal	GetIntegratorDSolTol (void) const

virtual void	OutputTypes (const bool fpred)

virtual void	SetDriveHandler (const DriveHandler *pDH)

virtual doublereal	Advance (Solver pS, const doublereal TStep, const doublereal dAlph, const StepChange StType, std::deque< MyVectorHandler > &qX, std::deque< MyVectorHandler * > &qXPrime, MyVectorHandler const pX, MyVectorHandler const pXPrime, integer &EffIter, doublereal &Err, doublereal &SolErr)=0

Public Member Functions inherited from NonlinearProblem
virtual	~NonlinearProblem (void)

Public Attributes
doublereal	db0Differential

doublereal	db0Algebraic

Protected Member Functions
void	UpdateDof (const int DCount, const DofOrder::Order Order, const VectorHandler *const pSol=0) const

virtual void	SetCoef (doublereal dT, doublereal dAlpha, enum StepChange NewStep)=0

Protected Member Functions inherited from StepIntegrator
template<class T >
void	UpdateLoop (const T const t, void(T::pUpd)(const int DCount, const DofOrder::Order Order, const VectorHandler const pSol) const, const VectorHandler const pSol=0) const

Additional Inherited Members
Public Types inherited from StepIntegrator
enum	{ DIFFERENTIAL = 0, ALGEBRAIC = 1 }

enum	StepChange { NEWSTEP, REPEATSTEP }

Protected Attributes inherited from ImplicitStepIntegrator
VectorHandler *	pXCurr

VectorHandler *	pXPrimeCurr

bool	bModResTest

Protected Attributes inherited from StepIntegrator
DataManager *	pDM

const DataManager::DofVecType *	pDofs

bool	outputPred

integer	MaxIters

doublereal	dTol

doublereal	dSolTol

integer	steps

integer	unkstates

Detailed Description

Definition at line 231 of file stepsol.h.

Constructor & Destructor Documentation

StepNIntegrator::StepNIntegrator	(	const integer	MaxIt,
		const doublereal	dT,
		const doublereal	dSolutionTol,
		const integer	stp,
		const bool	bmod_res_test
	)

Definition at line 446 of file stepsol.cc.

References NO_OP.

 : ImplicitStepIntegrator(MaxIt, dT, dSolutionTol, stp, 1, bmod_res_test),
 db0Differential(0.),
 db0Algebraic(0.)
 {
         NO_OP;
 }

StepNIntegrator::~StepNIntegrator ( void )

virtual

Definition at line 458 of file stepsol.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

void StepNIntegrator::Jacobian ( MatrixHandler * pJac ) const

virtual

Implements NonlinearProblem.

Definition at line 472 of file stepsol.cc.

References ASSERT, DataManager::AssJac(), DataManager::AssRes(), db0Differential, MatrixHandler::iGetNumCols(), MatrixHandler::iGetNumRows(), StepIntegrator::pDM, NaiveMatrixHandler::Reset(), Residual(), and Update().

 {
         ASSERT(pDM != NULL);
         pDM->AssJac(*pJac, db0Differential);
 
 #ifdef MBDYN_FDJAC
         // Finite difference check of Jacobian matrix
         // Uncomment this whenever you need to debug your new Jacobian
         // NOTE: might not be safe!
         if (pDM->bFDJac()) {
                 NaiveMatrixHandler fdjac(pJac->iGetNumRows());
                 fdjac.Reset();
                 MyVectorHandler basesol(pJac->iGetNumRows());
                 MyVectorHandler incsol(pJac->iGetNumRows());
                 MyVectorHandler inc(pJac->iGetNumRows());
                 Residual(&basesol);
                 doublereal ddd = 0.001;
                 for (integer i = 1; i <= pJac->iGetNumRows(); i++) {
                         incsol.Reset();
                         inc.Reset();
                         inc.PutCoef(i, ddd);
                         Update(&inc);
                         // std::cerr << pXPrimeCurr->operator()(30) << std::endl;
                         pDM->AssRes(incsol, db0Differential);
                         inc.Reset();
                         inc.PutCoef(i, -ddd);
                         Update(&inc);
                         incsol -= basesol;
                         incsol *= (1./(-ddd));
                         for (integer j = 1; j <= pJac->iGetNumCols(); j++) {
                                 fdjac.PutCoef(j, i, std::abs(incsol(j)) > 1.E-100 ? incsol(j) : 0.);
                         }
                 }
 
                 std::cerr << "\nxxxxxxxxxxxxxxx\n" << std::endl;
                 std::cerr << *pJac << std::endl;
                 std::cerr << "\n---------------\n" << std::endl;
                 std::cerr << fdjac << std::endl;
                 std::cerr << "\n===============\n" << std::endl;
         }
 #endif // MBDYN_FDJAC
 }

Here is the call graph for this function:

void StepNIntegrator::Residual ( VectorHandler * pRes ) const

virtual

Implements NonlinearProblem.

Definition at line 464 of file stepsol.cc.

References ASSERT, DataManager::AssRes(), db0Differential, and StepIntegrator::pDM.

Referenced by Jacobian().

 {
         ASSERT(pDM != NULL);
         pDM->AssRes(*pRes, db0Differential);
 }

Here is the call graph for this function:

virtual void StepNIntegrator::SetCoef	(	doublereal	dT,
		doublereal	dAlpha,
		enum StepChange	NewStep
	)

protectedpure virtual

Implemented in HopeSolver, MultistepSolver, Step2Integrator, ImplicitEulerIntegrator, and CrankNicolsonIntegrator.

Referenced by Step1Integrator::Advance().

doublereal StepNIntegrator::TestScale	(	const NonlinearSolverTest *	pTest,
		doublereal &	dAlgebraicEqu
	)		const

virtual

Reimplemented from ImplicitStepIntegrator.

Definition at line 547 of file stepsol.cc.

References db0Differential, and ImplicitStepIntegrator::TestScale().

 {
         const doublereal dDiffEqu = ImplicitStepIntegrator::TestScale(pTest, dAlgebraicEqu);
         dAlgebraicEqu = db0Differential;
 
         return dDiffEqu;
 }

Here is the call graph for this function:

void StepNIntegrator::Update ( const VectorHandler * pSol ) const

virtual

Implements NonlinearProblem.

Definition at line 538 of file stepsol.cc.

References ASSERT, DEBUGCOUTFNAME, StepIntegrator::pDM, DataManager::Update(), UpdateDof(), and StepIntegrator::UpdateLoop().

Referenced by Jacobian().

 {
         DEBUGCOUTFNAME("StepNIntegrator::Update");
         ASSERT(pDM != NULL);
 
         UpdateLoop(this, &StepNIntegrator::UpdateDof, pSol);
         pDM->Update();
 }

Here is the call graph for this function:

void StepNIntegrator::UpdateDof	(	const int	DCount,
		const DofOrder::Order	Order,
		const VectorHandler *const	pSol = `0`
	)		const

protected

Definition at line 516 of file stepsol.cc.

References db0Algebraic, db0Differential, DofOrder::DIFFERENTIAL, VectorHandler::IncCoef(), ImplicitStepIntegrator::pXCurr, and ImplicitStepIntegrator::pXPrimeCurr.

Referenced by Update().

 {
         doublereal d = pSol->operator()(DCount);
         if (Order == DofOrder::DIFFERENTIAL) {
                 pXPrimeCurr->IncCoef(DCount, d);
 
                 /* Nota: b0Differential e b0Algebraic
                  * possono essere distinti;
                  * in ogni caso sono calcolati
                  * dalle funzioni di predizione
                  * e sono dati globali */
                 pXCurr->IncCoef(DCount, db0Differential*d);
 
         } else {
                 pXCurr->IncCoef(DCount, d);
                 pXPrimeCurr->IncCoef(DCount, db0Algebraic*d);
         }
 }

Here is the call graph for this function:

Member Data Documentation

doublereal StepNIntegrator::db0Algebraic

Definition at line 236 of file stepsol.h.

Referenced by CrankNicolsonIntegrator::SetCoef(), ImplicitEulerIntegrator::SetCoef(), MultistepSolver::SetCoef(), HopeSolver::SetCoef(), and UpdateDof().

doublereal StepNIntegrator::db0Differential

Definition at line 235 of file stepsol.h.

Referenced by CrankNicolsonIntegrator::dPredictState(), ImplicitEulerIntegrator::dPredictState(), CrankNicolsonIntegrator::dPredState(), ImplicitEulerIntegrator::dPredState(), CrankNicolsonIntegrator::dPredStateAlg(), ImplicitEulerIntegrator::dPredStateAlg(), Jacobian(), Residual(), CrankNicolsonIntegrator::SetCoef(), ImplicitEulerIntegrator::SetCoef(), MultistepSolver::SetCoef(), HopeSolver::SetCoef(), TestScale(), and UpdateDof().

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

Public Member Functions

Public Attributes

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation