Changeset 1997 for GTP/trunk/Lib/Vis/Preprocessing/src/Mutation.cpp

Timestamp:

01/18/07 21:09:55 (17 years ago)

Author:

bittner

Message:

sil mutation

File:

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

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

@@ -18 +18 @@
 #define MUTATION_USE_CDF 0
 #define USE_SILHOUETTE_MUTATIONS 0
+
+#define USE_SIL_TERMINATION_MUTATION 1
+
 #define EVALUATE_MUTATION_STATS 1
 
@@ -257 +260 @@
 
 
+Vector3
+MutationBasedDistribution::ComputeSilhouetteTerminationMutation(const VssRay &ray,
+                                                                                                                                const Vector3 &origin,
+                                                                                                                                const AxisAlignedBox3 &box,
+                                                                                                                                const Vector3 &U,
+                                                                                                                                const Vector3 &V,
+                                                                                                                                const Vector2 vr2,
+                                                                                                                                const float radius
+                                                                                                                                )
+{
+  const int packetSize = 16;
+
+  static int hit_triangles[packetSize];
+  static float dist[packetSize];
+  static Vector3 dirs[packetSize];
+  // mutate the 
+  float alpha = RandomValue(0.0f, 2.0f*M_PI);
+  
+  // direction along which we will mutate the ray
+  Vector3 line = sin(alpha)*U + cos(alpha)*V;
+  
+  //  cout<<line<<endl;
+  // create 16 rays along the selected dir
+  
+  // cast rays to find silhouette ray
+  for (int i=0; i < packetSize; i++) {
+        dirs[i] = Normalize(ray.mTermination + ((radius/(packetSize - i))*line) - origin );
+        mlrtaStoreRayAS16(&origin.x,
+                                          &dirs[i].x,
+                                          i);
+  }
+  
+  mlrtaTraverseGroupAS16(&box.Min().x,
+                                                 &box.Max().x,
+                                                 hit_triangles,
+                                                 dist);
+  
+  for (int i=0; i < packetSize; i++) {
+        if (hit_triangles[i] == -1 || !box.IsInside(origin + dist[i]*dirs[i])) {
+          // break on first passing ray
+          break;
+        }
+  }
+  
+  if (i == packetSize) {
+        //      cerr<<"Warning: hit the same box here should never happen!"<<endl;
+        // shift the ray even a bit more
+        //      cout<<"W"<<i<<endl;
+        return ray.mTermination + (RandomValue(1.0f, 2.0f)*radius)*line;
+  }
+
+  //  cout<<i<<endl;
+  return dirs[i];
+}
 
 
@@ -385 +442 @@
                                                                   mutationRadius);
   
+#if USE_SIL_TERMINATION_MUTATION
+  termination += ComputeSilhouetteTerminationMutation(*ray,
+                                                                                                          origin,
+                                                                                                          box,
+                                                                                                          U, V,
+                                                                                                          Vector2(rr[2], rr[3]),
+                                                                                                          3.0f*objectRadius);
+#else
   termination += ComputeTerminationMutation(*ray, U, V,
                                                                                         Vector2(rr[2], rr[3]),
                                                                                         mutationRadius);
-  
+#endif
   Vector3 direction = termination - origin;