PeanoHilbert.h

Go to the documentation of this file.

// -*- C++ -*-

#ifndef _included_PeanoHilbert_h
#define _included_PeanoHilbert_h

#include "sfcIndex.h"

#include <cassert>

#ifndef PeanoHibertBase
#define PeanoHibertBase (2)
#endif

class PeanoHilbert : public sfcIndex
  {

public:

   inline PeanoHilbert()
        : sfcIndex() {}

   inline PeanoHilbert(const int dim, const int levs)
        : sfcIndex(dim,levs,PeanoHibertBase)
        {assert(dim<=3);}

   inline PeanoHilbert(PeanoHilbert const &other)
        : sfcIndex((sfcIndex const &)other) {}

private:
   void Invert(int *coords) const;

public:
   inline void SetIndex(const int index, const int idxlev)
     { sfcSetIndex(index,idxlev); }
   inline void SetIndex(const int index, const int idxlev, 
                        const int nlev)
     { sfcSetIndex(index,idxlev,nlev); }

   void Map(const int *coords);
   void Invert(int *coords, const int index);
   void Invert(int *coords, const int index, const int idxlev);
   void Invert(int *coords, const int index, 
               const int idxlev, const int nlev);

   PeanoHilbert *GetBox(const int lev, const int nlev) const;
   void SetBox(const int lev, const int nlev);

   PeanoHilbert *GetMax(const int lev, const int nlev) const;
   void SetMax(const int lev, const int nlev);

   void ResetBase(const int lev, const int nlev);
   
   /* Query routines */
   void GetCoords(int *coords, const int lev = 1) const;
   int *GetCoords(const int lev = 1) const;

   inline int GetCardinality(const int dim, const int lev) const
     {return (int)1<<(lev*dim);}
   inline int GetCardinality(const int lev) const
     {return (int)1<<(lev*sfcdim);}

   inline int GetDimCardinality(const int lev) const
     {return (int)1<<(sfclevs-lev);}

   inline int GetDimMax(const int lev) const
     {return (((int)1<<lev)-1)*((int)1<<(sfclevs-lev));}
   
private:

/*
// The following operations are used to orient the seed
// curve as required. The seed curve is "->"
// Axis are numbered as follows:
// "->" x: 1
// "|"  y: 2
// "/"  z: 3
// Directions are numbered as follows:
// "->" (positive x): 1
// "|" (positive y): 2
// "/" (positive z): 3
// phSwapDir ->(PH) = <-(PH)
// phTranslate ->(PH) = (PH)->
// phRotate ->(PH) (1,2) = |(PH)
// phRotate ->(PH) (1,3) = /(PH)
// phRotate |(PH) (2,1) = ->(PH)
// phRotate /(PH) (3,1) = ->(PH)
// Note:
// (1) "phRotate" is valid to/from the seed orientation only !
//     i.e. "to" or "from" should be "->" or 1.
// (2) Bit are numbered from 0 to DIM-1 such that the msb is 0 and
//     the lsb is Dim-1.
//     i.e. x = xb0 xb1 xb2 xb3 ... xbDIM-1
*/

   inline void phSwapDir(const int level, const int axis)
     {sfcToggleDigitWd(level,axis-1);}
   inline void phTranslate(const int level, const int axis)
     {sfcToggleDigitWd(level,axis-1);}
   inline void phRotate(const int level, const int from_axis, 
                        const int to_axis)
     {sfcSwapDigitWd(level,from_axis-1,to_axis-1);}

   void phOrient(const int level, const int pos);
   int phGray(const int level);
   int phGinv(const int level);

  };

#endif

---