47 #include "Teuchos_ParameterList.hpp" 48 #include "Teuchos_RefCountPtr.hpp" 52 using Teuchos::RefCountPtr;
56 #ifdef HAVE_IFPACK_EPETRAEXT 57 #include "EpetraExt_MatrixMatrix.h" 61 IsInitialized_(
false),
73 TEUCHOS_TEST_FOR_EXCEPT(!Acrs)
77 void Ifpack_IHSS::Destroy(){
82 int Ifpack_IHSS::Initialize(){
83 EigMaxIters_ = List_.get(
"ihss: eigenvalue max iterations",EigMaxIters_);
84 EigRatio_ = List_.get(
"ihss: ratio eigenvalue", EigRatio_);
85 NumSweeps_ = List_.get(
"ihss: sweeps",NumSweeps_);
93 int Ifpack_IHSS::SetParameters(Teuchos::ParameterList& parameterlist){
99 int Ifpack_IHSS::Compute(){
100 if(!IsInitialized_) Initialize();
106 Time_.ResetStartTime();
109 rv=EpetraExt::MatrixMatrix::Add(*A_,
false,.5,*A_,
true,.5,Aherm);
110 Aherm->FillComplete();
123 avec.MaxValue(&Alpha_);
126 for(
int i=0;i<Aherm->NumMyRows();i++) avec[i]+=Alpha_;
132 rv=EpetraExt::MatrixMatrix::Add(*A_,
false,.5,*A_,
true,-.5,Askew);
135 #ifndef EPETRA_NO_32BIT_GLOBAL_INDICES 138 int gid=Askew->
GRID(i);
144 #ifndef EPETRA_NO_64BIT_GLOBAL_INDICES 147 long long gid=Askew->
GRID64(i);
153 throw "Ifpack_IHSS::Compute: Unable to deduce GlobalIndices type";
159 Teuchos::ParameterList PLh=List_.sublist(
"ihss: hermetian list");
160 std::string htype=List_.get(
"ihss: hermetian type",
"ILU");
161 Pherm= Factory.
Create(htype, Aherm);
167 Teuchos::ParameterList PLs=List_.sublist(
"ihss: skew hermetian list");
168 std::string stype=List_.get(
"ihss: skew hermetian type",
"ILU");
169 Pskew= Factory.
Create(stype, Askew);
175 sprintf(Label_,
"IFPACK IHSS (H,S)=(%s/%s)",htype.c_str(),stype.c_str());
180 ComputeTime_ += Time_.ElapsedTime();
186 if(!IsComputed_)
return -1;
187 Time_.ResetStartTime();
188 bool initial_guess_is_zero=
false;
192 Teuchos::RefCountPtr<const Epetra_MultiVector> Xcopy;
198 initial_guess_is_zero=
true;
201 Xcopy = Teuchos::rcp( &X,
false );
215 for(
int i=0;i<NumSweeps_;i++){
217 if(!initial_guess_is_zero || i >0 ){
219 T2.Update(2*Alpha_,Y,-1,T1,1);
221 Pherm_->ApplyInverse(T2,Y);
225 T2.Scale(1.0,*Xcopy);
226 T2.Update(2*Alpha_,Y,-1,T1,1.0);
227 Pskew_->ApplyInverse(T2,Y);
232 ApplyInverseTime_ += Time_.ElapsedTime();
237 std::ostream& Ifpack_IHSS::Print(std::ostream& os)
const{
240 os<<
"Ifpack_IHSS"<<endl;
241 os<<
"-Hermetian preconditioner"<<endl;
242 os<<
"-Skew Hermetian preconditioner"<<endl;
256 int Ifpack_IHSS::PowerMethod(
Epetra_Operator * Op,
const int MaximumIterations,
double& lambda_max)
261 double RQ_top, RQ_bottom, norm;
270 for (
int iter = 0; iter < MaximumIterations; ++iter)
275 lambda_max = RQ_top / RQ_bottom;
const Epetra_Map & RowMap() const
virtual int InsertGlobalValues(int GlobalRow, int NumEntries, const double *Values, const int *Indices)
virtual const Epetra_Map & OperatorDomainMap() const=0
int FillComplete(bool OptimizeDataStorage=true)
Ifpack_CondestType
Ifpack_CondestType: enum to define the type of condition number estimate.
virtual int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
virtual const Epetra_Map & OperatorRangeMap() const=0
int Scale(double ScalarValue)
virtual int SetParameters(Teuchos::ParameterList &List)=0
Sets all parameters for the preconditioner.
int GRID(int LRID_in) const
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.
bool GlobalIndicesInt() const
long long GRID64(int LRID_in) const
static Ifpack_Preconditioner * Create(EPrecType PrecType, Epetra_RowMatrix *Matrix, const int overlap=0, bool overrideSerialDefault=false)
Creates an instance of Ifpack_Preconditioner given the enum value of the preconditioner type (can not...
virtual int Compute()=0
Computes all it is necessary to apply the preconditioner.
#define IFPACK_CHK_ERR(ifpack_err)
const Epetra_RowMatrix & A_
Reference to the matrix to be preconditioned, supposed symmetric.
bool GlobalIndicesLongLong() const
double ** Pointers() const