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

