00001 /*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*/ 00002 /* ******** *** SparseLib++ */ 00003 /* ******* ** *** *** *** */ 00004 /* ***** *** ******** ******** */ 00005 /* ***** *** ******** ******** R. Pozo */ 00006 /* ** ******* *** ** *** *** K. Remington */ 00007 /* ******** ******** A. Lumsdaine */ 00008 /*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*/ 00009 /* */ 00010 /* */ 00011 /* SparseLib++ : Sparse Matrix Library */ 00012 /* */ 00013 /* National Institute of Standards and Technology */ 00014 /* University of Notre Dame */ 00015 /* Authors: R. Pozo, K. Remington, A. Lumsdaine */ 00016 /* */ 00017 /* NOTICE */ 00018 /* */ 00019 /* Permission to use, copy, modify, and distribute this software and */ 00020 /* its documentation for any purpose and without fee is hereby granted */ 00021 /* provided that the above notice appear in all copies and supporting */ 00022 /* documentation. */ 00023 /* */ 00024 /* Neither the Institutions (National Institute of Standards and Technology, */ 00025 /* University of Notre Dame) nor the Authors make any representations about */ 00026 /* the suitability of this software for any purpose. This software is */ 00027 /* provided ``as is'' without expressed or implied warranty. */ 00028 /* */ 00029 /*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*/ 00030 00031 #ifndef DIAGPRE_H 00032 #define DIAGPRE_H 00033 00034 #include "vecdefs.h" 00035 #include VECTOR_H 00036 00037 #include "comprow_double.h" 00038 #include "compcol_double.h" 00039 00040 class DiagPreconditioner_double { 00041 00042 private: 00043 VECTOR_double diag_; 00044 00045 public: 00046 DiagPreconditioner_double (const CompCol_Mat_double &); 00047 DiagPreconditioner_double (const CompRow_Mat_double &); 00048 ~DiagPreconditioner_double (void) { }; 00049 VECTOR_double solve (const VECTOR_double &x) const; 00050 VECTOR_double trans_solve (const VECTOR_double &x) const; 00051 00052 const double& diag(int i) const { return diag_(i); } 00053 double& diag(int i) { return diag_(i); } 00054 }; 00055 00056 #endif