Changeset 3329 for GTP/trunk/App/Demos/Vis/FriendlyCulling/src/SampleGenerator.cpp

Timestamp:

02/26/09 18:39:04 (15 years ago)

Author:

mattausch

Message:

 

File:

• GTP/trunk/App/Demos/Vis/FriendlyCulling/src/SampleGenerator.cpp (modified) (1 diff)

GTP/trunk/App/Demos/Vis/FriendlyCulling/src/SampleGenerator.cpp

@@ -218 +218 @@
 void QuadraticDiscSampleGenerator2D::Generate(float *samples) const
 {
-#if 0
-        float r[2];
-        Sample2 *s = (Sample2 *)samples;
-
-        for (int i = 0; i < mNumSamples; ++ i)
-        {
-                mHalton->GetNext(r);
-
-                // create samples over disc: the sample density
-                // decreases quadratically with the distance to the origin
-                s[i].x = mRadius * r[1] * sin(2.0f * M_PI * r[0]);
-                s[i].y = mRadius * r[1] * cos(2.0f * M_PI * r[0]);
-        }
-#else
-
-        //PoissonDiscSampleGenerator2D::Generate(samples);
-        RandomSampleGenerator2D::Generate(samples);
-
-        Sample2 *s = (Sample2 *)samples;
-
-        // multiply with lenght to get quadratic dependence on the distance
-        for (int i = 0; i < mNumSamples; ++ i)
-        {
-                Sample2 &spl = s[i];
-
-                float len = sqrt(spl.x * spl.x + spl.y * spl.y);
-                spl.x *= len * mRadius;
-                spl.y *= len * mRadius;
-        }
-#endif
-}
+        Sample2 *s = (Sample2 *)samples;
+
+        int numSamples = 0;
+        float x[2];
+
+        static float total1 = 0;
+        static float total2 = 0;
+        static int totalSamples = 0;
+
+        for (int i = 0; i < mNumSamples; ++ i)
+        {
+                //x[0] = RandomValue(0, 1); x[1] = RandomValue(0, 1);
+                mHalton->GetNext(x);
+                
+                const float a = 2.0f * M_PI * x[0];
+                const float r = sqrt(x[1]);
+                //const float rad = mRadius * r * r * r * r;
+                const float rad = mRadius * r;// * r;
+
+                s[i].x = rad * cos(a);
+                s[i].y = rad * sin(a);
+        }
+}