Changeset 1996 for GTP/trunk/Lib/Vis/Preprocessing/src/SamplingStrategy.cpp

Timestamp:

01/18/07 20:20:48 (17 years ago)

Author:

mattausch

Message:

 

File:

• GTP/trunk/Lib/Vis/Preprocessing/src/SamplingStrategy.cpp (modified) (3 diffs)

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

@@ -20 +20 @@
 HaltonSequence HwGlobalLinesDistribution::sHalton;
 
+HaltonSequence ViewCellBasedDistribution::sHalton;
 
 SamplingStrategy::SamplingStrategy(Preprocessor &preprocessor): 
@@ -651 +652 @@
                                   } else
                                         if (strcmp(curr, "mutation")==0) {
-                                          mDistributions.push_back(new MutationBasedDistribution(mPreprocessor));
+                                                // temp matt: still no mutationstrategy!
+                                          //mDistributions.push_back(new MutationBasedDistribution(mPreprocessor));
                                         } 
         
@@ -738 +740 @@
 
 
-}
-
-
+bool ViewCellBasedDistribution::GenerateSample(SimpleRay &ray)
+{
+        Vector3 origin, direction; 
+
+        ViewCellContainer &viewCells = 
+                mPreprocessor.mViewCellsManager->GetViewCells();
+
+        Vector3 point;
+        Vector3 normal, normal2;
+        
+        float r[1];
+        sHalton.GetNext(1, r);
+
+        const int objIdx = (int)(r[0] * (float)mPreprocessor.mObjects.size() - 1.0f);
+                                //      = (int)RandomValue(0, (float)mPreprocessor.mObjects.size() - 0.5f);
+        Intersectable *object = mPreprocessor.mObjects[objIdx];
+
+        //cout << "obj: " << objIdx << endl;
+
+        object->GetRandomSurfacePoint(point, normal);
+        mViewCell->GetRandomSurfacePoint(origin, normal2);
+
+        direction = point - origin;
+
+        // move a little bit back to avoid piercing through walls
+        // that bound the view cell
+        origin -= 0.01f * normal2;
+
+        // $$ jb the pdf is yet not correct for all sampling methods!
+        const float c = Magnitude(direction);
+
+        if ((c <= Limits::Small) /*|| (DotProd(direction, normal) < 0)*/)
+        {
+                return false;
+        }
+
+        // $$ jb the pdf is yet not correct for all sampling methods!
+        const float pdf = 1.0f;
+
+        direction *= 1.0f / c;
+        ray = SimpleRay(origin, direction, VIEWCELL_BASED_DISTRIBUTION, pdf);
+
+        //cout << "ray: " << ray.mOrigin << " " << ray.mDirection << endl;
+
+        return true;
+}
+
+
+}
+
+