46 #include "Ifpack_ConfigDefs.h" 47 #include "Ifpack_Preconditioner.h" 49 #ifdef HAVE_IFPACK_EPETRAEXT 50 #include "Ifpack_Condest.h" 52 #include "Epetra_CompObject.h" 53 #include "Epetra_MultiVector.h" 54 #include "Epetra_Vector.h" 55 #include "Epetra_CrsGraph.h" 56 #include "Epetra_CrsMatrix.h" 57 #include "Epetra_BlockMap.h" 58 #include "Epetra_Map.h" 59 #include "Epetra_Object.h" 60 #include "Epetra_Comm.h" 61 #include "Epetra_CrsMatrix.h" 62 #include "Epetra_Time.h" 63 #include "Teuchos_RefCountPtr.hpp" 86 Ifpack_IHSS(Epetra_RowMatrix* A);
101 bool IsInitialized()
const 103 return(IsInitialized_);
112 bool IsComputed()
const 125 int SetParameters(Teuchos::ParameterList& parameterlist);
137 int SetUseTranspose(
bool UseTranspose_in) {UseTranspose_ = UseTranspose_in;
return(0);};
142 int Apply(
const Epetra_MultiVector& X,
143 Epetra_MultiVector& Y)
const 145 return(Multiply(
false,X,Y));
148 int Multiply(
bool Trans,
const Epetra_MultiVector& X,
149 Epetra_MultiVector& Y)
const{
return A_->Multiply(Trans,X,Y);}
165 int ApplyInverse(
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const;
168 double Condest(
const Ifpack_CondestType CT = Ifpack_Cheap,
169 const int MaxIters = 1550,
170 const double Tol = 1e-9,
171 Epetra_RowMatrix* Matrix_in = 0);
174 double Condest()
const 183 const char* Label()
const {
return(Label_);}
186 int SetLabel(
const char* Label_in)
188 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;
220 virtual int NumInitialize()
const 222 return(NumInitialize_);
226 virtual int NumCompute()
const 232 virtual int NumApplyInverse()
const 234 return(NumApplyInverse_);
238 virtual double InitializeTime()
const 240 return(InitializeTime_);
244 virtual double ComputeTime()
const 246 return(ComputeTime_);
250 virtual double ApplyInverseTime()
const 252 return(ApplyInverseTime_);
256 virtual double InitializeFlops()
const 261 virtual double ComputeFlops()
const 263 return(ComputeFlops_);
266 virtual double ApplyInverseFlops()
const 268 return(ApplyInverseFlops_);
277 Ifpack_IHSS(
const Ifpack_IHSS& RHS) :
282 Ifpack_IHSS& operator=(
const Ifpack_IHSS& RHS)
301 int Solve(
bool Trans,
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const;
304 #ifndef EPETRA_NO_32BIT_GLOBAL_INDICES 305 int NumGlobalRows()
const {
return(A_->NumGlobalRows());};
309 int NumGlobalCols()
const {
return(A_->NumGlobalCols());};
311 long long NumGlobalRows64()
const {
return(A_->NumGlobalRows64());};
312 long long NumGlobalCols64()
const {
return(A_->NumGlobalCols64());};
315 int NumMyRows()
const {
return(A_->NumMyRows());};
318 int NumMyCols()
const {
return(A_->NumMyCols());};
321 int PowerMethod(Epetra_Operator * Op,
const int MaximumIterations,
double& lambda_max);
333 Teuchos::RefCountPtr<Epetra_CrsMatrix> A_;
334 Teuchos::RefCountPtr<Epetra_RowMatrix> Aherm_;
335 Teuchos::RefCountPtr<Epetra_RowMatrix> Askew_;
336 Teuchos::RefCountPtr<Ifpack_Preconditioner> Pherm_;
337 Teuchos::RefCountPtr<Ifpack_Preconditioner> Pskew_;
339 Teuchos::ParameterList List_;
370 mutable int NumApplyInverse_;
372 double InitializeTime_;
376 mutable double ApplyInverseTime_;
378 double ComputeFlops_;
380 mutable double ApplyInverseFlops_;
382 mutable Epetra_Time Time_;
388 #pragma message("No EpetraEXT") 390 #warning "No EpetraEXT" Ifpack_ScalingType enumerable type.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.