Changeset 474 for trunk/VUT/GtpVisibilityPreprocessor

Timestamp:

12/21/05 10:53:47 (19 years ago)

Author:

mattausch

Message:

 

Location:

trunk/VUT/GtpVisibilityPreprocessor

Files:

• scripts/default.env (modified) (3 diffs)
• src/ViewCellBsp.cpp (modified) (1 diff)
• src/ViewCellsManager.cpp (modified) (6 diffs)
• src/VspBspTree.cpp (modified) (11 diffs)
• src/VspBspTree.h (modified) (2 diffs)
• src/VspKdTree.h (modified) (1 diff)
• src/VssPreprocessor.cpp (modified) (4 diffs)

trunk/VUT/GtpVisibilityPreprocessor/scripts/default.env

@@ -153 +153 @@
         
         height 5.0
-        maxViewCells 5000
+        maxViewCells 3000
         
         PostProcessing {
@@ -236 +236 @@
         
         # maximal candidates for split planes
-        maxPolyCandidates 50
-        maxRayCandidates 50
+        maxPolyCandidates 100
+        maxRayCandidates 100
         
         # maximal tested rays for split cost heuristics
@@ -251 +251 @@
         Termination {
                 # parameters used for autopartition
-                minRays                 100
+                minRays                 10
                 minPolygons             -1
                 maxDepth                30
-                minPvs                  100
-                minArea                 0.01
+                minPvs                  20
+                minArea                 0.0001
                 maxRayContribution      0.005
                 maxCostRatio            0.9
-                missTolerance           1
+                missTolerance           3
                 #maxAccRayLength        100
                 

trunk/VUT/GtpVisibilityPreprocessor/src/ViewCellBsp.cpp

@@ -344 +344 @@
 
         app << "#N_SPLITS ( Number of splits )\n" << splits << "\n";
-
-        app << "#N_PMAXDEPTHLEAVES ( Percentage of leaves at maximum depth )\n" 
-                <<      maxDepthNodes * 100 / (double)Leaves() << endl;
 
         app << "#N_PMAXDEPTHLEAVES ( Percentage of leaves at maximum depth )\n" 

trunk/VUT/GtpVisibilityPreprocessor/src/ViewCellsManager.cpp

@@ -268 +268 @@
         const int limit = min(mConstructionSamples, (int)sourceRays.size());
 
-        Debug << "size: " << sourceRays.size() << " limit " << limit << endl;
         VssRayContainer::const_iterator it, it_end = sourceRays.end();
         
@@ -478 +477 @@
 
         }
-
 
         cout << "starting post processing using " << mPostProcessSamples << " samples ... ";
@@ -1422 +1420 @@
 
         GetRaySets(rays, constructionRays, savedRays);
-        Debug << "rays " << rays.size() << " construction rays " << constructionRays.size() << endl;
+        
         mVspBspTree->Construct(constructionRays);
-        
         mVspBspTree->CollectViewCells(mViewCells);
         
@@ -1631 +1628 @@
         bool exportGeometry = false;
 
-        Debug << "here344" << endl;
         environment->GetBoolValue("VspBspTree.Visualization.exportRays", exportRays);
         environment->GetBoolValue("VspBspTree.Visualization.exportGeometry", exportGeometry);
@@ -1687 +1683 @@
                         ObjectPvsMap::iterator it = vc->GetPvs().mEntries.begin();
 
-                        exporter->SetWireframe();
-                        //exporter->SetFilled();
+                        //exporter->SetWireframe();
+                        exporter->SetFilled();
 
                         Material m;//= RandomMaterial();
@@ -1843 +1839 @@
                 }
         }
-        Debug << "here222211" << endl;
 }
 

trunk/VUT/GtpVisibilityPreprocessor/src/VspBspTree.cpp

@@ -38 +38 @@
 mRoot(NULL),
 mPvsUseArea(true),
-mNumCriteria(0)
+mCostNormalizer(Limits::Small)
 {
         mRootCell = new BspViewCell();
@@ -91 +91 @@
         Debug << "max polygon candidates: " << mMaxPolyCandidates << endl;
         Debug << "max plane candidates: " << mMaxRayCandidates << endl;
-
+        
         Debug << "Split plane strategy: ";
         if (mSplitPlaneStrategy & RANDOM_POLYGON)
+        {
                 Debug << "random polygon ";
-
+        }
         if (mSplitPlaneStrategy & AXIS_ALIGNED)
         {
-                ++ mNumCriteria;
                 Debug << "axis aligned ";
         }
         if (mSplitPlaneStrategy & LEAST_RAY_SPLITS)
         {
-                ++ mNumCriteria;
+                mCostNormalizer += mLeastRaySplitsFactor;
                 Debug << "least ray splits ";
         }
         if (mSplitPlaneStrategy & BALANCED_RAYS)
         {
-                ++ mNumCriteria;
+                mCostNormalizer += mBalancedRaysFactor;
                 Debug << "balanced rays ";
         }
         if (mSplitPlaneStrategy & PVS)
         {
-                ++ mNumCriteria;
+                mCostNormalizer += mPvsFactor;
                 Debug << "pvs";
         }
@@ -273 +273 @@
         }
 
