00001 #ifndef __F_BiCgStab_2LS_H__ 00002 #define __F_BiCgStab_2LS_H__ 00003 00034 #include <stdio.h> 00035 #include "F_IterativeLS.hpp" 00036 00037 00038 class F_BiCgStab_2LS : public F_IterativeLS 00039 { 00040 private: 00041 00042 double alpha; 00043 double beta; 00044 double gamma; 00045 00046 double rho0, rho1; 00047 double omega1, omega2, omega3; 00048 00049 FloatVector* r; // = ri 00050 FloatVector* u; 00051 FloatVector* v; 00052 FloatVector* w; 00053 00054 FloatVector* r0bar; 00055 FloatVector* s; 00056 FloatVector* t; 00057 00058 00059 protected: 00060 00061 virtual FloatVector* singleSolve(AFloatMatrix* A, FloatVector* b, FloatVector* x=0); 00062 00063 public: 00064 F_BiCgStab_2LS(int maxIterations, double epsilon=1.0e-4):F_IterativeLS(maxIterations, epsilon) 00065 { r=u=v=w= r0bar=s=t= 0;} 00066 virtual ~F_BiCgStab_2LS(){} 00067 00068 virtual FloatVector* solve(AFloatMatrix* m, FloatVector* b, FloatVector* dest); 00069 }; 00070 00071 #endif 00072 00073 00074 00075