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

#include <LOCA_MultiContinuation_ArcLengthConstraint.H>

Inheritance diagram for LOCA::MultiContinuation::ArcLengthConstraint:

[legend]

Collaboration diagram for LOCA::MultiContinuation::ArcLengthConstraint:

[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.

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.

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.

Public Member Functions inherited from LOCA::MultiContinuation::ConstraintInterface
	ConstraintInterface ()
	Constructor.

virtual	~ConstraintInterface ()
	Destructor.

virtual void	preProcessContinuationStep (LOCA::Abstract::Iterator::StepStatus)
	Perform any preprocessing before a continuation step starts.

virtual void	postProcessContinuationStep (LOCA::Abstract::Iterator::StepStatus)
	Perform any postprocessing after a continuation step finishes.

Implementation of LOCA::MultiContinuation::ConstraintInterface
virtual methods
Teuchos::RCP< LOCA::GlobalData >	globalData
	LOCA global data object.

Teuchos::RCP< LOCA::MultiContinuation::ArcLengthGroup >	arcLengthGroup
	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.

virtual void	copy (const 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 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.

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.

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 , are the solution and parameter components of the predictor direction respectively. Since the derivative of with respect to is just , the predictor tangent, this class implements the MVDX version of the constraint interface.

Member Function Documentation

◆ clone()

Teuchos::RCP< LOCA::MultiContinuation::ConstraintInterface > LOCA::MultiContinuation::ArcLengthConstraint::clone ( NOX::CopyType type = NOX::DeepCopy ) const

virtual

Cloning function.

Implements LOCA::MultiContinuation::ConstraintInterface.

◆ computeConstraints()

NOX::Abstract::Group::ReturnType LOCA::MultiContinuation::ArcLengthConstraint::computeConstraints ( )

virtual

Compute continuation constraint equations.

Implements LOCA::MultiContinuation::ConstraintInterface.

References LOCA::MultiContinuation::ExtendedMultiVector::clone(), and NOX::Abstract::Group::Ok.

◆ 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 if isValidG is false. If isValidG is true, then the dgdp contains on input.

Implements LOCA::MultiContinuation::ConstraintInterface.

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

◆ computeDX()

NOX::Abstract::Group::ReturnType LOCA::MultiContinuation::ArcLengthConstraint::computeDX ( )

virtual

Compute derivative of constraints w.r.t. solution vector x.

Implements LOCA::MultiContinuation::ConstraintInterface.

References NOX::Abstract::Group::Ok.

◆ copy()

void LOCA::MultiContinuation::ArcLengthConstraint::copy ( const ConstraintInterface & source )

virtual

Copy.

Implements LOCA::MultiContinuation::ConstraintInterface.

References conParamIDs, constraints, globalData, and isValidConstraints.

◆ getConstraints()

const NOX::Abstract::MultiVector::DenseMatrix & LOCA::MultiContinuation::ArcLengthConstraint::getConstraints ( ) const

virtual

Return constraint residuals.

Implements LOCA::MultiContinuation::ConstraintInterface.

◆ getDX()

const NOX::Abstract::MultiVector * LOCA::MultiContinuation::ArcLengthConstraint::getDX ( ) const

virtual

Return solution component of constraint derivatives.

Implements LOCA::MultiContinuation::ConstraintInterfaceMVDX.

References LOCA::MultiContinuation::ExtendedMultiVector::getXMultiVec().

◆ isConstraints()

bool LOCA::MultiContinuation::ArcLengthConstraint::isConstraints ( ) const

virtual

Return true if constraint residuals are valid.

Implements LOCA::MultiContinuation::ConstraintInterface.

◆ isDX()

bool LOCA::MultiContinuation::ArcLengthConstraint::isDX ( ) const

virtual

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

Implements LOCA::MultiContinuation::ConstraintInterface.

◆ isDXZero()

bool LOCA::MultiContinuation::ArcLengthConstraint::isDXZero ( ) const

virtual

Return true if solution component of constraint derivatives is zero.

Implements LOCA::MultiContinuation::ConstraintInterface.

◆ numConstraints()

int LOCA::MultiContinuation::ArcLengthConstraint::numConstraints ( ) const

virtual

Return number of constraints.

Implements LOCA::MultiContinuation::ConstraintInterface.

◆ setParam()

void LOCA::MultiContinuation::ArcLengthConstraint::setParam	(	int	paramID,
		double	val
	)

virtual

Sets parameter indexed by paramID.

Implements LOCA::MultiContinuation::ConstraintInterface.

◆ setParams()

void LOCA::MultiContinuation::ArcLengthConstraint::setParams	(	const std::vector< int > &	paramIDs,
		const NOX::Abstract::MultiVector::DenseMatrix &	vals
	)

virtual

Sets parameters indexed by paramIDs.

Implements LOCA::MultiContinuation::ConstraintInterface.

◆ setX()

void LOCA::MultiContinuation::ArcLengthConstraint::setX ( const NOX::Abstract::Vector & y )

virtual

Set the solution vector to y.

Implements LOCA::MultiContinuation::ConstraintInterface.

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

LOCA_MultiContinuation_ArcLengthConstraint.H
LOCA_MultiContinuation_ArcLengthConstraint.C

Public Member Functions

Implementation of LOCA::MultiContinuation::ConstraintInterface

Detailed Description

Member Function Documentation

◆ clone()

◆ computeConstraints()

◆ computeDP()

◆ computeDX()

◆ copy()

◆ getConstraints()

◆ getDX()

◆ isConstraints()

◆ isDX()

◆ isDXZero()

◆ numConstraints()

◆ setParam()

◆ setParams()

◆ setX()