43 #ifndef IFPACK_SPARSKIT_H 44 #define IFPACK_SPARSKIT_H 47 #ifdef HAVE_IFPACK_SPARSKIT 51 #include "Epetra_Vector.h" 52 #include "Epetra_RowMatrix.h" 55 class Epetra_MultiVector;
66 Ifpack_SPARSKIT(Epetra_RowMatrix* A);
69 virtual ~Ifpack_SPARSKIT();
81 int SetParameter(
const std::string Name,
const int Value)
85 int SetParameter(
const std::string Name,
const double Value)
90 const Epetra_RowMatrix&
Matrix()
const 102 return(IsInitialized_);
130 int ApplyInverse(
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const;
132 int Apply(
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const 147 const int MaxIters = 1550,
148 const double Tol = 1e-9,
149 Epetra_RowMatrix* Matrix = 0);
170 int SetUseTranspose(
bool UseTranspose) {UseTranspose_ = UseTranspose;
return(0);};
173 double NormInf()
const {
return(0.0);};
176 bool HasNormInf()
const {
return(
false);};
179 bool UseTranspose()
const {
return(UseTranspose_);};
182 const Epetra_Map & OperatorDomainMap()
const {
return(A_.OperatorDomainMap());};
185 const Epetra_Map & OperatorRangeMap()
const{
return(A_.OperatorRangeMap());};
188 const Epetra_Comm & Comm()
const{
return(Comm_);};
191 const char* Label()
const 193 return(Label_.c_str());
196 int SetLabel(
const char* Label)
203 virtual std::ostream&
Print(std::ostream& os)
const;
208 return(NumInitialize_);
220 return(NumApplyInverse_);
226 return(InitializeTime_);
232 return(ComputeTime_);
238 return(ApplyInverseTime_);
258 Epetra_RowMatrix& A_;
259 const Epetra_Comm& Comm_;
270 std::vector<double> alu_;
271 std::vector<int> jlu_;
272 std::vector<int> ju_;
276 std::vector<int> iperm_;
288 mutable int NumApplyInverse_;
291 double InitializeTime_;
295 mutable double ApplyInverseTime_;
298 double ComputeFlops_;
300 mutable double ApplyInverseFlops_;
304 #endif // HAVE_IFPACK_SPARSKIT virtual int NumInitialize() const =0
Returns the number of calls to Initialize().
virtual double ComputeTime() const =0
Returns the time spent in Compute().
virtual double ComputeFlops() const =0
Returns the number of flops in the computation phase.
virtual double ApplyInverseTime() const =0
Returns the time spent in ApplyInverse().
virtual double ApplyInverseFlops() const =0
Returns the number of flops in the application of the preconditioner.
virtual const Epetra_RowMatrix & Matrix() const =0
Returns a pointer to the matrix to be preconditioned.
virtual bool IsInitialized() const =0
Returns true if the preconditioner has been successfully initialized, false otherwise.
virtual std::ostream & Print(std::ostream &os) const =0
Prints basic information on iostream. This function is used by operator<<.
virtual int Initialize()=0
Computes all it is necessary to initialize the preconditioner.
virtual double InitializeTime() const =0
Returns the time spent in Initialize().
Ifpack_CondestType
Ifpack_CondestType: enum to define the type of condition number estimate.
virtual int SetParameters(Teuchos::ParameterList &List)=0
Sets all parameters for the preconditioner.
virtual double Condest() const =0
Returns the computed condition number estimate, or -1.0 if not computed.
Ifpack_ScalingType enumerable type.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.
virtual int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const =0
Applies the preconditioner to vector X, returns the result in Y.
virtual double InitializeFlops() const =0
Returns the number of flops in the initialization phase.
virtual int NumCompute() const =0
Returns the number of calls to Compute().
virtual bool IsComputed() const =0
Returns true if the preconditioner has been successfully computed, false otherwise.
virtual int Compute()=0
Computes all it is necessary to apply the preconditioner.
#define IFPACK_CHK_ERR(ifpack_err)
virtual int NumApplyInverse() const =0
Returns the number of calls to ApplyInverse().