Changeset 1919 for GTP/trunk/Lib/Vis/Preprocessing/src/ViewCellsManager.cpp

Timestamp:

12/21/06 23:52:10 (18 years ago)

Author:

mattausch

Message:

added mechanism for histogram on certain MB in hierarchymanager
no more bug in undersampling estimation
added sampling strategy to spatial box based (also for eval)
added testing scripts

File:

• GTP/trunk/Lib/Vis/Preprocessing/src/ViewCellsManager.cpp (modified) (11 diffs)

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

@@ -147 +147 @@
 
         // mix of sampling strategies
-        mStrategies.push_back(SamplingStrategy::OBJECT_DIRECTION_BASED_DISTRIBUTION);
-        
-        /*mStrategies.push_back(SamplingStrategy::OBJECT_BASED_DISTRIBUTION);
-        mStrategies.push_back(SamplingStrategy::REVERSE_VIEWSPACE_BORDER_BASED_DISTRIBUTION);
-        mStrategies.push_back(SamplingStrategy::REVERSE_OBJECT_BASED_DISTRIBUTION);
-        */
-        //mStrategies.push_back(SamplingStrategy::SPATIAL_BOX_BASED_DISTRIBUTION);
+        if (0)
+        {
+                mStrategies.push_back(SamplingStrategy::OBJECT_DIRECTION_BASED_DISTRIBUTION);
+        }
+        else
+        {
+                /*mStrategies.push_back(SamplingStrategy::OBJECT_BASED_DISTRIBUTION);
+                mStrategies.push_back(SamplingStrategy::REVERSE_VIEWSPACE_BORDER_BASED_DISTRIBUTION);
+                mStrategies.push_back(SamplingStrategy::REVERSE_OBJECT_BASED_DISTRIBUTION);
+                */
+        
+                mStrategies.push_back(SamplingStrategy::SPATIAL_BOX_BASED_DISTRIBUTION);
+        }
                 
     Debug << "casting strategies: ";
@@ -955 +961 @@
 
         it = viewCells.begin();
