Ifpack_TriDiContainer: a class to define containers for dense matrices. More...

#include <Ifpack_TriDiContainer.h>

Inheritance diagram for Ifpack_TriDiContainer:

Public Member Functions
virtual double	InitializeFlops () const
	Returns the flops in Initialize(). More...

virtual double	ComputeFlops () const
	Returns the flops in Compute(). More...

virtual double	ApplyFlops () const
	Returns the flops in Apply(). More...

virtual double	ApplyInverseFlops () const
	Returns the flops in ApplyInverse(). More...

virtual std::ostream &	Print (std::ostream &os) const
	Prints basic information on iostream. This function is used by operator<<. More...

Public Member Functions inherited from Ifpack_Container
virtual	~Ifpack_Container ()
	Destructor. More...

Private Member Functions
virtual int	Extract (const Epetra_RowMatrix &Matrix_in)
	Extract the submatrices identified by the ID set int ID(). More...

Private Attributes
int	NumRows_
	Number of rows in the container. More...

int	NumVectors_
	Number of vectors in the container. More...

Ifpack_SerialTriDiMatrix	NonFactoredMatrix_
	TriDi matrix, that contains the non-factored matrix. More...

Ifpack_SerialTriDiMatrix	Matrix_
	TriDi matrix. More...

Epetra_SerialDenseMatrix	LHS_
	SerialDense vector representing the LHS. More...

Epetra_SerialDenseMatrix	RHS_
	SerialDense vector representing the RHS. More...

Ifpack_SerialTriDiSolver	Solver_
	TriDi solver (solution will be get using LAPACK). More...

Epetra_IntSerialDenseVector	ID_
	Sets of local rows. More...

bool	KeepNonFactoredMatrix_
	If `true`, keeps a copy of the non-factored matrix. More...

bool	IsInitialized_
	If `true`, the container has been successfully initialized. More...

bool	IsComputed_
	If `true`, the container has been successfully computed. More...

std::string	Label_
	Label for `this` object. More...

double	ComputeFlops_
	Flops in Compute(). More...

double	ApplyFlops_
	Flops in Apply(). More...

double	ApplyInverseFlops_
	Flops in ApplyInverse(). More...

	Ifpack_TriDiContainer (const int NumRows_in, const int NumVectors_in=1)
	Default constructor. More...

	Ifpack_TriDiContainer (const Ifpack_TriDiContainer &rhs)
	Copy constructor. More...

virtual	~Ifpack_TriDiContainer ()
	Destructor. More...

Ifpack_TriDiContainer &	operator= (const Ifpack_TriDiContainer &rhs)
	Operator=. More...

virtual int	NumRows () const
	Returns the number of rows of the matrix and LHS/RHS. More...

virtual int	NumVectors () const
	Returns the number of vectors in LHS/RHS. More...

virtual int	SetNumVectors (const int NumVectors_in)
	Sets the number of vectors for LHS/RHS. More...

virtual double &	LHS (const int i, const int Vector=0)
	Returns the i-th component of the vector Vector of LHS. More...

virtual double &	RHS (const int i, const int Vector=0)
	Returns the i-th component of the vector Vector of RHS. More...

virtual int &	ID (const int i)
	Returns the ID associated to local row i. More...

virtual int	SetMatrixElement (const int row, const int col, const double value)
	Set the matrix element (row,col) to `value`. More...

virtual int	SetParameters (Teuchos::ParameterList &)
	Sets all necessary parameters. More...

virtual bool	IsInitialized () const
	Returns `true` is the container has been successfully initialized. More...

virtual bool	IsComputed () const
	Returns `true` is the container has been successfully computed. More...

virtual const char *	Label () const
	Returns the label of this container. More...

virtual int	SetKeepNonFactoredMatrix (const bool flag)
	If `flag` is `true`, keeps a copy of the non-factored matrix. More...

virtual bool	KeepNonFactoredMatrix () const
	Returns KeepNonFactoredMatrix_. More...

virtual const Epetra_SerialDenseMatrix &	LHS () const
	Returns the dense vector containing the LHS. More...

