---

extent.h

//
// extent.h
//

#ifndef _chemistry_qc_basis_extent_h
#define _chemistry_qc_basis_extent_h

#ifdef __GNUC__
#pragma interface
#endif

#include <float.h>
#include <chemistry/qc/basis/basis.h>

struct ExtentData {
    int shell;
    double bound;
    ExtentData() {}
    ExtentData(int s, double b): shell(s), bound(b) {}
};

ARRAY_dec(ExtentData);

class ShellExtent: public RefCount {
    double lower_[3];
    double resolution_;
    int n_[3];
    ArrayExtentData *contributing_shells_;
    ArrayExtentData null_;

    ArrayExtentData &data(int *b);
    double distance(double loc, int axis, int origin, int point);
    ArrayExtentData &data(int x, int y, int z);
  public:
    ShellExtent();
    ~ShellExtent();
    void init(const Ref<GaussianBasisSet>&,
              double resolution = 1.0, double tolerance = DBL_EPSILON);
    const ArrayExtentData &contributing_shells(int x, int y, int z)
        { return data(x,y,z); }
    const ArrayExtentData &contributing_shells(double x, double y, double z);
    void print(std::ostream &o = ExEnv::out());
    const int *n() const { return n_; }
    int n(int ixyz) const { return n_[ixyz]; }
    double lower(int ixyz) const { return lower_[ixyz]; }
    double resolution() const { return resolution_; }
};


#endif

// Local Variables:
// mode: c++
// c-file-style: "CLJ"
// End:

---