LOCA::BorderedSystem::EpetraAugmented Class Reference

Bordered system solver strategy based on augmenting the Jacobian operator. More...

#include <LOCA_BorderedSystem_EpetraAugmented.H>

Inheritance diagram for LOCA::BorderedSystem::EpetraAugmented:

[legend]Collaboration diagram for LOCA::BorderedSystem::EpetraAugmented:[legend]List of all members.


Public Member Functions
	EpetraAugmented (const Teuchos::RefCountPtr< LOCA::GlobalData > &global_data, const Teuchos::RefCountPtr< LOCA::Parameter::SublistParser > &topParams, const Teuchos::RefCountPtr< NOX::Parameter::List > &solverParams)
	Constructor.
virtual	~EpetraAugmented ()
	Destructor.
virtual void	setIsContiguous (bool flag)
	Set flag indicating whether F and A are stored in a continguous multivector.
virtual void	setMatrixBlocks (const Teuchos::RefCountPtr< const NOX::Abstract::Group > &group, const Teuchos::RefCountPtr< const NOX::Abstract::MultiVector > &blockA, const Teuchos::RefCountPtr< const LOCA::MultiContinuation::ConstraintInterface > &blockB, const Teuchos::RefCountPtr< const NOX::Abstract::MultiVector::DenseMatrix > &blockC)
	Set blocks.
virtual NOX::Abstract::Group::ReturnType	apply (const NOX::Abstract::MultiVector &X, const NOX::Abstract::MultiVector::DenseMatrix &Y, NOX::Abstract::MultiVector &U, NOX::Abstract::MultiVector::DenseMatrix &V) const
	Computed extended matrix-multivector product.
virtual NOX::Abstract::Group::ReturnType	applyTranspose (const NOX::Abstract::MultiVector &X, const NOX::Abstract::MultiVector::DenseMatrix &Y, NOX::Abstract::MultiVector &U, NOX::Abstract::MultiVector::DenseMatrix &V) const
	Computed extended matrix transpose-multivector product.
virtual NOX::Abstract::Group::ReturnType	applyInverse (NOX::Parameter::List &params, const NOX::Abstract::MultiVector F, const NOX::Abstract::MultiVector::DenseMatrix G, NOX::Abstract::MultiVector &X, NOX::Abstract::MultiVector::DenseMatrix &Y) const
	Solves the extended system using the technique described above.
Protected Member Functions
NOX::Abstract::Group::ReturnType	solveAZero (NOX::Parameter::List &params, const LOCA::MultiContinuation::ConstraintInterface BB, const NOX::Abstract::MultiVector::DenseMatrix CC, const NOX::Abstract::MultiVector F, const NOX::Abstract::MultiVector::DenseMatrix G, NOX::Abstract::MultiVector &X, NOX::Abstract::MultiVector::DenseMatrix &Y) const
	Solves the extended system of equations when the A block is zero using a simple block elimination.
NOX::Abstract::Group::ReturnType	solveBZero (NOX::Parameter::List &params, const NOX::Abstract::MultiVector AA, const NOX::Abstract::MultiVector::DenseMatrix CC, const NOX::Abstract::MultiVector F, const NOX::Abstract::MultiVector::DenseMatrix G, NOX::Abstract::MultiVector &X, NOX::Abstract::MultiVector::DenseMatrix &Y) const
	Solves the extended system of equations when the B block is zero using a simple block elimination.
Protected Attributes
Teuchos::RefCountPtr< LOCA::GlobalData >	globalData
	Global data object.
Teuchos::RefCountPtr< NOX::Parameter::List >	solverParams
	Solver parameters.
Teuchos::RefCountPtr< LOCA::EpetraNew::Group >	grp
	Pointer to group storing J.
Teuchos::RefCountPtr< const NOX::Abstract::MultiVector >	A
	Pointer to A block.
Teuchos::RefCountPtr< const NOX::Abstract::MultiVector >	B
	Pointer to B block.
Teuchos::RefCountPtr< const NOX::Abstract::MultiVector::DenseMatrix >	C
	Pointer to C block.
Teuchos::RefCountPtr< const LOCA::MultiContinuation::ConstraintInterfaceMVDX >	constraints
	Pointer to constraint interface.
int	numConstraints
	Number of constraint equations.
bool	isZeroA
	flag indicating whether A block is zero
bool	isZeroB
	flag indicating whether B block is zero
bool	isZeroC
	flag indicating whether C block is zero
