#include <stepsol.h>

Inheritance diagram for Step1Integrator:

Collaboration diagram for Step1Integrator:

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

virtual	~Step1Integrator (void)

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)

Public Member Functions inherited from StepNIntegrator
	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)

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

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

virtual void	Predict (void)

virtual doublereal	dPredictDerivative (const doublereal &dXm1, const doublereal &dXPm1, DofOrder::Order o=DofOrder::DIFFERENTIAL) const =0

virtual doublereal	dPredictState (const doublereal &dXm1, const doublereal &dXP, const doublereal &dXPm1, DofOrder::Order o=DofOrder::DIFFERENTIAL) const =0

virtual doublereal	dPredDer (const doublereal &dXm1, const doublereal &dXPm1) const =0

virtual doublereal	dPredState (const doublereal &dXm1, const doublereal &dXP, const doublereal &dXPm1) const =0

virtual doublereal	dPredDerAlg (const doublereal &dXm1, const doublereal &dXPm1) const =0

virtual doublereal	dPredStateAlg (const doublereal &dXm1, const doublereal &dXP, const doublereal &dXPm1) const =0

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

Protected Attributes
VectorHandler *	pXPrev

VectorHandler *	pXPrimePrev

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

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

enum	StepChange { NEWSTEP, REPEATSTEP }

Public Attributes inherited from StepNIntegrator
doublereal	db0Differential

doublereal	db0Algebraic

Detailed Description

Definition at line 266 of file stepsol.h.

Constructor & Destructor Documentation

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

Definition at line 560 of file stepsol.cc.

References NO_OP.

 : StepNIntegrator(MaxIt, dT, dSolutionTol, 1, bmod_res_test),
 pXPrev(NULL),
 pXPrimePrev(NULL)
 {
         NO_OP;
 }

Step1Integrator::~Step1Integrator ( void )

virtual

Definition at line 571 of file stepsol.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

doublereal Step1Integrator::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
	)

virtual

Implements StepIntegrator.

Definition at line 623 of file stepsol.cc.

References DataManager::AfterConvergence(), DataManager::AfterPredict(), ASSERT, DataManager::dGetTime(), StepIntegrator::dSolTol, StepIntegrator::dTol, DataManager::iGetNumDofs(), DataManager::LinkToSolution(), StepIntegrator::MaxIters, StepIntegrator::outputPred, StepIntegrator::pDM, Solver::pGetNonlinearSolver(), Predict(), ImplicitStepIntegrator::pXCurr, pXPrev, ImplicitStepIntegrator::pXPrimeCurr, pXPrimePrev, StepNIntegrator::SetCoef(), and NonlinearSolver::Solve().

 {
         ASSERT(pDM != NULL);
         pXCurr  = pX;
         pXPrev  = qX[0];
 
         pXPrimeCurr  = pXPrime;
         pXPrimePrev  = qXPrime[0];
 
         /* predizione */
         SetCoef(TStep, dAph, StType);
         Predict();
         pDM->LinkToSolution(*pXCurr, *pXPrimeCurr);
         pDM->AfterPredict();
 
 #ifdef DEBUG
         integer iNumDofs = pDM->iGetNumDofs();
         if (outputPred) {
                 std::cout << "After prediction, time=" << pDM->dGetTime() << std::endl;
                 std::cout << "Dof:      |    XCurr  ";
                 for (unsigned idx = 0; idx < qX.size(); idx++) {
                         std::cout << "|  XPrev[" << idx << "] ";
                 }
                 std::cout << "|   XPrime  ";
                 for (unsigned idx = 0; idx < qXPrime.size(); idx++) {
                         std::cout << "| XPPrev[" << idx << "] ";
                 }
                 std::cout << "|" << std::endl;
                 for (int iTmpCnt = 1; iTmpCnt <= iNumDofs; iTmpCnt++) {
                         std::cout << std::setw(8) << iTmpCnt << ": ";
                         std::cout << std::setw(12) << pX->operator()(iTmpCnt);
                         for (unsigned int ivec = 0; ivec < qX.size(); ivec++) {
                                 std::cout << std::setw(12)
                                 << (qX[ivec])->operator()(iTmpCnt);
                         }
                         std::cout << std::setw(12) << pXPrime->operator()(iTmpCnt);
                         for (unsigned int ivec = 0; ivec < qXPrime.size(); ivec++) {
                                 std::cout << std::setw(12)
                                 << (qXPrime[ivec])->operator()(iTmpCnt);
                         }
                         std::cout << " " << pDM->DataManager::GetDofDescription(iTmpCnt) << std::endl;
                 }
         }
 #endif /* DEBUG */
 
         Err = 0.;
         pS->pGetNonlinearSolver()->Solve(this, pS, MaxIters, dTol,
                         EffIter, Err, dSolTol, SolErr);
 
         /* if it gets here, it surely converged */
         pDM->AfterConvergence();
 
         return Err;
 }

