From: Brian Gough <bjg@network-theory.co.uk>
To: Juuso Heinisuo <heinisuo@pcu.helsinki.fi>
Cc: gsl-discuss@sources.redhat.com
Subject: Re: Is there a bug in gsl_odeiv?
Date: Wed, 19 Dec 2001 13:20:00 -0000 [thread overview]
Message-ID: <15330.54698.91242.909330@debian> (raw)
In-Reply-To: <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
----------------------------------------------------------------------
parent reply other threads:[~2001-12-19 13:20 UTC|newest]
Thread overview: expand[flat|nested] mbox.gz Atom feed
[parent not found: <200111021611.SAA27001@rock.helsinki.fi>]
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=15330.54698.91242.909330@debian \
--to=bjg@network-theory.co.uk \
--cc=gsl-discuss@sources.redhat.com \
--cc=heinisuo@pcu.helsinki.fi \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
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).