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