MPQC_IntegralEvaluator4_Impl.hh

// 
// File:          MPQC_IntegralEvaluator4_Impl.hh
// Symbol:        MPQC.IntegralEvaluator4-v0.2
// Symbol Type:   class
// Babel Version: 0.10.2
// Description:   Server-side implementation for MPQC.IntegralEvaluator4
// 
// WARNING: Automatically generated; only changes within splicers preserved
// 
// babel-version = 0.10.2
// 
 
#ifndef included_MPQC_IntegralEvaluator4_Impl_hh
#define included_MPQC_IntegralEvaluator4_Impl_hh
 
#ifndef included_sidl_cxx_hh
#include "sidl_cxx.hh"
#endif
#ifndef included_MPQC_IntegralEvaluator4_IOR_h
#include "MPQC_IntegralEvaluator4_IOR.h"
#endif
// 
// Includes for all method dependencies.
// 
#ifndef included_Chemistry_QC_GaussianBasis_DerivCenters_hh
#include "Chemistry_QC_GaussianBasis_DerivCenters.hh"
#endif
#ifndef included_Chemistry_QC_GaussianBasis_Molecular_hh
#include "Chemistry_QC_GaussianBasis_Molecular.hh"
#endif
#ifndef included_MPQC_IntegralEvaluator4_hh
#include "MPQC_IntegralEvaluator4.hh"
#endif
#ifndef included_sidl_BaseInterface_hh
#include "sidl_BaseInterface.hh"
#endif
#ifndef included_sidl_ClassInfo_hh
#include "sidl_ClassInfo.hh"
#endif
 
 
// DO-NOT-DELETE splicer.begin(MPQC.IntegralEvaluator4._includes)
#include <chemistry/qc/basis/gaussbas.h>
#include <chemistry/qc/intv3/cartitv3.h>
#include <chemistry/qc/intv3/intv3.h>
#ifdef HAVE_CINTS
  #include <chemistry/qc/cints/cints.h>
#endif
using namespace sc;
// DO-NOT-DELETE splicer.end(MPQC.IntegralEvaluator4._includes)
 
namespace MPQC { 
 
  class IntegralEvaluator4_impl
  // DO-NOT-DELETE splicer.begin(MPQC.IntegralEvaluator4._inherits)
 
  // Put additional inheritance here...
  // DO-NOT-DELETE splicer.end(MPQC.IntegralEvaluator4._inherits)
  {
 
  private:
    // Pointer back to IOR.
    // Use this to dispatch back through IOR vtable.
    IntegralEvaluator4 self;
 
    // DO-NOT-DELETE splicer.begin(MPQC.IntegralEvaluator4._implementation)
    Chemistry::Molecule molecule_;
    std::string evaluator_label_;
    Ref<sc::GaussianBasisSet> bs1_;
    Ref<sc::GaussianBasisSet> bs2_;
    Ref<sc::GaussianBasisSet> bs3_;
    Ref<sc::GaussianBasisSet> bs4_;
    Ref<sc::Integral> integral_;
    Ref<sc::TwoBodyInt> eval_;
    Ref<sc::TwoBodyDerivInt> deriv_eval_;
    int max_nshell4_;
    sidl::array<double> sidl_buffer_;
    const double* sc_buffer_;
    enum { two_body, two_body_deriv};
    int int_type_;
    int deriv_level_;
    std::string package_;
 
    // reorder stuff
    int bufn_;
    int **reorder_;
    double *buf_;
    double *temp_buffer_;
    int index_, con2_offset_, con3_offset_, con4_offset_, con_offset_,
        local2_offset_, local3_offset_, local4_offset_,
        c1_base_, c2_base_, c3_base_, c4_base_,
        s1_is_cart_, s2_is_cart_, s3_is_cart_, s4_is_cart_,
        s1_nfunc_, s2_nfunc_, s3_nfunc_, s4_nfunc_,
        nc1_, nc2_, nc3_, nc4_;
    sc::GaussianShell *s1_, *s2_, *s3_, *s4_;
 
    void reorder_intv3(int64_t,int64_t,int64_t,int64_t);
    void reorder_quartet( sc::GaussianShell*, sc::GaussianShell*,
                          sc::GaussianShell*, sc::GaussianShell*,
                          int, int, int, int, int );
    void reorder_intv3_inline(int64_t,int64_t,int64_t,int64_t);
    void initialize_reorder_intv3();
    void reorder_c4(int,int,int,int,int,int);
    // DO-NOT-DELETE splicer.end(MPQC.IntegralEvaluator4._implementation)
 
  private:
    // private default constructor (required)
    IntegralEvaluator4_impl() 
    {} 
 
  public:
    // sidl constructor (required)
    // Note: alternate Skel constructor doesn't call addref()
    // (fixes bug #275)
    IntegralEvaluator4_impl( struct MPQC_IntegralEvaluator4__object * s ) : 
      self(s,true) { _ctor(); }
 
    // user defined construction
    void _ctor();
 
    // virtual destructor (required)
    virtual ~IntegralEvaluator4_impl() { _dtor(); }
 
    // user defined destruction
    void _dtor();
 
    // static class initializer
    static void _load();
 
  public:
 
    void
    set_integral_package (
      /* in */ const ::std::string& label
    )
    throw () 
    ;
 
 
    void
    initialize (
      /* in */ ::Chemistry::QC::GaussianBasis::Molecular bs1,
      /* in */ ::Chemistry::QC::GaussianBasis::Molecular bs2,
      /* in */ ::Chemistry::QC::GaussianBasis::Molecular bs3,
      /* in */ ::Chemistry::QC::GaussianBasis::Molecular bs4,
      /* in */ const ::std::string& label,
      /* in */ int64_t max_deriv
    )
    throw () 
    ;
 
 
    void*
    get_buffer() throw () 
    ;
 
    void
    compute (
      /* in */ int64_t shellnum1,
      /* in */ int64_t shellnum2,
      /* in */ int64_t shellnum3,
      /* in */ int64_t shellnum4,
      /* in */ int64_t deriv_level,
      /* in */ ::Chemistry::QC::GaussianBasis::DerivCenters deriv_ctr
    )
    throw () 
    ;
 
 
    ::sidl::array<double>
    compute_array (
      /* in */ int64_t shellnum1,
      /* in */ int64_t shellnum2,
      /* in */ int64_t shellnum3,
      /* in */ int64_t shellnum4,
      /* in */ int64_t deriv_level,
      /* in */ ::Chemistry::QC::GaussianBasis::DerivCenters deriv_ctr
    )
    throw () 
    ;
 
  };  // end class IntegralEvaluator4_impl
 
} // end namespace MPQC
 
// DO-NOT-DELETE splicer.begin(MPQC.IntegralEvaluator4._misc)
// Put miscellaneous things here...
// DO-NOT-DELETE splicer.end(MPQC.IntegralEvaluator4._misc)
 
#endif