Changeset 611 for trunk/VUT/GtpVisibilityPreprocessor/src/ViewCellBsp.cpp

Timestamp:

02/10/06 13:04:39 (18 years ago)

Author:

mattausch

Message:

 

File:

• trunk/VUT/GtpVisibilityPreprocessor/src/ViewCellBsp.cpp (modified) (11 diffs)

trunk/VUT/GtpVisibilityPreprocessor/src/ViewCellBsp.cpp

@@ -262 +262 @@
 
         environment->GetFloatValue("BspTree.Construction.epsilon", mEpsilon);
-        
+
+        mSubdivisionStats.open("subdivisionStats.log");
+
     Debug << "BSP max depth: " << mTermMaxDepth << endl;
         Debug << "BSP min PVS: " << mTermMinPvs << endl;
@@ -841 +843 @@
                                                    geom);
 
+        mTotalCost = tData.mPvs * tData.mProbability / mBox.GetVolume();
+        mTotalPvsSize = tData.mPvs;
+        
+        mSubdivisionStats 
+                        << "#ViewCells\n1\n" <<  endl
+                        << "#RenderCostDecrease\n0\n" << endl 
+                        << "#TotalRenderCost\n" << mTotalCost << endl
+                        << "#AvgRenderCost\n" << mTotalPvsSize << endl;
+
         tStack.push(tData);
 
+        // used for intermediate time measurements and progress
+        long interTime = GetTime();     
+        int nleaves = 500;
+
         mStat.Start();
-        cout << "Contructing bsp tree ... ";
+        cout << "Contructing bsp tree ...\n";
         long startTime = GetTime();
         while (!tStack.empty()) 
@@ -858 +873 @@
             Debug << "BSP tree construction time spent at root: " 
                                   << TimeDiff(startTime, GetTime())*1e-3 << " secs" << endl;
+
+                if (mStat.Leaves() >= nleaves)
+                {
+                        nleaves += 500;
+                        
+                        cout << "leaves=" << mStat.Leaves() << endl;
+                        Debug << "needed "
+                                  << TimeDiff(interTime, GetTime())*1e-3 
+                                  << " secs to create 500 leaves" << endl;
+                        interTime = GetTime();
+                }
         }
 
@@ -937 +963 @@
                                                            new BoundedRayContainer(), 0, 0, new BspNodeGeometry());
 
+        int pvsData = tData.mPvs;
+        float cData = (float)tData.mPvs * tData.mProbability;
+
         // create new interior node and two leaf nodes
         BspInterior *interior = 
                 SubdivideNode(tData, tFrontData, tBackData, coincident);
 
+
+        if (1)
+        {
+                float cFront = (float)tFrontData.mPvs * tFrontData.mProbability;
+                float cBack = (float)tBackData.mPvs * tBackData.mProbability;
+                
+                float costDecr = 
+                        (cFront + cBack - cData) / mBox.GetVolume();
+
+                
+                mTotalCost += costDecr;
+                mTotalPvsSize += tFrontData.mPvs + tBackData.mPvs - pvsData;
+
+                mSubdivisionStats 
+                        << "#ViewCells\n" << mStat.Leaves() << endl
+                        << "#RenderCostDecrease\n" << -costDecr << endl 
+                        << "#TotalRenderCost\n" << mTotalCost << endl
+                        << "#AvgRenderCost\n" << mTotalPvsSize / mStat.Leaves() << endl;
+        }
 
         // extract view cells from coincident polygons according to plane normal
@@ -1006 +1054 @@
         BspLeaf *leaf = dynamic_cast<BspLeaf *>(tData.mNode);
 
-        long startTime;
-        if (0)
-        {
-                Debug << "*********************" << endl;
-                startTime = GetTime();
-        }
+        
         
         // select subdivision plane
@@ -1017 +1060 @@
                 new BspInterior(SelectPlane(leaf, tData)); 
         
-        if (0)
-        {
-                Debug << "time used for split plane selection: " 
-                          << TimeDiff(startTime, GetTime()) * 1e-3 << "s" << endl;
-        }
+        
 #ifdef _DEBUG
         Debug << interior << endl;
@@ -1027 +1066 @@
         
 
-        if (0)
-        {
-                Debug << "number of rays: " << (int)tData.mRays->size() << endl;
-                Debug << "number of polys: " << (int)tData.mPolygons->size() << endl;
-
-                startTime = GetTime();
-        }       
-        
         // subdivide rays into front and back rays
         SplitRays(interior->mPlane, *tData.mRays, *frontData.mRays, *backData.mRays);
         
-        if (0)
-        {
-                Debug << "time used for rays splitting: " << TimeDiff(startTime, GetTime())*1e-3 << "s" << endl;
-                startTime = GetTime();
-        }
+        
 
         // subdivide polygons with plane
@@ -1051 +1078 @@
                                                                           coincident);
 
-        if (0)
-        {
-                Debug << "time used for polygon splitting: " << TimeDiff(startTime, GetTime())*1e-3 << "s" << endl;
-        }
+        
 
     // compute pvs
@@ -1070 +1094 @@
         
                 
-                frontData.mProbability = frontData.mGeometry->GetVolume();
-                backData.mProbability = backData.mGeometry->GetVolume();
-        }
-
-        // compute accumulated ray length
-        //frontData.mAccRayLength = AccumulatedRayLength(*frontData.mRays);
-        //backData.mAccRayLength = AccumulatedRayLength(*backData.mRays);
+                if (mUseAreaForPvs)
+                {
+                        frontData.mProbability = frontData.mGeometry->GetArea();
+                        backData.mProbability = backData.mGeometry->GetArea();
+                }
+                else
+                {
+                        frontData.mProbability = frontData.mGeometry->GetVolume();
+                        backData.mProbability = tData.mProbability - frontData.mProbability;
+                }
+        }
 
         //-- create front and back leaf
@@ -1582 +1610 @@
 inline void BspTree::GenerateUniqueIdsForPvs()
 {
-        Intersectable::NewMail(); sBackId = ViewCell::sMailId;
-        Intersectable::NewMail(); sFrontId = ViewCell::sMailId;
-        Intersectable::NewMail(); sFrontAndBackId = ViewCell::sMailId;
+        ViewCell::NewMail(); sBackId = ViewCell::sMailId;
+        ViewCell::NewMail(); sFrontId = ViewCell::sMailId;
+        ViewCell::NewMail(); sFrontAndBackId = ViewCell::sMailId;
 }
 
@@ -1895 +1923 @@
                 return 0;
 
-        Intersectable::NewMail();
+        ViewCell::NewMail();
 
         Vector3 entp = ray.Extrap(mint);