Correcting two double 1e-15 constant calculations.

In recent commits started setting two 1e-15 values with the calculation:

= 1/std::numeric_limits<DBL>::digits10;

when I meant to code:

= 1/pow(10,std::numeric_limits<DBL>::digits10);

Both forms bad practice in not being specific as to the types for the
explicit values. The first returns integer 0, later cast to double and I got
away with it in all but my full set of test cases. Corrected:

= (DBL)1.0/pow((DBL)10.0,std::numeric_limits<DBL>::digits10);

Author: wfpokorny

source/core/math/polynomialsolver.cpp

@@ -744,7 +744,7 @@ static bool regula_falsa(const int order, const DBL *coef, DBL a, DBL b, DBL *va
 
                 break;
             }
-            else if (fabs(fx) < (1/std::numeric_limits<DBL>::digits10))
+            else if (fabs(fx) < ((DBL)1.0/pow((DBL)10.0,std::numeric_limits<DBL>::digits10)))
             {
                 *val = x;
 

source/core/shape/blob.cpp

@@ -143,10 +143,10 @@ const DBL DEPTH_TOLERANCE = MIN_ISECT_DEPTH;
 /// @note
 ///     Value for @ref INSIDE_TOLERANCE was since v1.0 set to 1.0e-6 which is a good
 ///     value if running single floats. With double floats a value of 1e-15 is better.
-///     This is an offset to the inside of the 0.0 density surface after the blob
+///     This is an offset to the outside of the 0.0 density surface after the blob
 ///     threshold is subtracted. Now using std::numeric_limits to set.
 ///
-const DBL INSIDE_TOLERANCE = 1/std::numeric_limits<DBL>::digits10;
+const DBL INSIDE_TOLERANCE = (DBL)1.0/pow((DBL)10.0,std::numeric_limits<DBL>::digits10);
 
 /* Ray enters/exits a component. */
 const int ENTERING = 0;