🚀 go-pugleaf

Ken Laberteaux <klaberte@trc.tellabs.com> writes:

> I am reading the 2nd Edition of _Numerical Recipes in C_.  I am
> specifically looking at laguer.c.
> 
> On p. 373 the line of code which implements Eq 9.5.7 (p. 372) is 
> 
>                h=Csub(g2,RCmul(2.0,Cdiv(f,b)));
> 
> However Eq 9.5.7 (p. 372) makes no reference of the 2, indicating that
> the correct version should be
> 
>                h=Csub(g2,Cdiv(f,b));
> 
> Can anyone shed some light into this inconsistency?

Not I, but for those viewing via DejaNews, there is a second
(performance) bug in another NR code.  The arbitrary-precision
multiplication code should call fmod(), but instead uses a loop,
producing a drop in performance of at least thirty times (I measured
it).

This is just to keep all the grumbles in one place, you understand.

ObOnTopic: ... and the liccense for the NR code is too restrictive.

-- 
ACTUALLY reachable as @free-lunch.demon.(whitehouse)co.uk:james+usenet