This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Absolute error estimates
- From: "J.Underwood" <J dot Underwood at open dot ac dot uk>
- To: gsl-discuss <gsl-discuss at sources dot redhat dot com>
- Date: Sun, 22 Feb 2004 00:55:21 -0000
- Subject: Absolute error estimates
Dear all,
I'm just finally getting around to porting some code to GSL for the calculation of wigner rotation matrices, and am modifying existing code to return the answers as gsl_sf_result structures. However, I am unable to determine how the error estimate should be calculated in terms of GSL_DBL_EPSILON for the various arithmetic operations - is this prescribed anywhere for GSL (I've rtfm and the Design Document to no avail). Please could you direct me to how the error estimate should be calculated in general (naively I would use normal error analysis, however looking through various parts of the GSL library code, this seems to be incorrect).
Thanks
Jonathan.