ImplicitStepIntegrator Class Reference

#include <stepsol.h>

Inheritance diagram for ImplicitStepIntegrator:

Collaboration diagram for ImplicitStepIntegrator:

Public Member Functions
	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
virtual doublereal	TestScale (const NonlinearSolverTest *pTest, doublereal &dCoef) 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)
virtual void	Residual (VectorHandler *pRes) const =0
virtual void	Jacobian (MatrixHandler *pJac) const =0
virtual void	Update (const VectorHandler *pSol) const =0

Protected Attributes
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

Private Attributes
MyVectorHandler	XTau
MyVectorHandler	SavedState
MyVectorHandler	SavedDerState
bool	bEvalProdCalledFirstTime

Additional Inherited Members
Public Types inherited from StepIntegrator
enum	{ DIFFERENTIAL = 0, ALGEBRAIC = 1 }
enum	StepChange { NEWSTEP, REPEATSTEP }
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

Detailed Description

Definition at line 151 of file stepsol.h.

Constructor & Destructor Documentation

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

Definition at line 124 of file stepsol.cc.

References NO_OP.

: StepIntegrator(MaxIt, dT, dSolutionTol, stp, sts),
bEvalProdCalledFirstTime(true),
pXCurr(0),
pXPrimeCurr(0),
bModResTest(bmod_res_test)
{
        NO_OP;
}

ImplicitStepIntegrator::~ImplicitStepIntegrator ( void  )

virtual

Definition at line 139 of file stepsol.cc.

References NO_OP.

{
        NO_OP;
}

Member Function Documentation

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

virtual

Implements NonlinearProblem.

Definition at line 145 of file stepsol.cc.

References ASSERT, bEvalProdCalledFirstTime, copysign(), f0, grad::fabs(), VectorHandler::iGetSize(), VectorHandler::InnerProd(), VectorHandler::Norm(), StepIntegrator::pDM, pXCurr, pXPrimeCurr, VectorHandler::Reset(), MyVectorHandler::Reset(), NonlinearProblem::Residual(), MyVectorHandler::Resize(), SavedDerState, SavedState, VectorHandler::ScalarMul(), MyVectorHandler::ScalarMul(), NonlinearProblem::Update(), DataManager::Update(), and XTau.

{
        /* matrix-free product
         *
         * J(XCurr) * w = -||w|| * (Res(XCurr + sigma * Tau * w/||w||) - f0) / (sigma * Tau)
         *
         */
        if (bEvalProdCalledFirstTime) {
                XTau.Resize(w.iGetSize());
                SavedState.Resize(w.iGetSize());
                SavedDerState.Resize(w.iGetSize());
                bEvalProdCalledFirstTime = false;
        }

        SavedState = *pXCurr;
        SavedDerState = *pXPrimeCurr;

        /* if w = 0; J * w = 0 */
        ASSERT(pDM != NULL);

        doublereal nw = w.Norm();
        if (nw < std::numeric_limits<doublereal>::epsilon()) {
                z.Reset();
                return;
        }
        doublereal sigma = pXCurr->InnerProd(w);
        sigma /=  nw;
        if (fabs(sigma) > std::numeric_limits<doublereal>::epsilon()) {
                doublereal xx = (fabs( sigma) <= 1.) ? 1. : fabs(sigma);
                Tau = copysign(Tau*xx, sigma);
        }
        Tau /= nw;
#ifdef DEBUG_ITERATIVE
        std::cout << "Tau " << Tau << std::endl;
#endif /* DEBUG_ITERATIVE */

        XTau.Reset();
        z.Reset();
        XTau.ScalarMul(w, Tau);
        Update(&XTau);
#ifdef  USE_EXTERNAL
        External::SendFreeze();
#endif /* USE_EXTERNAL */
        /* deal with throwing elements: do not honor their requests while perfoming matrix free update */
        try {
                Residual(&z);
        }
        catch (DataManager::ChangedEquationStructure) {
        }
        XTau.ScalarMul(XTau, -1.);

        /* riporta tutto nelle condizioni inziali */
#if 0
        Update(&XTau);
#endif
        *pXCurr = SavedState;
        *pXPrimeCurr = SavedDerState;
        pDM->Update();
        z -= f0;
        z.ScalarMul(z, -1./Tau);
}

