/Users/jab/Documents/Teaching/02585/GEL2_and_demos/GEL/src/Geometry/AncestorGrid.h

#ifndef __GEOMETRY_ANCESTORGRID_H
#define __GEOMETRY_ANCESTORGRID_H

#include "CGLA/Vec3i.h"

namespace Geometry 
{

        template<typename T, class ChildT>
        class AncestorGrid
        {
        public:
                typedef T DataType;

        private:
                CGLA::Vec3i dims;
  
        public:
  
                AncestorGrid(int _x_dim, int _y_dim, int _z_dim): 
                        dims(_x_dim,_y_dim,_z_dim) {}
  
                AncestorGrid(const CGLA::Vec3i& _dims): dims(_dims) {}
                
                bool in_domain(const CGLA::Vec3i& p) const
                {
                        for(int i=0; i<3; i++)
                                if (p[i]<0 || p[i] >= dims[i])
                                        return false;
                        return true;
                }
  
                const CGLA::Vec3i& get_dims() const {return dims;}

                const CGLA::Vec3i get_lo_corner() const {return CGLA::Vec3i(0);}

                const CGLA::Vec3i& get_hi_corner() const {return dims;}

                const T& operator[](const CGLA::Vec3i& p) const 
                {
                        return static_cast<const ChildT&>(*this)[p];
                }


                void store(const CGLA::Vec3i& p, const T& t)
                {
                        return static_cast<ChildT&>(*this).store(p,t);
                }
        };
}

#endif