Peter Rozovski writes: > result->err += GSL_DBL_EPSILON * (fabs(t1) + fabs(t2) + fabs(t3)); > result->val += 2.0 * GSL_DBL_EPSILON * fabs(result->val); > > It seems to me that in red line it should be result->err += ... > instead of result->val > += .
Yes, thanks for the bug report. -- Brian Gough Network Theory Ltd, Commercial support for GSL --- http://www.network-theory.co.uk/gsl/ _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
