mat.h

/* ... */

template<size_t Dimen>
class NDMatrix
{
public:
  /* this is imaginary, as of right now: */
  $for(size_t i = 0; i < Dimen; ++i) {
    /* Create one function definition for each dimen access:
     * float &At(int x0, int x1, etc...) { return m_array[_At(...)]; }
     */
  }
  
private:
  size_t 
  _At(int indices[Dimen])
  {
    /* calculate the correct index for access into linear array here . */
    return idx;
  }
};


/* ... */