simpleilu.h

Go to the documentation of this file.
00001 #ifndef  __SimpleIlu_h
00002 #define  __SimpleIlu_h
00003 
00004 #include  "simplematrix.h"
00005 #include  "matrixinterface.h"
00006 
00007 
00008 namespace Gascoigne
00009 {
00010 
00015 
00019 
00020 class SimpleIlu : public SimpleMatrix
00021 {
00022 protected:
00023 
00024   IntVector              p,q;
00025   mutable DoubleVector   yp;
00026 
00027   void hin(const DoubleVector& y) const;
00028   void her(DoubleVector& y) const;
00029   void backward() const;
00030   void forward () const;
00031   void backward_transpose() const;
00032   void forward_transpose () const;
00033 
00034 public:
00035 
00036 //
00038 //
00039 
00040     SimpleIlu() : SimpleMatrix() {}
00041     
00042       int   n()          const { return GetStencil()->n();};
00043     void zero() {SimpleMatrix::zero();}
00044     void ReInit(int n, int nentries);
00045     void copy_entries(const MatrixInterface*  A);
00046     void compute_ilu();
00047     void solve(DoubleVector& x) const;
00048     void solve_transpose(DoubleVector& x) const;
00049 };
00050 }
00051 
00052 #endif

Generated on Thu Sep 14 10:34:37 2006 for Gascoigne by  doxygen 1.4.7