NOX Development
Loading...
Searching...
No Matches
LOCA::AdaptiveStepper Class Reference

Implementation of LOCA::Abstract::Iterator for computing points along a continuation curve. More...

#include <LOCA_AdaptiveStepper.H>

Inheritance diagram for LOCA::AdaptiveStepper:
Collaboration diagram for LOCA::AdaptiveStepper:

Public Member Functions

 AdaptiveStepper (const Teuchos::RCP< Teuchos::ParameterList > &pList, const Teuchos::RCP< ::Thyra::AdaptiveSolutionManager > &mgr, const Teuchos::RCP< LOCA::GlobalData > &global_data, const Teuchos::RCP< NOX::StatusTest::Generic > &nt)
virtual ~AdaptiveStepper ()
 Destructor.
bool eigensolverReset (Teuchos::RCP< Teuchos::ParameterList > &newEigensolverList)
 Replaces the eigensolver parameter list.
virtual Teuchos::RCP< const LOCA::MultiContinuation::AbstractGroupgetSolutionGroup () const
 Return the current solution group.
virtual Teuchos::RCP< const LOCA::MultiContinuation::AbstractGroupgetBifurcationGroup () const
 Return the current bifurcation group.
virtual Teuchos::RCP< const Teuchos::ParameterList > getList () const
 Return the output parameters from the stepper algorithm.
virtual Teuchos::RCP< const NOX::Solver::GenericgetSolver () const
 Return the current nonlinear solver pointer.
virtual double getContinuationParameter () const
 Return the current continuation parameter from the underlying LOCA::MultiContinuation::AbstractStrategy.
virtual LOCA::Abstract::Iterator::IteratorStatus run ()
 Run the iterator.
Public Member Functions inherited from LOCA::Abstract::Iterator
 Iterator (Teuchos::ParameterList &p)
 Constructor.
 Iterator (const Iterator &it)
 Copy Constructor.
virtual ~Iterator ()
 Destructor.
virtual bool resetIterator (Teuchos::ParameterList &p)
 Reset the iterator to start a new iteration.
virtual IteratorStatus getIteratorStatus () const
 Return the status of the iterator.
virtual int getStepNumber () const
 Returns the number of accepted steps.
virtual int getNumFailedSteps () const
 Returns the number of failed steps.
virtual int getNumTotalSteps () const
 Returns the total number of steps attempted.

Protected Member Functions

void setSolutionGroup (const Teuchos::RCP< LOCA::MultiContinuation::AbstractGroup > &initialGuess, const double value)
void buildLOCAFactory ()
 Build LOCA Factory needed to support stepping strategy.
virtual LOCA::Abstract::Iterator::IteratorStatus start ()
 Initialize iterator.
virtual LOCA::Abstract::Iterator::IteratorStatus finish (LOCA::Abstract::Iterator::IteratorStatus iteratorStatus)
 Finalize iterator.
