Ifpack Package Browser (Single Doxygen Collection) Development

Ifpack_CrsRiluk: A class for constructing and using an incomplete lower/upper (ILU) factorization of a given Epetra_RowMatrix. More...
#include <Ifpack_CrsRiluk.h>
Public Member Functions  
Ifpack_CrsRiluk (const Ifpack_IlukGraph &Graph_in)  
Ifpack_CrsRiluk constuctor with variable number of indices per row.  
Ifpack_CrsRiluk (const Ifpack_CrsRiluk &Matrix)  
Copy constructor.  
virtual  ~Ifpack_CrsRiluk () 
Ifpack_CrsRiluk Destructor.  
int  InitValues (const Epetra_CrsMatrix &A) 
Initialize L and U with values from user matrix A.  
int  InitValues (const Epetra_VbrMatrix &A) 
Initialize L and U with values from user matrix A.  
bool  ValuesInitialized () const 
If values have been initialized, this query returns true, otherwise it returns false.  
void  SetRelaxValue (double RelaxValue) 
Set RILU(k) relaxation parameter.  
void  SetAbsoluteThreshold (double Athresh) 
Set absolute threshold value.  
void  SetRelativeThreshold (double Rthresh) 
Set relative threshold value.  
void  SetOverlapMode (Epetra_CombineMode OverlapMode) 
Set overlap mode type.  
int  SetParameters (const Teuchos::ParameterList ¶meterlist, bool cerr_warning_if_unused=false) 
Set parameters using a Teuchos::ParameterList object.  
int  Factor () 
Compute ILU factors L and U using the specified graph, diagonal perturbation thresholds and relaxation parameters.  
bool  Factored () const 
If factor is completed, this query returns true, otherwise it returns false.  
int  Solve (bool Trans, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const 
Returns the result of a Ifpack_CrsRiluk forward/back solve on a Epetra_MultiVector X in Y (works for Epetra_Vectors also).  
int  Multiply (bool Trans, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const 
Returns the result of multiplying U, D and L in that order on an Epetra_MultiVector X in Y.  
int  Condest (bool Trans, double &ConditionNumberEstimate) const 
Returns the maximum over all the condition number estimate for each local ILU set of factors.  
double  GetRelaxValue () 
Get RILU(k) relaxation parameter.  
double  GetAbsoluteThreshold () 
Get absolute threshold value.  
double  GetRelativeThreshold () 
Get relative threshold value.  
Epetra_CombineMode  GetOverlapMode () 
Get overlap mode type.  
int  NumGlobalRows () const 
Returns the number of global matrix rows.  
int  NumGlobalCols () const 
Returns the number of global matrix columns.  
int  NumGlobalNonzeros () const 
Returns the number of nonzero entries in the global graph.  
virtual int  NumGlobalBlockDiagonals () const 
Returns the number of diagonal entries found in the global input graph.  
int  NumMyRows () const 
Returns the number of local matrix rows.  
int  NumMyCols () const 
Returns the number of local matrix columns.  
int  NumMyNonzeros () const 
Returns the number of nonzero entries in the local graph.  
virtual int  NumMyBlockDiagonals () const 
Returns the number of diagonal entries found in the local input graph.  
virtual int  NumMyDiagonals () const 
Returns the number of nonzero diagonal values found in matrix.  
int  IndexBase () const 
Returns the index base for row and column indices for this graph.  
const Ifpack_IlukGraph &  Graph () const 
Returns the address of the Ifpack_IlukGraph associated with this factored matrix.  
const Epetra_CrsMatrix &  L () const 
Returns the address of the L factor associated with this factored matrix.  
const Epetra_Vector &  D () const 
Returns the address of the D factor associated with this factored matrix.  
const Epetra_CrsMatrix &  U () const 
Returns the address of the L factor associated with this factored matrix.  
Protected Member Functions  
void  SetFactored (bool Flag) 
void  SetValuesInitialized (bool Flag) 
bool  Allocated () const 
int  SetAllocated (bool Flag) 
int  BlockGraph2PointGraph (const Epetra_CrsGraph &BG, Epetra_CrsGraph &PG, bool Upper) 
Private Member Functions  
int  AllocateCrs () 
int  AllocateVbr () 
int  InitAllValues (const Epetra_RowMatrix &A, int MaxNumEntries) 
int  BlockMap2PointMap (const Epetra_BlockMap &BlockMap, Teuchos::RefCountPtr< Epetra_Map > *PointMap) 
int  GenerateXY (bool Trans, const Epetra_MultiVector &Xin, const Epetra_MultiVector &Yin, Teuchos::RefCountPtr< Epetra_MultiVector > *Xout, Teuchos::RefCountPtr< Epetra_MultiVector > *Yout) const 
Private Attributes  
bool  UserMatrixIsVbr_ 
bool  UserMatrixIsCrs_ 
bool  IsOverlapped_ 
const Ifpack_IlukGraph &  Graph_ 
Teuchos::RefCountPtr< Epetra_Map >  IlukRowMap_ 
Teuchos::RefCountPtr< Epetra_Map >  IlukDomainMap_ 
Teuchos::RefCountPtr< Epetra_Map >  IlukRangeMap_ 
Teuchos::RefCountPtr< const Epetra_Map >  U_DomainMap_ 
Teuchos::RefCountPtr< const Epetra_Map >  L_RangeMap_ 
const Epetra_Comm &  Comm_ 
Teuchos::RefCountPtr < Epetra_CrsMatrix >  L_ 
Teuchos::RefCountPtr < Epetra_CrsMatrix >  U_ 
Teuchos::RefCountPtr < Epetra_CrsGraph >  L_Graph_ 
Teuchos::RefCountPtr < Epetra_CrsGraph >  U_Graph_ 
Teuchos::RefCountPtr < Epetra_Vector >  D_ 
bool  UseTranspose_ 
int  NumMyDiagonals_ 
bool  Allocated_ 
bool  ValuesInitialized_ 
bool  Factored_ 
double  RelaxValue_ 
double  Athresh_ 
double  Rthresh_ 
double  Condest_ 
Teuchos::RefCountPtr < Epetra_MultiVector >  OverlapX_ 
Teuchos::RefCountPtr < Epetra_MultiVector >  OverlapY_ 
Teuchos::RefCountPtr < Epetra_MultiVector >  VbrX_ 
Teuchos::RefCountPtr < Epetra_MultiVector >  VbrY_ 
Epetra_CombineMode  OverlapMode_ 
Friends  
ostream &  operator<< (ostream &os, const Ifpack_CrsRiluk &A) 
<< operator will work for Ifpack_CrsRiluk.  
const char *  Label () const 
Returns a character string describing the operator.  
int  SetUseTranspose (bool UseTranspose_in) 
If set true, transpose of this operator will be applied.  
int  Apply (const Epetra_MultiVector &X, Epetra_MultiVector &Y) const 
Returns the result of a Epetra_Operator applied to a Epetra_MultiVector X in Y.  
int  ApplyInverse (const Epetra_MultiVector &X, Epetra_MultiVector &Y) const 
Returns the result of a Epetra_Operator inverse applied to an Epetra_MultiVector X in Y.  
double  NormInf () const 
Returns 0.0 because this class cannot compute Infnorm.  
bool  HasNormInf () const 
Returns false because this class cannot compute an Infnorm.  
bool  UseTranspose () const 
Returns the current UseTranspose setting.  
const Epetra_Map &  OperatorDomainMap () const 
Returns the Epetra_Map object associated with the domain of this operator.  
const Epetra_Map &  OperatorRangeMap () const 
Returns the Epetra_Map object associated with the range of this operator.  
const Epetra_Comm &  Comm () const 
Returns the Epetra_BlockMap object associated with the range of this matrix operator. 
Ifpack_CrsRiluk: A class for constructing and using an incomplete lower/upper (ILU) factorization of a given Epetra_RowMatrix.
The Ifpack_CrsRiluk class computes a "Relaxed" ILU factorization with level k fill of a given Epetra_CrsMatrix. The factorization that is produced is a function of several parameters:
The pattern of the matrix  All fill is derived from the original matrix nonzero structure. Level zero fill is defined as the original matrix pattern (nonzero structure), even if the matrix value at an entry is stored as a zero. (Thus it is possible to add entries to the ILU factors by adding zero entries the original matrix.)
Level of fill  Starting with the original matrix pattern as level fill of zero, the next level of fill is determined by analyzing the graph of the previous level and determining nonzero fill that is a result of combining entries that were from previous level only (not the current level). This rule limits fill to entries that are direct decendents from the previous level graph. Fill for level k is determined by applying this rule recursively. For sufficiently large values of k, the fill would eventually be complete and an exact LU factorization would be computed. Level of fill is defined during the construction of the Ifpack_IlukGraph object.
Level of overlap  All Ifpack preconditioners work on parallel distributed memory computers by using the row partitioning the user input matrix to determine the partitioning for local ILU factors. If the level of overlap is set to zero, the rows of the user matrix that are stored on a given processor are treated as a selfcontained local matrix and all column entries that reach to offprocessor entries are ignored. Setting the level of overlap to one tells Ifpack to increase the size of the local matrix by adding rows that are reached to by rows owned by this processor. Increasing levels of overlap are defined recursively in the same way. For sufficiently large levels of overlap, the entire matrix would be part of each processor's local ILU factorization process. Level of overlap is defined during the construction of the Ifpack_IlukGraph object.
Once the factorization is computed, applying the factorization \(LUy = x\) results in redundant approximations for any elements of y that correspond to rows that are part of more than one local ILU factor. The OverlapMode (changed by calling SetOverlapMode()) defines how these redundancies are handled using the Epetra_CombineMode enum. The default is to zero out all values of y for rows that were not part of the original matrix row distribution.
Fraction of relaxation  Ifpack_CrsRiluk computes the ILU factorization rowbyrow. As entries at a given row are computed, some number of them will be dropped because they do match the prescribed sparsity pattern. The relaxation factor determines how these dropped values will be handled. If the RelaxValue (changed by calling SetRelaxValue()) is zero, then these extra entries will by dropped. This is a classical ILU approach. If the RelaxValue is 1, then the sum of the extra entries will be added to the diagonal. This is a classical Modified ILU (MILU) approach. If RelaxValue is between 0 and 1, then RelaxValue times the sum of extra entries will be added to the diagonal.
For most situations, RelaxValue should be set to zero. For certain kinds of problems, e.g., reservoir modeling, there is a conservation principle involved such that any operator should obey a zero rowsum property. MILU was designed for these cases and you should set the RelaxValue to 1. For other situations, setting RelaxValue to some nonzero value may improve the stability of factorization, and can be used if the computed ILU factors are poorly conditioned.
Estimating Preconditioner Condition Numbers
For illconditioned matrices, we often have difficulty computing usable incomplete factorizations. The most common source of problems is that the factorization may encounter a small or zero pivot, in which case the factorization can fail, or even if the factorization succeeds, the factors may be so poorly conditioned that use of them in the iterative phase produces meaningless results. Before we can fix this problem, we must be able to detect it. To this end, we use a simple but effective condition number estimate for .
The condition of a matrix , called , is defined as in some appropriate norm . gives some indication of how many accurate floating point digits can be expected from operations involving the matrix and its inverse. A condition number approaching the accuracy of a given floating point number system, about 15 decimal digits in IEEE double precision, means that any results involving or may be meaningless.
The norm of a vector is defined as the maximum of the absolute values of the vector entries, and the norm of a matrix C is defined as . A crude lower bound for the is where . It is a lower bound because .
For our purposes, we want to estimate , where and are our incomplete factors. Edmond in his Ph.D. thesis demonstrates that provides an effective estimate for . Furthermore, since finding such that is a basic kernel for applying the preconditioner, computing this estimate of is performed by setting , calling the solve kernel to compute and then computing .
A priori Diagonal Perturbations
Given the above method to estimate the conditioning of the incomplete factors, if we detect that our factorization is too illconditioned we can improve the conditioning by perturbing the matrix diagonal and restarting the factorization using this more diagonally dominant matrix. In order to apply perturbation, prior to starting the factorization, we compute a diagonal perturbation of our matrix and perform the factorization on this perturbed matrix. The overhead cost of perturbing the diagonal is minimal since the first step in computing the incomplete factors is to copy the matrix into the memory space for the incomplete factors. We simply compute the perturbed diagonal at this point.
The actual perturbation values we use are the diagonal values with , , where is the matrix dimension and returns the sign of the diagonal entry. This has the effect of forcing the diagonal values to have minimal magnitude of and to increase each by an amount proportional to , and still keep the sign of the original diagonal entry.
Constructing Ifpack_CrsRiluk objects
Constructing Ifpack_CrsRiluk objects is a multistep process. The basic steps are as follows:
Note that, even after a matrix is constructed, it is possible to update existing matrix entries. It is not possible to create new entries.
Counting Floating Point Operations
Each Ifpack_CrsRiluk object keep track of the number of serial floating point operations performed using the specified object as the this argument to the function. The Flops() function returns this number as a double precision number. Using this information, in conjunction with the Epetra_Time class, one can get accurate parallel performance numbers. The ResetFlops() function resets the floating point counter.
Definition at line 196 of file Ifpack_CrsRiluk.h.
Ifpack_CrsRiluk::Ifpack_CrsRiluk  (  const Ifpack_IlukGraph &  Graph_in  ) 
Ifpack_CrsRiluk constuctor with variable number of indices per row.
Creates a Ifpack_CrsRiluk object and allocates storage.
In  Graph_in  Graph generated by Ifpack_IlukGraph. 
Definition at line 41 of file Ifpack_CrsRiluk.cpp.
Ifpack_CrsRiluk::Ifpack_CrsRiluk  (  const Ifpack_CrsRiluk &  Matrix  ) 
Copy constructor.
Definition at line 62 of file Ifpack_CrsRiluk.cpp.
Ifpack_CrsRiluk::~Ifpack_CrsRiluk  (  )  [virtual] 
Ifpack_CrsRiluk Destructor.
Definition at line 91 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::InitValues  (  const Epetra_CrsMatrix &  A  ) 
Initialize L and U with values from user matrix A.
Copies values from the user's matrix into the nonzero pattern of L and U.
In  A  User matrix to be factored. 
Definition at line 168 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::InitValues  (  const Epetra_VbrMatrix &  A  ) 
Initialize L and U with values from user matrix A.
Copies values from the user's matrix into the nonzero pattern of L and U.
In  A  User matrix to be factored. 
Definition at line 197 of file Ifpack_CrsRiluk.cpp.
bool Ifpack_CrsRiluk::ValuesInitialized  (  )  const [inline] 
If values have been initialized, this query returns true, otherwise it returns false.
Definition at line 236 of file Ifpack_CrsRiluk.h.
void Ifpack_CrsRiluk::SetRelaxValue  (  double  RelaxValue  )  [inline] 
Set RILU(k) relaxation parameter.
Definition at line 239 of file Ifpack_CrsRiluk.h.
void Ifpack_CrsRiluk::SetAbsoluteThreshold  (  double  Athresh  )  [inline] 
Set absolute threshold value.
Definition at line 242 of file Ifpack_CrsRiluk.h.
void Ifpack_CrsRiluk::SetRelativeThreshold  (  double  Rthresh  )  [inline] 
Set relative threshold value.
Definition at line 245 of file Ifpack_CrsRiluk.h.
void Ifpack_CrsRiluk::SetOverlapMode  (  Epetra_CombineMode  OverlapMode  )  [inline] 
Set overlap mode type.
Definition at line 248 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::SetParameters  (  const Teuchos::ParameterList &  parameterlist, 
bool  cerr_warning_if_unused = false 

) 
Set parameters using a Teuchos::ParameterList object.
Definition at line 148 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::Factor  (  void  ) 
Compute ILU factors L and U using the specified graph, diagonal perturbation thresholds and relaxation parameters.
This function computes the RILU(k) factors L and U using the current:
InitValues() must be called before the factorization can proceed.
Definition at line 345 of file Ifpack_CrsRiluk.cpp.
bool Ifpack_CrsRiluk::Factored  (  )  const [inline] 
If factor is completed, this query returns true, otherwise it returns false.
Definition at line 273 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::Solve  (  bool  Trans, 
const Epetra_MultiVector &  X,  
Epetra_MultiVector &  Y  
)  const 
Returns the result of a Ifpack_CrsRiluk forward/back solve on a Epetra_MultiVector X in Y (works for Epetra_Vectors also).
In  Trans If true, solve transpose problem. 
In  X  A Epetra_MultiVector of dimension NumVectors to solve for. 
Out  Y A Epetra_MultiVector of dimension NumVectorscontaining result. 
Definition at line 497 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::Multiply  (  bool  Trans, 
const Epetra_MultiVector &  X,  
Epetra_MultiVector &  Y  
)  const 
Returns the result of multiplying U, D and L in that order on an Epetra_MultiVector X in Y.
In  Trans If true, multiply by L^T, D and U^T in that order. 
In  X  A Epetra_MultiVector of dimension NumVectors to solve for. 
Out  Y A Epetra_MultiVector of dimension NumVectorscontaining result. 
Definition at line 538 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::Condest  (  bool  Trans, 
double &  ConditionNumberEstimate  
)  const 
Returns the maximum over all the condition number estimate for each local ILU set of factors.
This functions computes a local condition number estimate on each processor and return the maximum over all processor of the estimate.
In  Trans If true, solve transpose problem. 
Out  ConditionNumberEstimate  The maximum across all processors of the infinitynorm estimate of the condition number of the inverse of LDU. 
Definition at line 580 of file Ifpack_CrsRiluk.cpp.
double Ifpack_CrsRiluk::GetRelaxValue  (  )  [inline] 
Get RILU(k) relaxation parameter.
Definition at line 318 of file Ifpack_CrsRiluk.h.
double Ifpack_CrsRiluk::GetAbsoluteThreshold  (  )  [inline] 
Get absolute threshold value.
Definition at line 321 of file Ifpack_CrsRiluk.h.
double Ifpack_CrsRiluk::GetRelativeThreshold  (  )  [inline] 
Get relative threshold value.
Definition at line 324 of file Ifpack_CrsRiluk.h.
Epetra_CombineMode Ifpack_CrsRiluk::GetOverlapMode  (  )  [inline] 
Get overlap mode type.
Definition at line 327 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::NumGlobalRows  (  )  const [inline] 
Returns the number of global matrix rows.
Definition at line 331 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::NumGlobalCols  (  )  const [inline] 
Returns the number of global matrix columns.
Definition at line 334 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::NumGlobalNonzeros  (  )  const [inline] 
Returns the number of nonzero entries in the global graph.
Definition at line 337 of file Ifpack_CrsRiluk.h.
virtual int Ifpack_CrsRiluk::NumGlobalBlockDiagonals  (  )  const [inline, virtual] 
Returns the number of diagonal entries found in the global input graph.
Definition at line 340 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::NumMyRows  (  )  const [inline] 
Returns the number of local matrix rows.
Definition at line 343 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::NumMyCols  (  )  const [inline] 
Returns the number of local matrix columns.
Definition at line 346 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::NumMyNonzeros  (  )  const [inline] 
Returns the number of nonzero entries in the local graph.
Definition at line 349 of file Ifpack_CrsRiluk.h.
virtual int Ifpack_CrsRiluk::NumMyBlockDiagonals  (  )  const [inline, virtual] 
Returns the number of diagonal entries found in the local input graph.
Definition at line 352 of file Ifpack_CrsRiluk.h.
virtual int Ifpack_CrsRiluk::NumMyDiagonals  (  )  const [inline, virtual] 
Returns the number of nonzero diagonal values found in matrix.
Definition at line 355 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::IndexBase  (  )  const [inline] 
Returns the index base for row and column indices for this graph.
Definition at line 358 of file Ifpack_CrsRiluk.h.
const Ifpack_IlukGraph& Ifpack_CrsRiluk::Graph  (  )  const [inline] 
Returns the address of the Ifpack_IlukGraph associated with this factored matrix.
Definition at line 361 of file Ifpack_CrsRiluk.h.
const Epetra_CrsMatrix& Ifpack_CrsRiluk::L  (  )  const [inline] 
Returns the address of the L factor associated with this factored matrix.
Definition at line 364 of file Ifpack_CrsRiluk.h.
const Epetra_Vector& Ifpack_CrsRiluk::D  (  )  const [inline] 
Returns the address of the D factor associated with this factored matrix.
Definition at line 367 of file Ifpack_CrsRiluk.h.
const Epetra_CrsMatrix& Ifpack_CrsRiluk::U  (  )  const [inline] 
Returns the address of the L factor associated with this factored matrix.
Definition at line 370 of file Ifpack_CrsRiluk.h.
const char* Ifpack_CrsRiluk::Label  (  )  const [inline, virtual] 
Returns a character string describing the operator.
Reimplemented from Epetra_Object.
Definition at line 375 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::SetUseTranspose  (  bool  UseTranspose_in  )  [inline, virtual] 
If set true, transpose of this operator will be applied.
This flag allows the transpose of the given operator to be used implicitly. Setting this flag affects only the Apply() and ApplyInverse() methods. If the implementation of this interface does not support transpose use, this method should return a value of 1.
In  UseTranspose_in If true, multiply by the transpose of operator, otherwise just use operator. 
Implements Epetra_Operator.
Definition at line 387 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::Apply  (  const Epetra_MultiVector &  X, 
Epetra_MultiVector &  Y  
)  const [inline, virtual] 
Returns the result of a Epetra_Operator applied to a Epetra_MultiVector X in Y.
Note that this implementation of Apply does NOT perform a forward back solve with the LDU factorization. Instead it applies these operators via multiplication with U, D and L respectively. The ApplyInverse() method performs a solve.
In  X  A Epetra_MultiVector of dimension NumVectors to multiply with matrix. 
Out  Y A Epetra_MultiVector of dimension NumVectors containing result. 
Implements Epetra_Operator.
Definition at line 401 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::ApplyInverse  (  const Epetra_MultiVector &  X, 
Epetra_MultiVector &  Y  
)  const [inline, virtual] 
Returns the result of a Epetra_Operator inverse applied to an Epetra_MultiVector X in Y.
In this implementation, we use several existing attributes to determine how virtual method ApplyInverse() should call the concrete method Solve(). We pass in the UpperTriangular(), the Epetra_CrsMatrix::UseTranspose(), and NoDiagonal() methods. The most notable warning is that if a matrix has no diagonal values we assume that there is an implicit unit diagonal that should be accounted for when doing a triangular solve.
In  X  A Epetra_MultiVector of dimension NumVectors to solve for. 
Out  Y A Epetra_MultiVector of dimension NumVectors containing result. 
Implements Epetra_Operator.
Definition at line 418 of file Ifpack_CrsRiluk.h.
double Ifpack_CrsRiluk::NormInf  (  )  const [inline, virtual] 
Returns 0.0 because this class cannot compute Infnorm.
Implements Epetra_Operator.
Definition at line 422 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::HasNormInf  (  )  const [inline, virtual] 
Returns false because this class cannot compute an Infnorm.
Implements Epetra_Operator.
Definition at line 425 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::UseTranspose  (  )  const [inline, virtual] 
Returns the current UseTranspose setting.
Implements Epetra_Operator.
Definition at line 428 of file Ifpack_CrsRiluk.h.
const Epetra_Map& Ifpack_CrsRiluk::OperatorDomainMap  (  )  const [inline, virtual] 
Returns the Epetra_Map object associated with the domain of this operator.
Implements Epetra_Operator.
Definition at line 431 of file Ifpack_CrsRiluk.h.
const Epetra_Map& Ifpack_CrsRiluk::OperatorRangeMap  (  )  const [inline, virtual] 
Returns the Epetra_Map object associated with the range of this operator.
Implements Epetra_Operator.
Definition at line 434 of file Ifpack_CrsRiluk.h.
const Epetra_Comm& Ifpack_CrsRiluk::Comm  (  )  const [inline, virtual] 
Returns the Epetra_BlockMap object associated with the range of this matrix operator.
Implements Epetra_Operator.
Definition at line 437 of file Ifpack_CrsRiluk.h.
void Ifpack_CrsRiluk::SetFactored  (  bool  Flag  )  [inline, protected] 
Definition at line 441 of file Ifpack_CrsRiluk.h.
void Ifpack_CrsRiluk::SetValuesInitialized  (  bool  Flag  )  [inline, protected] 
Definition at line 442 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::Allocated  (  )  const [inline, protected] 
Definition at line 443 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::SetAllocated  (  bool  Flag  )  [inline, protected] 
Definition at line 444 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::BlockGraph2PointGraph  (  const Epetra_CrsGraph &  BG, 
Epetra_CrsGraph &  PG,  
bool  Upper  
)  [protected] 
Definition at line 598 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::AllocateCrs  (  )  [private] 
Definition at line 98 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::AllocateVbr  (  )  [private] 
Definition at line 110 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::InitAllValues  (  const Epetra_RowMatrix &  A, 
int  MaxNumEntries  
)  [private] 
Definition at line 235 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::BlockMap2PointMap  (  const Epetra_BlockMap &  BlockMap, 
Teuchos::RefCountPtr< Epetra_Map > *  PointMap  
)  [private] 
Definition at line 659 of file Ifpack_CrsRiluk.cpp.
int Ifpack_CrsRiluk::GenerateXY  (  bool  Trans, 
const Epetra_MultiVector &  Xin,  
const Epetra_MultiVector &  Yin,  
Teuchos::RefCountPtr< Epetra_MultiVector > *  Xout,  
Teuchos::RefCountPtr< Epetra_MultiVector > *  Yout  
)  const [private] 
Definition at line 686 of file Ifpack_CrsRiluk.cpp.
ostream& operator<<  (  ostream &  os, 
const Ifpack_CrsRiluk &  A  
)  [friend] 
<< operator will work for Ifpack_CrsRiluk.
Definition at line 747 of file Ifpack_CrsRiluk.cpp.
bool Ifpack_CrsRiluk::UserMatrixIsVbr_ [private] 
Definition at line 458 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::UserMatrixIsCrs_ [private] 
Definition at line 459 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::IsOverlapped_ [private] 
Definition at line 460 of file Ifpack_CrsRiluk.h.
const Ifpack_IlukGraph& Ifpack_CrsRiluk::Graph_ [private] 
Definition at line 461 of file Ifpack_CrsRiluk.h.
Definition at line 462 of file Ifpack_CrsRiluk.h.
Definition at line 463 of file Ifpack_CrsRiluk.h.
Definition at line 464 of file Ifpack_CrsRiluk.h.
Teuchos::RefCountPtr<const Epetra_Map> Ifpack_CrsRiluk::U_DomainMap_ [private] 
Definition at line 465 of file Ifpack_CrsRiluk.h.
Teuchos::RefCountPtr<const Epetra_Map> Ifpack_CrsRiluk::L_RangeMap_ [private] 
Definition at line 466 of file Ifpack_CrsRiluk.h.
const Epetra_Comm& Ifpack_CrsRiluk::Comm_ [private] 
Definition at line 467 of file Ifpack_CrsRiluk.h.
Definition at line 468 of file Ifpack_CrsRiluk.h.
Definition at line 469 of file Ifpack_CrsRiluk.h.
Definition at line 470 of file Ifpack_CrsRiluk.h.
Definition at line 471 of file Ifpack_CrsRiluk.h.
Teuchos::RefCountPtr<Epetra_Vector> Ifpack_CrsRiluk::D_ [private] 
Definition at line 472 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::UseTranspose_ [private] 
Definition at line 473 of file Ifpack_CrsRiluk.h.
int Ifpack_CrsRiluk::NumMyDiagonals_ [private] 
Definition at line 475 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::Allocated_ [private] 
Definition at line 476 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::ValuesInitialized_ [private] 
Definition at line 477 of file Ifpack_CrsRiluk.h.
bool Ifpack_CrsRiluk::Factored_ [private] 
Definition at line 478 of file Ifpack_CrsRiluk.h.
double Ifpack_CrsRiluk::RelaxValue_ [private] 
Definition at line 479 of file Ifpack_CrsRiluk.h.
double Ifpack_CrsRiluk::Athresh_ [private] 
Definition at line 480 of file Ifpack_CrsRiluk.h.
double Ifpack_CrsRiluk::Rthresh_ [private] 
Definition at line 481 of file Ifpack_CrsRiluk.h.
double Ifpack_CrsRiluk::Condest_ [mutable, private] 
Definition at line 482 of file Ifpack_CrsRiluk.h.
Teuchos::RefCountPtr<Epetra_MultiVector> Ifpack_CrsRiluk::OverlapX_ [mutable, private] 
Definition at line 484 of file Ifpack_CrsRiluk.h.
Teuchos::RefCountPtr<Epetra_MultiVector> Ifpack_CrsRiluk::OverlapY_ [mutable, private] 
Definition at line 485 of file Ifpack_CrsRiluk.h.
Teuchos::RefCountPtr<Epetra_MultiVector> Ifpack_CrsRiluk::VbrX_ [mutable, private] 
Definition at line 486 of file Ifpack_CrsRiluk.h.
Teuchos::RefCountPtr<Epetra_MultiVector> Ifpack_CrsRiluk::VbrY_ [mutable, private] 
Definition at line 487 of file Ifpack_CrsRiluk.h.
Definition at line 488 of file Ifpack_CrsRiluk.h.