Re: Division using FMAC, reciprocal estimates and Newton-Raphson - eg ia64, rs6000, SSE, ARM MaverickCrunch?

2008-05-12 Thread Andrew Haley
Paolo Bonzini wrote: I'd like to implement something similar for MaverickCrunch, using the integer 32-bit MAC functions, but there is no reciprocal estimate function on the MaverickCrunch. I guess a lookup table could be implemented, but how many entries will need to be generated, and how

Re: Division using FMAC, reciprocal estimates and Newton-Raphson - eg ia64, rs6000, SSE, ARM MaverickCrunch?

2008-05-11 Thread Hasjim Williams
On Sat, 10 May 2008 11:59:16 +0100, Martin Guy [EMAIL PROTECTED] said: On 5/9/08, Paolo Bonzini [EMAIL PROTECTED] wrote: The idea is to use integer arithmetic to compute the right exponent, and the lookup table to estimate the mantissa. The current soft-fp routine in libm seems to use a

Re: Division using FMAC, reciprocal estimates and Newton-Raphson - eg ia64, rs6000, SSE, ARM MaverickCrunch?

2008-05-10 Thread Martin Guy
On 5/9/08, Paolo Bonzini [EMAIL PROTECTED] wrote: The idea is to use integer arithmetic to compute the right exponent, and the lookup table to estimate the mantissa. I used something like this for square root: 1) shift the entire FP number by 1 to the right (logical right shift) 2) sum

Re: Division using FMAC, reciprocal estimates and Newton-Raphson - eg ia64, rs6000, SSE, ARM MaverickCrunch?

2008-05-10 Thread Andrew Haley
Paolo Bonzini wrote: I'd like to implement something similar for MaverickCrunch, using the integer 32-bit MAC functions, but there is no reciprocal estimate function on the MaverickCrunch. I guess a lookup table could be implemented, but how many entries will need to be generated, and how

Re: Division using FMAC, reciprocal estimates and Newton-Raphson - eg ia64, rs6000, SSE, ARM MaverickCrunch?

2008-05-10 Thread Andrew Haley
Andrew Haley wrote: Paolo Bonzini wrote: I'd like to implement something similar for MaverickCrunch, using the integer 32-bit MAC functions, but there is no reciprocal estimate function on the MaverickCrunch. I guess a lookup table could be implemented, but how many entries will need to be

Division using FMAC, reciprocal estimates and Newton-Raphson - eg ia64, rs6000, SSE, ARM MaverickCrunch?

2008-05-08 Thread Hasjim Williams
Hi all, I was looking for ways to improve the MaverickCrunch division routine on ARM ep93xx, and noticed that there are few other architectures that don't have a hardware divide. IA-64 has a frcpa instruction that returns an estimate of the reciprocal of a float or double. Likewise, RS-6000 has