46#include "Ifpack_ConfigDefs.h"
47#include "Ifpack_CondestType.h"
49#include "Ifpack_Preconditioner.h"
50#include "Ifpack_IKLU_Utils.h"
51#include "Epetra_Vector.h"
52#include "Epetra_CrsMatrix.h"
53#include "Epetra_Time.h"
54#include "Teuchos_RefCountPtr.hpp"
113 return(IsInitialized_);
146 double Condest(
const Ifpack_CondestType CT = Ifpack_Cheap,
147 const int MaxIters = 1550,
148 const double Tol = 1e-9,
167 int SetUseTranspose(
bool UseTranspose_in) {UseTranspose_ = UseTranspose_in;
return(0);};
202 return(Label_.c_str());
213 virtual std::ostream&
Print(std::ostream& os)
const;
218 return(NumInitialize_);
230 return(NumApplyInverse_);
236 return(InitializeTime_);
242 return(ComputeTime_);
248 return(ApplyInverseTime_);
259 return(ComputeFlops_);
264 return(ApplyInverseFlops_);
267 inline double LevelOfFill()
const {
268 return(LevelOfFill_);
291 return(DropTolerance_);
295#ifndef EPETRA_NO_32BIT_GLOBAL_INDICES
301 long long NumGlobalNonzeros64()
const {
303 return(
L().NumGlobalNonzeros64() +
U().NumGlobalNonzeros64() -
L().NumGlobalRows64());
340 Teuchos::RefCountPtr<Epetra_CrsMatrix> L_;
342 Teuchos::RefCountPtr<Epetra_CrsMatrix> U_;
354 double DropTolerance_;
372 mutable int NumApplyInverse_;
374 double InitializeTime_;
378 mutable double ApplyInverseTime_;
380 double ComputeFlops_;
382 mutable double ApplyInverseFlops_;
386 long long GlobalNonzeros_;
387 Teuchos::RefCountPtr<Epetra_SerialComm> SerialComm_;
388 Teuchos::RefCountPtr<Epetra_Map> SerialMap_;
Ifpack_ScalingType enumerable type.
virtual const Epetra_Map & OperatorDomainMap() const=0
virtual const Epetra_Map & OperatorRangeMap() const=0
Ifpack_IKLU: A class for constructing and using an incomplete LU factorization of a given Epetra_RowM...
virtual double ApplyInverseTime() const
Returns the time spent in ApplyInverse().
bool IsInitialized() const
Returns true if the preconditioner has been successfully initialized.
virtual int NumCompute() const
Returns the number of calls to Compute().
double DropTolerance() const
Gets the dropping tolerance.
virtual ~Ifpack_IKLU()
Ifpack_IKLU Destructor.
int SetUseTranspose(bool UseTranspose_in)
If set true, transpose of this operator will be applied.
bool IsComputed() const
If factor is completed, this query returns true, otherwise it returns false.
int Compute()
Compute IC factor U using the specified graph, diagonal perturbation thresholds and relaxation parame...
int SetParameters(Teuchos::ParameterList ¶meterlis)
Set parameters using a Teuchos::ParameterList object.
const Epetra_RowMatrix & Matrix() const
Returns a reference to the matrix to be preconditioned.
const char * Label() const
Returns the label of this object.
int NumGlobalNonzeros() const
Returns the number of nonzero entries in the global graph.
virtual int NumInitialize() const
Returns the number of calls to Initialize().
bool HasNormInf() const
Returns false because this class cannot compute an Inf-norm.
double AbsoluteThreshold() const
Get absolute threshold value.
int NumMyNonzeros() const
Returns the number of nonzero entries in the local graph.
Ifpack_IKLU(const Epetra_RowMatrix *A)
Ifpack_IKLU constuctor with variable number of indices per row.
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.
double RelativeThreshold() const
Get relative threshold value.
virtual double InitializeTime() const
Returns the time spent in Initialize().
bool UseTranspose() const
Returns the current UseTranspose setting.
double Condest() const
Returns the computed estimated condition number, or -1.0 if no computed.
virtual double ComputeTime() const
Returns the time spent in Compute().
int Initialize()
Initialize L and U with values from user matrix A.
const Epetra_CrsMatrix & L() const
Returns a reference to the L factor.
const Epetra_CrsMatrix & U() const
Returns a reference to the U factor.
double RelaxValue() const
Set relative threshold value.
const Epetra_Map & OperatorRangeMap() const
Returns the Epetra_Map object associated with the range of this operator.
virtual int NumApplyInverse() const
Returns the number of calls to ApplyInverse().
const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this operator.
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.
virtual double InitializeFlops() const
Returns the number of flops in the initialization phase.
const Epetra_Comm & Comm() const
Returns the Epetra_BlockMap object associated with the range of this matrix operator.
virtual std::ostream & Print(std::ostream &os) const
Prints basic information on iostream. This function is used by operator<<.
int SetLabel(const char *Label_in)
Sets the label for this object.
virtual double ApplyInverseFlops() const
Returns the number of flops in the application of the preconditioner.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.