Changeset 645 for GTP/trunk/Lib/Vis/Preprocessing/src/ViewCellBsp.cpp

Timestamp:

02/15/06 18:03:34 (18 years ago)

Author:

mattausch

Message:

 

File:

• GTP/trunk/Lib/Vis/Preprocessing/src/ViewCellBsp.cpp (modified) (9 diffs)

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

@@ -963 +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
@@ -973 +971 @@
         if (1)
         {
+                int pvsData = tData.mPvs;
+                float cData = (float)pvsData * tData.mProbability;
+
                 float cFront = (float)tFrontData.mPvs * tFrontData.mProbability;
                 float cBack = (float)tBackData.mPvs * tBackData.mProbability;
                 
-                float costDecr = 
-                        (cFront + cBack - cData) / mBox.GetVolume();
+                Debug << "front pvs: " << tFrontData.mPvs << " p: " << tFrontData.mProbability << endl;
+                Debug << "back pvs: " << tBackData.mPvs << " p: " << tBackData.mProbability << endl;
+                Debug << "pvs: " << pvsData << " p: " << tData.mProbability << endl;
+
+                Debug << "data: " << cData << " front: " << cFront << " back: " << cBack << endl;
+                Debug << "cost ratio: " << (cFront + cBack) / cData << endl;
+
+                if ((cFront + cBack) / cData > 1)
+                        Debug << "ERROR!!" << endl;
+
+                float costDecr = (cFront + cBack - cData) / mBox.GetVolume();
 
                 
@@ -1017 +1027 @@
 }
 
+
 void BspTree::ExtractViewCells(BspTraversalData &frontData,
                                                            BspTraversalData &backData, 
@@ -1044 +1055 @@
         }
 }
+
 
 BspInterior *BspTree::SubdivideNode(BspTraversalData &tData,
@@ -1091 +1103 @@
                                                                            interior->mPlane,
                                                                            mBox,
-                                                                           mEpsilon);
+                                                                           0.000000000001);
+                                                                           //mEpsilon);
         
                 
@@ -1103 +1116 @@
                         frontData.mProbability = frontData.mGeometry->GetVolume();
                         backData.mProbability = tData.mProbability - frontData.mProbability;
+
+                        Debug << "************" << endl;
+                        Debug << "plane: " << interior->mPlane << endl;
+                        Debug << "box: " << mBox << endl;
+                        Debug << "all: " << tData.mProbability << " f: " << frontData.mProbability << " b: " << backData.mProbability << endl;
+                        Debug << "real back: " << backData.mGeometry->GetVolume() << endl;
+                        Debug << "data geom: " << tData.mGeometry->mPolys.size() << " f: " << frontData.mGeometry->mPolys.size() << " b: " << backData.mGeometry->mPolys.size() << endl;
                 }
         }
@@ -1443 +1463 @@
                 const float candidateCost = 
                         SplitPlaneCost(poly->GetSupportingPlane(), data);
+
+                Debug << "supporting plane: " << poly->GetSupportingPlane() << endl;
 
                 if (candidateCost < lowestCost)
@@ -1982 +2004 @@
                         if (!leaf->mViewCell->Mailed())
                         {
-                          //                            ray.bspIntersections.push_back(Ray::BspIntersection(maxt, leaf));
+                                //ray.bspIntersections.push_back(Ray::BspIntersection(maxt, leaf));
                                 leaf->mViewCell->Mail();
                                 ++ hits;
@@ -2061 +2083 @@
                         else // ray and plane are coincident
                         {
-                                // WHAT TO DO IN THIS CASE ?
-                                //break;
-                                node = in->GetFront();
-                                continue;
+                                // NOTE: what to do if ray is coincident with plane?
+                                if (extSide < 0)
+                                        node = in->GetBack();
+                                else 
+                                        node = in->GetFront();
+                                                                
+                                continue; // no far child
                         }