reset default tolerance to smaller value

Post/OctreePost.cpp

@@ -393,22 +393,8 @@ bool OctreePost::_getValue(void *in, int nbComp, double P[3], int timestep,
   return true;
 } 
 
-bool OctreePost::searchScalar(double x, double y, double z, double *values, 
-                              int step, double *size)
-{
-  bool a =  searchScalar_(x,y,z,values,step,size);
-  if (!a){
-    element::setTolerance(10.);
-    a =  searchScalar_(x,y,z,values,step,size);
-    element::setTolerance(1.e-3);
-  }    
-  if (!a)printf("cannot find %g %g %g\n",x,y,z);
-  return a;
-}
-
-
-bool OctreePost::searchScalar_(double x, double y, double z, double *values, 
-                              int step, double *size)
+bool OctreePost::_searchScalar(double x, double y, double z, double *values, 
+                               int step, double *size)
 {
   double P[3] = {x, y, z};
 
@@ -438,6 +424,20 @@ bool OctreePost::searchScalar_(double x, double y, double z, double *values,
   return false;
 }
 
+bool OctreePost::searchScalar(double x, double y, double z, double *values, 
+                              int step, double *size)
+{
+  bool a = _searchScalar(x, y, z, values, step, size);
+  if (!a){
+    double oldeps = element::getTolerance();
+    element::setTolerance(10.);
+    a = _searchScalar(x, y, z, values, step, size);
+    element::setTolerance(oldeps);
+  }    
+  if (!a) Msg::Debug("No element found containing point (%g,%g,%g)", x, y, z);
+  return a;
+}
+
 bool OctreePost::searchVector(double x, double y, double z, double *values, 
                               int step, double *size)
 {

Post/OctreePost.h

@@ -30,8 +30,8 @@ class OctreePost
                  double *elementSize);
   bool _getValue(void *in, int nbComp, double P[3], int step, 
                  double *values, double *elementSize);
-  bool searchScalar_(double x, double y, double z, double *values, 
-                    int step = -1, double *size = 0);
+  bool _searchScalar(double x, double y, double z, double *values, 
+                     int step = -1, double *size = 0);
  public :
   OctreePost(PView *);
   ~OctreePost();

Post/shapeFunctions.cpp

+// Gmsh - Copyright (C) 1997-2009 C. Geuzaine, J.-F. Remacle
+//
+// See the LICENSE.txt file for license information. Please report all
+// bugs and problems to <gmsh@geuz.org>.
+
 #include "shapeFunctions.h"
-double element::ONE = 1.+1.e-2;
-double element::ZERO = -1.e-2;
+
+double element::ONE = 1. + 1.e-6;
+double element::ZERO = -1.e-6;

Post/shapeFunctions.h

@@ -17,13 +17,8 @@ protected:
   static double ONE, ZERO;
 public:
   element(double *x, double *y, double *z) : _x(x), _y(y), _z(z) {}
-  static void setTolerance (const double tol){
-    ONE = 1.+tol;
-    ZERO = -tol;
-  }
-  static double getTolerance () {
-    return -ZERO;
-  }
+  static void setTolerance (const double tol){ ONE = 1. + tol; ZERO = -tol; }
+  static double getTolerance () { return -ZERO; }
   virtual ~element(){}
   virtual int getDimension() = 0;
   virtual int getNumNodes() = 0;