/*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*/
/*             ********   ***                                 SparseLib++    */
/*          *******  **  ***       ***      ***               v. 1.5c        */
/*           *****      ***     ******** ********                            */
/*            *****    ***     ******** ********              R. Pozo        */
/*       **  *******  ***   **   ***      ***                 K. Remington   */
/*        ********   ********                                 A. Lumsdaine   */
/*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*/
/*                                                                           */
/*                                                                           */
/*                     SparseLib++ : Sparse Matrix Library                   */
/*                                                                           */
/*               National Institute of Standards and Technology              */
/*                        University of Notre Dame                           */
/*              Authors: R. Pozo, K. Remington, A. Lumsdaine                 */
/*                                                                           */
/*                                 NOTICE                                    */
/*                                                                           */
/* Permission to use, copy, modify, and distribute this software and         */
/* its documentation for any purpose and without fee is hereby granted       */
/* provided that the above notice appear in all copies and supporting        */
/* documentation.                                                            */
/*                                                                           */
/* Neither the Institutions (National Institute of Standards and Technology, */
/* University of Notre Dame) nor the Authors make any representations about  */
/* the suitability of this software for any purpose.  This software is       */
/* provided ``as is'' without expressed or implied warranty.                 */
/*                                                                           */
/*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*/

#ifndef DIAGPRE_H
#define DIAGPRE_H

//#include "vecdefs.h"
#include "vecdefs_GR.h"   // modif GR
//#include VECTOR_H

#include "comprow_double.h"
#include "compcol_double.h"

class DiagPreconditioner_double {

 private:
  VECTOR_double diag_;

 public:
  DiagPreconditioner_double (const CompCol_Mat_double &);
  DiagPreconditioner_double (const CompRow_Mat_double &);
  ~DiagPreconditioner_double (void) { };
  VECTOR_double solve (const VECTOR_double &x) const;
  VECTOR_double trans_solve (const VECTOR_double &x) const;
  
 protected : // modif GR : les deux fonctions suivantes ne sont pas gˇnˇriques
             // d'o¯ le status protected 
  const double&         diag(int i) const { return diag_(i); }
  double&           diag(int i) { return diag_(i); }
};

#endif
