46 #include "Ifpack_ConfigDefs.h" 47 #ifdef HAVE_IFPACK_HIPS 50 #include "Ifpack_Preconditioner.h" 51 #include "Ifpack_Condest.h" 53 #include "Epetra_CompObject.h" 54 #include "Epetra_MultiVector.h" 55 #include "Epetra_Vector.h" 56 #include "Epetra_CrsGraph.h" 57 #include "Epetra_CrsMatrix.h" 58 #include "Epetra_BlockMap.h" 59 #include "Epetra_Map.h" 60 #include "Epetra_Object.h" 61 #include "Epetra_Comm.h" 62 #include "Epetra_CrsMatrix.h" 63 #include "Epetra_Time.h" 64 #include "Teuchos_RefCountPtr.hpp" 85 Ifpack_HIPS(Epetra_RowMatrix* A);
102 return(IsInitialized_);
136 int SetUseTranspose(
bool UseTranspose_in) {UseTranspose_ = UseTranspose_in;
return(0);};
141 int Apply(
const Epetra_MultiVector& X,
142 Epetra_MultiVector& Y)
const 144 return(Multiply(
false,X,Y));
147 int Multiply(
bool Trans,
const Epetra_MultiVector& X,
148 Epetra_MultiVector& Y)
const{
return A_->Multiply(Trans,X,Y);}
164 int ApplyInverse(
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const;
167 double Condest(
const Ifpack_CondestType CT = Ifpack_Cheap,
168 const int MaxIters = 1550,
169 const double Tol = 1e-9,
170 Epetra_RowMatrix* Matrix_in = 0);
182 const char* Label()
const {
return(Label_);}
185 int SetLabel(
const char* Label_in)
187 strcpy(Label_,Label_in);
193 double NormInf()
const {
return(0.0);};
196 bool HasNormInf()
const {
return(
false);};
199 bool UseTranspose()
const {
return(UseTranspose_);};
202 const Epetra_Map & OperatorDomainMap()
const {
return(A_->OperatorDomainMap());};
205 const Epetra_Map & OperatorRangeMap()
const{
return(A_->OperatorRangeMap());};
208 const Epetra_Comm & Comm()
const{
return(A_->Comm());};
211 const Epetra_RowMatrix&
Matrix()
const 217 virtual std::ostream&
Print(std::ostream& os)
const;
222 return(NumInitialize_);
234 return(NumApplyInverse_);
240 return(InitializeTime_);
246 return(ComputeTime_);
252 return(ApplyInverseTime_);
263 return(ComputeFlops_);
268 return(ApplyInverseFlops_);
277 Ifpack_HIPS(
const Ifpack_HIPS& RHS) :
282 Ifpack_HIPS& operator=(
const Ifpack_HIPS& RHS)
301 int Solve(
bool Trans,
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const;
305 int NumGlobalRows()
const {
return(A_->NumGlobalRows());};
308 int NumGlobalCols()
const {
return(A_->NumGlobalCols());};
311 int NumMyRows()
const {
return(A_->NumMyRows());};
314 int NumMyCols()
const {
return(A_->NumMyCols());};
326 Teuchos::RefCountPtr<Epetra_RowMatrix> A_;
327 Teuchos::ParameterList List_;
328 Teuchos::RefCountPtr<Epetra_Map> RowMap0_;
329 Teuchos::RefCountPtr<Epetra_Map> ColMap0_;
336 bool ValuesInitialized_;
354 mutable int NumApplyInverse_;
356 double InitializeTime_;
360 mutable double ApplyInverseTime_;
362 double ComputeFlops_;
364 mutable double ApplyInverseFlops_;
366 mutable Epetra_Time Time_;
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().
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.
virtual int NumApplyInverse() const =0
Returns the number of calls to ApplyInverse().