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