Index: /issm/trunk/src/c/shared/Numerics/BrentSearch.cpp
===================================================================
--- /issm/trunk/src/c/shared/Numerics/BrentSearch.cpp	(revision 476)
+++ /issm/trunk/src/c/shared/Numerics/BrentSearch.cpp	(revision 477)
@@ -64,5 +64,5 @@
 
 	//initialize the other variables
-	seps=sqrt(DBL_EPSILON); //precision of a double
+	seps=sqrt(DBL_EPSILON);    //precision of a double
 	distance=0.0;              //new_x=old_x + distance
 	gold=0.5*(3.0-sqrt(5.0));  //gold = 1 - golden ratio
@@ -115,8 +115,9 @@
 			intervalgold=distance;
 
-			// Is the parabola acceptable
+			// Is the parabola acceptable (we use seps here because in some configuration parab_num==parab_den*(xmax-xbest)
+			// and the result is not repeatable anymore
 			if (( sqrt(pow(parab_num,2)) < sqrt(pow(0.5*parab_den*oldintervalgold,2))) &&
-						(parab_num>parab_den*(xmin-xbest)) &&
-						(parab_num<parab_den*(xmax-xbest))){
+						(parab_num>parab_den*(xmin-xbest)+seps) && 
+						(parab_num<parab_den*(xmax-xbest)-seps)){
 
 				// Yes, parabolic interpolation step