On Wed, 10 Nov 2010, Bruce Evans wrote:
[...unrelated]
Just noticed a minor bug in the float version: from the PR: % diff --git a/lib/msun/src/e_jnf.c b/lib/msun/src/e_jnf.c % index 3bbf7b7..d045bb05 100644 % --- a/lib/msun/src/e_jnf.c % +++ b/lib/msun/src/e_jnf.c % @@ -152,7 +152,12 @@ __ieee754_jnf(int n, float x) % } % } % } % - b = (t*__ieee754_j0f(x)/b); % + z = __ieee754_j0f(x); % + w = __ieee754_j1f(x); % + if (fabs(z) >= fabs(w)) This should use fabsf(). % + b = (t*z/b); % + else % + b = (t*w/a); % } % } % if(sgn==1) return -b; else return b; % Bruce _______________________________________________ freebsd-bugs@freebsd.org mailing list http://lists.freebsd.org/mailman/listinfo/freebsd-bugs To unsubscribe, send any mail to "freebsd-bugs-unsubscr...@freebsd.org"
