LOCA::Hopf::MinimallyAugmented::Constraint Class Reference

Implementation of LOCA::MultiContinuation::ConstraintInterfaceMVDX for computing Hopf bifurcations for the minimally augmented Hopf formulation. More...

#include <LOCA_Hopf_MinimallyAugmented_Constraint.H>

Inheritance diagram for LOCA::Hopf::MinimallyAugmented::Constraint:

Inheritance graph
[legend]
Collaboration diagram for LOCA::Hopf::MinimallyAugmented::Constraint:

Collaboration graph
[legend]

List of all members.

Public Member Functions

 Constraint (const Teuchos::RCP< LOCA::GlobalData > &global_data, const Teuchos::RCP< LOCA::Parameter::SublistParser > &topParams, const Teuchos::RCP< Teuchos::ParameterList > &hpfParams, const Teuchos::RCP< LOCA::Hopf::MinimallyAugmented::AbstractGroup > &g, bool is_symmetric, const NOX::Abstract::Vector &a_real, const NOX::Abstract::Vector &a_imag, const NOX::Abstract::Vector *b_real, const NOX::Abstract::Vector *b_imag, int bif_param, double freq)
 Constructor.
 Constraint (const Constraint &source, NOX::CopyType type=NOX::DeepCopy)
 Copy constructor.
virtual ~Constraint ()
 Destructor.
virtual void setGroup (const Teuchos::RCP< LOCA::Hopf::MinimallyAugmented::AbstractGroup > &g)
 Set the group pointer.
virtual void setFrequency (double freq)
 Set Hopf frequency.
virtual Teuchos::RCP< const
NOX::Abstract::Vector
getLeftNullVecReal () const
 Returns real component of left null vector w.
virtual Teuchos::RCP< const
NOX::Abstract::Vector
getLeftNullVecImag () const
 Returns imaginary component of left null vector w.
virtual Teuchos::RCP< const
NOX::Abstract::Vector
getRightNullVecReal () const
 Returns real component of right null vector v.
virtual Teuchos::RCP< const
NOX::Abstract::Vector
getRightNullVecImag () const
 Returns imaginary component of right null vector v.
virtual double getSigmaReal () const
 Returns real component of sigma.
virtual double getSigmaImag () const
 Returns imaginary component of sigma.
virtual
NOX::Abstract::Group::ReturnType 
computeDOmega (NOX::Abstract::MultiVector::DenseMatrix &domega)
 Compute derivative of sigma w.r.t. frequency omega.
Implementation of LOCA::MultiContinuation::ConstraintInterface
virtual methods

virtual void copy (const LOCA::MultiContinuation::ConstraintInterface &source)
 Copy.
virtual Teuchos::RCP
< LOCA::MultiContinuation::ConstraintInterface
clone (NOX::CopyType type=NOX::DeepCopy) const
 Cloning function.
virtual int numConstraints () const
 Return number of constraints.
virtual void setX (const NOX::Abstract::Vector &y)
 Set the solution vector to y.
virtual void setParam (int paramID, double val)
 Sets parameter indexed by paramID.
virtual void setParams (const vector< int > &paramIDs, const NOX::Abstract::MultiVector::DenseMatrix &vals)
 Sets parameters indexed by paramIDs.
virtual
NOX::Abstract::Group::ReturnType 
computeConstraints ()
 Compute continuation constraint equations.
virtual
NOX::Abstract::Group::ReturnType 
computeDX ()
 Compute derivative of constraints w.r.t. solution vector x.
virtual
NOX::Abstract::Group::ReturnType 
computeDP (const vector< int > &paramIDs, NOX::Abstract::MultiVector::DenseMatrix &dgdp, bool isValidG)
 Compute derivative of constraints w.r.t. supplied parameters.
virtual bool isConstraints () const
 Return true if constraint residuals are valid.
virtual bool isDX () const
 Return true if derivatives of constraints w.r.t. x are valid.
virtual const
NOX::Abstract::MultiVector::DenseMatrix
getConstraints () const
 Return constraint residuals.
virtual const
NOX::Abstract::MultiVector
getDX () const
 Return solution component of constraint derivatives.
virtual bool isDXZero () const
 Return true if solution component of constraint derivatives is zero.
