|
Packit |
67cb25 |
2011-04-25 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* This module is now superseded by ode-initval2
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-08-12 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rk4.c (rk4_alloc): fix order of free statements to avoid memory
|
|
Packit |
67cb25 |
leak
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-07-09 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* step.c (gsl_odeiv_step_free): handle NULL argument in free
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c (gsl_odeiv_evolve_free): handle NULL argument in free
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* control.c (gsl_odeiv_control_free): handle NULL argument in free
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-04-21 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gear2.c (gear2_apply): save initial values of y for possible
|
|
Packit |
67cb25 |
failures
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-04-20 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c (gsl_odeiv_evolve_apply): restore original t value when
|
|
Packit |
67cb25 |
step fails
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-07-03 Brian Gough <bjg@hp2.network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* Makefile.am (INCLUDES): use top_srcdir instead of top_builddir
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-01-28 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test.c (test_evolve_negative_h): test evolution with negative
|
|
Packit |
67cb25 |
steps as well as positive
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-01-03 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c (gsl_odeiv_evolve_apply): handle the case where
|
|
Packit |
67cb25 |
GSL_ODEIV_HADJ_DEC does not represent a real decrease (e.g. where
|
|
Packit |
67cb25 |
x+h = x due to finite precision)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-07-29 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c (gsl_odeiv_evolve_apply): notify user of step-size
|
|
Packit |
67cb25 |
which caused any failure by returning it
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-06-13 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_odeiv.h, control.c (gsl_odeiv_control_hadjust): change y0 to
|
|
Packit |
67cb25 |
y in prototype to avoid conflict with bessel y0() in math.h
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test.c (test_odeiv_stepper): change y0 to ystart
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-04-08 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* Fixed all functions to consistently return status instead of
|
|
Packit |
67cb25 |
GSL_EBADFUNC.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-12-29 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test.c (test_compare_vanderpol): avoid variable size array
|
|
Packit |
67cb25 |
(test_compare_oregonator): as above
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-12-23 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test.c (test_evolve_xsin): changed description from "failures"
|
|
Packit |
67cb25 |
to "errors"
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-12-22 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rk2simp.c (rk2simp_step): added missing static declaration
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gear1.c (gear1_step): added missing static declaration
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gear2.c (gear2_step): added missing static declaration
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rk4imp.c (rk4imp_step): added missing static declaration
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rk4.c (rk4_step): added missing static declaration
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rk2imp.c (rk2imp_step): added missing static declaration
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-12-01 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* updated all solvers to ensure correct handling of error codes
|
|
Packit |
67cb25 |
from user defined functions
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test.c (rhs_xsin): added a test which includes failures of the
|
|
Packit |
67cb25 |
function itself and jacobian
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-11-25 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gear2.c (gear2_alloc): added check for non-allocation of primer,
|
|
Packit |
67cb25 |
and two-step error estimation
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-11-23 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gear1.c (gear1_apply): use 90% error bound
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rk2imp.c (rk2imp_apply): use 90% error bound
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rk2simp.c (rk2simp_apply): use 90% error bound
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-11-18 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test.c: new test program (Tuomo Keskitalo)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-11-12 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c (gsl_odeiv_evolve_apply): check for internal stepper
|
|
Packit |
67cb25 |
failure (Tuomo Keskitalo)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bsimp.c: gives exact dydt_out (Tuomo Keskitalo)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rk2simp.c: new semi-implicit solver (Tuomo Keskitalo)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gear1.c rkf45.c rkck.c rk8pd.c rk4imp.c rk2imp.c rk2.c rk4.c:
|
|
Packit |
67cb25 |
fix error estimate, exact derivatives on output (Tuomo Keskitalo)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-05-28 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bsimp.c (bsimp_apply): fix prototype of function to match
|
|
Packit |
67cb25 |
definition in typedef by removing const on double args.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sat Aug 3 19:02:37 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* cscal.c (struct): use matlab style error control, with scalar
|
|
Packit |
67cb25 |
relative error and vector absolute error
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Wed Oct 3 13:19:28 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c (gsl_odeiv_evolve_apply): report error if step
|
|
Packit |
67cb25 |
direction does not match interval direction
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* cstd.c (std_control_hadjust): allow for negative step-sizes
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c (gsl_odeiv_evolve_apply): allow for integrating
|
|
Packit |
67cb25 |
backwards in time
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sat Sep 29 21:02:36 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bsimp.c: use the condition |D^-1 Delta| >> 1 to detect
|
|
Packit |
67cb25 |
singularities, as described in the paper and implemented in METAN1
|
|
Packit |
67cb25 |
by Deuflhard et al.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c (gsl_odeiv_evolve_apply): reset the final_step flag
|
|
Packit |
67cb25 |
whenever the step size is decreased
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sun Jul 1 22:42:02 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bsimp.c: modified to use new-style vector views
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Fri Jun 22 11:59:24 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* evolve.c: keep track of failed step count
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Jun 12 11:34:23 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* rkf45.c: added Runge-Kutta Fehlberg method, RKF45
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* reorganized to use standard gsl conventions
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Mon Apr 23 10:26:22 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* unified error handling conventions to _e for error handling
|
|
Packit |
67cb25 |
functions and no suffix for plain functions, so _impl functions
|
|
Packit |
67cb25 |
are no longer needed.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* removed tests for EFAULT, since EFAULT should only apply to
|
|
Packit |
67cb25 |
invalid non-null pointers.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2000-05-14 Steve Robbins <steve@nyongwa.montreal.qc.ca>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* Makefile.am (test_LDADD): put libgslpermutation ahead of
|
|
Packit |
67cb25 |
libgslvector in the link beccause the former uses
|
|
Packit |
67cb25 |
`gsl_check_range' which is defined in the latter.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Apr 11 19:59:01 2000 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bsimp.c (bsimp_step_local): changed gsl_la prefix to gsl_linalg,
|
|
Packit |
67cb25 |
new naming convention
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Fri Oct 1 15:46:13 1999 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bsimp.c: converted to use new-style block/vector code
|
|
Packit |
67cb25 |
|