virtual LOCA::Abstract::Iterator::StepStatus preprocess (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Preprocess step.
virtual LOCA::Abstract::Iterator::StepStatus adapt (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Adaptation step.
virtual LOCA::Abstract::Iterator::StepStatus compute (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Compute step.
virtual LOCA::Abstract::Iterator::StepStatus postprocess (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Postprocess step.
LOCA::Abstract::Iterator::IteratorStatus stop (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Check stopping criteria.
virtual LOCA::Abstract::Iterator::IteratorStatus iterate ()
 Perform iteration.
Teuchos::RCP< LOCA::MultiContinuation::AbstractGroupbuildConstrainedGroup (const Teuchos::RCP< LOCA::MultiContinuation::AbstractGroup > &grp)
 Builds a constrained group is specified through parameter lists.
LOCA::Abstract::Iterator::StepStatus computeStepSize (LOCA::Abstract::Iterator::StepStatus stepStatus, double &stepSize)
 Compute a new step size based on results from the previous step.
virtual void printInitializationInfo ()
 Print to the screen the initialization information.
virtual void printStartStep ()
 Print to the screen information for the beginning of a new continuation step.
virtual void printRelaxationStep ()
 Print to the screen information for the beginning of a restart relaxation step.
virtual void printEndStep (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Print to the screen information for the end of a successful continuation step.
virtual void printRelaxationEndStep (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Print to the screen information for the end of a relaxation step.
virtual void printEndInfo ()
 Print to the screen the concluding information.
virtual bool withinThreshold ()
 Checks to see if continuation parameter is within threshold of bound.
Protected Member Functions inherited from LOCA::Abstract::Iterator
 Iterator ()
 Constructor.
virtual void setLastIteration ()
 Set iterator status as LastIteration.
virtual bool isLastIteration ()
 Returns true of the iteratorStatus is LastIteration.
virtual StepStatus computeStepStatus (StepStatus preStatus, StepStatus compStatus, StepStatus postStatus)
 Compute step, iterator status from status of each component of step.

Protected Attributes

Teuchos::RCP< ::Thyra::AdaptiveSolutionManagermgr
 Adaptive solution manager provides interface to remesh/resize functions.
Teuchos::RCP< LOCA::GlobalDataglobalData
 Global data.
Teuchos::RCP< LOCA::Parameter::SublistParserparsedParams
 Parsed parameter list.
Teuchos::RCP< LOCA::MultiPredictor::AbstractStrategypredictor
 Predictor strategy.
Teuchos::RCP< LOCA::MultiContinuation::AbstractStrategycurGroupPtr
 Current continuation group.
Teuchos::RCP< LOCA::MultiContinuation::AbstractStrategyprevGroupPtr
 Continuation group from previous step.
Teuchos::RCP< LOCA::Eigensolver::AbstractStrategyeigensolver
 Eigensolver for computing eigenvalues after each continuation step.
Teuchos::RCP< LOCA::SaveEigenData::AbstractStrategysaveEigenData
 Strategy for saving eigenvectors/values.
Teuchos::RCP< LOCA::MultiContinuation::AbstractGroupbifGroupPtr
 Pointer to bifurcation group.
Teuchos::RCP< NOX::StatusTest::GenericnoxStatusTestPtr
 Pointer to the NOX status test.
Teuchos::RCP< Teuchos::ParameterList > paramListPtr
 Pointer to parameter list passed in constructor/reset method.
Teuchos::RCP< Teuchos::ParameterList > stepperList
 Pointer to stepper parameters.
Teuchos::RCP< NOX::Solver::GenericsolverPtr
 Solver (a.k.a corrector).
Teuchos::RCP< LOCA::MultiContinuation::ExtendedVectorcurPredictorPtr
 Current predictor direction.
Teuchos::RCP< LOCA::MultiContinuation::ExtendedVectorprevPredictorPtr
 previous predictor direction
Teuchos::RCP< LOCA::StepSize::AbstractStrategystepSizeStrategyPtr
 Step size control strategy.
Teuchos::RCP< Teuchos::ParameterList > bifurcationParams
 Other parameter lists.
Teuchos::RCP< Teuchos::ParameterList > firstStepperParams
std::string conParamName
 Name of continuation parameter.
std::vector< int > conParamIDs
 Index of continuation parameter.
double startValue
 Starting value of continuation parameter, $\lambda_b $.
double maxValue
 Largest allowable value of continuation parameter.
double minValue
 Smallest allowable value of continuation parameter.
double stepSize
 Current step size (change in the continuation parameter), $\Delta \lambda_i $.
int maxNonlinearSteps
 Maximum number of Newton iterations per continuation step, $ N_{max} $. Defaults to 15.
double targetValue
 Target parameter value for last step (either maxValue or minValue).
bool isTargetStep
 Flag indicating if this is an extra step to hit target value.
bool doTangentFactorScaling
 Flag indicating whether to do tangent factor step size scaling.
double tangentFactor
 Tangent factor.
double minTangentFactor
 Minimum tangent factor.
double tangentFactorExponent
 Tangent factor exponent.
bool calcEigenvalues
 Flag indicating whether to compute eigenvalues after each step.
bool return_failed_on_max_steps
 Flag indicating whether to return failed upon reaching max steps.
bool max_steps_exceeded
 Flag indicating that the failed step is only due to the maximum number of steps being exceeded.
Protected Attributes inherited from LOCA::Abstract::Iterator
int stepNumber
 Current step number of continuation algorithm (does not include failed steps).
int numFailedSteps
 Number of continuation steps that have failed.
int numTotalSteps
 Total number of steps attempeted (includes failed and successful steps).
int maxSteps
 Maximum number of continuation steps to take. Defaults to 100.
IteratorStatus iteratorStatus
 Current status of the iterator.

Additional Inherited Members

Public Types inherited from LOCA::Abstract::Iterator
enum  IteratorStatus { LastIteration = 2 , Finished = 1 , Failed = 0 , NotFinished = -1 }
 Enumerated type for status of the iterator. More...
enum  StepStatus { Successful = 1 , Unsuccessful = 0 , Provisional = -1 }
 Enumerated type for status of each step of iterator. More...

Detailed Description

Implementation of LOCA::Abstract::Iterator for computing points along a continuation curve.

The AdaptiveStepper class implements the pure virtual methods of the LOCA::Abstract::Iterator for iteratively computing points along a continuation curve.

Member Function Documentation

◆ adapt()

◆ buildConstrainedGroup()

Teuchos::RCP< LOCA::MultiContinuation::AbstractGroup > LOCA::AdaptiveStepper::buildConstrainedGroup ( const Teuchos::RCP< LOCA::MultiContinuation::AbstractGroup > & grp)
protected

Builds a constrained group is specified through parameter lists.

References LOCA::ParameterVector::getIndex(), globalData, and parsedParams.

◆ buildLOCAFactory()

void LOCA::AdaptiveStepper::buildLOCAFactory ( )
protected

Build LOCA Factory needed to support stepping strategy.

References eigensolver, globalData, parsedParams, predictor, saveEigenData, and stepSizeStrategyPtr.

Referenced by adapt().

◆ compute()

◆ computeStepSize()

◆ eigensolverReset()

bool LOCA::AdaptiveStepper::eigensolverReset ( Teuchos::RCP< Teuchos::ParameterList > & newEigensolverList)

Replaces the eigensolver parameter list.

References eigensolver, globalData, and parsedParams.

◆ finish()

◆ getBifurcationGroup()

Teuchos::RCP< const LOCA::MultiContinuation::AbstractGroup > LOCA::AdaptiveStepper::getBifurcationGroup ( ) const
virtual

Return the current bifurcation group.

If the current bifurcation method is "None", then the returned group is the same as getSolutionGroup(), otherwise this method returns the current bifurcation group (e.g., a turning point group).

References curGroupPtr.

◆ getContinuationParameter()

double LOCA::AdaptiveStepper::getContinuationParameter ( ) const
virtual

Return the current continuation parameter from the underlying LOCA::MultiContinuation::AbstractStrategy.

References curGroupPtr.

Referenced by iterate(), and run().

◆ getList()

Teuchos::RCP< const Teuchos::ParameterList > LOCA::AdaptiveStepper::getList ( ) const
virtual

Return the output parameters from the stepper algorithm.

References paramListPtr.

◆ getSolutionGroup()

Teuchos::RCP< const LOCA::MultiContinuation::AbstractGroup > LOCA::AdaptiveStepper::getSolutionGroup ( ) const
virtual

Return the current solution group.

References curGroupPtr.

◆ getSolver()

Teuchos::RCP< const NOX::Solver::Generic > LOCA::AdaptiveStepper::getSolver ( ) const
virtual

Return the current nonlinear solver pointer.

Will throw an error if the solver does not exist yet.

References globalData, and solverPtr.

◆ iterate()

◆ postprocess()

◆ preprocess()

◆ printEndStep()

void LOCA::AdaptiveStepper::printEndStep ( LOCA::Abstract::Iterator::StepStatus stepStatus)
protectedvirtual

Print to the screen information for the end of a successful continuation step.

References conParamName, curGroupPtr, globalData, prevGroupPtr, solverPtr, LOCA::Abstract::Iterator::stepNumber, NOX::Utils::StepperIteration, and LOCA::Abstract::Iterator::Successful.

Referenced by compute(), finish(), and start().

◆ printInitializationInfo()

void LOCA::AdaptiveStepper::printInitializationInfo ( )
protectedvirtual

Print to the screen the initialization information.

References globalData, LOCA::Abstract::Iterator::maxSteps, maxValue, minValue, startValue, NOX::Utils::StepperIteration, and stepperList.

◆ printRelaxationEndStep()

void LOCA::AdaptiveStepper::printRelaxationEndStep ( LOCA::Abstract::Iterator::StepStatus stepStatus)
protectedvirtual

Print to the screen information for the end of a relaxation step.

References conParamName, curGroupPtr, globalData, solverPtr, LOCA::Abstract::Iterator::stepNumber, NOX::Utils::StepperIteration, and LOCA::Abstract::Iterator::Successful.

Referenced by adapt().

◆ printRelaxationStep()

void LOCA::AdaptiveStepper::printRelaxationStep ( )
protectedvirtual

Print to the screen information for the beginning of a restart relaxation step.

References globalData, LOCA::Abstract::Iterator::stepNumber, and NOX::Utils::StepperIteration.

Referenced by adapt().

◆ printStartStep()

void LOCA::AdaptiveStepper::printStartStep ( )
protectedvirtual

Print to the screen information for the beginning of a new continuation step.

References conParamName, curGroupPtr, globalData, isTargetStep, prevGroupPtr, LOCA::Abstract::Iterator::stepNumber, NOX::Utils::StepperIteration, stepperList, stepSize, stepSizeStrategyPtr, and targetValue.

Referenced by compute(), finish(), and start().

◆ run()

◆ start()

◆ stop()

◆ withinThreshold()

bool LOCA::AdaptiveStepper::withinThreshold ( )
protectedvirtual

Checks to see if continuation parameter is within threshold of bound.

References curGroupPtr, parsedParams, stepperList, and targetValue.

Referenced by stop().

Member Data Documentation

◆ calcEigenvalues

bool LOCA::AdaptiveStepper::calcEigenvalues
protected

Flag indicating whether to compute eigenvalues after each step.

Referenced by adapt(), postprocess(), and start().

◆ conParamIDs

std::vector<int> LOCA::AdaptiveStepper::conParamIDs
protected

Index of continuation parameter.

Referenced by adapt(), finish(), and start().

◆ conParamName

std::string LOCA::AdaptiveStepper::conParamName
protected

Name of continuation parameter.

Referenced by printEndStep(), printRelaxationEndStep(), and printStartStep().

◆ curGroupPtr

◆ curPredictorPtr

Teuchos::RCP<LOCA::MultiContinuation::ExtendedVector> LOCA::AdaptiveStepper::curPredictorPtr
protected

Current predictor direction.

Referenced by adapt(), computeStepSize(), finish(), postprocess(), preprocess(), and start().

◆ doTangentFactorScaling

bool LOCA::AdaptiveStepper::doTangentFactorScaling
protected

Flag indicating whether to do tangent factor step size scaling.

Referenced by computeStepSize(), and postprocess().

◆ eigensolver

Teuchos::RCP<LOCA::Eigensolver::AbstractStrategy> LOCA::AdaptiveStepper::eigensolver
protected

Eigensolver for computing eigenvalues after each continuation step.

Referenced by buildLOCAFactory(), and eigensolverReset().

◆ globalData

◆ isTargetStep

bool LOCA::AdaptiveStepper::isTargetStep
protected

Flag indicating if this is an extra step to hit target value.

Referenced by finish(), and printStartStep().

◆ max_steps_exceeded

bool LOCA::AdaptiveStepper::max_steps_exceeded
protected

Flag indicating that the failed step is only due to the maximum number of steps being exceeded.

Referenced by run(), and stop().

◆ maxValue

double LOCA::AdaptiveStepper::maxValue
protected

Largest allowable value of continuation parameter.

Referenced by computeStepSize(), printInitializationInfo(), and stop().

◆ mgr

Teuchos::RCP< ::Thyra::AdaptiveSolutionManager > LOCA::AdaptiveStepper::mgr
protected

Adaptive solution manager provides interface to remesh/resize functions.

Referenced by adapt(), iterate(), and run().

◆ minTangentFactor

double LOCA::AdaptiveStepper::minTangentFactor
protected

Minimum tangent factor.

Referenced by postprocess().

◆ minValue

double LOCA::AdaptiveStepper::minValue
protected

Smallest allowable value of continuation parameter.

Referenced by computeStepSize(), printInitializationInfo(), and stop().

◆ noxStatusTestPtr

Teuchos::RCP<NOX::StatusTest::Generic> LOCA::AdaptiveStepper::noxStatusTestPtr
protected

Pointer to the NOX status test.

Referenced by adapt(), finish(), preprocess(), and start().

◆ paramListPtr

Teuchos::RCP<Teuchos::ParameterList> LOCA::AdaptiveStepper::paramListPtr
protected

Pointer to parameter list passed in constructor/reset method.

Referenced by getList().

◆ parsedParams

Teuchos::RCP<LOCA::Parameter::SublistParser> LOCA::AdaptiveStepper::parsedParams
protected

◆ predictor

Teuchos::RCP<LOCA::MultiPredictor::AbstractStrategy> LOCA::AdaptiveStepper::predictor
protected

Predictor strategy.

Referenced by adapt(), buildLOCAFactory(), finish(), and start().

◆ prevGroupPtr

Teuchos::RCP<LOCA::MultiContinuation::AbstractStrategy> LOCA::AdaptiveStepper::prevGroupPtr
protected

Continuation group from previous step.

Referenced by adapt(), finish(), preprocess(), printEndStep(), printStartStep(), start(), and stop().

◆ prevPredictorPtr

Teuchos::RCP<LOCA::MultiContinuation::ExtendedVector> LOCA::AdaptiveStepper::prevPredictorPtr
protected

previous predictor direction

Referenced by adapt(), postprocess(), and start().

◆ return_failed_on_max_steps

bool LOCA::AdaptiveStepper::return_failed_on_max_steps
protected

Flag indicating whether to return failed upon reaching max steps.

Referenced by stop().

◆ saveEigenData

Teuchos::RCP<LOCA::SaveEigenData::AbstractStrategy> LOCA::AdaptiveStepper::saveEigenData
protected

Strategy for saving eigenvectors/values.

Referenced by buildLOCAFactory().

◆ solverPtr

Teuchos::RCP<NOX::Solver::Generic> LOCA::AdaptiveStepper::solverPtr
protected

◆ startValue

double LOCA::AdaptiveStepper::startValue
protected

Starting value of continuation parameter, $\lambda_b $.

Referenced by printInitializationInfo().

◆ stepperList

Teuchos::RCP<Teuchos::ParameterList> LOCA::AdaptiveStepper::stepperList
protected

Pointer to stepper parameters.

Referenced by adapt(), finish(), printInitializationInfo(), printStartStep(), start(), and withinThreshold().

◆ stepSize

double LOCA::AdaptiveStepper::stepSize
protected

Current step size (change in the continuation parameter), $\Delta \lambda_i $.

Referenced by adapt(), finish(), preprocess(), printStartStep(), and start().

◆ stepSizeStrategyPtr

Teuchos::RCP<LOCA::StepSize::AbstractStrategy> LOCA::AdaptiveStepper::stepSizeStrategyPtr
protected

Step size control strategy.

Referenced by buildLOCAFactory(), computeStepSize(), and printStartStep().

◆ tangentFactor

double LOCA::AdaptiveStepper::tangentFactor
protected

Tangent factor.

Referenced by computeStepSize(), and postprocess().

◆ tangentFactorExponent

double LOCA::AdaptiveStepper::tangentFactorExponent
protected

Tangent factor exponent.

Referenced by computeStepSize().

◆ targetValue

double LOCA::AdaptiveStepper::targetValue
protected

Target parameter value for last step (either maxValue or minValue).

Referenced by computeStepSize(), finish(), printStartStep(), stop(), and withinThreshold().


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