virtual void postProcessContinuationStep (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Perform any postprocessing after a continuation step finishes.

Protected Attributes

Teuchos::RCP< LOCA::GlobalDataglobalData
 Pointer LOCA global data object.
Teuchos::RCP
< LOCA::Parameter::SublistParser
parsedParams
 Parsed top-level parameters.
Teuchos::RCP
< Teuchos::ParameterList > 
hopfParams
 Bifurcation parameter list.
Teuchos::RCP
< LOCA::Hopf::MinimallyAugmented::AbstractGroup
grpPtr
 Pointer to base group that defines $F$.
Teuchos::RCP
< NOX::Abstract::MultiVector
a_vector
 Vector for $a$.
Teuchos::RCP
< NOX::Abstract::MultiVector
b_vector
 Vector for $b$.
Teuchos::RCP
< NOX::Abstract::MultiVector
w_vector
 Stores left null vector.
Teuchos::RCP
< NOX::Abstract::MultiVector
v_vector
 Stores right null vector.
Teuchos::RCP
< NOX::Abstract::MultiVector
Cv_vector
 Stores C*v.
Teuchos::RCP
< NOX::Abstract::MultiVector
sigma_x
 Stores sigma_x.
NOX::Abstract::MultiVector::DenseMatrix constraints
 Constraint values.
Teuchos::RCP
< LOCA::BorderedSolver::AbstractStrategy
borderedSolver
 Stores bordered solver strategy.
double dn
 Stores vector length as a double.
double sigma_scale
 Stores scale factor on sigma.
bool isSymmetric
 Flag indicating whether Jacobian is symmetric.
bool isValidConstraints
 Flag indicating whether constraints are valid.
bool isValidDX
 Flag indicating whether sigma_x is valid.
vector< int > bifParamID
 Stores the bifurcation parameter index.
double omega
 Stores Hopf frequency.
bool updateVectorsEveryContinuationStep
 Flag indicating whether to update $a$ and $b$ every continuation step.
bool updateVectorsEveryIteration
 Flag indicating whether to update $a$ and $b$ every nonlinear iteration.

Private Member Functions

Constraintoperator= (const Constraint &source)
 Prohibit generation and use of operator=().


Detailed Description

Implementation of LOCA::MultiContinuation::ConstraintInterfaceMVDX for computing Hopf bifurcations for the minimally augmented Hopf formulation.

This class implements the turning point constraint equation $\sigma(x,p,\omega) = 0$ for the minimally augmented Hopf formulation where $\sigma$ is defined via

\[ \begin{bmatrix} J+i\omega M & a \\ b^H & 0 \end{bmatrix} \begin{bmatrix} v \\ \sigma_1 \end{bmatrix} = \begin{bmatrix} 0 \\ n \end{bmatrix}, \]

\[ \begin{bmatrix} J^T-i\omega M^T & b \\ a^H & 0 \end{bmatrix} \begin{bmatrix} w \\ \sigma_2 \end{bmatrix} = \begin{bmatrix} 0 \\ n \end{bmatrix}, \]

\[ \sigma = w^H J+i\omega M v/n \]

for any vectors $a$ and $b$ in $\mathbb{C}^n$. Using these relationships, it is easy to show

\[ \begin{split} \sigma_x &= (w^H(J+i\omega M)v)_x/n = w^H(J+i\omega M)_x v/n \\ \sigma_p &= (w^H(J+i\omega M)v)_p/n = w^H(J+i\omega M)_p v/n \end{split} \]

The class is intialized via the hpfParams parameter list argument to the constructor. The parameters this class recognizes are:

Definition at line 132 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.


Constructor & Destructor Documentation

LOCA::Hopf::MinimallyAugmented::Constraint::Constraint ( const Teuchos::RCP< LOCA::GlobalData > &  global_data,
const Teuchos::RCP< LOCA::Parameter::SublistParser > &  topParams,
const Teuchos::RCP< Teuchos::ParameterList > &  hpfParams,
const Teuchos::RCP< LOCA::Hopf::MinimallyAugmented::AbstractGroup > &  g,
bool  is_symmetric,
const NOX::Abstract::Vector a_real,
const NOX::Abstract::Vector a_imag,
const NOX::Abstract::Vector b_real,
const NOX::Abstract::Vector b_imag,
int  bif_param,
double  freq 
)

LOCA::Hopf::MinimallyAugmented::Constraint::Constraint ( const Constraint source,
NOX::CopyType  type = NOX::DeepCopy 
)

LOCA::Hopf::MinimallyAugmented::Constraint::~Constraint (  )  [virtual]

Destructor.

Definition at line 154 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.


Member Function Documentation

void LOCA::Hopf::MinimallyAugmented::Constraint::setGroup ( const Teuchos::RCP< LOCA::Hopf::MinimallyAugmented::AbstractGroup > &  g  )  [virtual]

Set the group pointer.

This method should be called when ever the constrained group is copied, since we don't explicitly copy the underlying group here.

Definition at line 160 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References grpPtr.

void LOCA::Hopf::MinimallyAugmented::Constraint::setFrequency ( double  freq  )  [virtual]

Set Hopf frequency.

Definition at line 167 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References isValidConstraints, isValidDX, and omega.

Teuchos::RCP< const NOX::Abstract::Vector > LOCA::Hopf::MinimallyAugmented::Constraint::getLeftNullVecReal (  )  const [virtual]

Returns real component of left null vector w.

Definition at line 176 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References w_vector.

Teuchos::RCP< const NOX::Abstract::Vector > LOCA::Hopf::MinimallyAugmented::Constraint::getLeftNullVecImag (  )  const [virtual]

Returns imaginary component of left null vector w.

Definition at line 183 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References w_vector.

Teuchos::RCP< const NOX::Abstract::Vector > LOCA::Hopf::MinimallyAugmented::Constraint::getRightNullVecReal (  )  const [virtual]

Returns real component of right null vector v.

Definition at line 190 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References v_vector.

Teuchos::RCP< const NOX::Abstract::Vector > LOCA::Hopf::MinimallyAugmented::Constraint::getRightNullVecImag (  )  const [virtual]

Returns imaginary component of right null vector v.

Definition at line 197 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References v_vector.

double LOCA::Hopf::MinimallyAugmented::Constraint::getSigmaReal (  )  const [virtual]

Returns real component of sigma.

Definition at line 204 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References constraints.

double LOCA::Hopf::MinimallyAugmented::Constraint::getSigmaImag (  )  const [virtual]

Returns imaginary component of sigma.

Definition at line 211 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References constraints.

void LOCA::Hopf::MinimallyAugmented::Constraint::copy ( const LOCA::MultiContinuation::ConstraintInterface source  )  [virtual]

Teuchos::RCP< LOCA::MultiContinuation::ConstraintInterface > LOCA::Hopf::MinimallyAugmented::Constraint::clone ( NOX::CopyType  type = NOX::DeepCopy  )  const [virtual]

Cloning function.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 258 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References Constraint().

int LOCA::Hopf::MinimallyAugmented::Constraint::numConstraints (  )  const [virtual]

Return number of constraints.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 265 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

void LOCA::Hopf::MinimallyAugmented::Constraint::setX ( const NOX::Abstract::Vector y  )  [virtual]

Set the solution vector to y.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 272 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References grpPtr, isValidConstraints, and isValidDX.

void LOCA::Hopf::MinimallyAugmented::Constraint::setParam ( int  paramID,
double  val 
) [virtual]

Sets parameter indexed by paramID.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 281 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References grpPtr, isValidConstraints, and isValidDX.

void LOCA::Hopf::MinimallyAugmented::Constraint::setParams ( const vector< int > &  paramIDs,
const NOX::Abstract::MultiVector::DenseMatrix vals 
) [virtual]

Sets parameters indexed by paramIDs.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 290 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References grpPtr, isValidConstraints, and isValidDX.

NOX::Abstract::Group::ReturnType LOCA::Hopf::MinimallyAugmented::Constraint::computeConstraints (  )  [virtual]

NOX::Abstract::Group::ReturnType LOCA::Hopf::MinimallyAugmented::Constraint::computeDX (  )  [virtual]

NOX::Abstract::Group::ReturnType LOCA::Hopf::MinimallyAugmented::Constraint::computeDP ( const vector< int > &  paramIDs,
NOX::Abstract::MultiVector::DenseMatrix dgdp,
bool  isValidG 
) [virtual]

Compute derivative of constraints w.r.t. supplied parameters.

The first column of dgdp should be filled with the constraint residuals $g$ if isValidG is false. If isValidG is true, then the dgdp contains $g$ on input.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 577 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References computeConstraints(), constraints, globalData, grpPtr, isValidConstraints, NOX::Abstract::Group::Ok, omega, sigma_scale, v_vector, and w_vector.

bool LOCA::Hopf::MinimallyAugmented::Constraint::isConstraints (  )  const [virtual]

Return true if constraint residuals are valid.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 622 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References isValidConstraints.

bool LOCA::Hopf::MinimallyAugmented::Constraint::isDX (  )  const [virtual]

Return true if derivatives of constraints w.r.t. x are valid.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 629 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References isValidDX.

const NOX::Abstract::MultiVector::DenseMatrix & LOCA::Hopf::MinimallyAugmented::Constraint::getConstraints (  )  const [virtual]

Return constraint residuals.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 636 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References constraints.

const NOX::Abstract::MultiVector * LOCA::Hopf::MinimallyAugmented::Constraint::getDX (  )  const [virtual]

Return solution component of constraint derivatives.

Implements LOCA::MultiContinuation::ConstraintInterfaceMVDX.

Definition at line 643 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References sigma_x.

bool LOCA::Hopf::MinimallyAugmented::Constraint::isDXZero (  )  const [virtual]

Return true if solution component of constraint derivatives is zero.

Implements LOCA::MultiContinuation::ConstraintInterface.

Definition at line 650 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

void LOCA::Hopf::MinimallyAugmented::Constraint::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 update the $a$ and $b$ vectors to $w$ and $v$ respectively if requested.

Reimplemented from LOCA::MultiContinuation::ConstraintInterface.

Definition at line 657 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References a_vector, b_vector, dn, globalData, NOX::Utils::StepperDetails, LOCA::Abstract::Iterator::Successful, updateVectorsEveryContinuationStep, v_vector, and w_vector.

NOX::Abstract::Group::ReturnType LOCA::Hopf::MinimallyAugmented::Constraint::computeDOmega ( NOX::Abstract::MultiVector::DenseMatrix domega  )  [virtual]

Compute derivative of sigma w.r.t. frequency omega.

Definition at line 679 of file LOCA_Hopf_MinimallyAugmented_Constraint.C.

References computeConstraints(), globalData, grpPtr, isValidConstraints, NOX::Abstract::Group::Ok, NOX::ShapeCopy, sigma_scale, v_vector, and w_vector.

Constraint& LOCA::Hopf::MinimallyAugmented::Constraint::operator= ( const Constraint source  )  [private]

Prohibit generation and use of operator=().


Member Data Documentation

Parsed top-level parameters.

Definition at line 286 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), Constraint(), and copy().

