43 #ifndef IFPACK_POLYNOMIAL_H 44 #define IFPACK_POLYNOMIAL_H 48 #include "Teuchos_RefCountPtr.hpp" 49 #include "Teuchos_LAPACK.hpp" 50 #include "Teuchos_SerialDenseMatrix.hpp" 65 #ifdef HAVE_IFPACK_EPETRAEXT 225 const int MaxIters = 1550,
226 const double Tol = 1e-9,
239 virtual std::ostream&
Print(std::ostream & os)
const;
305 const int MaximumIterations,
311 const int MaximumIterations,
312 double& lambda_min,
double& lambda_max);
314 #ifdef HAVE_IFPACK_EPETRAEXT 317 int CG(
const int MaximumIterations,
318 double& lambda_min,
double& lambda_max);
321 int PowerMethod(
const int MaximumIterations,
double& lambda_max);
327 const int MaximumIterations,
328 double& lambda_real_min,
double& lambda_real_max,
329 double& lambda_imag_min,
double& lambda_imag_max);
389 bool ComputeCondest_;
417 Teuchos::RefCountPtr<const Epetra_RowMatrix>
Matrix_;
422 #ifdef HAVE_IFPACK_EPETRAEXT 423 Teuchos::ParameterList BlockList_;
425 Teuchos::RefCountPtr<EpetraExt_PointToBlockDiagPermute> InvBlockDiagonal_;
435 Teuchos::RefCountPtr<Epetra_Time>
Time_;
444 #endif // IFPACK_POLYNOMIAL_H
int GMRES(const Epetra_Operator &Operator, const Epetra_Vector &InvPointDiagonal, const int MaximumIterations, double &lambda_real_min, double &lambda_real_max, double &lambda_imag_min, double &lambda_imag_max)
Uses AztecOO's GMRES to estimate the height and width of the spectrum.
bool IsInitialized_
If true, the preconditioner has been computed successfully.
Ifpack_Polynomial(const Ifpack_Polynomial &)
Copy constructor (PRIVATE, should not be used)
bool ZeroStartingSolution_
If true, the starting solution is always the zero vector.
Teuchos::RefCountPtr< Epetra_Vector > InvDiagonal_
Contains the inverse of diagonal elements of Matrix.
virtual const char * Label() const
int NumCompute_
Contains the number of successful call to Compute().
virtual bool IsComputed() const
Returns true if the preconditioner has been successfully computed.
double ApplyInverseFlops_
Contain sthe number of flops for ApplyInverse().
virtual double ApplyInverseFlops() const
Returns the number of flops for the application of the preconditioner.
virtual double ApplyInverseTime() const
Returns the time spent in ApplyInverse().
double Condest_
Contains the estimated condition number.
virtual std::ostream & Print(std::ostream &os) const
Prints object to an output stream.
virtual int Initialize()
Computes all it is necessary to initialize the preconditioner.
double ApplyInverseTime_
Contains the time for all successful calls to ApplyInverse().
virtual double InitializeTime() const
Returns the time spent in Initialize().
long long NumGlobalNonzeros_
Number of global nonzeros.
int LSPointsReal_
Contains the number of discretization points of the least squares problem.
double RealEigRatio_
Contains the ratio such that the rectangular domain in the complex plane is [-LambdaRealMax_ / EigRat...
virtual int SetUseTranspose(bool UseTranspose_in)
virtual int ApplyInverse(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Applies the preconditioner to X, returns the result in Y.
long long NumGlobalRows_
Number of global rows.
virtual bool IsInitialized() const
Returns true if the preconditioner has been successfully initialized, false otherwise.
virtual double ComputeFlops() const
Returns the number of flops in the computation phase.
virtual bool HasNormInf() const
Returns true if the this object can provide an approximate Inf-norm, false otherwise.
double ComputeFlops_
Contains the number of flops for Compute().
bool UseTranspose_
If true, use the tranpose of Matrix_.
Teuchos::RefCountPtr< const Epetra_RowMatrix > Matrix_
Pointers to the matrix to be preconditioned as an Epetra_RowMatrix.
std::string Label_
Contains the label of this object.
virtual int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const
Applies the matrix to an Epetra_MultiVector.
int NumApplyInverse_
Contains the number of successful call to ApplyInverse().
Teuchos::RefCountPtr< Epetra_Time > Time_
Time object to track timing.
Ifpack_CondestType
Ifpack_CondestType: enum to define the type of condition number estimate.
static int CG(const Epetra_Operator &Operator, const Epetra_Vector &InvPointDiagonal, const int MaximumIterations, double &lambda_min, double &lambda_max)
Uses AztecOO's CG to estimate lambda_min and lambda_max.
Ifpack_Polynomial & operator=(const Ifpack_Polynomial &)
operator = (PRIVATE, should not be used)
virtual ~Ifpack_Polynomial()
Destructor.
virtual int Compute()
Computes the preconditioners.
virtual int NumApplyInverse() const
Returns the number of calls to ApplyInverse().
std::vector< double > coeff_
coefficients of the polynomial
int PolyDegree_
Contains the degree of the least squares polynomial.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.
virtual const Epetra_Map & OperatorRangeMap() const
Returns the Epetra_Map object associated with the range of this operator.
Teuchos::RefCountPtr< const Epetra_Operator > Operator_
Pointers to the matrix to be preconditioned as an Epetra_Operator.
double InitializeTime_
Contains the time for all successful calls to Initialize().
double MinDiagonalValue_
Contains the minimum value on the diagonal.
virtual double Condest() const
Returns the condition number estimate, or -1.0 if not computed.
Ifpack_Polynomial(const Epetra_Operator *Matrix)
Ifpack_Polynomial constructor with given Epetra_Operator/Epetra_RowMatrix.
virtual const Epetra_Comm & Comm() const
Returns a pointer to the Epetra_Comm communicator associated with this operator.
int NumMyNonzeros_
Number of local nonzeros.
double ComputeTime_
Contains the time for all successful calls to Compute().
virtual bool UseTranspose() const
Returns the current UseTranspose setting.
virtual int SetParameters(Teuchos::ParameterList &List)
Sets all the parameters for the preconditioner.
virtual double InitializeFlops() const
Returns the number of flops in the initialization phase.
virtual double ComputeTime() const
Returns the time spent in Compute().
Ifpack_Polynomial: class for preconditioning with least squares polynomials in Ifpack.
bool IsComputed_
If true, the preconditioner has been computed successfully.
int NumMyRows_
Number of local rows.
virtual const Epetra_Map & OperatorDomainMap() const
Returns the Epetra_Map object associated with the domain of this operator.
bool IsIndefinite_
If true, have to compute polynomial for a spectrum with negative eigenvalues.
virtual int NumInitialize() const
Returns the number of calls to Initialize().
int EigMaxIters_
Max number of iterations to use in eigenvalue estimation (if automatic).
bool SolveNormalEquations_
Run on the normal equations.
static int PowerMethod(const Epetra_Operator &Operator, const Epetra_Vector &InvPointDiagonal, const int MaximumIterations, double &LambdaMax)
Simple power method to compute lambda_max.
bool UseBlockMode_
Use Block Preconditioning.
bool IsRowMatrix_
If true, the Operator_ is an Epetra_RowMatrix.
double LambdaRealMin_
Bounds on the spectrum.
int NumInitialize_
Contains the number of successful calls to Initialize().
bool IsComplex_
If true, have to compute polynomial for a spectrum with nonzero imaginary part.
virtual int NumCompute() const
Returns the number of calls to Compute().
virtual void SetLabel()
Sets the label.
virtual double NormInf() const
Returns the infinity norm of the global matrix (not implemented)
virtual const Epetra_RowMatrix & Matrix() const
Returns a pointer to the matrix to be preconditioned.