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 The main O<span style='position: relative; top: 0.3em; font-size: 0.8em'>2</span>scl O$_2$scl names...
double rescale_error(double err, const double result_abs, const double result_asc)
QUADPACK's nonlinear rescaling of the absolute-error estimate.