+        mTermMinArea *= mBox.SurfaceArea(); // normalize
         mStat.polys = (int)polys.size();
 
@@ -315 +316 @@
 
             tStack.pop();
-
                 // subdivide leaf node
                 BspNode *r = Subdivide(tStack, tData);
@@ -412 +412 @@
                 ++ maxCostMisses;
 
-                /*if (maxCostMisses >= mTermMissTolerance)
+                if (maxCostMisses >= mTermMissTolerance)
                 {
                         // terminate branch because of max cost
                         ++ mStat.maxCostNodes; 
             return leaf;
-                }*/
+                }
         }
 
@@ -471 +471 @@
                                                                            mEpsilon);
         
-                // area is normalized with view space area
-                frontData.mArea = frontData.mGeometry->GetArea() / mBox.SurfaceArea();
-                backData.mArea = backData.mGeometry->GetArea() / mBox.SurfaceArea();
+                frontData.mArea = frontData.mGeometry->GetArea();
+                backData.mArea = backData.mGeometry->GetArea();
         }
 
@@ -810 +809 @@
                 const int candidateIdx = (int)RandomValue(0, (Real)(-- maxIdx));
                 //Debug << "current Idx: " << maxIdx << " cand idx " << candidateIdx << endl;
-
+                
                 Polygon3 *poly = (*data.mPolygons)[candidateIdx];
-                
+
                 // swap candidate to the end to avoid testing same plane
                 std::swap((*data.mPolygons)[maxIdx], (*data.mPolygons)[candidateIdx]);
-
+        
                 //Polygon3 *poly = (*data.mPolygons)[(int)RandomValue(0, (int)polys.size() - 1)];
 
@@ -829 +828 @@
         }
         
-        //Debug << "lowest: " << lowestCost << endl;
-
         //-- choose candidate planes extracted from rays
-        // we use different methods chosen with
-        // equal probability
+        //-- different methods are available
         for (int i = 0; i < mMaxRayCandidates; ++ i)
         {
                 plane = ChooseCandidatePlane3(*data.mRays);
-
                 const float candidateCost = SplitPlaneCost(plane, data);
 
@@ -846 +841 @@
                 }
         }
+
+#ifdef _DEBUG
+        Debug << "plane lowest cost: " << lowestCost << endl;
+#endif
 
         // cost ratio miss
         if (lowestCost > mTermMaxCostRatio)
                 return false;
-
-#ifdef _DEBUG
-        Debug << "plane lowest cost: " << lowestCost << endl;
-#endif
 
         return true;
@@ -1011 +1006 @@
         {
                 const float oldCost = pOverall * (float)data.mPvs + Limits::Small;
-
-                cost += mPvsFactor * (frontPvs * pFront + (backPvs * pBack)) / oldCost;
-
-                // give penalty to unbalanced split
-                if (0)
-                if (((pFront * 0.2 + Limits::Small) > pBack) || 
-                        (pFront < (pBack * 0.2 + Limits::Small)))
-                        cost += 0.5;
+                cost += mPvsFactor * (frontPvs * pFront + backPvs * pBack) / oldCost;
+
+                //cost += mPvsFactor * 0.5 * (frontPvs * pFront + backPvs * pBack) / oldCost;
+                //Debug << "new cost: " << cost << " over" << frontPvs * pFront + backPvs * pBack << " old cost " << oldCost << endl;
+        
+                if (0) // give penalty to unbalanced split
+                        if (((pFront * 0.2 + Limits::Small) > pBack) || 
+                                (pFront < (pBack * 0.2 + Limits::Small)))
+                                        cost += 0.5;
         }
 
 #ifdef _DEBUG
-//      Debug << "totalpvs: " << pvs << " ptotal: " << pOverall
-//                << " frontpvs: " << frontPvs << " pFront: " << pFront 
-//                << " backpvs: " << backPvs << " pBack: " << pBack << endl << endl;
+        Debug << "totalpvs: " << data.mPvs << " ptotal: " << pOverall
+                  << " frontpvs: " << frontPvs << " pFront: " << pFront 
+                  << " backpvs: " << backPvs << " pBack: " << pBack << endl << endl;
 #endif
-
-        // normalize by number of criteria
-        return cost / (float)(mNumCriteria + Limits::Small);
+        
+        // normalize cost by sum of linear factors
+        return cost / (float)mCostNormalizer;
 }
 
@@ -1118 +1114 @@
         if (data.mDepth < mStat.minDepth)
                 mStat.minDepth = data.mDepth;
