LOCA::TurningPoint::MinimallyAugmented::ModifiedConstraint Class Reference

Implementation of LOCA::MultiContinuation::ConstraintInterfaceMVDX for computing turning points for the minimally augmented turning point formulation. More...

#include <LOCA_TurningPoint_MinimallyAugmented_ModifiedConstraint.H>

Inheritance diagram for LOCA::TurningPoint::MinimallyAugmented::ModifiedConstraint:

[legend]Collaboration diagram for LOCA::TurningPoint::MinimallyAugmented::ModifiedConstraint:[legend]List of all members.


Public Member Functions
	ModifiedConstraint (const Teuchos::RefCountPtr< LOCA::GlobalData > &global_data, const Teuchos::RefCountPtr< LOCA::Parameter::SublistParser > &topParams, const Teuchos::RefCountPtr< Teuchos::ParameterList > &tpParams, const Teuchos::RefCountPtr< LOCA::TurningPoint::MinimallyAugmented::AbstractGroup > &g, bool is_symmetric, const NOX::Abstract::Vector &a, const NOX::Abstract::Vector *b, int bif_param)
	Constructor.
	ModifiedConstraint (const ModifiedConstraint &source, NOX::CopyType type=NOX::DeepCopy)
	Copy constructor.
virtual	~ModifiedConstraint ()
	Destructor.
void	setNewtonUpdates (const NOX::Abstract::Vector &dx, double dp, double step)
Implementation of LOCA::MultiContinuation::ConstraintInterface
virtual methods
virtual void	copy (const LOCA::MultiContinuation::ConstraintInterface &source)
	Copy.
virtual Teuchos::RefCountPtr< LOCA::MultiContinuation::ConstraintInterface >	clone (NOX::CopyType type=NOX::DeepCopy) const
	Cloning function.
virtual NOX::Abstract::Group::ReturnType	computeConstraints ()
	Compute continuation constraint equations.
virtual void	preProcessContinuationStep (LOCA::Abstract::Iterator::StepStatus stepStatus)
	Perform any preprocessing before a continuation step starts.
virtual void	postProcessContinuationStep (LOCA::Abstract::Iterator::StepStatus stepStatus)
	Perform any postprocessing after a continuation step finishes.
Protected Attributes
Teuchos::RefCountPtr< NOX::Abstract::MultiVector >	w_vector_update
	Stores update to left null vector.
Teuchos::RefCountPtr< NOX::Abstract::MultiVector >	v_vector_update
	Stores update to right null vector.
Teuchos::RefCountPtr< NOX::Abstract::MultiVector >	w_residual
	Stores left null vector residual.
Teuchos::RefCountPtr< NOX::Abstract::MultiVector >	v_residual
	Stores right null vector residual.
Teuchos::RefCountPtr< NOX::Abstract::MultiVector >	deltaX
	Stores solution update.
NOX::Abstract::MultiVector::DenseMatrix	sigma1
	Stores sigma_1.
NOX::Abstract::MultiVector::DenseMatrix	sigma2
	Stores sigma_1.
double	deltaP
	Stores parameter update.
bool	isFirstSolve
bool	includeNewtonTerms
	Flag indicating whether to include the newton update terms.

Detailed Description

Implementation of LOCA::MultiContinuation::ConstraintInterfaceMVDX for computing turning points for the minimally augmented turning point formulation.

This class is a modification of LOCA::TurningPoint::MinimallyAugmented::Constraint where updates are computed to the left and right null vectors $w$ and $v$ every nonlinear iteration instead of solving for them directly:

$\begin{bmatrix} J & a \\ b^T & 0 \end{bmatrix} \begin{bmatrix} \Delta v \\ \Delta \sigma_1 \end{bmatrix} = - \begin{bmatrix} J v + \sigma_1 a + (Jv)_x\Delta x + (Jv)_p \Delta p\\ b^T a - n \end{bmatrix},$

$\begin{bmatrix} J^T & b \\ a^T & 0 \end{bmatrix} \begin{bmatrix} \Delta w \\ \Delta \sigma_2 \end{bmatrix} = - \begin{bmatrix} J^T w + \sigma_2 b + (J^T w)_x \Delta x + (J^T w)_p \Delta p \\ a^T w - n \end{bmatrix},$

The class is intialized via the tpParams parameter list argument to the constructor. This class recognizes all paramters for LOCA::TurningPoint::MinimallyAugmented::Constraint plus the following:

"Include Newton Terms" -- [bool] (default: false) - Flag indicating whether to include the $\Delta x$ and $\Delta p$ terms above when computing the null vector updates

Member Function Documentation

void LOCA::TurningPoint::MinimallyAugmented::ModifiedConstraint::postProcessContinuationStep ( LOCA::Abstract::Iterator::StepStatus stepStatus ) [virtual]

Perform any postprocessing after a continuation step finishes.
The stepStatus argument indicates whether the step was successful. Here we set up the constraint class to solve for $w$ and $v$ for the first nonlinear iteration.
Reimplemented from LOCA::TurningPoint::MinimallyAugmented::Constraint.

void LOCA::TurningPoint::MinimallyAugmented::ModifiedConstraint::preProcessContinuationStep ( LOCA::Abstract::Iterator::StepStatus stepStatus ) [virtual]

Perform any preprocessing before a continuation step starts.
The stepStatus argument indicates whether the previous step was successful. Here we set up the constraint class to solve for $w$ and $v$ for the first nonlinear iteration.
Reimplemented from LOCA::MultiContinuation::ConstraintInterface.

Member Data Documentation

bool LOCA::TurningPoint::MinimallyAugmented::ModifiedConstraint::isFirstSolve [protected]

Flag that indicates whether we're in the first solve per continuation step.

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

LOCA_TurningPoint_MinimallyAugmented_ModifiedConstraint.H
LOCA_TurningPoint_MinimallyAugmented_ModifiedConstraint.C

Generated on Thu Sep 18 12:38:36 2008 for NOX by

1.3.9.1

LOCA::TurningPoint::MinimallyAugmented::ModifiedConstraint Class Reference

Public Member Functions

Protected Attributes

Detailed Description

Member Function Documentation

Member Data Documentation