Blame integration/positivity.c
|
Packit |
67cb25 |
/* Compare the integral of f(x) with the integral of |f(x)|
|
|
Packit |
67cb25 |
to determine if f(x) covers both positive and negative values */
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static inline int
|
|
Packit |
67cb25 |
test_positivity (double result, double resabs);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static inline int
|
|
Packit |
67cb25 |
test_positivity (double result, double resabs)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
int status = (fabs (result) >= (1 - 50 * GSL_DBL_EPSILON) * resabs);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return status;
|
|
Packit |
67cb25 |
}
|