public inbox for gsl-discuss@sourceware.org
 help / color / mirror / Atom feed
From: "Robert G. Brown" <rgb@phy.duke.edu>
To: Gert Van den Eynde <gvdeynde@sckcen.be>
Cc: gsl-discuss@sources.redhat.com
Subject: Re: Laguerre Polynomial overflow...
Date: Wed, 02 Oct 2002 04:32:00 -0000	[thread overview]
Message-ID: <Pine.LNX.4.44.0210020711550.1527-100000@lilith.rgb.private.net> (raw)
In-Reply-To: <20021002090340.413f6a83.gvdeynde@sckcen.be>

On Wed, 2 Oct 2002, Gert Van den Eynde wrote:

> 
> Dear Robert,

> I cannot reproduce this using gsl-1.2 (my own compilation) on a SuSE
> 8.0 linux machine. Could you try to compile and run the attached program
> and see if that fails too?

Before starting to mess with the internal code I upgraded to the gsl 1.1
I had available in the RH 7.3 RPM set.  It already worked at that point.
Obviously little brownies were at work removing odd convergence/overflow
bugs in between...;-)

The moral, I suppose, is that I should make sure I'm using the latest
version before filing a bug report.  I'm a Bad Dog...

I appreciate it, though, and apologize for wasting your time.

I'll definitely compile and work through 1.2 before worrying about my
NEXT gsl problem -- the lack of (visible) adaptive stepsizes in the gsl
ode solvers.  Thus far I've had to solve on a very fine grid to get
decent accuracy (bound state eigensolutions integrated out from the
origin spend a bit of time crossing the axis where they oscillate and
eventually become stiff), but from the manual it looks like some
adaptive stepsize routines have already made their way into the later
revisions.  I can't wait to try them out.

BTW, the GSL is a fabulous toolset and I don't mean to sound critical at
all.  NR sucks (yes, I've read the rant, and then there is the
strangulation noose -- I mean "copyright").  I'd like to contribute, if
I ever have the chance.  For example, I have (somewhere, I'd have to dig
it out) a C translation of Forsythe, Malcom and Moler's "QUANC8" (8
panel Newton-Cotes quadrature) if that is of interest.  I'm not really a
computer scientist, but I've done a lot of numerical computations over
the last decade or two.  The other thing I might be able to contribute
is at least simple genetic optimization algorithms and (if I could
figure out how to frame them) neural network algorithms, as this is one
of the things I do.

   rgb

-- 
Robert G. Brown	                       http://www.phy.duke.edu/~rgb/
Duke University Dept. of Physics, Box 90305
Durham, N.C. 27708-0305
Phone: 1-919-660-2567  Fax: 919-660-2525     email:rgb@phy.duke.edu



  reply	other threads:[~2002-10-02 11:32 UTC|newest]

Thread overview: 5+ messages / expand[flat|nested]  mbox.gz  Atom feed  top
2002-10-01 11:40 Robert G. Brown
2002-10-02  0:03 ` Gert Van den Eynde
2002-10-02  4:32   ` Robert G. Brown [this message]
2002-10-02 11:49     ` Brian Gough
2002-10-02 12:32     ` Brian Gough

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=Pine.LNX.4.44.0210020711550.1527-100000@lilith.rgb.private.net \
    --to=rgb@phy.duke.edu \
    --cc=gsl-discuss@sources.redhat.com \
    --cc=gvdeynde@sckcen.be \
    /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).