Block elimination strategy for solving a block lower-triangular system. More...

#include <LOCA_BorderedSolver_LowerTriangularBlockElimination.H>

Collaboration diagram for LOCA::BorderedSolver::LowerTriangularBlockElimination:

Public Member Functions
	LowerTriangularBlockElimination (const Teuchos::RCP< LOCA::GlobalData > &global_data)
	Constructor. More...

virtual	~LowerTriangularBlockElimination ()
	Destructor.

NOX::Abstract::Group::ReturnType	solve (Teuchos::ParameterList &params, const LOCA::BorderedSolver::AbstractOperator &op, const LOCA::MultiContinuation::ConstraintInterface &B, const NOX::Abstract::MultiVector::DenseMatrix &C, 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 as described above with B specified as a LOCA::MultiContinuation::ConstraintInterface object. More...

NOX::Abstract::Group::ReturnType	solve (Teuchos::ParameterList &params, const LOCA::BorderedSolver::AbstractOperator &op, const NOX::Abstract::MultiVector &B, const NOX::Abstract::MultiVector::DenseMatrix &C, 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 as described above with B specified as a NOX::Abstract::MultiVector. More...

NOX::Abstract::Group::ReturnType	solveTranspose (Teuchos::ParameterList &params, const LOCA::BorderedSolver::AbstractOperator &op, const LOCA::MultiContinuation::ConstraintInterface &B, const NOX::Abstract::MultiVector::DenseMatrix &C, 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 tranpose of J and C as described above with B specified as a LOCA::MultiContinuation::ConstraintInterface object. More...

NOX::Abstract::Group::ReturnType	solveTranspose (Teuchos::ParameterList &params, const LOCA::BorderedSolver::AbstractOperator &op, const NOX::Abstract::MultiVector &B, const NOX::Abstract::MultiVector::DenseMatrix &C, 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 tranpose of J and C as described above with B specified as a NOX::Abstract::MultiVector object. More...

Protected Attributes
Teuchos::RCP< LOCA::GlobalData >	globalData
	Global data object.

Detailed Description

Block elimination strategy for solving a block lower-triangular system.

This class solves the extended system of equations

$\begin{bmatrix} op(J) & 0 \\ B^T & op(C) \end{bmatrix} \begin{bmatrix} X \\ Y \end{bmatrix} = \begin{bmatrix} F \\ G \end{bmatrix}$

via block elimination:

$\begin{aligned} X &= op(J)^{-1} F \\ Y &= op(C)^{-1}(G-B^T X) \end{aligned}$

where $op$ represents either the identity operation or the transpose. $C$ must be nonzero, while $B$ , $F$ or $G$ may be zero. $B$ may be specified either as a NOX::Abstract::MultiVector or a LOCA::MultiContinuation::ConstraintInterface object. The solve for the non-transposed system is implemented by the solve() method, while the solve for the transposed system is implemented by the solveTranspose() method.

Constructor & Destructor Documentation

LOCA::BorderedSolver::LowerTriangularBlockElimination::LowerTriangularBlockElimination ( const Teuchos::RCP< LOCA::GlobalData > & global_data )

Constructor.

Parameters

global_data [in] Global data object

Member Function Documentation

NOX::Abstract::Group::ReturnType LOCA::BorderedSolver::LowerTriangularBlockElimination::solve	(	Teuchos::ParameterList &	params,
		const LOCA::BorderedSolver::AbstractOperator &	op,
		const LOCA::MultiContinuation::ConstraintInterface &	B,
		const NOX::Abstract::MultiVector::DenseMatrix &	C,
		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 as described above with B specified as a LOCA::MultiContinuation::ConstraintInterface object.

Either F, or G may be zero by passing NULL.

References LOCA::BorderedSolver::AbstractOperator::applyInverse(), Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::assign(), Teuchos::Copy, NOX::Abstract::Group::Failed, Teuchos::LAPACK< OrdinalType, ScalarType >::GETRF(), Teuchos::LAPACK< OrdinalType, ScalarType >::GETRS(), NOX::Abstract::MultiVector::init(), LOCA::MultiContinuation::ConstraintInterface::isDXZero(), LOCA::MultiContinuation::ConstraintInterface::multiplyDX(), Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::numCols(), Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::numRows(), NOX::Abstract::Group::Ok, Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::putScalar(), Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::stride(), and Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::values().

Referenced by LOCA::BorderedSolver::Bordering::applyInverse(), LOCA::BorderedSolver::EpetraAugmented::applyInverse(), LOCA::BorderedSolver::TpetraHouseholder::applyInverse(), and LOCA::BorderedSolver::EpetraHouseholder::applyInverse().

NOX::Abstract::Group::ReturnType LOCA::BorderedSolver::LowerTriangularBlockElimination::solve	(	Teuchos::ParameterList &	params,
		const LOCA::BorderedSolver::AbstractOperator &	op,
		const NOX::Abstract::MultiVector &	B,
		const NOX::Abstract::MultiVector::DenseMatrix &	C,
		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 as described above with B specified as a NOX::Abstract::MultiVector.

Either F, or G may be zero by passing NULL.

References Teuchos::rcp().

NOX::Abstract::Group::ReturnType LOCA::BorderedSolver::LowerTriangularBlockElimination::solveTranspose	(	Teuchos::ParameterList &	params,
		const LOCA::BorderedSolver::AbstractOperator &	op,
		const LOCA::MultiContinuation::ConstraintInterface &	B,
		const NOX::Abstract::MultiVector::DenseMatrix &	C,
		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 tranpose of J and C as described above with B specified as a LOCA::MultiContinuation::ConstraintInterface object.

Either F, or G may be zero by passing NULL.

References LOCA::BorderedSolver::AbstractOperator::applyInverseTranspose(), Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::assign(), Teuchos::Copy, NOX::Abstract::Group::Failed, Teuchos::LAPACK< OrdinalType, ScalarType >::GETRF(), Teuchos::LAPACK< OrdinalType, ScalarType >::GETRS(), NOX::Abstract::MultiVector::init(), LOCA::MultiContinuation::ConstraintInterface::isDXZero(), LOCA::MultiContinuation::ConstraintInterface::multiplyDX(), Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::numCols(), Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::numRows(), NOX::Abstract::Group::Ok, Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::putScalar(), Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::stride(), and Teuchos::SerialDenseMatrix< OrdinalType, ScalarType >::values().

Referenced by LOCA::BorderedSolver::Bordering::applyInverseTranspose(), LOCA::BorderedSolver::TpetraHouseholder::applyInverseTranspose(), and LOCA::BorderedSolver::EpetraHouseholder::applyInverseTranspose().

NOX::Abstract::Group::ReturnType LOCA::BorderedSolver::LowerTriangularBlockElimination::solveTranspose	(	Teuchos::ParameterList &	params,
		const LOCA::BorderedSolver::AbstractOperator &	op,
		const NOX::Abstract::MultiVector &	B,
		const NOX::Abstract::MultiVector::DenseMatrix &	C,
		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 tranpose of J and C as described above with B specified as a NOX::Abstract::MultiVector object.

Either F, or G may be zero by passing NULL.

References Teuchos::rcp().

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

LOCA_BorderedSolver_LowerTriangularBlockElimination.H
LOCA_BorderedSolver_LowerTriangularBlockElimination.C

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation