NOX Development
Loading...
Searching...
No Matches
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
LOCA::Epetra::LowRankUpdateRowMatrix Class Reference

An Epetra row matrix for implementing the operator $P = J + U V^T$. More...

#include <LOCA_Epetra_LowRankUpdateRowMatrix.H>

Inheritance diagram for LOCA::Epetra::LowRankUpdateRowMatrix:
Inheritance graph
[legend]
Collaboration diagram for LOCA::Epetra::LowRankUpdateRowMatrix:
Collaboration graph
[legend]

Public Member Functions

 LowRankUpdateRowMatrix (const Teuchos::RCP< LOCA::GlobalData > &global_data, const Teuchos::RCP< Epetra_RowMatrix > &jacRowMatrix, const Teuchos::RCP< Epetra_MultiVector > &U_multiVec, const Teuchos::RCP< Epetra_MultiVector > &V_multiVec, bool setup_for_solve, bool include_UV_terms)
 Constructor.
 
virtual ~LowRankUpdateRowMatrix ()
 Destructor.
 
virtual const Epetra_BlockMapMap () const
 Returns a reference to the Epetra_BlockMap for this object.
 
virtual int NumMyRowEntries (int MyRow, int &NumEntries) const
 Returns the number of nonzero entries in MyRow.
 
virtual int MaxNumEntries () const
 Returns the maximum of NumMyRowEntries() over all rows.
 
virtual int ExtractMyRowCopy (int MyRow, int Length, int &NumEntries, double *Values, int *Indices) const
 Returns a copy of the specified local row in user-provided arrays.
 
virtual int ExtractDiagonalCopy (Epetra_Vector &Diagonal) const
 Returns a copy of the main diagonal in a user-provided vector.
 