bool	isContiguous
	flag indicating whether F and A are contiguous

Detailed Description

Bordered system solver strategy based on augmenting the Jacobian operator.

This class solves the extended system of equations

$\begin{bmatrix} J & A \\ B^T & C \end{bmatrix} \begin{bmatrix} X \\ Y \end{bmatrix} = \begin{bmatrix} F \\ G \end{bmatrix}$

by forming an augmented Epetra_Operator representing

$\begin{bmatrix} J & A \\ B^T & C \end{bmatrix}$

by creating a new Epetra_Map for the additional equations.

Constructor & Destructor Documentation

LOCA::BorderedSystem::EpetraAugmented::EpetraAugmented ( const Teuchos::RefCountPtr< LOCA::GlobalData > & global_data,

const Teuchos::RefCountPtr< LOCA::Parameter::SublistParser > & topParams,

const Teuchos::RefCountPtr< NOX::Parameter::List > & solverParams

)

Constructor.

Parameters:

global_data [in] Global data object

topParams [in] Parsed top-level parameter list

solverParams [in] Bordered solver parameters. Currently none are referenced.

Member Function Documentation

NOX::Abstract::Group::ReturnType LOCA::BorderedSystem::EpetraAugmented::apply ( const NOX::Abstract::MultiVector & X,

const NOX::Abstract::MultiVector::DenseMatrix & Y,

NOX::Abstract::MultiVector & U,

NOX::Abstract::MultiVector::DenseMatrix & V

) const [virtual]

Computed extended matrix-multivector product.
Computes
$\begin{bmatrix} U \\ V \end{bmatrix} = \begin{bmatrix} J & A \\ B^T & C \end{bmatrix} \begin{bmatrix} X \\ Y \end{bmatrix} = \begin{bmatrix} J*X + A*Y \\ B^T*X + C*Y \end{bmatrix}.$

Implements LOCA::BorderedSystem::AbstractStrategy.

NOX::Abstract::Group::ReturnType LOCA::BorderedSystem::EpetraAugmented::applyInverse ( NOX::Parameter::List & params,

const NOX::Abstract::MultiVector * F,

const NOX::Abstract::MultiVector::DenseMatrix * G,

NOX::Abstract::MultiVector & X,

NOX::Abstract::MultiVector::DenseMatrix & Y

) const [virtual]

Solves the extended system using the technique described above.
The params argument is the linear solver parameters. If isZeroF or isZeroG is true, than the corresponding F or G pointers may be NULL.
Note that if either the A or B blocks are zero, the system is solved using a simple block elimination scheme instead of the Householder scheme.
Implements LOCA::BorderedSystem::AbstractStrategy.

NOX::Abstract::Group::ReturnType LOCA::BorderedSystem::EpetraAugmented::applyTranspose ( const NOX::Abstract::MultiVector & X,

const NOX::Abstract::MultiVector::DenseMatrix & Y,

NOX::Abstract::MultiVector & U,

NOX::Abstract::MultiVector::DenseMatrix & V

) const [virtual]

Computed extended matrix transpose-multivector product.
Computes
$\begin{bmatrix} U \\ V \end{bmatrix} = \begin{bmatrix} J^T & B \\ A^T & C \end{bmatrix} \begin{bmatrix} X \\ Y \end{bmatrix} = \begin{bmatrix} J^T*X + B*Y \\ A^T*X + C^T*Y \end{bmatrix}.$

Implements LOCA::BorderedSystem::AbstractStrategy.

void LOCA::BorderedSystem::EpetraAugmented::setMatrixBlocks ( const Teuchos::RefCountPtr< const NOX::Abstract::Group > & group,

const Teuchos::RefCountPtr< const NOX::Abstract::MultiVector > & blockA,

const Teuchos::RefCountPtr< const LOCA::MultiContinuation::ConstraintInterface > & blockB,

const Teuchos::RefCountPtr< const NOX::Abstract::MultiVector::DenseMatrix > & blockC

) [virtual]

Set blocks.
The blockA or blockC pointer may be null if either is zero. Whether block B is zero will be determined by querying blockB via ConstraintInterface::isConstraintDerivativesXZero.
Implements LOCA::BorderedSystem::AbstractStrategy.

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

LOCA_BorderedSystem_EpetraAugmented.H
LOCA_BorderedSystem_EpetraAugmented.C

Generated on Thu Sep 18 12:40:43 2008 for NOX by

1.3.9.1

LOCA::BorderedSystem::EpetraAugmented Class Reference

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation