43 #ifndef IFPACK_POINTRELAXATION_H 44 #define IFPACK_POINTRELAXATION_H 46 #include "Ifpack_ConfigDefs.h" 47 #include "Ifpack_Preconditioner.h" 49 #include "Epetra_Vector.h" 50 #include "Epetra_Time.h" 51 #include "Epetra_RowMatrix.h" 52 #include "Epetra_Import.h" 54 #include "Teuchos_RefCountPtr.hpp" 59 class Epetra_MultiVector;
63 class Epetra_CrsMatrix;
156 UseTranspose_ = UseTranspose_in;
173 virtual inline int Apply(
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const 175 if (IsComputed() ==
false)
178 if (X.NumVectors() != Y.NumVectors())
181 IFPACK_CHK_ERR(Matrix_->Multiply(UseTranspose(),X,Y));
196 virtual int ApplyInverse(
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const;
207 virtual const char * Label()
const 209 return(Label_.c_str());
215 return(UseTranspose_);
225 virtual const Epetra_Comm & Comm()
const;
228 virtual const Epetra_Map & OperatorDomainMap()
const;
231 virtual const Epetra_Map & OperatorRangeMap()
const;
233 virtual int Initialize();
237 return(IsInitialized_);
247 virtual int Compute();
253 virtual const Epetra_RowMatrix&
Matrix()
const 259 virtual double Condest(
const Ifpack_CondestType CT = Ifpack_Cheap,
260 const int MaxIters = 1550,
261 const double Tol = 1e-9,
262 Epetra_RowMatrix* Matrix = 0);
271 virtual int SetParameters(Teuchos::ParameterList& List);
274 virtual std::ostream& Print(std::ostream & os)
const;
283 return(NumInitialize_);
295 return(NumApplyInverse_);
301 return(InitializeTime_);
307 return(ComputeTime_);
313 return(ApplyInverseTime_);
325 return(ComputeFlops_);
331 return(ApplyInverseFlops_);
341 virtual int ApplyInverseJacobi(
const Epetra_MultiVector& X,
342 Epetra_MultiVector& Y)
const;
345 virtual int ApplyInverseGS(
const Epetra_MultiVector& X,
346 Epetra_MultiVector& Y)
const;
348 virtual int ApplyInverseGS_RowMatrix(
const Epetra_MultiVector& X,
349 Epetra_MultiVector& Y)
const;
351 virtual int ApplyInverseGS_CrsMatrix(
const Epetra_CrsMatrix* A,
352 const Epetra_MultiVector& X,
353 Epetra_MultiVector& Y)
const;
355 virtual int ApplyInverseGS_FastCrsMatrix(
const Epetra_CrsMatrix* A,
356 const Epetra_MultiVector& X,
357 Epetra_MultiVector& Y)
const;
358 virtual int ApplyInverseGS_LocalFastCrsMatrix(
const Epetra_CrsMatrix* A,
359 const Epetra_MultiVector& X,
360 Epetra_MultiVector& Y)
const;
363 virtual int ApplyInverseSGS(
const Epetra_MultiVector& X,
364 Epetra_MultiVector& Y)
const;
366 virtual int ApplyInverseSGS_RowMatrix(
const Epetra_MultiVector& X,
367 Epetra_MultiVector& Y)
const;
369 virtual int ApplyInverseSGS_CrsMatrix(
const Epetra_CrsMatrix* A,
370 const Epetra_MultiVector& X,
371 Epetra_MultiVector& Y)
const;
373 virtual int ApplyInverseSGS_FastCrsMatrix(
const Epetra_CrsMatrix* A,
374 const Epetra_MultiVector& X,
375 Epetra_MultiVector& Y)
const;
377 virtual int ApplyInverseSGS_LocalFastCrsMatrix(
const Epetra_CrsMatrix* A,
378 const Epetra_MultiVector& X,
379 Epetra_MultiVector& Y)
const;
387 virtual void SetLabel();
409 mutable int NumApplyInverse_;
411 double InitializeTime_;
415 mutable double ApplyInverseTime_;
417 double ComputeFlops_;
419 mutable double ApplyInverseFlops_;
426 double DampingFactor_;
435 bool ComputeCondest_;
440 double MinDiagonalValue_;
449 long long NumGlobalRows_;
451 long long NumGlobalNonzeros_;
453 Teuchos::RefCountPtr<const Epetra_RowMatrix> Matrix_;
455 Teuchos::RefCountPtr<Epetra_Import> Importer_;
457 mutable Teuchos::RefCountPtr<Epetra_Vector> Diagonal_;
459 Teuchos::RefCountPtr<Epetra_Time> Time_;
463 bool ZeroStartingSolution_;
472 int NumLocalSmoothingIndices_;
474 int * LocalSmoothingIndices_;
482 #endif // IFPACK_POINTRELAXATION_H virtual int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Applies the matrix to an Epetra_MultiVector.
virtual int NumApplyInverse() const
Returns the number of calls to ApplyInverse().
virtual double ComputeFlops() const
Returns the number of flops in the computation phase.
virtual int SetUseTranspose(bool UseTranspose_in)
virtual double ApplyInverseFlops() const
Returns the number of flops for the application of the preconditioner.
virtual int NumCompute() const
Returns the number of calls to Compute().
virtual double Condest() const
Returns the condition number estimate, or -1.0 if not computed.
virtual double ComputeTime() const
Returns the time spent in Compute().
virtual double NormInf() const
Returns the infinity norm of the global matrix (not implemented)
virtual bool UseTranspose() const
Returns the current UseTranspose setting.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.
virtual bool IsInitialized() const
Returns true if the preconditioner has been successfully initialized, false otherwise.
Ifpack_PointRelaxation: a class to define point relaxation preconditioners of for Epetra_RowMatrix's...
virtual const Epetra_RowMatrix & Matrix() const
Returns a pointer to the matrix to be preconditioned.
virtual bool IsComputed() const
Returns true if the preconditioner has been successfully computed.
virtual int NumInitialize() const
Returns the number of calls to Initialize().
virtual double InitializeFlops() const
Returns the number of flops in the initialization phase.
virtual bool HasNormInf() const
Returns true if the this object can provide an approximate Inf-norm, false otherwise.
virtual double InitializeTime() const
Returns the time spent in Initialize().
virtual double ApplyInverseTime() const
Returns the time spent in ApplyInverse().
virtual ~Ifpack_PointRelaxation()
Destructor.