ClosePackIterator.h

Go to the documentation of this file.

//                                                         //
// Copyright (c) 2003-2014 by The University of Queensland //
// Centre for Geoscience Computing                         //
// http://earth.uq.edu.au/centre-geoscience-computing      //
//                                                         //
// Primary Business: Brisbane, Queensland, Australia       //
// Licensed under the Open Software License version 3.0    //
// http://www.opensource.org/licenses/osl-3.0.php          //
//                                                         //


#ifndef ESYS_LSMCLOSEPACKITERATOR_H
#define ESYS_LSMCLOSEPACKITERATOR_H

#include "Foundation/BoundingBox.h"
#include "Foundation/vec3.h"
#include "Geometry/Vec3L.h"
#include "Geometry/ClosePackOrientation.h"

namespace esys
{
  namespace lsm
  {
    template <int NI, int NJ, int NK>
    class TmplMatrix
    {
    public:
      TmplMatrix();

      TmplMatrix(const TmplMatrix &m);

      TmplMatrix &operator=(const TmplMatrix &m);

      const double &operator()(int i, int j, int k) const;

      double &operator()(int i, int j, int k);

      int getNumI() const;

      int getNumJ() const;

      int getNumK() const;

    private:
      double m_matrix[NI][NJ][NK];
    };

    class ClosePackIterator
    {
    public:
      static const double SQRT_1_OVER_3;
      static const double SQRT_8_OVER_3;
      static const double SQRT_3;

      inline ClosePackIterator();

      inline ClosePackIterator(
        int numI,
        int numJ,
        int numK,
        double sphereRadius,
        ClosePackOrientation orientation = DEFAULT_ORIENT
      );

      inline bool hasNext() const;

      inline Vec3 next();

      inline double getRadius() const;

    protected:
      typedef TmplMatrix<3,6,6> OffsetMatrix;

      inline void incrementDimIndex();

      inline double getOffset(int i) const;

      inline const Vec3 &getMinPt() const;

      inline void setMinPt(const Vec3 &pt) const;

      inline void setDimRepeat(const Vec3L &dimRepeat);

      inline void setOffsetMatrix(const OffsetMatrix &offsetMatrix);

    private:
      static  Vec3L s_orientationDimMap[NUM_ORIENTATIONS];
      double        m_radius;
      Vec3          m_minPt;
      OffsetMatrix  m_offsetMatrix;
      Vec3L         m_dimRepeat;
      Vec3L         m_dimCount;
      Vec3L         m_dimIdx;
      Vec3L         m_dim;
    };
  }
}

#include "Geometry/ClosePackIterator.hpp"

#endif