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 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,UTF8 Path: g2news2.google.com!news2.google.com!border1.nntp.dca.giganews.com!nntp.giganews.com!newsfeed00.sul.t-online.de!newsfeed01.sul.t-online.de!t-online.de!newsfeed.straub-nv.de!de-l.enfer-du-nord.net!gegeweb.org!aioe.org!nospam From: "John B. Matthews" Newsgroups: comp.lang.ada Subject: Re: ANN: miscellaneous Math routines, GPL'd Date: Wed, 13 May 2009 09:40:02 -0400 Organization: The Wasteland Message-ID: References: <9dbcdcd0-3c15-4165-9336-3d5778df2ca3@e20g2000vbc.googlegroups.com> NNTP-Posting-Host: ib4TTflHUauJidfWP/+Rjw.user.aioe.org Mime-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit X-Complaints-To: abuse@aioe.org X-Notice: Filtered by postfilter v. 0.7.9 Cancel-Lock: sha1:xMZOiucpC09uY7B+FlkGNhzBQU0= User-Agent: MT-NewsWatcher/3.5.3b3 (Intel Mac OS X) Xref: g2news2.google.com comp.lang.ada:5808 Date: 2009-05-13T09:40:02-04:00 List-Id: In article , johnscpg@googlemail.com wrote: > On May 12, 7:45 pm, "John B. Matthews" 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: > > > > [...] > 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