51 #include "Epetra_Vector.h"
52 #include "Epetra_CrsMatrix.h"
53 #include "Epetra_Time.h"
54 #include "Teuchos_RefCountPtr.hpp"
147 const int MaxIters = 1550,
148 const double Tol = 1e-9,
213 virtual std::ostream&
Print(std::ostream& os)
const;
295 #ifndef EPETRA_NO_32BIT_GLOBAL_INDICES
340 Teuchos::RefCountPtr<Epetra_CrsMatrix>
L_;
342 Teuchos::RefCountPtr<Epetra_CrsMatrix>
U_;
const Epetra_CrsMatrix & U() const
Returns a reference to the U factor.
Ifpack_IKLU(const Epetra_RowMatrix *A)
Ifpack_IKLU constuctor with variable number of indices per row.
const Epetra_RowMatrix & Matrix() const
Returns a reference to the matrix to be preconditioned.
bool IsComputed_
true if this object has been computed
const Epetra_RowMatrix & A_
reference to the matrix to be preconditioned.
double ApplyInverseTime_
Contains the time for all successful calls to ApplyInverse().
const Epetra_CrsMatrix & L() const
Returns a reference to the L factor.
virtual std::ostream & Print(std::ostream &os) const
Prints basic information on iostream. This function is used by operator<<.
bool IsInitialized_
true if this object has been initialized
Teuchos::RefCountPtr< Epetra_SerialComm > SerialComm_
int SetParameters(Teuchos::ParameterList ¶meterlis)
Set parameters using a Teuchos::ParameterList object.
long long NumGlobalNonzeros64() const
int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
virtual const Epetra_Map & OperatorDomainMap() const =0
virtual int NumApplyInverse() const
Returns the number of calls to ApplyInverse().
int NumMyNonzeros_
Number of local nonzeros.
int NumApplyInverse_
Contains the number of successful call to ApplyInverse().
const Epetra_Map & OperatorRangeMap() const
Returns the Epetra_Map object associated with the range of this operator.
Epetra_Time Time_
Used for timing purposed.
virtual double InitializeTime() const
Returns the time spent in Initialize().
int NumGlobalNonzeros() const
Returns the number of nonzero entries in the global graph.
double DropTolerance_
Discards all elements below this tolerance.
const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this operator.
virtual double ComputeTime() const
Returns the time spent in Compute().
int NumMyNonzeros() const
Returns the number of nonzero entries in the local graph.
Ifpack_CondestType
Ifpack_CondestType: enum to define the type of condition number estimate.
csr * csrA_
Containers for the matrix storage and permutation.
virtual const Epetra_Map & OperatorRangeMap() const =0
long long GlobalNonzeros_
Global number of nonzeros in L and U factors.
double ApplyInverseFlops_
Contain sthe number of flops for ApplyInverse().
Ifpack_ScalingType enumerable type.
double LevelOfFill_
Level-of-fill.
Teuchos::RefCountPtr< Epetra_CrsMatrix > U_
U factor.
virtual double ApplyInverseFlops() const
Returns the number of flops in the application of the preconditioner.
const char * Label() const
Returns the label of this object.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.
Ifpack_IKLU & operator=(const Ifpack_IKLU &)
operator= (should never be used)
bool IsInitialized() const
Returns true if the preconditioner has been successfully initialized.
int NumCompute_
Contains the number of successful call to Compute().
virtual double ApplyInverseTime() const
Returns the time spent in ApplyInverse().
bool IsComputed() const
If factor is completed, this query returns true, otherwise it returns false.
double Athresh_
Absolute threshold.
int NumInitialize_
Contains the number of successful calls to Initialize().
bool UseTranspose() const
Returns the current UseTranspose setting.
std::string Label_
Label for this object.
virtual int NumCompute() const
Returns the number of calls to Compute().
double Condest_
Condition number estimate.
double LevelOfFill() const
bool HasNormInf() const
Returns false because this class cannot compute an Inf-norm.
virtual ~Ifpack_IKLU()
Ifpack_IKLU Destructor.
double DropTolerance() const
Gets the dropping tolerance.
double Rthresh_
Relative threshold.
virtual double InitializeFlops() const
Returns the number of flops in the initialization phase.
double InitializeTime_
Contains the time for all successful calls to Initialize().
int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Returns the result of a Ifpack_IKLU forward/back solve on a Epetra_MultiVector X in Y...
double Relax_
relaxation value
Ifpack_IKLU: A class for constructing and using an incomplete LU factorization of a given Epetra_RowM...
double RelativeThreshold() const
Get relative threshold value.
double RelaxValue() const
Set relative threshold value.
Teuchos::RefCountPtr< Epetra_CrsMatrix > L_
L factor.
bool UseTranspose_
true if transpose has to be used.
const Epetra_Comm & Comm() const
Returns the Epetra_BlockMap object associated with the range of this matrix operator.
virtual int NumInitialize() const
Returns the number of calls to Initialize().
csrn * csrnN_
Container for the L and U factor.
double NormInf() const
Returns 0.0 because this class cannot compute Inf-norm.
Teuchos::RefCountPtr< Epetra_Map > SerialMap_
const Epetra_Comm & Comm_
Reference to the communicator object.
int SetUseTranspose(bool UseTranspose_in)
If set true, transpose of this operator will be applied.
int SetLabel(const char *Label_in)
Sets the label for this object.
double ComputeTime_
Contains the time for all successful calls to Compute().
double AbsoluteThreshold() const
Get absolute threshold value.
int NumMyRows_
Number of local rows.
int Initialize()
Initialize L and U with values from user matrix A.
double ComputeFlops_
Contains the number of flops for Compute().
double Condest() const
Returns the computed estimated condition number, or -1.0 if no computed.
Ifpack_IKLU(const Ifpack_IKLU &RHS)
Copy constructor (should never be used)
void Destroy()
Releases all allocated memory.
int Compute()
Compute IC factor U using the specified graph, diagonal perturbation thresholds and relaxation parame...
virtual double ComputeFlops() const
Returns the number of flops in the computation phase.