50#include "Epetra_Vector.h"
51#include "Epetra_CrsMatrix.h"
52#include "Epetra_Time.h"
53#include "Teuchos_RefCountPtr.hpp"
171 const int MaxIters = 1550,
172 const double Tol = 1e-9,
183#ifndef EPETRA_NO_32BIT_GLOBAL_INDICES
239 virtual std::ostream&
Print(std::ostream& os)
const;
351 Teuchos::RefCountPtr<Epetra_CrsMatrix>
H_;
397 template<
typename int_type>
Ifpack_CondestType
Ifpack_CondestType: enum to define the type of condition number estimate.
@ Ifpack_Cheap
cheap estimate
Ifpack_ScalingType enumerable type.
virtual const Epetra_Map & OperatorDomainMap() const=0
virtual const Epetra_Map & OperatorRangeMap() const=0
Ifpack_ICT: A class for constructing and using an incomplete Cholesky factorization of a given Epetra...
const Epetra_Comm & Comm() const
Returns the Epetra_BlockMap object associated with the range of this matrix operator.
double Rthresh_
Relative threshold.
double LevelOfFill() const
Returns the level-of-fill.
const Epetra_RowMatrix & A_
Reference to the matrix to be preconditioned, supposed symmetric.
double ComputeFlops_
Contains the number of flops for Compute().
int NumGlobalNonzeros() const
Returns the number of nonzero entries in the global graph.
const Epetra_Comm & Comm_
Reference to the communicator.
double DropTolerance_
During factorization, drop all values below this.
double RelaxValue() const
Returns the relaxation value.
double Relax_
Relaxation value.
bool UseTranspose_
If true, use the transpose of the matrix.
int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Returns the result of a Ifpack_ICT forward/back solve on a Epetra_MultiVector X in Y.
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.
virtual int NumCompute() const
Returns the number of calls to Compute().
Ifpack_ICT & operator=(const Ifpack_ICT &)
Should not be used.
virtual double InitializeTime() const
Returns the time spent in Initialize().
virtual std::ostream & Print(std::ostream &os) const
Prints basic information on iostream. This function is used by operator<<.
bool UseTranspose() const
Returns the current UseTranspose setting.
double DropTolerance() const
Returns the drop threshold.
double Condest_
Contains the estimate of the condition number, if -1.0 if not computed.
virtual double ApplyInverseTime() const
Returns the time spent in ApplyInverse().
double Condest() const
Returns the computed condition number estimate, or -1.0 if not computed.
int NumInitialize_
Contains the number of successful calls to Initialize().
virtual double ComputeTime() const
Returns the time spent in Compute().
void Destroy()
Destroys all data associated to the preconditioner.
double Athresh_
Absolute threshold.
bool HasNormInf() const
Returns false because this class cannot compute an Inf-norm.
Teuchos::RefCountPtr< Epetra_SerialComm > SerialComm_
double ComputeTime_
Contains the time for all successful calls to Compute().
long long NumGlobalNonzeros64() const
const Epetra_CrsMatrix & H() const
Returns the address of the D factor associated with this factored matrix.
virtual double InitializeFlops() const
Returns the number of flops in the initialization phase.
double ApplyInverseTime_
Contains the time for all successful calls to ApplyInverse().
const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this operator.
Ifpack_ICT(const Ifpack_ICT &rhs)
Should not be used.
long long GlobalNonzeros_
Global number of nonzeros in L and U factors.
double RelativeThreshold() const
Returns the relative threshold.
double InitializeTime_
Contains the time for all successful calls to Initialize().
virtual int NumApplyInverse() const
Returns the number of calls to ApplyInverse().
int SetParameters(Teuchos::ParameterList ¶meterlis)
Set parameters using a Teuchos::ParameterList object.
bool IsComputed() const
If factor is completed, this query returns true, otherwise it returns false.
double AbsoluteThreshold() const
Returns the absolute threshold.
virtual double ApplyInverseFlops() const
Returns the number of flops in all applications of ApplyInverse().
int Initialize()
Initialize L and U with values from user matrix A.
virtual ~Ifpack_ICT()
Ifpack_ICT Destructor.
bool IsComputed_
If true, the preconditioner has been successfully computed.
Teuchos::RefCountPtr< Epetra_Map > SerialMap_
double LevelOfFill_
Level of fill.
int NumMyRows_
Number of local rows in the matrix.
int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
bool IsInitialized() const
Returns true is the preconditioner has been successfully initialized.
const char * Label() const
int NumApplyInverse_
Contains the number of successful call to ApplyInverse().
Epetra_Time Time_
Used for timing purposes.
int NumCompute_
Contains the number of successful call to Compute().
virtual double ComputeFlops() const
Returns the number of flops in all applications of Compute().
int Compute()
Compute IC factor U using the specified graph, diagonal perturbation thresholds and relaxation parame...
int SetLabel(const char *Label_in)
Teuchos::RefCountPtr< Epetra_CrsMatrix > H_
Contains the Cholesky factorization.
const Epetra_Map & OperatorRangeMap() const
Returns the Epetra_Map object associated with the range of this operator.
std::string Label_
Label of this object.
virtual int NumInitialize() const
Returns the number of calls to Initialize().
double NormInf() const
Returns 0.0 because this class cannot compute Inf-norm.
double ApplyInverseFlops_
Contain sthe number of flops for ApplyInverse().
const Epetra_RowMatrix & Matrix() const
Returns a reference to the matrix to be preconditioned.
bool IsInitialized_
If true, the preconditioner has been successfully initialized.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.