source: GTP/trunk/Lib/Vis/Preprocessing/src/Halton.cpp @ 2575

#include "Halton.h"

#ifdef PERFTIMER  
#include "Timer/PerfTimer.h"
#endif

#define PREGENERATE_HALTON 0

namespace GtpVisibilityPreprocessor {

  
#ifdef PERFTIMER  
PerfTimer haltonTimer;
#endif
  
// global halton generator

  Halton2 halton2;
  float Halton2::_invBases[2];

  //float Halton<6>::_invBases[6];
  
  template<int T> float Halton<T>::_invBases[T];
  
  Halton<1> dummmyHalton1(true);
  Halton<2> dummmyHalton2(true);
  Halton<3> dummmyHalton3(true);
  Halton<4> dummmyHalton4(true);
  Halton<5> dummmyHalton5(true);
  Halton<6> dummmyHalton6(true);


  int HaltonSequence::sPregeneratedDim = 0;
  int HaltonSequence::sPregeneratedNumber = 0;
  float *HaltonSequence::sPregeneratedValues = NULL;
  
#if PREGENERATE_HALTON
  // pregenerate first 100k numbers for 6 dimensions
  static HaltonSequence staticHalton(6, 200000);
#endif
  
  // special construtor for pregenerating static halton sequences
  HaltonSequence::HaltonSequence(const int dim,
                                                                 const int number) {

        sPregeneratedDim = 0;
        sPregeneratedNumber = 0;
        sPregeneratedValues = new float[number*dim];

        float *p = sPregeneratedValues;
        int i;
        for (i=0; i < number; i++, p+=dim) {
          GetNext(dim, p);
        }
        
        sPregeneratedDim = dim;
        sPregeneratedNumber = number;
  }

void
HaltonSequence::GetNext(const int dimensions, float *p)
{
  //  haltonTimer.Entry();
#if PREGENERATE_HALTON
  if (index <= sPregeneratedNumber && dimensions <= sPregeneratedDim) {
        float *pp = sPregeneratedValues + (sPregeneratedDim*(index-1));
        for (int i=0; i < dimensions; i++)
          p[i] = pp[i];
  } else
#endif
        {
          for (int i=0; i < dimensions; i++)
                p[i] = (float)GetNumber(i+1);
        }
  GenerateNext();
  //  haltonTimer.Exit();
}

}