Here is the call graph for this function:

virtual doublereal Step1Integrator::dPredDer	(	const doublereal &	dXm1,
		const doublereal &	dXPm1
	)		const

protectedpure virtual

Implemented in ImplicitEulerIntegrator, and CrankNicolsonIntegrator.

Referenced by PredictDof().

virtual doublereal Step1Integrator::dPredDerAlg	(	const doublereal &	dXm1,
		const doublereal &	dXPm1
	)		const

protectedpure virtual

Implemented in ImplicitEulerIntegrator, and CrankNicolsonIntegrator.

Referenced by PredictDof().

virtual doublereal Step1Integrator::dPredictDerivative	(	const doublereal &	dXm1,
		const doublereal &	dXPm1,
		DofOrder::Order	o = `DofOrder::DIFFERENTIAL`
	)		const

protectedpure virtual

Implemented in ImplicitEulerIntegrator, and CrankNicolsonIntegrator.

virtual doublereal Step1Integrator::dPredictState	(	const doublereal &	dXm1,
		const doublereal &	dXP,
		const doublereal &	dXPm1,
		DofOrder::Order	o = `DofOrder::DIFFERENTIAL`
	)		const

protectedpure virtual

Implemented in ImplicitEulerIntegrator, and CrankNicolsonIntegrator.

virtual doublereal Step1Integrator::dPredState	(	const doublereal &	dXm1,
		const doublereal &	dXP,
		const doublereal &	dXPm1
	)		const

protectedpure virtual

Implemented in ImplicitEulerIntegrator, and CrankNicolsonIntegrator.

Referenced by PredictDof().

virtual doublereal Step1Integrator::dPredStateAlg	(	const doublereal &	dXm1,
		const doublereal &	dXP,
		const doublereal &	dXPm1
	)		const

protectedpure virtual

Implemented in ImplicitEulerIntegrator, and CrankNicolsonIntegrator.

Referenced by PredictDof().

void Step1Integrator::Predict ( void )

protectedvirtual

Definition at line 615 of file stepsol.cc.

References ASSERT, DEBUGCOUTFNAME, StepIntegrator::pDM, PredictDof(), and StepIntegrator::UpdateLoop().

Referenced by Advance().

 {
         DEBUGCOUTFNAME("Step1Integrator::Predict");
         ASSERT(pDM != NULL);
         UpdateLoop(this, &Step1Integrator::PredictDof);
 }

Here is the call graph for this function:

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

protected

Definition at line 581 of file stepsol.cc.

References DofOrder::ALGEBRAIC, DofOrder::DIFFERENTIAL, dPredDer(), dPredDerAlg(), dPredState(), dPredStateAlg(), MBDYN_EXCEPT_ARGS, VectorHandler::PutCoef(), ImplicitStepIntegrator::pXCurr, pXPrev, ImplicitStepIntegrator::pXPrimeCurr, and pXPrimePrev.

Referenced by Predict().

 {
         if (Order == DofOrder::DIFFERENTIAL) {
                 doublereal dXnm1 = pXPrev->operator()(DCount);
                 doublereal dXPnm1 = pXPrimePrev->operator()(DCount);
                 doublereal dXPn = dPredDer(dXnm1, dXPnm1);
                 doublereal dXn = dPredState(dXnm1, dXPn, dXPnm1);
                 pXPrimeCurr->PutCoef(DCount, dXPn);
                 pXCurr->PutCoef(DCount, dXn);
 
         } else if (Order == DofOrder::ALGEBRAIC) {
                 doublereal dXnm1 = pXPrev->operator()(DCount);
                 doublereal dXInm1 =
                         pXPrimePrev->operator()(DCount);
 
                 doublereal dXn = dPredDerAlg(dXInm1, dXnm1);
                 doublereal dXIn = dPredStateAlg(dXInm1, dXn, dXnm1);
 
                 pXCurr->PutCoef(DCount, dXn);
                 pXPrimeCurr->PutCoef(DCount, dXIn);
 
         } else {
                 silent_cerr("Step1Integrator::"
                         "PredictDof(): "
                         "unknown order for local dof "
                         << DCount << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 }