00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #ifndef _math_optimize_gdiis_h
00029 #define _math_optimize_gdiis_h
00030
00031 #ifdef __GNUC__
00032 #pragma interface
00033 #endif
00034
00035 #include <util/state/state.h>
00036 #include <util/class/class.h>
00037 #include <math/scmat/matrix.h>
00038 #include <math/optimize/function.h>
00039 #include <math/optimize/opt.h>
00040 #include <math/optimize/update.h>
00041
00042 namespace sc {
00043
00044
00045
00046
00047 class GDIISOpt: public Optimize {
00048 protected:
00049 int nsave;
00050 int diis_iter;
00051
00052 double maxabs_gradient;
00053 double convergence_;
00054 double accuracy_;
00055
00056 RefSCVector *coords_;
00057 RefSCVector *grad_;
00058 RefSCVector *error_;
00059
00060 RefSymmSCMatrix ihessian_;
00061 Ref<HessianUpdate> update_;
00062
00063 public:
00064 GDIISOpt(const Ref<KeyVal>&);
00065 GDIISOpt(StateIn&);
00066 ~GDIISOpt();
00067 void save_data_state(StateOut&);
00068
00069 void init();
00070 int update();
00071 };
00072
00073 }
00074
00075 #endif
00076
00077
00078
00079
00080