Re: A perfect power, and then?

2012-10-27 Thread Torbjorn Granlund
ni...@lysator.liu.se (Niels Möller) writes: Perhaps you could add some of your remarks as a comment to the file? Makes sense. Not tonight, though. I saw your edits. And I edited them sligtly, adding some more FIXMEs and clarifying the operation of binv_root. Please read and tell me

Re: A perfect power, and then?

2012-10-27 Thread Niels Möller
Torbjorn Granlund t...@gmplib.org writes: I realise that Use a small table to get starting value might not be easy to implement for root since one might need k tables. Other possibilities would be: I think we discussed some months ago. IIRC, to get a starting value for a^{1/k}, it should

Re: A perfect power, and then?

2012-10-27 Thread Niels Möller
ni...@lysator.liu.se (Niels Möller) writes: broot.c:mpn_broot and your mpn_xxx that computes a^{1/n-1}, perhaps call the latter mpn_brootinvm1 I'll start with this one, then. Here's an initial patch, integrating my code from a few months ago. Some TODO items: 0. Support in