| [0b990d] | 1 | //
 | 
|---|
 | 2 | // int2e.cc
 | 
|---|
 | 3 | //
 | 
|---|
 | 4 | // Copyright (C) 1996 Limit Point Systems, Inc.
 | 
|---|
 | 5 | //
 | 
|---|
 | 6 | // Author: Curtis Janssen <cljanss@limitpt.com>
 | 
|---|
 | 7 | // Maintainer: LPS
 | 
|---|
 | 8 | //
 | 
|---|
 | 9 | // This file is part of the SC Toolkit.
 | 
|---|
 | 10 | //
 | 
|---|
 | 11 | // The SC Toolkit is free software; you can redistribute it and/or modify
 | 
|---|
 | 12 | // it under the terms of the GNU Library General Public License as published by
 | 
|---|
 | 13 | // the Free Software Foundation; either version 2, or (at your option)
 | 
|---|
 | 14 | // any later version.
 | 
|---|
 | 15 | //
 | 
|---|
 | 16 | // The SC Toolkit is distributed in the hope that it will be useful,
 | 
|---|
 | 17 | // but WITHOUT ANY WARRANTY; without even the implied warranty of
 | 
|---|
 | 18 | // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 | 
|---|
 | 19 | // GNU Library General Public License for more details.
 | 
|---|
 | 20 | //
 | 
|---|
 | 21 | // You should have received a copy of the GNU Library General Public License
 | 
|---|
 | 22 | // along with the SC Toolkit; see the file COPYING.LIB.  If not, write to
 | 
|---|
 | 23 | // the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
 | 
|---|
 | 24 | //
 | 
|---|
 | 25 | // The U.S. Government is granted a limited license as per AL 91-7.
 | 
|---|
 | 26 | //
 | 
|---|
 | 27 | 
 | 
|---|
 | 28 | #ifdef __GNUG__
 | 
|---|
 | 29 | #pragma implementation
 | 
|---|
 | 30 | #endif
 | 
|---|
 | 31 | 
 | 
|---|
 | 32 | #include <chemistry/qc/basis/integral.h>
 | 
|---|
 | 33 | #include <chemistry/qc/intv3/int2e.h>
 | 
|---|
 | 34 | #include <chemistry/qc/intv3/utils.h>
 | 
|---|
 | 35 | 
 | 
|---|
 | 36 | using namespace sc;
 | 
|---|
 | 37 | 
 | 
|---|
 | 38 | Int2eV3::Int2eV3(Integral *integral,
 | 
|---|
 | 39 |                  const Ref<GaussianBasisSet>& b1,
 | 
|---|
 | 40 |                  const Ref<GaussianBasisSet>& b2,
 | 
|---|
 | 41 |                  const Ref<GaussianBasisSet>& b3,
 | 
|---|
 | 42 |                  const Ref<GaussianBasisSet>& b4,
 | 
|---|
 | 43 |                  int order, size_t storage) :
 | 
|---|
 | 44 |   integral_(integral),
 | 
|---|
 | 45 |   grp_(integral->messagegrp()),
 | 
|---|
 | 46 |   store(0),
 | 
|---|
 | 47 |   int_Qvec(0),
 | 
|---|
 | 48 |   int_Rvec(0)
 | 
|---|
 | 49 | {
 | 
|---|
 | 50 |   bs1_ = b1;
 | 
|---|
 | 51 |   bs2_ = b2;
 | 
|---|
 | 52 |   bs3_ = b3;
 | 
|---|
 | 53 |   bs4_ = b4;
 | 
|---|
 | 54 | 
 | 
|---|
 | 55 |   int_unit2 = bs2_.null();
 | 
|---|
 | 56 |   int_unit4 = bs4_.null();
 | 
|---|
 | 57 | 
 | 
|---|
 | 58 |   transform_init();
 | 
|---|
 | 59 |   int_initialize_offsets2();
 | 
|---|
 | 60 |   int_initialize_erep(storage,order,bs1_,bs2_,bs3_,bs4_);
 | 
|---|
 | 61 |   if (!(int_unit2 || int_unit4)) {
 | 
|---|
 | 62 |     if (order==0) {
 | 
|---|
 | 63 |       init_bounds();
 | 
|---|
 | 64 |       }
 | 
|---|
 | 65 |     else if (order==1) {
 | 
|---|
 | 66 |       init_bounds_1der();
 | 
|---|
 | 67 |     }
 | 
|---|
 | 68 |   }
 | 
|---|
 | 69 | }
 | 
|---|
 | 70 | 
 | 
|---|
 | 71 | Int2eV3::~Int2eV3()
 | 
|---|
 | 72 | {
 | 
|---|
 | 73 |   transform_done();
 | 
|---|
 | 74 |   int_done_offsets2();
 | 
|---|
 | 75 |   int_done_erep();
 | 
|---|
 | 76 |   if (int_integral_storage) {
 | 
|---|
 | 77 |     done_storage();
 | 
|---|
 | 78 |   }
 | 
|---|
 | 79 |   done_bounds();
 | 
|---|
 | 80 |   done_bounds_1der();
 | 
|---|
 | 81 | }
 | 
|---|
 | 82 | 
 | 
|---|
 | 83 | /////////////////////////////////////////////////////////////////////////////
 | 
|---|
 | 84 | 
 | 
|---|
 | 85 | // Local Variables:
 | 
|---|
 | 86 | // mode: c++
 | 
|---|
 | 87 | // c-file-style: "CLJ-CONDENSED"
 | 
|---|
 | 88 | // End:
 | 
|---|