whats happening ?
Torbjorn Granlund
tege@swox.com
23 Oct 2002 10:04:18 +0200
Jason Moxham <J.L.Moxham@maths.soton.ac.uk> writes:
On Tuesday 22 Oct 2002 8:25 pm, Torbjorn Granlund wrote:
> Jason Moxham <J.L.Moxham@maths.soton.ac.uk> writes:
>
> I done some preliminary work on a few algorithms ,and even in a
> completely un-optimized form have got some good results for some of the
> following functions , other results , not so good :(
>
> mpz_root using varible precision , division free newton method
>
> What is needed is mpn_rootrem with these properties. mpz_root
> should not change.
The varible precision root method I implemented has proven bounds , so cannot
ever be in error , I assume this is your objection , if not , what is it ?
Well, we're clearly in misunderstanding.
No, I don't object to code that is correct. As a matter of fact,
I much prefer correct code to incorrect code. It is a very
reasonable misunderstanding that you thought I liked incorrect
code. ;-)
My point was that replacing mpz_root isn't the right change since
the real root extraction is performed at the mpn level. In the
development sources, there are several call sites for
mpn_rootrem.
--
Torbjörn