Bill Allombert on Fri, 07 Nov 2003 12:29:09 +0100 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: Experimenting with the Alpha build. |
On Fri, Nov 07, 2003 at 12:09:56PM +0100, Karim BELABAS wrote: > On Thu, 6 Nov 2003, Phil Carmody wrote: > > --- Karim BELABAS <Karim.Belabas@math.u-psud.fr> wrote: > >> On Thu, 6 Nov 2003, Phil Carmody wrote: > >>> I was wondering if -ffast-math would do to Pari/GP/GMP on the alpha. > >>> The speed-ups weren't significant on the whole (I should have tried > >>> -mno-ieee too perhaps) > >> > >> This is expected: there's almost no floating point code in PARI (everything > >> is emulated using 'unsigned long' arithmetic). Might be considered a > >> shortcoming. > >> > >> There are lots of variables of type 'double', but they are used in parts of > >> the code which is negligible timewise. And never where stability / accuracy > >> is important. I.e --fast-math should be a no-op with respect to performances > >> or results. > > > > I also recompiled GMP with -ffast-math too for the test. In theory GMP might > > perform some FP maths (FFT multiplication? or is it NTT?). > > NTT from a quick glance at the code. I don't think GMP uses doubles much, > either. mp_limb_t is always a variant of some unsigned integer (int, long, > long long). What is NTT ? GMP FFT use modular arithmetic, not double. Cheers, Bill.