--- a/src/meshTools/triSurface/surfaceFeatures/surfaceFeatures.C	2009-09-04 12:51:54.000000000 +0200
+++ a/src/meshTools/triSurface/surfaceFeatures/surfaceFeatures.C	2009-09-04 12:54:25.000000000 +0200
@@ -762,10 +762,11 @@ 
     const scalarField& maxDist
 ) const
 {
+    treeBoundBox tbb(samples);
     // Build tree out of all samples.
     octree<octreeDataPoint> ppTree
     (
-        treeBoundBox(samples),      // overall search domain
+        tbb,      // overall search domain
         octreeDataPoint(samples),   // all information needed to do checks
         1,          // min levels
         20.0,       // maximum ratio of cubes v.s. cells
@@ -863,10 +864,11 @@ 
     scalar maxSearch = max(maxDist);
     vector span(maxSearch, maxSearch, maxSearch);
 
+    treeBoundBox tbb(samples);
     // octree.shapes holds reference!
     octree<octreeDataPoint> ppTree
     (
-        treeBoundBox(samples),      // overall search domain
+        tbb,      // overall search domain
         octreeDataPoint(samples),   // all information needed to do checks
         1,          // min levels
         20.0,       // maximum ratio of cubes v.s. cells