Changeset 2105 for GTP/trunk/Lib/Vis/Preprocessing/src/Halton.cpp

Timestamp:

02/08/07 15:46:23 (17 years ago)

Author:

bittner

Message:

simple ray separated

File:

• GTP/trunk/Lib/Vis/Preprocessing/src/Halton.cpp (modified) (2 diffs)

GTP/trunk/Lib/Vis/Preprocessing/src/Halton.cpp

@@ -1 +1 @@
 #include "Halton.h"
 #include "PerfTimer.h"
+
+#define PREGENERATE_HALTON 0
 
 namespace GtpVisibilityPreprocessor {
@@ -25 +27 @@
   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();
-  for (int i=0; i < dimensions; i++)
-        p[i] = (float)GetNumber(i+1);
+  //  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();
+  //  haltonTimer.Exit();
 }