1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
|
2011-04-25 Brian Gough <bjg@network-theory.co.uk>
* This module is now superseded by ode-initval2
2009-08-12 Brian Gough <bjg@network-theory.co.uk>
* rk4.c (rk4_alloc): fix order of free statements to avoid memory
leak
2009-07-09 Brian Gough <bjg@network-theory.co.uk>
* step.c (gsl_odeiv_step_free): handle NULL argument in free
* evolve.c (gsl_odeiv_evolve_free): handle NULL argument in free
* control.c (gsl_odeiv_control_free): handle NULL argument in free
2009-04-21 Brian Gough <bjg@network-theory.co.uk>
* gear2.c (gear2_apply): save initial values of y for possible
failures
2009-04-20 Brian Gough <bjg@network-theory.co.uk>
* evolve.c (gsl_odeiv_evolve_apply): restore original t value when
step fails
2008-07-03 Brian Gough <bjg@hp2.network-theory.co.uk>
* Makefile.am (INCLUDES): use top_srcdir instead of top_builddir
2008-01-28 Brian Gough <bjg@network-theory.co.uk>
* test.c (test_evolve_negative_h): test evolution with negative
steps as well as positive
2008-01-03 Brian Gough <bjg@network-theory.co.uk>
* evolve.c (gsl_odeiv_evolve_apply): handle the case where
GSL_ODEIV_HADJ_DEC does not represent a real decrease (e.g. where
x+h = x due to finite precision)
2006-07-29 Brian Gough <bjg@network-theory.co.uk>
* evolve.c (gsl_odeiv_evolve_apply): notify user of step-size
which caused any failure by returning it
2006-06-13 Brian Gough <bjg@network-theory.co.uk>
* gsl_odeiv.h, control.c (gsl_odeiv_control_hadjust): change y0 to
y in prototype to avoid conflict with bessel y0() in math.h
* test.c (test_odeiv_stepper): change y0 to ystart
2005-04-08 Brian Gough <bjg@network-theory.co.uk>
* Fixed all functions to consistently return status instead of
GSL_EBADFUNC.
2004-12-29 Brian Gough <bjg@network-theory.co.uk>
* test.c (test_compare_vanderpol): avoid variable size array
(test_compare_oregonator): as above
2004-12-23 Brian Gough <bjg@network-theory.co.uk>
* test.c (test_evolve_xsin): changed description from "failures"
to "errors"
2004-12-22 Brian Gough <bjg@network-theory.co.uk>
* rk2simp.c (rk2simp_step): added missing static declaration
* gear1.c (gear1_step): added missing static declaration
* gear2.c (gear2_step): added missing static declaration
* rk4imp.c (rk4imp_step): added missing static declaration
* rk4.c (rk4_step): added missing static declaration
* rk2imp.c (rk2imp_step): added missing static declaration
2004-12-01 Brian Gough <bjg@network-theory.co.uk>
* updated all solvers to ensure correct handling of error codes
from user defined functions
* test.c (rhs_xsin): added a test which includes failures of the
function itself and jacobian
2004-11-25 Brian Gough <bjg@network-theory.co.uk>
* gear2.c (gear2_alloc): added check for non-allocation of primer,
and two-step error estimation
2004-11-23 Brian Gough <bjg@network-theory.co.uk>
* gear1.c (gear1_apply): use 90% error bound
* rk2imp.c (rk2imp_apply): use 90% error bound
* rk2simp.c (rk2simp_apply): use 90% error bound
2004-11-18 Brian Gough <bjg@network-theory.co.uk>
* test.c: new test program (Tuomo Keskitalo)
2004-11-12 Brian Gough <bjg@network-theory.co.uk>
* evolve.c (gsl_odeiv_evolve_apply): check for internal stepper
failure (Tuomo Keskitalo)
* bsimp.c: gives exact dydt_out (Tuomo Keskitalo)
* rk2simp.c: new semi-implicit solver (Tuomo Keskitalo)
* gear1.c rkf45.c rkck.c rk8pd.c rk4imp.c rk2imp.c rk2.c rk4.c:
fix error estimate, exact derivatives on output (Tuomo Keskitalo)
2004-05-28 Brian Gough <bjg@network-theory.co.uk>
* bsimp.c (bsimp_apply): fix prototype of function to match
definition in typedef by removing const on double args.
Sat Aug 3 19:02:37 2002 Brian Gough <bjg@network-theory.co.uk>
* cscal.c (struct): use matlab style error control, with scalar
relative error and vector absolute error
Wed Oct 3 13:19:28 2001 Brian Gough <bjg@network-theory.co.uk>
* evolve.c (gsl_odeiv_evolve_apply): report error if step
direction does not match interval direction
* cstd.c (std_control_hadjust): allow for negative step-sizes
* evolve.c (gsl_odeiv_evolve_apply): allow for integrating
backwards in time
Sat Sep 29 21:02:36 2001 Brian Gough <bjg@network-theory.co.uk>
* bsimp.c: use the condition |D^-1 Delta| >> 1 to detect
singularities, as described in the paper and implemented in METAN1
by Deuflhard et al.
* evolve.c (gsl_odeiv_evolve_apply): reset the final_step flag
whenever the step size is decreased
Sun Jul 1 22:42:02 2001 Brian Gough <bjg@network-theory.co.uk>
* bsimp.c: modified to use new-style vector views
Fri Jun 22 11:59:24 2001 Brian Gough <bjg@network-theory.co.uk>
* evolve.c: keep track of failed step count
Tue Jun 12 11:34:23 2001 Brian Gough <bjg@network-theory.co.uk>
* rkf45.c: added Runge-Kutta Fehlberg method, RKF45
* reorganized to use standard gsl conventions
Mon Apr 23 10:26:22 2001 Brian Gough <bjg@network-theory.co.uk>
* unified error handling conventions to _e for error handling
functions and no suffix for plain functions, so _impl functions
are no longer needed.
* removed tests for EFAULT, since EFAULT should only apply to
invalid non-null pointers.
2000-05-14 Steve Robbins <steve@nyongwa.montreal.qc.ca>
* Makefile.am (test_LDADD): put libgslpermutation ahead of
libgslvector in the link beccause the former uses
`gsl_check_range' which is defined in the latter.
Tue Apr 11 19:59:01 2000 Brian Gough <bjg@network-theory.co.uk>
* bsimp.c (bsimp_step_local): changed gsl_la prefix to gsl_linalg,
new naming convention
Fri Oct 1 15:46:13 1999 Brian Gough <bjg@network-theory.co.uk>
* bsimp.c: converted to use new-style block/vector code
|