virtual int Multiply (bool TransA, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
 Returns the result of a Epetra_RowMatrix multiplied by a Epetra_MultiVector X in Y.
 
virtual int Solve (bool Upper, bool Trans, bool UnitDiagonal, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
 Returns result of a local-only solve using a triangular Epetra_RowMatrix with Epetra_MultiVectors X and Y.
 
virtual int InvRowSums (Epetra_Vector &x) const
 Computes the sum of absolute values of the rows of the Epetra_RowMatrix, results returned in x.
 
virtual int LeftScale (const Epetra_Vector &x)
 Scales the Epetra_RowMatrix on the left with a Epetra_Vector x.
 
virtual int InvColSums (Epetra_Vector &x) const
 Computes the sum of absolute values of the columns of the Epetra_RowMatrix, results returned in x.
 
virtual int RightScale (const Epetra_Vector &x)
 Scales the Epetra_RowMatrix on the right with a Epetra_Vector x.
 
virtual bool Filled () const
 If FillComplete() has been called, this query returns true, otherwise it returns false.
 
virtual double NormInf () const
 Returns the infinity norm of the global matrix.
 
virtual double NormOne () const
 Returns the one norm of the global matrix.
 
virtual int NumGlobalNonzeros () const
 Returns the number of nonzero entries in the global matrix.
 
virtual long long NumGlobalNonzeros64 () const
 
virtual int NumGlobalRows () const
 Returns the number of global matrix rows.
 
virtual long long NumGlobalRows64 () const
 
virtual int NumGlobalCols () const
 Returns the number of global matrix columns.
 
virtual long long NumGlobalCols64 () const
 
virtual int NumGlobalDiagonals () const
 Returns the number of global nonzero diagonal entries, based on global row/column index comparisons.
 
virtual long long NumGlobalDiagonals64 () const
 
virtual int NumMyNonzeros () const
 Returns the number of nonzero entries in the calling processor's portion of the matrix.
 
virtual int NumMyRows () const
 Returns the number of matrix rows owned by the calling processor.
 
virtual int NumMyCols () const
 Returns the number of matrix columns owned by the calling processor.
 
virtual int NumMyDiagonals () const
 Returns the number of local nonzero diagonal entries, based on global row/column index comparisons.
 
virtual bool LowerTriangular () const
 If matrix is lower triangular in local index space, this query returns true, otherwise it returns false.
 
virtual bool UpperTriangular () const
 If matrix is upper triangular in local index space, this query returns true, otherwise it returns false.
 
virtual const Epetra_MapRowMatrixRowMap () const
 Returns the Epetra_Map object associated with the rows of this matrix.
 
virtual const Epetra_MapRowMatrixColMap () const
 Returns the Epetra_Map object associated with the columns of this matrix.
 
virtual const Epetra_ImportRowMatrixImporter () const
 Returns the Epetra_Import object that contains the import operations for distributed operations.
 
- Public Member Functions inherited from LOCA::Epetra::LowRankUpdateOp
 LowRankUpdateOp (const Teuchos::RCP< LOCA::GlobalData > &global_data, const Teuchos::RCP< Epetra_Operator > &jacOperator, const Teuchos::RCP< const Epetra_MultiVector > &U_multiVec, const Teuchos::RCP< const Epetra_MultiVector > &V_multiVec, bool setup_for_solve)
 Constructor.
 
virtual ~LowRankUpdateOp ()
 Destructor.
 
virtual int SetUseTranspose (bool UseTranspose)
 Set to true if the transpose of the operator is requested.
 
virtual int Apply (const Epetra_MultiVector &Input, Epetra_MultiVector &Result) const
 Returns the result of a Epetra_Operator applied to a Epetra_MultiVector Input in Result as described above.
 
virtual int ApplyInverse (const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
 This method does nothing.
 
virtual double NormInf () const
 Returns an approximate infinity norm of the operator matrix.
 
virtual const char * Label () const
 Returns a character std::string describing the operator.
 
virtual bool UseTranspose () const
 Returns the current UseTranspose setting. Always returns false.
 
virtual bool HasNormInf () const
 Returns true if the this object can provide an approximate Inf-norm, false otherwise.
 
virtual const Epetra_CommComm () const
 Returns a reference to the Epetra_Comm communicator associated with this operator.
 
virtual const Epetra_MapOperatorDomainMap () const
 Returns the Epetra_Map object associated with the domain of this matrix operator.
 
virtual const Epetra_MapOperatorRangeMap () const
 Returns the Epetra_Map object associated with the range of this matrix operator.
 
virtual int SetUseTranspose (bool UseTranspose)=0
 
virtual int Apply (const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
 
virtual int ApplyInverse (const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
 
virtual double NormInf () const=0
 
virtual const char * Label () const=0
 
virtual bool UseTranspose () const=0
 
virtual bool HasNormInf () const=0
 
virtual const Epetra_CommComm () const=0
 
virtual const Epetra_MapOperatorDomainMap () const=0
 
virtual const Epetra_MapOperatorRangeMap () const=0
 
virtual int SetUseTranspose (bool UseTranspose)=0
 
virtual int Apply (const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
 
virtual int ApplyInverse (const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
 
virtual double NormInf () const=0
 
virtual const char * Label () const=0
 
virtual bool UseTranspose () const=0
 
virtual bool HasNormInf () const=0
 
virtual const Epetra_CommComm () const=0
 
virtual const Epetra_MapOperatorDomainMap () const=0
 
virtual const Epetra_MapOperatorRangeMap () const=0
 
virtual int NumMyRowEntries (int MyRow, int &NumEntries) const=0
 
virtual int MaxNumEntries () const=0
 
virtual int ExtractMyRowCopy (int MyRow, int Length, int &NumEntries, double *Values, int *Indices) const=0
 
virtual int ExtractDiagonalCopy (Epetra_Vector &Diagonal) const=0
 
virtual int Multiply (bool TransA, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
 
virtual int Solve (bool Upper, bool Trans, bool UnitDiagonal, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
 
virtual int InvRowSums (Epetra_Vector &x) const=0
 
virtual int LeftScale (const Epetra_Vector &x)=0
 
virtual int InvColSums (Epetra_Vector &x) const=0
 
virtual int RightScale (const Epetra_Vector &x)=0
 
virtual bool Filled () const=0
 
virtual double NormInf () const=0
 
virtual double NormOne () const=0
 
virtual int NumGlobalNonzeros () const=0
 
virtual int NumGlobalRows () const=0
 
virtual int NumGlobalCols () const=0
 
virtual int NumGlobalDiagonals () const=0
 
virtual int NumMyNonzeros () const=0
 
virtual int NumMyRows () const=0
 
virtual int NumMyCols () const=0
 
virtual int NumMyDiagonals () const=0
 
virtual bool LowerTriangular () const=0
 
virtual bool UpperTriangular () const=0
 
virtual const Epetra_MapRowMatrixRowMap () const=0
 
virtual const Epetra_MapRowMatrixColMap () const=0
 
virtual const Epetra_ImportRowMatrixImporter () const=0
 
virtual int NumMyRowEntries (int MyRow, int &NumEntries) const=0
 
virtual int MaxNumEntries () const=0
 
virtual int ExtractMyRowCopy (int MyRow, int Length, int &NumEntries, double *Values, int *Indices) const=0
 
virtual int ExtractDiagonalCopy (Epetra_Vector &Diagonal) const=0
 
virtual int Multiply (bool TransA, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
 
virtual int Solve (bool Upper, bool Trans, bool UnitDiagonal, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
 
virtual int InvRowSums (Epetra_Vector &x) const=0
 
virtual int LeftScale (const Epetra_Vector &x)=0
 
virtual int InvColSums (Epetra_Vector &x) const=0
 
virtual int RightScale (const Epetra_Vector &x)=0
 
virtual bool Filled () const=0
 
virtual double NormInf () const=0
 
virtual double NormOne () const=0
 
virtual int NumGlobalNonzeros () const=0
 
virtual int NumGlobalRows () const=0
 
virtual int NumGlobalCols () const=0
 
virtual int NumGlobalDiagonals () const=0
 
virtual int NumMyNonzeros () const=0
 
virtual int NumMyRows () const=0
 
virtual int NumMyCols () const=0
 
virtual int NumMyDiagonals () const=0
 
virtual bool LowerTriangular () const=0
 
virtual bool UpperTriangular () const=0
 
virtual const Epetra_MapRowMatrixRowMap () const=0
 
virtual const Epetra_MapRowMatrixColMap () const=0
 
virtual const Epetra_ImportRowMatrixImporter () const=0
 
virtual const Epetra_BlockMapMap () const=0
 
virtual const Epetra_BlockMapMap () const=0
 

Protected Member Functions

double computeUV (int MyRow, int MyCol) const
 Compute MyRow, MyCol entry of $U V^T$.
 

Protected Attributes

Teuchos::RCP< Epetra_RowMatrixJ_rowMatrix
 Stores row matrix representing J.
 
Teuchos::RCP< Epetra_MultiVectornonconst_U
 Stores pointer to non-const U.
 
Teuchos::RCP< Epetra_MultiVectornonconst_V
 Stores pointer to non-const V.
 
bool includeUV
 Flag indicating whether to include U*V^T terms.
 
int m
 Number of columns in U and V.
 
const Epetra_BlockMapU_map
 Map for U.
 
const Epetra_BlockMapV_map
 Map for V.
 
const Epetra_BlockMaprow_map
 Row map for J.
 
- Protected Attributes inherited from LOCA::Epetra::LowRankUpdateOp
Teuchos::RCP< LOCA::GlobalDataglobalData
 Global data object.
 
std::string label
 Label for operator.
 
Epetra_LocalMap localMap
 Local map for generating temporary matrices.
 
Teuchos::RCP< Epetra_OperatorJ
 Stores operator representing J.
 
Teuchos::RCP< const Epetra_MultiVectorU
 Stores multivector representing U.
 
Teuchos::RCP< const Epetra_MultiVectorV
 Stores multivector representing V.
 
bool useTranspose
 Flag indicating whether to use the transpose.
 
Teuchos::RCP< Epetra_MultiVectortmpMat
 Temporary matrix.
 
Teuchos::RCP< Epetra_MultiVectorJinvU
 Stores J^{-1}*U.
 
Teuchos::RCP< Epetra_MultiVectorlu
 Stores LU factorization of I + V^T*J^{-1}*U.
 
std::vector< int > ipiv
 Stores pivots for LU factorization.
 
Teuchos::LAPACK< int, double > lapack
 Lapack wrappers.
 

Detailed Description

An Epetra row matrix for implementing the operator $P = J + U V^T$.

This class implements the Epetra_RowMatrix interface for $P = J + U V^T$ where $J$ is an Epetra_RowMatrix and $U$ and $V$ are Epetra_MultiVectors. It is derived from LOCA::Epetra::LowRankUpdateOp to implement the Epetra_Operator interface. The interface here implements the Epetra_RowMatrix interface when the matrix $J$ is itself a row matrix. This allows preconditioners to be computed and scaling in linear systems to be performed when using this operator. The implementation here merely adds the corresponding entries for $U V^T$ to the rows of $J$. Note however this is only an approximation to the true matrix $J + U V^T$.

This class assumes $U$ and $V$ have the same distribution as the rows of $J$.

Constructor & Destructor Documentation

◆ LowRankUpdateRowMatrix()

LOCA::Epetra::LowRankUpdateRowMatrix::LowRankUpdateRowMatrix ( const Teuchos::RCP< LOCA::GlobalData > &  global_data,
const Teuchos::RCP< Epetra_RowMatrix > &  jacRowMatrix,
const Teuchos::RCP< Epetra_MultiVector > &  U_multiVec,
const Teuchos::RCP< Epetra_MultiVector > &  V_multiVec,
bool  setup_for_solve,
bool  include_UV_terms 
)

Constructor.

Parameters
global_data[in] The global data object
jacRowMatrix[in] Jacobian operator J as a row matrix
U_multiVec[in] Multivector representing U
V_multiVec[in] Multivector representing V
setup_for_solve[in] Setup data structures for ApplyInverse()
include_UV_terms[in] Include $U V^T$ terms in RowMatrix routines ExtractRowCopy(), ExtactDiagonalCopy(), InvRowSums(), InvColSums(), NormInf() and NormOne().

Member Function Documentation

◆ ExtractDiagonalCopy()

int LOCA::Epetra::LowRankUpdateRowMatrix::ExtractDiagonalCopy ( Epetra_Vector Diagonal) const
virtual

Returns a copy of the main diagonal in a user-provided vector.

Implements Epetra_RowMatrix.

◆ ExtractMyRowCopy()

int LOCA::Epetra::LowRankUpdateRowMatrix::ExtractMyRowCopy ( int  MyRow,
int  Length,
int &  NumEntries,
double *  Values,
int *  Indices 
) const
virtual

Returns a copy of the specified local row in user-provided arrays.

Implements Epetra_RowMatrix.

◆ Filled()

bool LOCA::Epetra::LowRankUpdateRowMatrix::Filled ( ) const
virtual

If FillComplete() has been called, this query returns true, otherwise it returns false.

Implements Epetra_RowMatrix.

◆ InvColSums()

int LOCA::Epetra::LowRankUpdateRowMatrix::InvColSums ( Epetra_Vector x) const
virtual

Computes the sum of absolute values of the columns of the Epetra_RowMatrix, results returned in x.

Implements Epetra_RowMatrix.

References Epetra_MultiVector::PutScalar(), and Epetra_MultiVector::Reciprocal().

◆ InvRowSums()

int LOCA::Epetra::LowRankUpdateRowMatrix::InvRowSums ( Epetra_Vector x) const
virtual

Computes the sum of absolute values of the rows of the Epetra_RowMatrix, results returned in x.

Implements Epetra_RowMatrix.

References Epetra_MultiVector::PutScalar(), and Epetra_MultiVector::Reciprocal().

◆ LeftScale()

int LOCA::Epetra::LowRankUpdateRowMatrix::LeftScale ( const Epetra_Vector x)
virtual

Scales the Epetra_RowMatrix on the left with a Epetra_Vector x.

Implements Epetra_RowMatrix.

◆ LowerTriangular()

bool LOCA::Epetra::LowRankUpdateRowMatrix::LowerTriangular ( ) const
virtual

If matrix is lower triangular in local index space, this query returns true, otherwise it returns false.

Implements Epetra_RowMatrix.

◆ Map()

const Epetra_BlockMap & LOCA::Epetra::LowRankUpdateRowMatrix::Map ( ) const
virtual

Returns a reference to the Epetra_BlockMap for this object.

Implements Epetra_SrcDistObject.

◆ MaxNumEntries()

int LOCA::Epetra::LowRankUpdateRowMatrix::MaxNumEntries ( ) const
virtual

Returns the maximum of NumMyRowEntries() over all rows.

Implements Epetra_RowMatrix.

◆ Multiply()

int LOCA::Epetra::LowRankUpdateRowMatrix::Multiply ( bool  TransA,
const Epetra_MultiVector X,
Epetra_MultiVector Y 
) const
virtual

Returns the result of a Epetra_RowMatrix multiplied by a Epetra_MultiVector X in Y.

Implements Epetra_RowMatrix.

References Epetra_MultiVector::Multiply(), and Epetra_MultiVector::NumVectors().

◆ NormInf()

double LOCA::Epetra::LowRankUpdateRowMatrix::NormInf ( ) const
virtual

Returns the infinity norm of the global matrix.

Implements Epetra_RowMatrix.

References Epetra_MultiVector::MaxValue(), and Epetra_MultiVector::Reciprocal().

◆ NormOne()

double LOCA::Epetra::LowRankUpdateRowMatrix::NormOne ( ) const
virtual

Returns the one norm of the global matrix.

Implements Epetra_RowMatrix.

References Epetra_MultiVector::MaxValue(), and Epetra_MultiVector::Reciprocal().

◆ NumGlobalCols()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumGlobalCols ( ) const
virtual

Returns the number of global matrix columns.

Implements Epetra_RowMatrix.

◆ NumGlobalDiagonals()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumGlobalDiagonals ( ) const
virtual

Returns the number of global nonzero diagonal entries, based on global row/column index comparisons.

Implements Epetra_RowMatrix.

◆ NumGlobalNonzeros()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumGlobalNonzeros ( ) const
virtual

Returns the number of nonzero entries in the global matrix.

Implements Epetra_RowMatrix.

◆ NumGlobalRows()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumGlobalRows ( ) const
virtual

Returns the number of global matrix rows.

Implements Epetra_RowMatrix.

◆ NumMyCols()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumMyCols ( ) const
virtual

Returns the number of matrix columns owned by the calling processor.

Implements Epetra_RowMatrix.

◆ NumMyDiagonals()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumMyDiagonals ( ) const
virtual

Returns the number of local nonzero diagonal entries, based on global row/column index comparisons.

Implements Epetra_RowMatrix.

◆ NumMyNonzeros()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumMyNonzeros ( ) const
virtual

Returns the number of nonzero entries in the calling processor's portion of the matrix.

Implements Epetra_RowMatrix.

◆ NumMyRowEntries()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumMyRowEntries ( int  MyRow,
int &  NumEntries 
) const
virtual

Returns the number of nonzero entries in MyRow.

Implements Epetra_RowMatrix.

◆ NumMyRows()

int LOCA::Epetra::LowRankUpdateRowMatrix::NumMyRows ( ) const
virtual

Returns the number of matrix rows owned by the calling processor.

Implements Epetra_RowMatrix.

◆ RightScale()

int LOCA::Epetra::LowRankUpdateRowMatrix::RightScale ( const Epetra_Vector x)
virtual

Scales the Epetra_RowMatrix on the right with a Epetra_Vector x.

Implements Epetra_RowMatrix.

◆ RowMatrixColMap()

const Epetra_Map & LOCA::Epetra::LowRankUpdateRowMatrix::RowMatrixColMap ( ) const
virtual

Returns the Epetra_Map object associated with the columns of this matrix.

Implements Epetra_RowMatrix.

◆ RowMatrixImporter()

const Epetra_Import * LOCA::Epetra::LowRankUpdateRowMatrix::RowMatrixImporter ( ) const
virtual

Returns the Epetra_Import object that contains the import operations for distributed operations.

Implements Epetra_RowMatrix.

◆ RowMatrixRowMap()

const Epetra_Map & LOCA::Epetra::LowRankUpdateRowMatrix::RowMatrixRowMap ( ) const
virtual

Returns the Epetra_Map object associated with the rows of this matrix.

Implements Epetra_RowMatrix.

◆ Solve()

int LOCA::Epetra::LowRankUpdateRowMatrix::Solve ( bool  Upper,
bool  Trans,
bool  UnitDiagonal,
const Epetra_MultiVector X,
Epetra_MultiVector Y 
) const
virtual

Returns result of a local-only solve using a triangular Epetra_RowMatrix with Epetra_MultiVectors X and Y.

Implements Epetra_RowMatrix.

◆ UpperTriangular()

bool LOCA::Epetra::LowRankUpdateRowMatrix::UpperTriangular ( ) const
virtual

If matrix is upper triangular in local index space, this query returns true, otherwise it returns false.

Implements Epetra_RowMatrix.


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