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" 64 #include "EpetraExt_Transpose_RowMatrix.h" 87 Ifpack_SORa(Epetra_RowMatrix* A);
102 bool IsInitialized()
const 104 return(IsInitialized_);
113 bool IsComputed()
const 126 int SetParameters(Teuchos::ParameterList& parameterlist);
138 int SetUseTranspose(
bool UseTranspose_in) {UseTranspose_ = UseTranspose_in;
return(0);};
143 int Apply(
const Epetra_MultiVector& X,
144 Epetra_MultiVector& Y)
const 146 return(Multiply(
false,X,Y));
149 int Multiply(
bool Trans,
const Epetra_MultiVector& X,
150 Epetra_MultiVector& Y)
const{
return A_->Multiply(Trans,X,Y);}
166 int ApplyInverse(
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const;
169 double Condest(
const Ifpack_CondestType CT = Ifpack_Cheap,
170 const int MaxIters = 1550,
171 const double Tol = 1e-9,
172 Epetra_RowMatrix* Matrix_in = 0);
175 double Condest()
const 184 const char* Label()
const {
return(Label_);}
187 int SetLabel(
const char* Label_in)
189 strcpy(Label_,Label_in);
194 double NormInf()
const {
return(0.0);};
197 bool HasNormInf()
const {
return(
false);};
200 bool UseTranspose()
const {
return(UseTranspose_);};
203 const Epetra_Map & OperatorDomainMap()
const {
return(A_->OperatorDomainMap());};
206 const Epetra_Map & OperatorRangeMap()
const{
return(A_->OperatorRangeMap());};
209 const Epetra_Comm & Comm()
const{
return(A_->Comm());};
212 const Epetra_RowMatrix& Matrix()
const 218 virtual double GetLambdaMax()
const{
return LambdaMax_;}
221 virtual double GetOmega()
const{
if(UseGlobalDamping_)
return 12.0/(11.0*LambdaMax_);
else return 1.0;}
224 virtual std::ostream& Print(std::ostream& os)
const;
227 virtual int NumInitialize()
const 229 return(NumInitialize_);
233 virtual int NumCompute()
const 239 virtual int NumApplyInverse()
const 241 return(NumApplyInverse_);
245 virtual double InitializeTime()
const 247 return(InitializeTime_);
251 virtual double ComputeTime()
const 253 return(ComputeTime_);
257 virtual double ApplyInverseTime()
const 259 return(ApplyInverseTime_);
263 virtual double InitializeFlops()
const 268 virtual double ComputeFlops()
const 270 return(ComputeFlops_);
273 virtual double ApplyInverseFlops()
const 275 return(ApplyInverseFlops_);
284 Ifpack_SORa(
const Ifpack_SORa& RHS) :
289 Ifpack_SORa& operator=(
const Ifpack_SORa& RHS)
308 int Solve(
bool Trans,
const Epetra_MultiVector& X, Epetra_MultiVector& Y)
const;
311 #ifndef EPETRA_NO_32BIT_GLOBAL_INDICES 312 int NumGlobalRows()
const {
return(A_->NumGlobalRows());};
316 int NumGlobalCols()
const {
return(A_->NumGlobalCols());};
318 long long NumGlobalRows64()
const {
return(A_->NumGlobalRows64());};
319 long long NumGlobalCols64()
const {
return(A_->NumGlobalCols64());};
322 int NumMyRows()
const {
return(A_->NumMyRows());};
325 int NumMyCols()
const {
return(A_->NumMyCols());};
328 int PowerMethod(
const int MaximumIterations,
double& lambda_max);
336 template<
typename int_type>
343 Teuchos::RefCountPtr<Epetra_CrsMatrix> Acrs_;
344 Teuchos::RefCountPtr<Epetra_RowMatrix> A_;
345 Teuchos::RefCountPtr<Epetra_CrsMatrix> W_;
346 Teuchos::RefCountPtr<Epetra_Vector> Wdiag_;
347 Teuchos::ParameterList List_;
373 bool HaveOAZBoundaries_;
375 bool UseInterprocDamping_;
377 bool UseGlobalDamping_;
382 mutable int NumApplyInverse_;
384 double InitializeTime_;
388 mutable double ApplyInverseTime_;
390 double ComputeFlops_;
392 mutable double ApplyInverseFlops_;
394 mutable Epetra_Time Time_;
Ifpack_ScalingType enumerable type.
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.