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
00029 #ifndef _chemistry_qc_mbpt_csgrad34qb_h
00030 #define _chemistry_qc_mbpt_csgrad34qb_h
00031
00032 #ifdef __GNUC__
00033 #pragma interface
00034 #endif
00035
00036 #include <util/misc/regtime.h>
00037 #include <util/group/memory.h>
00038 #include <util/group/thread.h>
00039 #include <chemistry/qc/basis/integral.h>
00040
00041 class CSGrad34Qbtr: public Thread {
00042 private:
00043 Ref<MessageGrp> msg;
00044 Ref<MemoryGrp> mem;
00045 Ref<TwoBodyInt> tbint;
00046 Ref<TwoBodyDerivInt> tbintder;
00047 Ref<GaussianBasisSet> basis;
00048 Ref<ThreadLock> lock;
00049 Ref<RegionTimer> timer;
00050 int mythread;
00051 int nthread;
00052 int ni;
00053 int nocc;
00054 int nfzc;
00055 int i_offset;
00056 int aointder_computed;
00057 int me;
00058 int nproc;
00059 double tol;
00060 double **scf_vector;
00061 int debug;
00062 int dynamic_;
00063 int dograd;
00064 int natom;
00065 double *Lpi;
00066 double **ginter;
00067 public:
00068 CSGrad34Qbtr(int mythread_a, int nthread_a,
00069 int me_a, int nproc_a,
00070 const Ref<MemoryGrp> &mem_a,
00071 const Ref<MessageGrp> &msg_a,
00072 const Ref<ThreadLock> &lock_a,
00073 const Ref<GaussianBasisSet> &basis_a,
00074 const Ref<TwoBodyInt> &tbint_a,
00075 const Ref<TwoBodyDerivInt> &tbintder_a,
00076 int nocc_a, int nfzc_a,
00077 double **scf_vector_a,
00078 double tol_a, int debug_a,
00079 int dynamic_a, int dograd_a, int natom_a);
00080 ~CSGrad34Qbtr();
00081
00082 void set_i_offset(int ioff) { i_offset = ioff; }
00083 void set_ni(int nivalue) { ni = nivalue; }
00084 void run();
00085 double *get_Lpi() { return Lpi ;}
00086 double **get_ginter() { return ginter ;}
00087 int get_aointder_computed() { return aointder_computed ;}
00088 };
00089
00090 #endif
00091
00092
00093
00094
00095
00096
00097