Here is the call graph for this function:

doublereal ImplicitStepIntegrator::TestScale ( const NonlinearSolverTest *  pTest, doublereal &  dCoef  ) const

virtual

Implements NonlinearProblem.

Reimplemented in StepNIntegrator, and DerivativeSolver.

Definition at line 210 of file stepsol.cc.

References ASSERT, bModResTest, DofOrder::DIFFERENTIAL, NonlinearSolverTest::dScaleCoef(), VectorHandler::iGetSize(), StepIntegrator::pDofs, and pXPrimeCurr.

Referenced by StepNIntegrator::TestScale().

{
        dCoef = 1.;

        if (bModResTest) {
#ifdef USE_MPI
#warning "ImplicitStepIntegrator::TestScale() not available with Schur solution"
#endif /* USE_MPI */

                doublereal dXPr = 0.;

                DataManager::DofIterator_const CurrDof = pDofs->begin();

                for (int iCntp1 = 1; iCntp1 <= pXPrimeCurr->iGetSize();
                        iCntp1++, ++CurrDof)
                {
                        ASSERT(CurrDof != pDofs->end());

                        if (CurrDof->Order == DofOrder::DIFFERENTIAL) {
                                doublereal d = pXPrimeCurr->operator()(iCntp1);
                                doublereal d2 = d*d;

                                doublereal ds = pTest->dScaleCoef(iCntp1);
                                doublereal ds2 = ds*ds;
                                d2 *= ds2;

                                dXPr += d2;
                        }
                        /* else if ALGEBRAIC: non aggiunge nulla */
                }

                return 1./(1. + dXPr);

        } else {
                return 1.;
        }
}

Here is the call graph for this function:

Member Data Documentation

bool ImplicitStepIntegrator::bEvalProdCalledFirstTime

mutableprivate

Definition at line 160 of file stepsol.h.

Referenced by EvalProd().

bool ImplicitStepIntegrator::bModResTest

protected

Definition at line 165 of file stepsol.h.

Referenced by TestScale().

VectorHandler* ImplicitStepIntegrator::pXCurr

protected

Definition at line 163 of file stepsol.h.

Referenced by ThirdOrderIntegrator::Advance(), DerivativeSolver::Advance(), Step1Integrator::Advance(), Step2Integrator::Advance(), EvalProd(), ThirdOrderIntegrator::Jacobian(), ThirdOrderIntegrator::Predict(), Step1Integrator::PredictDof(), Step2Integrator::PredictDof(), ThirdOrderIntegrator::PredictDof_for_AfterPredict(), ThirdOrderIntegrator::RealPredictDof(), ThirdOrderIntegrator::Residual(), ThirdOrderIntegrator::UpdateDof(), DerivativeSolver::UpdateDof(), and StepNIntegrator::UpdateDof().

VectorHandler* ImplicitStepIntegrator::pXPrimeCurr

protected

Definition at line 164 of file stepsol.h.

Referenced by ThirdOrderIntegrator::Advance(), DerivativeSolver::Advance(), Step1Integrator::Advance(), Step2Integrator::Advance(), EvalProd(), ThirdOrderIntegrator::Jacobian(), ThirdOrderIntegrator::Predict(), Step1Integrator::PredictDof(), Step2Integrator::PredictDof(), ThirdOrderIntegrator::PredictDof_for_AfterPredict(), ThirdOrderIntegrator::RealPredictDof(), ThirdOrderIntegrator::Residual(), TestScale(), ThirdOrderIntegrator::UpdateDof(), DerivativeSolver::UpdateDof(), and StepNIntegrator::UpdateDof().

MyVectorHandler ImplicitStepIntegrator::SavedDerState

mutableprivate

Definition at line 159 of file stepsol.h.

Referenced by EvalProd().

MyVectorHandler ImplicitStepIntegrator::SavedState

mutableprivate

Definition at line 158 of file stepsol.h.

Referenced by EvalProd().

MyVectorHandler ImplicitStepIntegrator::XTau

mutableprivate

Definition at line 157 of file stepsol.h.

Referenced by EvalProd().

---

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

• stepsol.h
• stepsol.cc