Step2Integrator Class Reference

#include <stepsol.h>

Inheritance diagram for Step2Integrator:

Collaboration diagram for Step2Integrator:

Public Member Functions
	Step2Integrator (const integer MaxIt, const doublereal dT, const doublereal dSolutionTol, const bool bmod_res_test)
virtual	~Step2Integrator (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 &dXm2, const doublereal &dXPm1, const doublereal &dXPm2, DofOrder::Order o=DofOrder::DIFFERENTIAL) const =0
virtual doublereal	dPredictState (const doublereal &dXm1, const doublereal &dXm2, const doublereal &dXP, const doublereal &dXPm1, const doublereal &dXPm2, DofOrder::Order o=DofOrder::DIFFERENTIAL) const =0
virtual doublereal	dPredDer (const doublereal &dXm1, const doublereal &dXm2, const doublereal &dXPm1, const doublereal &dXPm2) const =0
virtual doublereal	dPredState (const doublereal &dXm1, const doublereal &dXm2, const doublereal &dXP, const doublereal &dXPm1, const doublereal &dXPm2) const =0
virtual doublereal	dPredDerAlg (const doublereal &dXm1, const doublereal &dXPm1, const doublereal &dXPm2) const =0
virtual doublereal	dPredStateAlg (const doublereal &dXm1, const doublereal &dXP, const doublereal &dXPm1, const doublereal &dXPm2) const =0
virtual void	SetCoef (doublereal dT, doublereal dAlpha, enum StepChange NewStep)=0
Protected Member Functions inherited from StepNIntegrator
void	UpdateDof (const int DCount, const DofOrder::Order Order, const VectorHandler *const pSol=0) const
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 *	pXPrev2
VectorHandler *	pXPrimePrev
VectorHandler *	pXPrimePrev2
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 430 of file stepsol.h.

Constructor & Destructor Documentation

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

Definition at line 692 of file stepsol.cc.

References NO_OP.

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

Step2Integrator::~Step2Integrator ( void  )

virtual

Definition at line 705 of file stepsol.cc.

References NO_OP.

{
        NO_OP;
}

Member Function Documentation

doublereal Step2Integrator::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 765 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, pXPrev2, ImplicitStepIntegrator::pXPrimeCurr, pXPrimePrev, pXPrimePrev2, SetCoef(), and NonlinearSolver::Solve().

{
        ASSERT(pDM != NULL);
        pXCurr  = pX;
        pXPrev  = qX[0];
        pXPrev2 = qX[1];

        pXPrimeCurr  = pXPrime;
        pXPrimePrev  = qXPrime[0];
        pXPrimePrev2 = qXPrime[1];

        /* 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 Step2Integrator::dPredDer ( const doublereal &  dXm1, const doublereal &  dXm2, const doublereal &  dXPm1, const doublereal &  dXPm2  ) const

protectedpure virtual

Implemented in HopeSolver, and MultistepSolver.

Referenced by PredictDof().

virtual doublereal Step2Integrator::dPredDerAlg ( const doublereal &  dXm1, const doublereal &  dXPm1, const doublereal &  dXPm2  ) const

protectedpure virtual

Implemented in HopeSolver, and MultistepSolver.

Referenced by PredictDof().

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

protectedpure virtual

Implemented in HopeSolver, and MultistepSolver.

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

protectedpure virtual

Implemented in HopeSolver, and MultistepSolver.

virtual doublereal Step2Integrator::dPredState ( const doublereal &  dXm1, const doublereal &  dXm2, const doublereal &  dXP, const doublereal &  dXPm1, const doublereal &  dXPm2  ) const

protectedpure virtual

Implemented in HopeSolver, and MultistepSolver.

Referenced by PredictDof().

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

protectedpure virtual

Implemented in HopeSolver, and MultistepSolver.

Referenced by PredictDof().

void Step2Integrator::Predict ( void  )

protectedvirtual

Definition at line 757 of file stepsol.cc.

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

Referenced by Advance().

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

Here is the call graph for this function:

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

protected

Definition at line 714 of file stepsol.cc.

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

Referenced by Predict().

{
        if (Order == DofOrder::DIFFERENTIAL) {
                doublereal dXnm1 = pXPrev->operator()(DCount);
                doublereal dXnm2 = pXPrev2->operator()(DCount);
                doublereal dXPnm1 = pXPrimePrev->operator()(DCount);

                doublereal dXPnm2 =
                        pXPrimePrev2->operator()(DCount);
                doublereal dXPn = dPredDer(dXnm1, dXnm2,
                                dXPnm1, dXPnm2);
                doublereal dXn = dPredState(dXnm1, dXnm2,
                                dXPn, dXPnm1, dXPnm2);

                pXPrimeCurr->PutCoef(DCount, dXPn);
                pXCurr->PutCoef(DCount, dXn);

        } else if (Order == DofOrder::ALGEBRAIC) {
                doublereal dXnm1 = pXPrev->operator()(DCount);
                doublereal dXnm2 = pXPrev2->operator()(DCount);
                doublereal dXInm1 =
                        pXPrimePrev->operator()(DCount);

                doublereal dXn = dPredDerAlg(dXInm1,
                                dXnm1, dXnm2);
                doublereal dXIn = dPredStateAlg(dXInm1,
                                dXn, dXnm1, dXnm2);

                pXCurr->PutCoef(DCount, dXn);
                pXPrimeCurr->PutCoef(DCount, dXIn);

        } else {
                silent_cerr("Step2Integrator::"
                        "PredictDof(): "
                        "unknown order for local dof "
                        << DCount << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
}

Here is the call graph for this function:

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

protectedpure virtual

Implements StepNIntegrator.

Implemented in HopeSolver, and MultistepSolver.

Referenced by Advance().

Member Data Documentation

VectorHandler* Step2Integrator::pXPrev

protected

Definition at line 434 of file stepsol.h.

Referenced by Advance(), and PredictDof().

VectorHandler * Step2Integrator::pXPrev2

protected

Definition at line 434 of file stepsol.h.

Referenced by Advance(), and PredictDof().

VectorHandler* Step2Integrator::pXPrimePrev

protected

Definition at line 435 of file stepsol.h.

Referenced by Advance(), and PredictDof().

VectorHandler * Step2Integrator::pXPrimePrev2

protected

Definition at line 435 of file stepsol.h.

Referenced by Advance(), and PredictDof().

---

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

• stepsol.h
• stepsol.cc