00001 #ifndef __StdTikhonovPhillipsSolver_H__
00002 #define __StdTikhonovPhillipsSolver_H__
00003 
00004 
00005 
00033 #include "F_CgLS_ShiftedLS.hpp"
00034 #include "ATikhonovPhillipsSolver.hpp"
00035 
00036 class StdTikhonovPhillipsSolver : public ATikhonovPhillipsSolver
00037 {
00038 protected:
00039 
00040 F_IterativeShiftedLS* linearSystem;
00041 
00042 public:
00043 StdTikhonovPhillipsSolver(int iter= 10, double adelta=0.0001):ATikhonovPhillipsSolver(10, adelta)
00044         { linearSystem= new F_CgLS_ShiftedLS(10, alpha0, adelta); }
00045 virtual ~StdTikhonovPhillipsSolver(){}  
00046 
00047 
00048 virtual FloatVector* solve(AFloatMatrix*        m, FloatVector* b, FloatVector* dest);
00049 virtual FloatVector* solve(AFloatMatrix*        m, FloatVector* b, bool factorisation=true, FloatVector* dest=0)
00050         { return solve(m, b, dest); }
00051 };
00052 
00053 #endif
00054 
00055