Re: ANN: miscellaneous Math routines, GPL'd

[johnscpg@googlemail.com]

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>
>
> --
> John B. Matthews
> trashgod at gmail dot com
> <http://sites.google.com/site/drjohnbmatthews>

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). If I ever find the time I'll hook
up your root finder to my Extended Precision package.

cheers
jonathan