Changeset 1524

Timestamp:

09/28/06 09:35:15 (18 years ago)

Author:

mattausch

Message:

 

Location:

GTP/trunk/Lib/Vis/Preprocessing/src

Files:

• GvsPreprocessor.cpp (modified) (9 diffs)
• Plane3.cpp (modified) (2 diffs)
• Preprocessor.cpp (modified) (2 diffs)

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

@@ -71 +71 @@
         mViewCellsManager->ComputeSampleContribution(*vssRay, true, false);
 
-        if ((vssRay->mPvsContribution - oldContri) > 0) // new contribution
+         // some pvs contribution for this ray?
+        if ((vssRay->mPvsContribution - oldContri) > 0)
         {
                 //cout << "h";
@@ -94 +95 @@
         const int indexL = (index == 0) ? 2 : index - 1;
 
-        const Vector3 a = hitTriangle.mVertices[index] - ray.GetDir();
+        const Vector3 a = hitTriangle.mVertices[index] - ray.GetOrigin();
         const Vector3 b = hitTriangle.mVertices[indexU] - hitTriangle.mVertices[index];
         const Vector3 c = hitTriangle.mVertices[index] - hitTriangle.mVertices[indexL];
@@ -102 +103 @@
         const Vector3 dir1 = Normalize(CrossProd(a, b)); //N((pi-xp)×(pi+1- pi));
         const Vector3 dir2 = Normalize(CrossProd(a, c)); // N((pi-xp)×(pi- pi-1))
-        const Vector3 dir3 = DotProd(dir1, dir2) > 0 ? // N((pi-xp)×di,i-1+di,i+1×(pi-xp))
-                Normalize(dir2 + dir1) : Normalize(CrossProd(a, dir2) + CrossProd(dir1, a)); 
+        const Vector3 dir3 = DotProd(dir2, dir1) > 0 ? // N((pi-xp)×di,i-1+di,i+1×(pi-xp))
+                Normalize(dir2 + dir1) : Normalize(CrossProd(a, dir1) + CrossProd(dir2, a)); 
 
         // compute the new three hit points
@@ -113 +114 @@
         const Vector3 pt3 = hitTriangle.mVertices[index] + eps * len * dir3;
         
+        vertices.push_back(pt2);
+        vertices.push_back(pt3);
         vertices.push_back(pt1);
-        //vertices.push_back(pt2);
-        vertices.push_back(pt3);
 }
 
@@ -196 +197 @@
         {
                 hitTriangle = dynamic_cast<TriangleIntersectable *>(tObj)->GetItem();
-                cout << "t: " << hitTriangle << endl << endl;
         }
         else
@@ -217 +217 @@
         {
                 const Vector3 rayDir = (*vit) - currentRay.GetOrigin();
-                simpleRays.push_back(SimpleRay(*vit, rayDir));
+                simpleRays.push_back(SimpleRay(currentRay.GetOrigin(), rayDir));
         }
 
@@ -226 +226 @@
         vizContainer.push_back(dummy);
 
-        // establish visibility
+        // cast rays to triangle vertices and determine visibility
         VssRayContainer vssRays;
         CastRays(simpleRays, vssRays, false);
-
+cout << "here211 " << simpleRays.size() << " " << vssRays.size() << endl;
         // add to ray queue
         EnqueueRays(vssRays);
-#if 0
+        const int n = (int)enlargedTriangle.size();
+        
+#if 1
     // recursivly subdivide each edge
-        for (int i = 0; i < 9; ++ i)
+        for (int i = 0; i < n; ++ i)
         {
                 SubdivideEdge(
                         hitTriangle,
                         enlargedTriangle[i], 
-                        enlargedTriangle[(i + 1) % 9], 
+                        enlargedTriangle[(i + 1) % n], 
                         *vssRays[i], 
-                        *vssRays[(i + 1) % 9],
+                        *vssRays[(i + 1) % n],
                         currentRay);
         }
@@ -343 +345 @@
                 castSamples += mInitialSamples;
 #endif
-                //cout << "\ncast " << castSamples << " of " << mSamplesPerPass << endl;
+                cout << "\nround finished, cast " << castSamples << " of " << mSamplesPerPass << endl;
         }
 
@@ -417 +419 @@
         for (vit = vizContainer.begin(); vit != vit_end; ++ vit)
         {
-                cout << "v";
+                exporter->SetWireframe();
                 exporter->ExportPolygon((*vit).enlargedTriangle);
                 //Material m;
-                //Polygon3 poly = Polygon3((*vit).originalTriangle);
-//              exporter->ExportPolygon(&poly);
+                exporter->SetFilled();
+                Polygon3 poly = Polygon3((*vit).originalTriangle);
+                exporter->ExportPolygon(&poly);
         }
 

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

@@ -75 +75 @@
         dir = CrossProd(plane1.mNormal, plane2.mNormal);
     
-        float abs;
         const int index = dir.DrivingAxis();
    
@@ -98 +97 @@
 
         Normalize(dir);
-
         return SimpleRay(point, dir);
 }

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

@@ -832 +832 @@
         else
         {
+#ifdef GTP_INTERNAL
                 mRayCaster = new IntelRayCaster(*this, externKdTree);
                 cout << "ray cast method: intel" << endl;
+#endif
         }
 
@@ -847 +849 @@
                                            )
 {
-        long t1 = GetTime();
+        const long t1 = GetTime();
+
+        for (int i = 0; i < (int)rays.size();)
+        {
+                if (i + 16 < (int)rays.size()) 
+                {
+                        mRayCaster->CastRays16(
+                                i,
+                                rays,                           
+                                vssRays,
+                                mViewCellsManager->GetViewSpaceBox(),
+                                castDoubleRays);
+                        i += 16;
+                }
+                else 
+                {
+                        mRayCaster->CastRay(
+                                rays[i].mOrigin,
+                                rays[i].mDirection,
+                                rays[i].mPdf,
+                                vssRays,
+                                mViewCellsManager->GetViewSpaceBox(),
+                                castDoubleRays);
+                        i ++;
+                }
+        }
+
+#if DEBUB_RAYCAST
+        long t2 = GetTime();
+        if (i % 10000 == 0)
+                cout<<".";
 
         if (castDoubleRays)
-        {
-                for (int i = 0; i < (int)rays.size(); ) 
-                {
-                        if (i + 16 < (int)rays.size()) 
-                        {
-                                mRayCaster->CastRays16(
-                                        i,
-                                        rays,                           
-                                        vssRays,
-                                        mViewCellsManager->GetViewSpaceBox(),
-                                        castDoubleRays);
-                                i += 16;
-                        } 
-                        else 
-                        {
-                                mRayCaster->CastRay(
-                                        rays[i].mOrigin,
-                                        rays[i].mDirection,
-                                        rays[i].mPdf,
-                                        vssRays,
-                                        mViewCellsManager->GetViewSpaceBox(),
-                                        castDoubleRays);
-                                i++;
-                        }
-                        
-                        if (i % 10000 == 0)
-                                cout<<".";
-                }
-
-                long t2 = GetTime();
-
-                if (castDoubleRays)
-                        cout << 2 * rays.size() / (1e3f * TimeDiff(t1, t2)) << "M rays/s" << endl;
-                else
-                        cout << rays.size() / (1e3f * TimeDiff(t1, t2)) << "M rays/s" << endl;
-        }
-        
+                cout << 2 * rays.size() / (1e3f * TimeDiff(t1, t2)) << "M rays/s" << endl;
+        else
+                cout << rays.size() / (1e3f * TimeDiff(t1, t2)) << "M rays/s" << endl;
+#endif  
 }