00001 #ifndef __F_MRII_LS_H__ 00002 #define __F_MRII_LS_H__ 00003 00004 00034 #include <stdio.h> 00035 #include "F_IterativeLS.hpp" 00036 00037 00038 class F_MRII_LS : public F_IterativeLS 00039 { 00040 private: 00041 00042 double alpha; 00043 double beta; 00044 double rho; 00045 00046 FloatVector* r; 00047 00048 FloatVector* v; 00049 FloatVector* v_1; 00050 00051 FloatVector* w; 00052 FloatVector* w_1; 00053 FloatVector* sw; 00054 00055 FloatVector* tmp; 00056 00057 protected: 00058 00059 virtual FloatVector* singleSolve(AFloatMatrix* A, FloatVector* b, FloatVector* x=0); 00060 00061 public: 00062 F_MRII_LS(int maxIterations, double epsilon=1.0e-4):F_IterativeLS(maxIterations, epsilon) 00063 { r=v=v_1=w=w_1=sw= 0;} 00064 virtual ~F_MRII_LS(){} 00065 00066 virtual FloatVector* solve(AFloatMatrix* m, FloatVector* b, FloatVector* dest); 00067 }; 00068 00069 #endif 00070 00071 00072
|
