51 #include "Epetra_CompObject.h" 52 #include "Epetra_MultiVector.h" 53 #include "Epetra_Vector.h" 54 #include "Epetra_CrsGraph.h" 55 #include "Epetra_CrsMatrix.h" 56 #include "Epetra_BlockMap.h" 57 #include "Epetra_Map.h" 58 #include "Epetra_Object.h" 59 #include "Epetra_Comm.h" 60 #include "Epetra_RowMatrix.h" 61 #include "Epetra_Time.h" 62 #include "Teuchos_RefCountPtr.hpp" 177 const int MaxIters = 1550,
178 const double Tol = 1e-9,
234 virtual std::ostream&
Print(std::ostream& os)
const;
336 #ifndef EPETRA_NO_32BIT_GLOBAL_INDICES 374 #ifndef EPETRA_NO_32BIT_GLOBAL_INDICES 392 Teuchos::RefCountPtr<Epetra_RowMatrix>
A_;
393 Teuchos::RefCountPtr<Ifpack_IlukGraph>
Graph_;
402 Teuchos::RefCountPtr<Epetra_CrsMatrix>
L_;
404 Teuchos::RefCountPtr<Epetra_CrsMatrix>
U_;
408 Teuchos::RefCountPtr<Epetra_Vector>
D_;
char Label_[160]
Label of this object.
Teuchos::RefCountPtr< Epetra_CrsGraph > CrsGraph_
double InitializeTime_
Contains the time for all successful calls to Initialize().
Teuchos::RefCountPtr< Epetra_CrsGraph > L_Graph_
const Epetra_Comm & Comm_
const Epetra_Comm & Comm() const
Returns the Epetra_BlockMap object associated with the range of this matrix operator.
double Athresh_
absolute threshold
Teuchos::RefCountPtr< Epetra_Vector > D_
Diagonal of factors.
virtual int NumInitialize() const
Returns the number of calls to Initialize().
int SetUseTranspose(bool UseTranspose_in)
If set true, transpose of this operator will be applied.
int NumMyNonzeros() const
Returns the number of nonzero entries in the local graph.
bool IsComputed() const
If factor is completed, this query returns true, otherwise it returns false.
Epetra_RowMatrix & Matrix()
Returns a reference to the matrix.
Ifpack_IlukGraph: A class for constructing level filled graphs for use with ILU(k) class precondition...
double RelativeThreshold() const
Get relative threshold value.
const Epetra_CrsMatrix & L() const
Returns the address of the L factor associated with this factored matrix.
Ifpack_ILU(const Ifpack_ILU &RHS)
Copy constructor (should never be used)
const Epetra_CrsMatrix & U() const
Returns the address of the L factor associated with this factored matrix.
virtual double InitializeTime() const
Returns the time spent in Initialize().
Teuchos::RefCountPtr< Ifpack_IlukGraph > Graph_
int Multiply(bool Trans, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
virtual int NumMyDiagonals() const
Returns the number of nonzero diagonal values found in matrix.
const Epetra_Map * U_DomainMap_
bool IsComputed_
If true, the preconditioner has been successfully computed.
Teuchos::RefCountPtr< Epetra_RowMatrix > A_
Pointer to the Epetra_RowMatrix to factorize.
const char * Label() const
Returns a character string describing the operator.
const Epetra_Map * L_RangeMap_
Ifpack_ILU: A class for constructing and using an incomplete lower/upper (ILU) factorization of a giv...
int NumGlobalRows() const
Returns the number of global matrix rows.
Ifpack_CondestType
Ifpack_CondestType: enum to define the type of condition number estimate.
const Epetra_Map & OperatorRangeMap() const
Returns the Epetra_Map object associated with the range of this operator.
int SetLabel(const char *Label_in)
Sets label for this object.
long long NumGlobalRows64() const
double ApplyInverseFlops_
Contain sthe number of flops for ApplyInverse().
virtual int NumApplyInverse() const
Returns the number of calls to ApplyInverse().
int NumCompute_
Contains the number of successful call to Compute().
int Compute()
Compute ILU factors L and U using the specified graph, diagonal perturbation thresholds and relaxatio...
Ifpack_ScalingType enumerable type.
int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
int InitAllValues(const Epetra_RowMatrix &A, int MaxNumEntries)
int SetParameters(Teuchos::ParameterList ¶meterlist)
Set parameters using a Teuchos::ParameterList object.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.
int LevelOfFill() const
Returns the level of fill.
Epetra_Time Time_
Used for timing issues.
Teuchos::RefCountPtr< Epetra_CrsMatrix > L_
Contains the L factors.
int NumGlobalCols() const
Returns the number of global matrix columns.
double RelaxValue() const
Get ILU(k) relaxation parameter.
virtual double ApplyInverseFlops() const
Returns the number of flops in the application of the preconditioner.
long long IndexBase64() const
virtual double ComputeTime() const
Returns the time spent in Compute().
Teuchos::RefCountPtr< Epetra_Map > IlukRangeMap_
int Initialize()
Initialize the preconditioner, does not touch matrix values.
double ApplyInverseTime_
Contains the time for all successful calls to ApplyInverse().
int IndexBase() const
Returns the index base for row and column indices for this graph.
Ifpack_ILU & operator=(const Ifpack_ILU &RHS)
operator= (should never be used)
virtual int NumGlobalBlockDiagonals() const
Returns the number of diagonal entries found in the global input graph.
virtual double InitializeFlops() const
Returns the number of flops in the initialization phase.
Teuchos::RefCountPtr< Epetra_CrsMatrix > U_
Contains the U factors.
double AbsoluteThreshold() const
Get absolute threshold value.
bool HasNormInf() const
Returns false because this class cannot compute an Inf-norm.
bool IsInitialized() const
Returns true if the preconditioner has been successfully initialized.
const Epetra_RowMatrix & Matrix() const
Returns a reference to the matrix to be preconditioned.
virtual long long NumGlobalBlockDiagonals64() const
Ifpack_ILU(Epetra_RowMatrix *A)
Constructor.
virtual double ApplyInverseTime() const
Returns the time spent in ApplyInverse().
long long NumGlobalCols64() const
double NormInf() const
Returns 0.0 because this class cannot compute Inf-norm.
virtual double ComputeFlops() const
Returns the number of flops in the computation phase.
virtual std::ostream & Print(std::ostream &os) const
Prints on stream basic information about this object.
int Solve(bool Trans, const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Returns the result of a Ifpack_ILU forward/back solve on a Epetra_MultiVector X in Y...
int NumApplyInverse_
Contains the number of successful call to ApplyInverse().
int NumGlobalNonzeros() const
Returns the number of nonzero entries in the global graph.
bool UseTranspose() const
Returns the current UseTranspose setting.
int NumInitialize_
Contains the number of successful calls to Initialize().
Teuchos::RefCountPtr< Epetra_Map > IlukRowMap_
void Destroy()
Destroys all internal data.
Teuchos::RefCountPtr< Epetra_CrsGraph > U_Graph_
bool IsInitialized_
If true, the preconditioner has been successfully initialized.
const Epetra_Vector & D() const
Returns the address of the D factor associated with this factored matrix.
const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this operator.
Teuchos::RefCountPtr< Epetra_Map > IlukDomainMap_
double ComputeFlops_
Contains the number of flops for Compute().
int NumMyCols() const
Returns the number of local matrix columns.
int NumMyRows() const
Returns the number of local matrix rows.
const Ifpack_IlukGraph & Graph() const
Returns the address of the Ifpack_IlukGraph associated with this factored matrix. ...
long long NumGlobalNonzeros64() const
virtual int NumMyBlockDiagonals() const
Returns the number of diagonal entries found in the local input graph.
double ComputeTime_
Contains the time for all successful calls to Compute().
double Condest() const
Returns the computed estimated condition number, or -1.0 if not computed.
virtual int NumCompute() const
Returns the number of calls to Compute().
double Rthresh_
relative threshold
double RelaxValue_
Relaxation value.
int LevelOfFill_
Level of fill.
double Condest_
condition number estimate
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...