From mboxrd@z Thu Jan  1 00:00:00 1970
X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on polar.synack.me
X-Spam-Level: 
X-Spam-Status: No, score=-1.9 required=5.0 tests=BAYES_00,FREEMAIL_FROM
	autolearn=ham autolearn_force=no version=3.4.4
X-Google-Thread: a07f3367d7,523c3848a9d03183
X-Google-Attributes: gida07f3367d7,public,usenet
X-Google-NewGroupId: yes
X-Google-Language: ENGLISH,ASCII
Path: 
 g2news2.google.com!postnews.google.com!o27g2000vbd.googlegroups.com!not-for-mail
From: johnscpg@googlemail.com
Newsgroups: comp.lang.ada
Subject: Re: ANN: miscellaneous Math routines, GPL'd
Date: Wed, 13 May 2009 02:39:50 -0700 (PDT)
Organization: http://groups.google.com
Message-ID: 
 <cb936fa3-df11-4563-b283-b0bfed0ba070@o27g2000vbd.googlegroups.com>
References: 
 <9dbcdcd0-3c15-4165-9336-3d5778df2ca3@e20g2000vbc.googlegroups.com>
	<nospam-51F0B4.14453112052009@news.aioe.org>
NNTP-Posting-Host: 143.117.23.126
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable
X-Trace: posting.google.com 1242207590 24598 127.0.0.1 (13 May 2009 09:39:50
 GMT)
X-Complaints-To: groups-abuse@google.com
NNTP-Posting-Date: Wed, 13 May 2009 09:39:50 +0000 (UTC)
Complaints-To: groups-abuse@google.com
Injection-Info: o27g2000vbd.googlegroups.com; posting-host=143.117.23.126;
	posting-account=Jzt5lQoAAAB4PhTgRLOPGuTLd_K1LY-C
User-Agent: G2/1.0
X-HTTP-UserAgent: Mozilla/5.0 (X11; U; Linux x86_64; en-GB; rv:1.9.0.9)
	Gecko/2009050519 Iceweasel/3.0.6 (Debian-3.0.6-1),gzip(gfe),gzip(gfe)
Xref: g2news2.google.com comp.lang.ada:5805
Date: 2009-05-13T02:39:50-07:00
List-Id: <comp.lang.ada>

On May 12, 7:45=A0pm, "John B. Matthews" <nos...@nospam.invalid> wrote:
> In article
> <9dbcdcd0-3c15-4165-9336-3d5778df2...@e20g2000vbc.googlegroups.com>,
>
>
>
> =A0johns...@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:
>
> > =A0http://web.am.qub.ac.uk/users/j.parker/miscellany
>
> > The full set is tarred in the file:
>
> > =A0 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