virtual const Epetra_SerialDenseMatrix &	RHS () const
	Returns the dense vector containing the RHS. More...

virtual const Ifpack_SerialTriDiMatrix &	Matrix () const
	Returns the dense matrix or its factors. More...

virtual const Ifpack_SerialTriDiMatrix &	NonFactoredMatrix () const
	Returns the non-factored dense matrix (only if stored). More...

virtual const Epetra_IntSerialDenseVector &	ID () const
	Returns the integer dense vector of IDs. More...

virtual int	Initialize ()
	Initialize the container. More...

virtual int	Compute (const Epetra_RowMatrix &Matrix_in)
	Finalizes the linear system matrix and prepares for the application of the inverse. More...

virtual int	Apply ()
	Apply the matrix to RHS, results are stored in LHS. More...

virtual int	ApplyInverse ()
	Apply the inverse of the matrix to RHS, results are stored in LHS. More...

Detailed Description

Ifpack_TriDiContainer: a class to define containers for dense matrices.

To understand what an IFPACK container is, please refer to the documentation of the pure virtual class Ifpack_Container. Currently, containers are used by class Ifpack_BlockRelaxation.

Using block methods, one needs to store all diagonal blocks and to be also to apply the inverse of each diagonal block. Using class Ifpack_TriDiContainer, one can store the blocks as dense matrices, which can be advantageous when the blocks are small. Otherwise, class Ifpack_SparseContainer is probably more appropriate.

A typical use of a container is as follows:

#include "Ifpack_TriDiContainer.h"
...
// local matrix of (5,5), with two vectors for solution and rhs.
Ifpack_Container* Container = new
  Ifpack_TriDiContainer(5,5);
// assign local rows 1, 5, 12, 13, 16 to this container
Container(0) = 1;
Container(1) = 5;
Container(2) = 12;
Container(3) = 13;
Container(4) = 16;
// Now extract the submatrix corresponding to rows and columns:
// 1. initialize the container.
Container.Initialize();
// 2. extract matrix values from an Epetra_RowMatrix A,
// and compute LU factors of the submatrix identified by rows
// and columns 1, 5, 12, 13 and 16 using LAPACK
Container.Compute(A);
// We can set the RHS as follows:
Container.RHS(0) = 1.0;
Container.RHS(1) = 2.0;
Container.RHS(2) = 3.0;
Container.RHS(3) = 4.0;
Container.RHS(4) = 5.0;
// The linear system with the submatrix is solved as follows:
Container.ApplyInverse().

A call to Compute() computes the LU factorization of the linear system matrix, using LAPACK (more precisely, by calling the corresponding routines in Ifpack_SerialTriDiSolver). The default behavior is to store the matrix factors by overwriting the linear system matrix itself. This way, method Apply() fails, as the original matrix does no longer exists. An alternative is to call KeepNonFactoredMatrix(true), which forces Ifpack_TriDiContainer to maintain in memory a copy of the non-factored matrix.

Author: Marzio Sala, SNL 9214.

Date: Last update Nov-04.

Definition at line 124 of file Ifpack_TriDiContainer.h.

Constructor & Destructor Documentation

Ifpack_TriDiContainer::Ifpack_TriDiContainer	(	const int	NumRows_in,
		const int	NumVectors_in = `1`
	)

inline

Default constructor.

Definition at line 131 of file Ifpack_TriDiContainer.h.

Ifpack_TriDiContainer::Ifpack_TriDiContainer ( const Ifpack_TriDiContainer & rhs )

inline

Copy constructor.

Definition at line 143 of file Ifpack_TriDiContainer.h.

virtual Ifpack_TriDiContainer::~Ifpack_TriDiContainer ( )

inlinevirtual

Destructor.

Definition at line 159 of file Ifpack_TriDiContainer.h.

Member Function Documentation

Ifpack_TriDiContainer& Ifpack_TriDiContainer::operator= ( const Ifpack_TriDiContainer & rhs )

inline

Operator=.

Definition at line 166 of file Ifpack_TriDiContainer.h.

int Ifpack_TriDiContainer::NumRows ( ) const

virtual

Returns the number of rows of the matrix and LHS/RHS.

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation