41 #ifndef O2SCL_INTE_GSL_H
42 #define O2SCL_INTE_GSL_H
50 #include <gsl/gsl_machine.h>
52 #ifndef DOXYGEN_NO_O2NS
107 const double result_asc) {
111 if (result_asc != 0 && err != 0) {
113 double scale=pow((200*err/result_asc),1.5);
116 err=result_asc*scale;
122 double dbl_eps=std::numeric_limits<double>::epsilon();
123 double dbl_min=std::numeric_limits<double>::min();
127 if (result_abs > dbl_min/(50*dbl_eps)) {
129 double min_err=50*dbl_eps*result_abs;
142 #ifndef DOXYGEN_NO_O2NS