Main Page | Namespace List | Class Hierarchy | Alphabetical List | Class List | File List | Namespace Members | Class Members | File Members | Related Pages | Search for 

LOCA::BorderedSolver::UpperTriangularBlockElimination Class Reference

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

#include <LOCA_BorderedSolver_UpperTriangularBlockElimination.H>

Collaboration diagram for LOCA::BorderedSolver::UpperTriangularBlockElimination:

[legend]
List of all members.

Public Member Functions

 UpperTriangularBlockElimination (const Teuchos::RefCountPtr< LOCA::GlobalData > &global_data)
 Constructor.
virtual ~UpperTriangularBlockElimination ()
 Destructor.
NOX::Abstract::Group::ReturnType solve (Teuchos::ParameterList &params, const NOX::Abstract::Group &grp, const NOX::Abstract::MultiVector *A, 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.
NOX::Abstract::Group::ReturnType solveTranspose (Teuchos::ParameterList &params, const LOCA::Abstract::TransposeSolveGroup &grp, const NOX::Abstract::MultiVector *A, 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.

Protected Attributes

Teuchos::RefCountPtr< LOCA::GlobalData > globalData
 Global data object.

Detailed Description

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

This class solves the extended system of equations

\[ \begin{bmatrix} op(J) & A \\ 0 & op(C) \end{bmatrix} \begin{bmatrix} X \\ Y \end{bmatrix} = \begin{bmatrix} F \\ G \end{bmatrix} \]

via block elimination:

\[ \begin{aligned} Y &= op(C)^{-1}G \\ X &= op(J)^{-1}(F - A Y) \end{aligned} \]

where $op$ represents either the identity operation or the transpose. $C$ must be nonzero, while $A$, $F$ or $G$ may be zero. 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. Note that for the transpose solve, the group representing $J$ must implement the LOCA::Abstract::TransposeSolveGroup interface.

Constructor & Destructor Documentation

LOCA::BorderedSolver::UpperTriangularBlockElimination::UpperTriangularBlockElimination (  const Teuchos::RefCountPtr< LOCA::GlobalData > &  global_data  ) 
 

Constructor.

Parameters:
global_data [in] Global data object

Member Function Documentation

NOX::Abstract::Group::ReturnType LOCA::BorderedSolver::UpperTriangularBlockElimination::solve (  Teuchos::ParameterList &  params,
const NOX::Abstract::Group &  grp,
const NOX::Abstract::MultiVector *  A,
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.

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

NOX::Abstract::Group::ReturnType LOCA::BorderedSolver::UpperTriangularBlockElimination::solveTranspose (  Teuchos::ParameterList &  params,
const LOCA::Abstract::TransposeSolveGroup &  grp,
const NOX::Abstract::MultiVector *  A,
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.

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

The documentation for this class was generated from the following files:
Generated on Thu Sep 18 12:38:30 2008 for NOX by doxygen 1.3.9.1