* Re: Is there a bug in gsl_odeiv?
[not found] <200111021611.SAA27001@rock.helsinki.fi>
@ 2001-12-19 13:20 ` Brian Gough
0 siblings, 0 replies; only message in thread
From: Brian Gough @ 2001-12-19 13:20 UTC (permalink / raw)
To: Juuso Heinisuo; +Cc: gsl-discuss
> 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
----------------------------------------------------------------------
^ permalink raw reply [flat|nested] only message in thread
only message in thread, other threads:[~2001-12-19 13:20 UTC | newest]
Thread overview: (only message) (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
[not found] <200111021611.SAA27001@rock.helsinki.fi>
2001-12-19 13:20 ` Is there a bug in gsl_odeiv? Brian Gough
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).