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 #ifdef __GNUG__
00029 #pragma interface
00030 #endif
00031
00032 #ifndef _chemistry_qc_cints_primpairs_h
00033 #define _chemistry_qc_cints_primpairs_h
00034
00035 #include <util/ref/ref.h>
00036 #include <chemistry/qc/basis/basis.h>
00037
00038 namespace sc {
00039
00040 typedef struct {
00041 double P[3];
00042 double gamma;
00043 double ovlp;
00044 } prim_pair_t;
00045
00047 class PrimPairsCints : public RefCount {
00048 Ref<GaussianBasisSet> bs1_;
00049 Ref<GaussianBasisSet> bs2_;
00050 int nprim1_;
00051 int nprim2_;
00052 prim_pair_t *prim_pair_;
00053
00054 public:
00055 PrimPairsCints(const Ref<GaussianBasisSet>&,
00056 const Ref<GaussianBasisSet>&);
00057 ~PrimPairsCints();
00058
00059 prim_pair_t* prim_pair(int p1, int p2) const { return prim_pair_ + p1*nprim2_ + p2; };
00060 double P(int p1, int p2, int xyz) const { return prim_pair_[p1*nprim2_ + p2].P[xyz]; };
00061 double gamma(int p1, int p2) const { return prim_pair_[p1*nprim2_ + p2].gamma; };
00062 double ovlp(int p1, int p2) const { return prim_pair_[p1*nprim2_ + p2].ovlp; };
00063
00064 };
00065
00066 }
00067
00068 #endif
00069
00070
00071
00072
00073