+        
+        for (int j = 0; j < intervals; ++ j)
+        {
+                float volDif = 0;
+                int smallerDif = 0;
+
+                while ((i < (int)viewCells.size()) && 
+                           (viewCells[i]->GetPvsCost() < currentPvs))
+                {
+                        volDif += viewCells[i]->GetVolume();
+                        volSum += viewCells[i]->GetVolume();
+
+                        ++ i;
+                        ++ smallerDif;
+                        ++ smallerSum;
+                }
+                
+                //              if (0 && (i < (int)viewCells.size()))
+                //                Debug << "new pvs cost increase: " << mViewCellsTree->GetPvsCost(viewCells[i]) 
+                //                              << " " << currentPvs << endl;
+                
+                const float volRatioDif = volDif / totalVol;
+                const float volRatioSum = volSum / totalVol;
+
+                outstream << "#Pass\n" << j << endl;
+                outstream << "#Pvs\n" << currentPvs << endl;
+                outstream << "#ViewCellsDif\n" << smallerDif << endl;
+                outstream << "#ViewCellsSum\n" << smallerSum << endl;   
+                outstream << "#VolumeDif\n" << volRatioDif << endl << endl;
+                outstream << "#VolumeSum\n" << volRatioSum << endl << endl;
+        
+                //-- increase current pvs size to define next interval
+                currentPvs += stepSize;
+        }
+
+        outstream.close();
+}
+
+
+void ViewCellsManager::EvalViewCellHistogramForPvsSize(const string filename, 
+                                                                                                           ViewCellContainer &viewCells)
+{
+        std::ofstream outstream;
+        outstream.open(filename.c_str());
+
+        float maxPvs, maxVal, minVal;
+        
+        // sort by pvs size
+        sort(viewCells.begin(), viewCells.end(), ViewCell::SmallerPvs);
+
+        maxPvs = viewCells.back()->GetPvsCost();
+        minVal = 0;
+
+        // hack: normalize pvs size
+        int histoMaxVal;
+        Environment::GetSingleton()->GetIntValue("Preprocessor.histogram.maxValue", histoMaxVal);
+        maxVal = max((float)histoMaxVal, maxPvs);
+                
+        Debug << "histogram minpvssize: " << minVal << " maxpvssize: " << maxVal 
+                << " real maxpvs " << maxPvs << endl;
+
+        int histoIntervals;
+        Environment::GetSingleton()->GetIntValue("Preprocessor.histogram.intervals", histoIntervals);
+        const int intervals = min(histoIntervals, (int)viewCells.size());
+
+        const float range = maxVal - minVal;
+        int stepSize = (int)(range / intervals);
+
+        // set step size to avoid endless loop
+        if (!stepSize) stepSize = 1;
+        
+        Debug << "intervals " << histoIntervals << endl;
+        Debug << "stepsize: " << stepSize << endl;
+        cout << "intervals " << histoIntervals << endl;
+        cout << "stepsize: " << stepSize << endl;
+
+        const float totalRenderCost = mViewCellsTree->GetRoot()->GetRenderCost();
+        const float totalVol = GetViewSpaceBox().GetVolume();
+
+        float currentPvs = minVal;
+        
+        int i = 0;
+        int j = 0;
+        float volSum = 0;
+        int smallerSum = 0;
+
+        //ViewCellContainer::const_iterator it = viewCells.begin(), it_end = viewCells.end();   
         
         for (int j = 0; j < intervals; ++ j)
@@ -6031 +6124 @@
         Environment::GetSingleton()->GetIntValue("ViewCells.Evaluation.stepSize", splitsStepSize);
         
+        bool useHisto;
+        int histoMem;
+
+        Environment::GetSingleton()->GetBoolValue("ViewCells.Evaluation.histogram", useHisto);
+        Environment::GetSingleton()->GetIntValue("ViewCells.Evaluation.histoMem", histoMem);
+
         Debug << "step size: " << splitsStepSize << endl;
         Debug << "view cell evaluation samples per pass: " << samplesPerPass << endl;
@@ -6047 +6146 @@
 
         // mix of sampling strategies
-        evalStrats.push_back(SamplingStrategy::OBJECT_DIRECTION_BASED_DISTRIBUTION);
-
-        /*evalStrats.push_back(SamplingStrategy::OBJECT_BASED_DISTRIBUTION);
-        evalStrats.push_back(SamplingStrategy::REVERSE_VIEWSPACE_BORDER_BASED_DISTRIBUTION);
-        evalStrats.push_back(SamplingStrategy::REVERSE_OBJECT_BASED_DISTRIBUTION);
-*/
-        //evalStrats.push_back(SamplingStrategy::SPATIAL_BOX_BASED_DISTRIBUTION);
+        if (0)
+        {
+                evalStrats.push_back(SamplingStrategy::OBJECT_DIRECTION_BASED_DISTRIBUTION);
+        }
+        else
+        {
+                /*evalStrats.push_back(SamplingStrategy::OBJECT_BASED_DISTRIBUTION);
+                evalStrats.push_back(SamplingStrategy::REVERSE_VIEWSPACE_BORDER_BASED_DISTRIBUTION);
+                evalStrats.push_back(SamplingStrategy::REVERSE_OBJECT_BASED_DISTRIBUTION);
+                */
+                evalStrats.push_back(SamplingStrategy::SPATIAL_BOX_BASED_DISTRIBUTION);
+        }
                         
     Debug << "casting eval strategies: ";
@@ -6064 +6168 @@
                 cout << evalStrats[i] << " ";
         cout << endl;
+
+        int pass = 0;
 
         while (castSamples < numSamples)
@@ -6122 +6228 @@
 
                         ofstream ofstr(filename.c_str());
-                        mHierarchyManager->EvaluateSubdivision2(ofstr, splitsStepSize, false);
+                        mHierarchyManager->EvaluateSubdivision2(ofstr, splitsStepSize, false, useHisto, histoMem, pass);
 
                         timeDiff = TimeDiff(startTime, GetTime());
@@ -6167 +6273 @@
                                          << (float)pvsSize * ObjectPvs::GetEntrySize() << endl;
                         }
+                        ++ pass;
                 }
 
@@ -6174 +6281 @@
         ////////////
         //-- histogram
-
+#if 1
         const int numLeaves = mViewCellsTree->GetNumInitialViewCells(mViewCellsTree->GetRoot());
-        bool useHisto;
         int histoStepSize;
 
@@ -6185 +6291 @@
         {
                 // evaluate view cells in a histogram           
-                char s[64];
+                char str[64];
 
                 // hack: just show final view cells
@@ -6195 +6301 @@
 
                         cout << "computing histogram for " << pass << " view cells" << endl;
-#if 0
-                        //-- evaluate histogram for render cost
-                        sprintf(s, "-%09d-histo.log", pass);
-                        filename = string(statsPrefix) + string(s);
-
-                        EvalViewCellHistogram(filename, pass);
-
-#endif
+
                         //////////////////////////////////////////
             //-- evaluate histogram for pvs size
@@ -6208 +6307 @@
                         cout << "computing pvs histogram for " << pass << " view cells" << endl;
 
-                        sprintf(s, "-%09d-histo-pvs.log", pass);
-                        filename = string(statsPrefix) + string(s);
+                        sprintf(str, "-%09d-histo-pvs2.log", pass);
+                        filename = string(statsPrefix) + string(str);
 
                         EvalViewCellHistogramForPvsSize(filename, pass);
                 }
         }
-}
-
 #endif
-
-
-
-}
+}
+
+#endif
+
+
+
+}