1 Thu Dec 21 21:46:54 2000 Brian Gough <bjg@network-theory.co.uk>
3 * changed err to abserr, to make it clear that it is an absolute
4 error in line with other modules.
6 Thu Nov 2 20:08:14 2000 Brian Gough <bjg@network-theory.co.uk>
8 * added support for a workspace so that the user does not have to
11 * made the names of the functions consistent as either levin_u,
12 for the full u transform with error estimate, or levin_utrunc for
13 the u transform with only a truncation error estimate.
15 Mon Apr 24 21:15:27 2000 Brian Gough <bjg@network-theory.co.uk>
17 * gsl_sum.h: added #include <stdlib.h> for size_t
19 Mon Nov 1 12:50:17 1999 Brian Gough <bjg@network-theory.co.uk>
21 * test.c (main): added tests using series for the Euler constant
24 * resolved problems with spurious failures by replacing the
25 directly computed truncation error by an estimate which varies
26 more smoothly. I have used the average of the previous two values,
27 which seems to give a reliable estimate of the truncation
28 error. The direct evaluation of the truncation error sometimes
29 fluctuated wildly, due to cancellation effects.
31 Thu Oct 28 12:05:47 1999 Brian Gough <bjg@network-theory.co.uk>
33 * test.c: cleaned up tests, now find that everything works in
34 double-precision but not extended-precision where there are two
35 failures depending on the optimization level
37 * levin_uerr.c (gsl_sum_levin_u_accel_minmax): changed loop
38 maximum from <=n to <n to avoid access of uninitialized memory,
39 need further correctness checks on algorithm though.
41 Thu Oct 7 12:03:36 1999 Brian Gough <bjg@network-theory.co.uk>
43 * levin_u.c levin_uerr.c: changed DBL_MAX to GSL_DBL_MAX since we
46 Sat Feb 6 20:35:26 1999 Brian Gough <bjg@netsci.freeserve.co.uk>
48 * test.c: adjusted the precision check to allow for "infinite
49 accuracy" which occurs when two results agree to machine precision
51 Thu Nov 19 13:10:19 1998 Brian Gough <bjg@vvv.lanl.gov>
53 * added an n_used parameter to all routines which gives the number
54 of terms actually used
56 Tue Nov 17 12:31:03 1998 Brian Gough <bjg@vvv.lanl.gov>
58 * test.c: added #include <config.h>
60 * renamed test_sum.c to test.c
62 * renamed all the functions so that _with_derivs is now the
63 default and _trunc is the case of no error estimates from the
66 * test_sum.c (main): cleaned up tests
68 Mon Nov 9 22:05:45 1998 Brian Gough <bjg@vvv.lanl.gov>
70 * levin_u.c (gsl_sum_levin_u_accel_minmax): got rid of noise
71 variables since they aren't used in the non-derivative case
73 1998-11-06 <bjg@ancho.lanl.gov>
75 * test_sum.c: replace variable N by macro to avoid variable length
78 Tue Oct 27 18:06:16 1998 Brian Gough <bjg@vvv.lanl.gov>
80 * levin_u.c: added in noise but it does not seem to be giving the
81 right answer for the error estimate. The actual value for the
82 accelerated sum is correct though. Check toms/602 for the original