[355] | 1 | #include "Halton.h"
|
---|
[2575] | 2 |
|
---|
| 3 | #ifdef PERFTIMER
|
---|
[2227] | 4 | #include "Timer/PerfTimer.h"
|
---|
[2575] | 5 | #endif
|
---|
[355] | 6 |
|
---|
[2105] | 7 | #define PREGENERATE_HALTON 0
|
---|
| 8 |
|
---|
[863] | 9 | namespace GtpVisibilityPreprocessor {
|
---|
[860] | 10 |
|
---|
[2575] | 11 |
|
---|
| 12 | #ifdef PERFTIMER
|
---|
[2076] | 13 | PerfTimer haltonTimer;
|
---|
[2575] | 14 | #endif
|
---|
| 15 |
|
---|
[355] | 16 | // global halton generator
|
---|
| 17 |
|
---|
[2076] | 18 | Halton2 halton2;
|
---|
| 19 | float Halton2::_invBases[2];
|
---|
[860] | 20 |
|
---|
[2575] | 21 | //float Halton<6>::_invBases[6];
|
---|
[2076] | 22 |
|
---|
[2575] | 23 | template<int T> float Halton<T>::_invBases[T];
|
---|
| 24 |
|
---|
[2076] | 25 | Halton<1> dummmyHalton1(true);
|
---|
| 26 | Halton<2> dummmyHalton2(true);
|
---|
| 27 | Halton<3> dummmyHalton3(true);
|
---|
| 28 | Halton<4> dummmyHalton4(true);
|
---|
| 29 | Halton<5> dummmyHalton5(true);
|
---|
| 30 | Halton<6> dummmyHalton6(true);
|
---|
[2105] | 31 |
|
---|
| 32 |
|
---|
| 33 | int HaltonSequence::sPregeneratedDim = 0;
|
---|
| 34 | int HaltonSequence::sPregeneratedNumber = 0;
|
---|
| 35 | float *HaltonSequence::sPregeneratedValues = NULL;
|
---|
[2076] | 36 |
|
---|
[2105] | 37 | #if PREGENERATE_HALTON
|
---|
| 38 | // pregenerate first 100k numbers for 6 dimensions
|
---|
| 39 | static HaltonSequence staticHalton(6, 200000);
|
---|
| 40 | #endif
|
---|
| 41 |
|
---|
| 42 | // special construtor for pregenerating static halton sequences
|
---|
| 43 | HaltonSequence::HaltonSequence(const int dim,
|
---|
| 44 | const int number) {
|
---|
| 45 |
|
---|
| 46 | sPregeneratedDim = 0;
|
---|
| 47 | sPregeneratedNumber = 0;
|
---|
| 48 | sPregeneratedValues = new float[number*dim];
|
---|
| 49 |
|
---|
| 50 | float *p = sPregeneratedValues;
|
---|
| 51 | int i;
|
---|
| 52 | for (i=0; i < number; i++, p+=dim) {
|
---|
| 53 | GetNext(dim, p);
|
---|
| 54 | }
|
---|
| 55 |
|
---|
| 56 | sPregeneratedDim = dim;
|
---|
| 57 | sPregeneratedNumber = number;
|
---|
| 58 | }
|
---|
| 59 |
|
---|
[2076] | 60 | void
|
---|
| 61 | HaltonSequence::GetNext(const int dimensions, float *p)
|
---|
| 62 | {
|
---|
[2105] | 63 | // haltonTimer.Entry();
|
---|
| 64 | #if PREGENERATE_HALTON
|
---|
| 65 | if (index <= sPregeneratedNumber && dimensions <= sPregeneratedDim) {
|
---|
| 66 | float *pp = sPregeneratedValues + (sPregeneratedDim*(index-1));
|
---|
| 67 | for (int i=0; i < dimensions; i++)
|
---|
| 68 | p[i] = pp[i];
|
---|
| 69 | } else
|
---|
| 70 | #endif
|
---|
| 71 | {
|
---|
| 72 | for (int i=0; i < dimensions; i++)
|
---|
| 73 | p[i] = (float)GetNumber(i+1);
|
---|
| 74 | }
|
---|
[2076] | 75 | GenerateNext();
|
---|
[2105] | 76 | // haltonTimer.Exit();
|
---|
[2076] | 77 | }
|
---|
| 78 |
|
---|
| 79 | }
|
---|