FgmresIte.h

// =======================================================================
// @(#)ItMethod.h        1.6 10/18/99
// 
// Class:       TFgmresIte
// Purpose:     defines the fixed point iteration
//
// Author:      Volker John
//
// History:     start of implementation 24.10.2000
//
// =======================================================================
#ifndef __FGMRESITE__
#define __FGMRESITE__

#ifdef _MPI   
   #ifdef __3D__
    #include <ParFECommunicator3D.h>
   #else
    #include <ParFECommunicator2D.h>
   #endif
#endif 
      
#include <ItMethod.h>

class TFgmresIte : public TItMethod
{
  protected:

  double *s;
  double *cosi;
  double *sn;

  double *d;
  double *z;
  
  double **H;
  double **v;
  double **zv;

#ifdef _MPI   
   #ifdef __3D__
      TParFECommunicator3D *ParComm;
   #else
      TParFECommunicator2D *ParComm;
   #endif
#endif  
      
  public:
    TFgmresIte(MatVecProc *MatVec, DefectProc *Defect, TItMethod *Prec,
               int n_aux, int N_Unknowns, int scalar);
 
#ifdef _MPI
    TFgmresIte(MatVecProc *MatVec, DefectProc *Defect, TItMethod *Prec,
               int n_aux, int N_Unknowns, int scalar,
  #ifdef  __3D__                   
                   TParFECommunicator3D *parcomm
  #else                
                               TParFECommunicator2D *parcomm
  #endif
          );
#endif

    virtual ~TFgmresIte();
    
    int Iterate(TSquareMatrix **A, TMatrix **B, double *sol, 
                double *rhs);    
};
#endif