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_csgrads2pdm_h
00030 #define _chemistry_qc_mbpt_csgrads2pdm_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 namespace sc {
00042
00043 #define PRINT_BIGGEST_INTS 0
00044
00045 class CSGradS2PDM: public Thread {
00046 private:
00047 int mythread;
00048 int nthread;
00049 int me;
00050 int nproc;
00051 Ref<ThreadLock> lock;
00052 Ref<GaussianBasisSet> basis;
00053 Ref<TwoBodyDerivInt> tbintder;
00054 const double *PHF;
00055 const double *P2AO;
00056 int tol;
00057 int debug;
00058 int dynamic;
00059
00060 double **ginter;
00061 double **hf_ginter;
00062
00063 void accum_contrib(double **sum, double **contribs);
00064 public:
00065 CSGradS2PDM(int mythread_a, int nthread_a,
00066 int me_a, int nproc_a,
00067 const Ref<ThreadLock> &lock_a,
00068 const Ref<GaussianBasisSet> &basis_a,
00069 const Ref<TwoBodyDerivInt> &tbintder_a,
00070 const double *PHF_a, const double *P2AO_a,
00071 int tol_a, int debug_a, int dynamic_a);
00072
00073 ~CSGradS2PDM();
00074 void accum_mp2_contrib(double **ginter);
00075 void accum_hf_contrib(double **hf_ginter);
00076 void run();
00077 };
00078
00079 }
00080
00081 #endif
00082
00083
00084
00085
00086
00087
00088