Cube.hpp

Go to the documentation of this file.

#include <vector>

//______________________________________________________________________________
template<typename T>  // See section on templates for more
class Cube {
public:
  Cube(unsigned nx, unsigned ny, unsigned nz);
  // Throws a BadSize object if either size is zero
  class BadSize { };

  void resize(unsigned nx, unsigned ny, unsigned nz);
 
  // No need for any of The Big Three!
 
  // Access methods to get the (i,j) element:
  T&       operator() (unsigned i, unsigned j, unsigned k);        // subscript operators often come in pairs
  const T& operator() (unsigned i, unsigned j, unsigned k) const;  // subscript operators often come in pairs
    // These throw a BoundsViolation object if i or j is too big
  class BoundsViolation { };
 
  unsigned nx() const;  // # x-row elements in this matrix
  unsigned ny() const;  // # y-row elements columns in this matrix
  unsigned nz() const;  // # z-row elements in this matrix
 
private:
  // note: need to put > > > instead of >>>: otherwise the compiler thinks we have the >> operator
  std::vector<std::vector<std::vector<T> > > _data;
};
 
template<typename T>
inline unsigned Cube<T>::nx() const
{ return _data.size(); }
 
template<typename T>
inline unsigned Cube<T>::ny() const
{ return _data[0].size(); }

template<typename T>
inline unsigned Cube<T>::nz() const
{ return _data[0][0].size(); }
 
template<typename T>
inline T& Cube<T>::operator() (unsigned i, unsigned j, unsigned k)
{
  if (i >= nx() || j >= ny() || k >= nz()) throw BoundsViolation();
  return _data[i][j][k];
}
 
template<typename T>
inline const T& Cube<T>::operator() (unsigned i, unsigned j, unsigned k) const
{
  if (i >= nx() || j >= ny() || k >= nz()) throw BoundsViolation();
  return _data[i][j][k];
}

//______________________________________________________________________________
template<typename T>
void
Cube<T>::resize(unsigned nx, unsigned ny, unsigned nz)
{
  if (nx == 0 || ny == 0 || nz == 0)
    throw BadSize();
  for (unsigned i = 0; i < nx; ++i){    
    _data[i].resize(ny);
    for (unsigned j = 0; j < ny; ++j)
      _data[i][j].resize(nz);  
  }
}

//______________________________________________________________________________
template<typename T>
Cube<T>::Cube(unsigned nx, unsigned ny, unsigned nz)
  : _data (nx)
{
  resize(nx, ny, nz);
} 

---