Teuchos::RCP<Teuchos::ParameterList> LOCA::Hopf::MinimallyAugmented::Constraint::hopfParams [protected]

Bifurcation parameter list.

Definition at line 289 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by Constraint(), and copy().

Pointer to base group that defines $F$.

Definition at line 292 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), computeDOmega(), computeDP(), computeDX(), setGroup(), setParam(), setParams(), and setX().

Stores C*v.

Definition at line 307 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), and copy().

Stores sigma_x.

Definition at line 310 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeDX(), copy(), and getDX().

Stores bordered solver strategy.

Definition at line 316 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), Constraint(), and copy().

Stores vector length as a double.

Definition at line 319 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), copy(), and postProcessContinuationStep().

Stores scale factor on sigma.

Definition at line 322 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), computeDOmega(), computeDP(), computeDX(), and copy().

Flag indicating whether Jacobian is symmetric.

Definition at line 327 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), Constraint(), and copy().

Flag indicating whether constraints are valid.

Definition at line 330 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), computeDOmega(), computeDP(), computeDX(), Constraint(), copy(), isConstraints(), setFrequency(), setParam(), setParams(), and setX().

Flag indicating whether sigma_x is valid.

Definition at line 333 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeDX(), Constraint(), copy(), isDX(), setFrequency(), setParam(), setParams(), and setX().

Stores the bifurcation parameter index.

Definition at line 336 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by copy().

Stores Hopf frequency.

Definition at line 339 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), computeDP(), computeDX(), copy(), and setFrequency().

Flag indicating whether to update $a$ and $b$ every continuation step.

Definition at line 345 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by Constraint(), copy(), and postProcessContinuationStep().

Flag indicating whether to update $a$ and $b$ every nonlinear iteration.

Definition at line 351 of file LOCA_Hopf_MinimallyAugmented_Constraint.H.

Referenced by computeConstraints(), Constraint(), and copy().


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

Generated on Thu Dec 17 11:02:51 2009 for Nonlinear Solver Project by  doxygen 1.5.9