Paul Zimmermann <Paul.Zimmermann@inria.fr> writes: > I guess it should be 1.0 instead. The whole expression should be evaluated > in double precision (z2 and the cp[] are doubles), and then only rounded to > float. Oh sure, and if it were any number other than ONE... ;-)