/* Compare the integral of f(x) with the integral of |f(x)|

   to determine if f(x) covers both positive and negative values */

static inline int

test_positivity (double result, double resabs);

static inline int

test_positivity (double result, double resabs)

{

  int status = (fabs (result) >= (1 - 50 * GSL_DBL_EPSILON) * resabs);

  return status;

}