-
+        
+        if (data.mDepth >= mTermMaxDepth)
+                ++ mStat.maxDepthNodes;
+
+        if (data.mPvs < mTermMinPvs)
+                ++ mStat.minPvsNodes;
+
+        if ((int)data.mRays->size() < mTermMinRays)
+                ++ mStat.minRaysNodes;
+
+        if (data.GetAvgRayContribution() > mTermMaxRayContribution)
+                ++ mStat.maxRayContribNodes;
+        
+        if (data.mArea <= mTermMinArea) 
+        {
+                //Debug << "area: " << data.mArea / mBox.SurfaceArea() << " min area: " << mTermMinArea / mBox.SurfaceArea() << endl;
+                ++ mStat.minAreaNodes;
+        }
         // accumulate depth to compute average depth
         mStat.accumDepth += data.mDepth;
-        
-
-        if (data.mDepth >= mTermMaxDepth)
-                ++ mStat.maxDepthNodes;
-
-        if (data.mPvs < mTermMinPvs)
-                ++ mStat.minPvsNodes;
-
-        if ((int)data.mRays->size() < mTermMinRays)
-                ++ mStat.minRaysNodes;
-
-        if (data.GetAvgRayContribution() > mTermMaxRayContribution)
-                ++ mStat.maxRayContribNodes;
-        
-        if (data.mGeometry->GetArea() <= mTermMinArea) 
-                ++ mStat.minAreaNodes;
 
 #ifdef _DEBUG

trunk/VUT/GtpVisibilityPreprocessor/src/VspBspTree.h

@@ -112 +112 @@
                 {}
 
+                /** Returns cost of the traversal data.
+                */
+                float GetCost() const
+                {
+#if 0
+                        return mPvs * mArea;
+#endif
+#if 1
+                        return (float)(mPvs * (int)mRays->size());
+#endif
+#if 0
+                        return (float)mPvs;
+#endif
+#if 0
+                        return mArea * (float)mRays->size();
+#endif
+                }
+
                 friend bool operator<(const VspBspTraversalData &a, const VspBspTraversalData &b)
                 {
-#if 0
-                        return a.mPvs * a.mArea < b.mPvs * b.mArea;
-#endif
-#if 1
-                        return a.mPvs * (int)a.mRays->size() < b.mPvs * (int)b.mRays->size();
-#endif
-#if 0
-                        return a.mPvs < b.mPvs;
-#endif
-#if 0
-                        return 
-                                a.mPvs / (float)(a.mRays->size() + Limits::Small())
-                                >
-                                b.mPvs / (float)(b.mRays->size() + Limits::Small());
-#endif
-#if 0
-                        return a.mPvs * (int)a.mRays->size() <  b.mPvs * (int)b.mRays->size();
-#endif
+                        return a.GetCost() < b.GetCost();
                 }
     };
@@ -544 +545 @@
         /// maximal number of test rays used to evaluate candidate split plane
         int mMaxTests;
-        /// number of different bsp split plane criteria
-        int mNumCriteria;
+        /// normalizes different bsp split plane criteria
+        float mCostNormalizer;
         /// maximal number of view cells
         int mMaxViewCells;

trunk/VUT/GtpVisibilityPreprocessor/src/VspKdTree.h

@@ -477 +477 @@
 #if 0
                         return
-                                leafa->rays.size()*a.mBox.GetVolume()
+                                leafa->rays.size() * a.mBox.GetVolume()
                                 <
-                                leafb->rays.size()*b.mBox.GetVolume();
+                                leafb->rays.size() * b.mBox.GetVolume();
 #endif
 #if 1
                         return
-                                leafa->GetPvsSize()*a.mBox.GetVolume()
+                                leafa->GetPvsSize() * a.mBox.GetVolume()
                                 <
-                                leafb->GetPvsSize()*b.mBox.GetVolume();
+                                leafb->GetPvsSize() * b.mBox.GetVolume();
 #endif
 #if 0

trunk/VUT/GtpVisibilityPreprocessor/src/VssPreprocessor.cpp

@@ -165 +165 @@
   if (!use2dSampling) {
         Vector3 normal;
-        int i = RandomValue(0, mObjects.size()-1);
+        int i = (int)RandomValue(0, (Real)((int)mObjects.size()-1));
         Intersectable *object = mObjects[i];
         object->GetRandomSurfacePoint(point, normal);
@@ -561 +561 @@
 
         mViewCellsManager->Visualize(mObjects, storedRays);
-
-        Debug << "here2" << endl;
   }
 
@@ -569 +567 @@
   
   mRenderSimulator->RenderScene();
-  Debug << "here3" << endl;
+  
   SimulationStatistics ss;
   mRenderSimulator->GetStatistics(ss);
-  Debug << "here244" << endl;
+  
   cout << " finished" << endl;
   cout << ss << endl;
@@ -579 +577 @@
   delete vssTree;
 
-  Debug << "here4" << endl;
   return true;
 }