47#include "Teuchos_ParameterList.hpp"
48#include "Teuchos_RefCountPtr.hpp"
52using Teuchos::RefCountPtr;
56#ifdef HAVE_IFPACK_EPETRAEXT
61 IsInitialized_(
false),
73 TEUCHOS_TEST_FOR_EXCEPT(!Acrs)
77void Ifpack_IHSS::Destroy(){
82int 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_);
93int Ifpack_IHSS::SetParameters(Teuchos::ParameterList& parameterlist){
99int Ifpack_IHSS::Compute(){
100 if(!IsInitialized_) Initialize();
106 Time_.ResetStartTime();
110 Aherm->FillComplete();
123 avec.MaxValue(&Alpha_);
126 for(
int i=0;i<Aherm->NumMyRows();i++) avec[i]+=Alpha_;
135#ifndef EPETRA_NO_32BIT_GLOBAL_INDICES
136 if(Askew->RowMap().GlobalIndicesInt()) {
137 for(
int i=0;i<Askew->NumMyRows();i++) {
138 int gid=Askew->GRID(i);
139 Askew->InsertGlobalValues(gid,1,&Alpha_,&gid);
144#ifndef EPETRA_NO_64BIT_GLOBAL_INDICES
145 if(Askew->RowMap().GlobalIndicesLongLong()) {
146 for(
int i=0;i<Askew->NumMyRows();i++) {
147 long long gid=Askew->GRID64(i);
148 Askew->InsertGlobalValues(gid,1,&Alpha_,&gid);
153 throw "Ifpack_IHSS::Compute: Unable to deduce GlobalIndices type";
155 Askew->FillComplete();
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);
162 Pherm->SetParameters(PLh);
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);
170 Pskew->SetParameters(PLs);
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;
194 if (X.
Pointers()[0] == Y.Pointers()[0]){
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();
237std::ostream& Ifpack_IHSS::Print(std::ostream& os)
const{
240 os<<
"Ifpack_IHSS"<<endl;
241 os<<
"-Hermetian preconditioner"<<endl;
242 os<<
"-Skew Hermetian preconditioner"<<endl;
256int 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;
Ifpack_CondestType
Ifpack_CondestType: enum to define the type of condition number estimate.
#define IFPACK_CHK_ERR(ifpack_err)
static int Add(const Epetra_CrsMatrix &A, bool transposeA, double scalarA, Epetra_CrsMatrix &B, double scalarB)
double ** Pointers() const
virtual int Apply(const Epetra_MultiVector &X, Epetra_MultiVector &Y) const=0
virtual const Epetra_Map & OperatorDomainMap() const=0
virtual const Epetra_Map & OperatorRangeMap() const=0
Ifpack_Preconditioner: basic class for preconditioning in Ifpack.