Bill Allombert on Tue, 05 Feb 2008 10:41:09 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: integer n-roots like sqrtint()


On Tue, Feb 05, 2008 at 08:31:27AM +0100, Jeroen Demeyer wrote:
> Karim Belabas wrote:
> >Yes, since perfect powers are the bad case, you could use
> >
> >  if (!ispower(x, n, &z), z = floor(sqrtn(x, n)));
> >  z;
> 
> Yes, I guess that works (at least if the realprecision is sufficiently 
> high).
> 
> I'm wondering if this can be considered a bug in PARI.  Do you think it 

Well, sqrtn could be slightly more accurate since 6. is exactly
representable.

> makes sense to change the implementation of sqrtn() to first call 
> ispower() when the argument is an exact integer?

If we improve sqrtn accuracy we should do it for all input, not just
perfect powers.

Cheers,
Bill.