1 2008-01-28 Brian Gough <bjg@network-theory.co.uk>
3 * test.c (test_evolve_negative_h): test evolution with negative
4 steps as well as positive
6 2008-01-03 Brian Gough <bjg@network-theory.co.uk>
8 * evolve.c (gsl_odeiv_evolve_apply): handle the case where
9 GSL_ODEIV_HADJ_DEC does not represent a real decrease (e.g. where
10 x+h = x due to finite precision)
12 2006-07-29 Brian Gough <bjg@network-theory.co.uk>
14 * evolve.c (gsl_odeiv_evolve_apply): notify user of step-size
15 which caused any failure by returning it
17 2006-06-13 Brian Gough <bjg@network-theory.co.uk>
19 * gsl_odeiv.h, control.c (gsl_odeiv_control_hadjust): change y0 to
20 y in prototype to avoid conflict with bessel y0() in math.h
22 * test.c (test_odeiv_stepper): change y0 to ystart
24 2005-04-08 Brian Gough <bjg@network-theory.co.uk>
26 * Fixed all functions to consistently return status instead of
29 2004-12-29 Brian Gough <bjg@network-theory.co.uk>
31 * test.c (test_compare_vanderpol): avoid variable size array
32 (test_compare_oregonator): as above
34 2004-12-23 Brian Gough <bjg@network-theory.co.uk>
36 * test.c (test_evolve_xsin): changed description from "failures"
39 2004-12-22 Brian Gough <bjg@network-theory.co.uk>
41 * rk2simp.c (rk2simp_step): added missing static declaration
43 * gear1.c (gear1_step): added missing static declaration
45 * gear2.c (gear2_step): added missing static declaration
47 * rk4imp.c (rk4imp_step): added missing static declaration
49 * rk4.c (rk4_step): added missing static declaration
51 * rk2imp.c (rk2imp_step): added missing static declaration
53 2004-12-01 Brian Gough <bjg@network-theory.co.uk>
55 * updated all solvers to ensure correct handling of error codes
56 from user defined functions
58 * test.c (rhs_xsin): added a test which includes failures of the
59 function itself and jacobian
61 2004-11-25 Brian Gough <bjg@network-theory.co.uk>
63 * gear2.c (gear2_alloc): added check for non-allocation of primer,
64 and two-step error estimation
66 2004-11-23 Brian Gough <bjg@network-theory.co.uk>
68 * gear1.c (gear1_apply): use 90% error bound
70 * rk2imp.c (rk2imp_apply): use 90% error bound
72 * rk2simp.c (rk2simp_apply): use 90% error bound
74 2004-11-18 Brian Gough <bjg@network-theory.co.uk>
76 * test.c: new test program (Tuomo Keskitalo)
78 2004-11-12 Brian Gough <bjg@network-theory.co.uk>
80 * evolve.c (gsl_odeiv_evolve_apply): check for internal stepper
81 failure (Tuomo Keskitalo)
83 * bsimp.c: gives exact dydt_out (Tuomo Keskitalo)
85 * rk2simp.c: new semi-implicit solver (Tuomo Keskitalo)
87 * gear1.c rkf45.c rkck.c rk8pd.c rk4imp.c rk2imp.c rk2.c rk4.c:
88 fix error estimate, exact derivatives on output (Tuomo Keskitalo)
90 2004-05-28 Brian Gough <bjg@network-theory.co.uk>
92 * bsimp.c (bsimp_apply): fix prototype of function to match
93 definition in typedef by removing const on double args.
95 Sat Aug 3 19:02:37 2002 Brian Gough <bjg@network-theory.co.uk>
97 * cscal.c (struct): use matlab style error control, with scalar
98 relative error and vector absolute error
100 Wed Oct 3 13:19:28 2001 Brian Gough <bjg@network-theory.co.uk>
102 * evolve.c (gsl_odeiv_evolve_apply): report error if step
103 direction does not match interval direction
105 * cstd.c (std_control_hadjust): allow for negative step-sizes
107 * evolve.c (gsl_odeiv_evolve_apply): allow for integrating
110 Sat Sep 29 21:02:36 2001 Brian Gough <bjg@network-theory.co.uk>
112 * bsimp.c: use the condition |D^-1 Delta| >> 1 to detect
113 singularities, as described in the paper and implemented in METAN1
116 * evolve.c (gsl_odeiv_evolve_apply): reset the final_step flag
117 whenever the step size is decreased
119 Sun Jul 1 22:42:02 2001 Brian Gough <bjg@network-theory.co.uk>
121 * bsimp.c: modified to use new-style vector views
123 Fri Jun 22 11:59:24 2001 Brian Gough <bjg@network-theory.co.uk>
125 * evolve.c: keep track of failed step count
127 Tue Jun 12 11:34:23 2001 Brian Gough <bjg@network-theory.co.uk>
129 * rkf45.c: added Runge-Kutta Fehlberg method, RKF45
131 * reorganized to use standard gsl conventions
133 Mon Apr 23 10:26:22 2001 Brian Gough <bjg@network-theory.co.uk>
135 * unified error handling conventions to _e for error handling
136 functions and no suffix for plain functions, so _impl functions
137 are no longer needed.
139 * removed tests for EFAULT, since EFAULT should only apply to
140 invalid non-null pointers.
142 2000-05-14 Steve Robbins <steve@nyongwa.montreal.qc.ca>
144 * Makefile.am (test_LDADD): put libgslpermutation ahead of
145 libgslvector in the link beccause the former uses
146 `gsl_check_range' which is defined in the latter.
148 Tue Apr 11 19:59:01 2000 Brian Gough <bjg@network-theory.co.uk>
150 * bsimp.c (bsimp_step_local): changed gsl_la prefix to gsl_linalg,
151 new naming convention
154 Fri Oct 1 15:46:13 1999 Brian Gough <bjg@network-theory.co.uk>
156 * bsimp.c: converted to use new-style block/vector code