[PATCH] Faster x^(2^m) for mpexp using __sqr

[Andreas Jaeger]

On Tuesday, February 12, 2013 20:04:42 Siddhesh Poyarekar wrote:
> Sorry, there was a stray bit of code that I forgot to remove in the
> previous patch.  Here's the updated patch; everything else remains the
> same.


Thanks, this is fine.

I noticed that there's some other code in dbl-64 that does squaring and 
could use __sqr instead of __mul, e.g.:
sysdeps/ieee754/dbl-64/sincos32.c:  __mul(x,x,&x2,p);
sysdeps/ieee754/dbl-64/mpatan.c:    __mul(x,x,&mpsm,p);
sysdeps/ieee754/dbl-64/mpsqrt.c:    __mul(&mpu,&mpu,&mpt1,p);

Could you check whether using __sqr in those - and perhaps other cases 
helps as well?

Andreas
-- 
 Andreas Jaeger aj@{suse.com,opensuse.org} Twitter/Identica: jaegerandi
  SUSE LINUX Products GmbH, Maxfeldstr. 5, 90409 Nürnberg, Germany
   GF: Jeff Hawn,Jennifer Guild,Felix Imendörffer,HRB16746 (AG Nürnberg)
    GPG fingerprint = 93A3 365E CE47 B889 DF7F  FED1 389A 563C C272 A126