NOX  Development
Public Member Functions | Protected Attributes | List of all members
LOCA::MultiContinuation::ArcLengthConstraint Class Reference

Implementation of LOCA::MultiContinuation::ConstraintInterfaceMVDX for arclength continuation. More...

#include <LOCA_MultiContinuation_ArcLengthConstraint.H>

Inheritance diagram for LOCA::MultiContinuation::ArcLengthConstraint:
Inheritance graph
[legend]
Collaboration diagram for LOCA::MultiContinuation::ArcLengthConstraint:
Collaboration graph
[legend]

Public Member Functions

 ArcLengthConstraint (const Teuchos::RCP< LOCA::GlobalData > &global_data, const Teuchos::RCP< LOCA::MultiContinuation::ArcLengthGroup > &grp)
 Constructor.
 
 ArcLengthConstraint (const ArcLengthConstraint &source, NOX::CopyType type=NOX::DeepCopy)
 Copy constructor.
 
 ~ArcLengthConstraint ()
 Destructor.
 
virtual void setArcLengthGroup (const Teuchos::RCP< LOCA::MultiContinuation::ArcLengthGroup > &grp)
 Set pointer to arclength group.
 
Implementation of LOCA::MultiContinuation::ConstraintInterface

virtual methods

virtual void copy (const ConstraintInterface &source)
 Copy.
 
virtual Teuchos::RCP< LOCA::MultiContinuation::ConstraintInterfaceclone (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 std::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 std::vector< int > &paramIDs, NOX::Abstract::MultiVector::DenseMatrix &dgdp, bool isValidG)
 Compute derivative of constraints w.r.t. supplied parameters. More...
 
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::DenseMatrixgetConstraints () const
 Return constraint residuals.
 
virtual const NOX::Abstract::MultiVectorgetDX () const
 Return solution component of constraint derivatives.
 
virtual bool isDXZero () const
 Return true if solution component of constraint derivatives is zero.
 
- Public Member Functions inherited from LOCA::MultiContinuation::ConstraintInterfaceMVDX
 ConstraintInterfaceMVDX ()
 Constructor.
 
virtual ~ConstraintInterfaceMVDX ()
 Destructor.
 
virtual NOX::Abstract::Group::ReturnType multiplyDX (double alpha, const NOX::Abstract::MultiVector &input_x, NOX::Abstract::MultiVector::DenseMatrix &result_p) const
 Compute result_p = alpha * dg/dx * input_x. More...
 
virtual NOX::Abstract::Group::ReturnType addDX (Teuchos::ETransp transb, double alpha, const NOX::Abstract::MultiVector::DenseMatrix &b, double beta, NOX::Abstract::MultiVector &result_x) const
 Compute result_x = alpha * dg/dx^T * op(b) + beta * result_x. More...
 
- Public Member Functions inherited from LOCA::MultiContinuation::ConstraintInterface
 ConstraintInterface ()
 Constructor.
 
virtual ~ConstraintInterface ()
 Destructor.
 
virtual void preProcessContinuationStep (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Perform any preprocessing before a continuation step starts. More...
 
virtual void postProcessContinuationStep (LOCA::Abstract::Iterator::StepStatus stepStatus)
 Perform any postprocessing after a continuation step finishes. More...
 

Protected Attributes

Teuchos::RCP< LOCA::GlobalDataglobalData
 LOCA global data object.
 
Teuchos::RCP< LOCA::MultiContinuation::ArcLengthGrouparcLengthGroup
 Pointer to arc-length group.
 
NOX::Abstract::MultiVector::DenseMatrix constraints
 Constraint values.
 
bool isValidConstraints
 Flag indicating whether constraints are valid.
 
std::vector< int > conParamIDs
 Continuation parameter IDs.
 

Detailed Description

Implementation of LOCA::MultiContinuation::ConstraintInterfaceMVDX for arclength continuation.

This class implements the arclength constraint equation for pseudo-arclength continuation:

\[ g(x,p,x_0,p_0,x^\ast,p^\ast,v,\Delta s)= (x-x^\ast)^Tv_x + (p-p^\ast) v_p - \Delta s \]

where $v_x$, $v_p$ are the solution and parameter components of the predictor direction $v$ respectively. Since the derivative of $g$ with respect to $x$ is just $v$, the predictor tangent, this class implements the MVDX version of the constraint interface.

Member Function Documentation

◆ computeDP()

NOX::Abstract::Group::ReturnType LOCA::MultiContinuation::ArcLengthConstraint::computeDP ( const std::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.

References LOCA::Extended::MultiVector::getScalar(), and NOX::Abstract::Group::Ok.


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