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 #if defined(__GNUC__)
00029 #pragma interface
00030 #endif
00031
00032 #ifndef _chemistry_qc_intv3_tranform_h
00033 #define _chemistry_qc_intv3_tranform_h
00034
00035 #include <chemistry/qc/basis/gaussshell.h>
00036 #include <chemistry/qc/basis/transform.h>
00037 #include <chemistry/qc/intv3/macros.h>
00038 #include <chemistry/qc/intv3/int1e.h>
00039 #include <chemistry/qc/intv3/int2e.h>
00040
00041 class Integral;
00042
00043 class SphericalTransformComponentV3 : public SphericalTransformComponent {
00044 public:
00045 void init(int a, int b, int c, double coef, int pureindex) {
00046 a_ = a;
00047 b_ = b;
00048 c_ = c;
00049 coef_ = coef;
00050 pureindex_ = pureindex;
00051 cartindex_ = INT_CARTINDEX(a+b+c,a,b);
00052 }
00053 };
00054
00055 class SphericalTransformV3 : public SphericalTransform {
00056 public:
00057 SphericalTransformV3(int l, int subl=-1):SphericalTransform(l,subl) {
00058 init();
00059 }
00060
00061 SphericalTransformComponent * new_components() {
00062 return new SphericalTransformComponentV3[n_+1];
00063 }
00064 };
00065
00066 class ISphericalTransformV3 : public ISphericalTransform {
00067 public:
00068 ISphericalTransformV3(int l, int subl=-1):ISphericalTransform(l,subl) {
00069 init();
00070 }
00071
00072 SphericalTransformComponent * new_components() {
00073 return new SphericalTransformComponentV3[n_+1];
00074 }
00075 };
00076
00077 #endif
00078
00079
00080
00081
00082