2011-04-25 Brian Gough * This module is now superseded by ode-initval2 2009-08-12 Brian Gough * rk4.c (rk4_alloc): fix order of free statements to avoid memory leak 2009-07-09 Brian Gough * 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 * gear2.c (gear2_apply): save initial values of y for possible failures 2009-04-20 Brian Gough * evolve.c (gsl_odeiv_evolve_apply): restore original t value when step fails 2008-07-03 Brian Gough * Makefile.am (INCLUDES): use top_srcdir instead of top_builddir 2008-01-28 Brian Gough * test.c (test_evolve_negative_h): test evolution with negative steps as well as positive 2008-01-03 Brian Gough * 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 * evolve.c (gsl_odeiv_evolve_apply): notify user of step-size which caused any failure by returning it 2006-06-13 Brian Gough * 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 * Fixed all functions to consistently return status instead of GSL_EBADFUNC. 2004-12-29 Brian Gough * test.c (test_compare_vanderpol): avoid variable size array (test_compare_oregonator): as above 2004-12-23 Brian Gough * test.c (test_evolve_xsin): changed description from "failures" to "errors" 2004-12-22 Brian Gough * 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 * 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 * gear2.c (gear2_alloc): added check for non-allocation of primer, and two-step error estimation 2004-11-23 Brian Gough * 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 * test.c: new test program (Tuomo Keskitalo) 2004-11-12 Brian Gough * 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 * 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 * 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 * 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 * 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 * bsimp.c: modified to use new-style vector views Fri Jun 22 11:59:24 2001 Brian Gough * evolve.c: keep track of failed step count Tue Jun 12 11:34:23 2001 Brian Gough * rkf45.c: added Runge-Kutta Fehlberg method, RKF45 * reorganized to use standard gsl conventions Mon Apr 23 10:26:22 2001 Brian Gough * 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 * 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 * bsimp.c (bsimp_step_local): changed gsl_la prefix to gsl_linalg, new naming convention Fri Oct 1 15:46:13 1999 Brian Gough * bsimp.c: converted to use new-style block/vector code