1 /* Compare the integral of f(x) with the integral of |f(x)|
2    to determine if f(x) covers both positive and negative values */
3 
4 static inline int
5 test_positivity (double result, double resabs);
6 
7 static inline int
8 test_positivity (double result, double resabs)
9 {
10   int status = (fabs (result) >= (1 - 50 * GSL_DBL_EPSILON) * resabs);
11 
12   return status;
13 }
14