comp.lang.ada
 help / color / mirror / Atom feed
From: "John B. Matthews" <nospam@nospam.invalid>
Subject: Re: ANN: miscellaneous Math routines, GPL'd
Date: Wed, 13 May 2009 09:40:02 -0400
Date: 2009-05-13T09:40:02-04:00	[thread overview]
Message-ID: <nospam-576E21.09400213052009@news.aioe.org> (raw)
In-Reply-To: cb936fa3-df11-4563-b283-b0bfed0ba070@o27g2000vbd.googlegroups.com

In article 
<cb936fa3-df11-4563-b283-b0bfed0ba070@o27g2000vbd.googlegroups.com>,
 johnscpg@googlemail.com wrote:

> On May 12, 7:45 pm, "John B. Matthews" <nos...@nospam.invalid> wrote:
> > In article
> > <9dbcdcd0-3c15-4165-9336-3d5778df2...@e20g2000vbc.googlegroups.com>,
> >
> >  johns...@googlemail.com wrote:
> > > For those of you who like Ada with their numerics (is there
> > > anyone who doesn't?!) I have released a collection of math
> > > routines under the GPL license. Find them at:
> >
> > >  http://web.am.qub.ac.uk/users/j.parker/miscellany
> >
> > > The full set is tarred in the file:
> >
> > >   miscellany.10may09.tar.gz
> >
> > > in the directory given above.
> >
> > > Most are old classics I have found useful over the years
> > > (SVD, QR, LU, Runge-Kutta, FFT, Arbitrary precision floating point).
> > > The random number generators are very new; the documentation
> > > should explain why I recommend them. The random number generators
> > > and the Arbitrary precision floating point are designed to
> > > make good use of the new 64-bit CPU's.
> >
> > Excellent. You might like to see my experimental implementation of
> > Generic_Roots, using the Durand-Kerner-Weierstrass method and
> > distributed under the GNAT modified GPL:
> >
> > <http://home.roadrunner.com/~jbmatthews/misc/groots.html>
[...]
> Thanks for the pointer.  Root finding BTW is one of those problems 
> that cries out for extended precision floating pt.  Its rarely 
> time-sensitive (rarely in inner loops), but very sensitive to 
> precision ...  last time I did it, it was on mathematica in extended 
> precision .. I prefer Ada!  If I did in Fortran I'ld always use 
> REAL*16 (Quad-precision, 32 digit floats) but only one Fortran (intel 
> ifort) supports this as far as I know (except IBM's fortran on 
> powerpc).

Ah, yes. I noticed the difference running examples on PPC & x86. IIUC, 
the latter uses an 80-bit format internally.

> If I ever find the time I'll hook up your root finder to my 
> Extended Precision package.

Outstanding. As both are GMGPL, please feel free to incorporate any of 
my work in your project.

-- 
John B. Matthews
trashgod at gmail dot com
<http://sites.google.com/site/drjohnbmatthews>



  reply	other threads:[~2009-05-13 13:40 UTC|newest]

Thread overview: 9+ messages / expand[flat|nested]  mbox.gz  Atom feed  top
2009-05-12 16:38 ANN: miscellaneous Math routines, GPL'd johnscpg
2009-05-12 18:45 ` John B. Matthews
2009-05-13  9:39   ` johnscpg
2009-05-13 13:40     ` John B. Matthews [this message]
2009-05-12 19:03 ` Dmitry A. Kazakov
2009-05-12 21:02 ` lanceboyle
2009-05-13  9:25   ` johnscpg
2009-05-13 11:52     ` gautier_niouzes
2009-05-13 13:43       ` johnscpg
replies disabled

This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox