This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Is there a bug in gsl_odeiv?
- To: Juuso Heinisuo <heinisuo at pcu dot helsinki dot fi>
- Subject: Re: Is there a bug in gsl_odeiv?
- From: Brian Gough <bjg at network-theory dot co dot uk>
- Date: Fri, 2 Nov 2001 17:19:38 +0000 (GMT)
- CC: jungman at lanl dot gov, vuorinen at csc dot fi, arasila at csc dot fi
- CC: gsl-discuss at sources dot redhat dot com
- References: <200111021611.SAA27001@rock.helsinki.fi>
> Is there a bug in gsl_odeiv?
No, there is a bug in your program:
bjg|debian> diff -c ode.c.~1~ ode.c
*** ode.c.~1~ Fri Nov 2 16:37:45 2001
--- ode.c Fri Nov 2 17:04:17 2001
***************
*** 76,82 ****
{
f[0] = y[1];
f[1] = y[2];
! f[3] = 1.0;
return GSL_SUCCESS;
}
--- 76,82 ----
{
f[0] = y[1];
f[1] = y[2];
! f[2] = 1.0;
return GSL_SUCCESS;
}
You had me worrried for a few minutes there.
I recommend using 'checkergcc' to catch problems like that.
For bsimp the jacobian is required, which explains the second problem
you encountered. In fact, in general it is required (although not
always used). If there is somewhere in the manual that says otherwise
let me know.
--
Brian Gough
----------------------------------------------------------------------
Network Theory Ltd Phone: +44 (0)117 3179309
15 Royal Park WWW: http://www.network-theory.co.uk/
Clifton Email: bjg@network-theory.co.uk
Bristol BS8 3AL
----------------------------------------------------------------------