Re: bin/144306: Nasty bug in jn(3)

The following reply was made to PR bin/144306; it has been noted by GNATS.

From: "Steven G. Kargl" <kargl@troutmask.apl.washington.edu>
To: FreeBSD-gnats-submit@FreeBSD.org, freebsd-bugs@FreeBSD.org
Cc:  
Subject: Re: bin/144306: Nasty bug in jn(3)
Date: Thu, 25 Feb 2010 20:15:33 -0800 (PST)

 This patch fixes the problem.

 Index: e_jn.c
 ===================================================================
 --- e_jn.c	(revision 204219)
 +++ e_jn.c	(working copy)
 @@ -200,7 +200,12 @@
  			}
  	     	    }
  		}
 -	    	b = (t*__ieee754_j0(x)/b);
 +		z = __ieee754_j0(x);
 +		w = __ieee754_j1(x);
 +		if (fabs(z) >= fabs(w))
 +	    		b = (t*z/b);
 +		else
 +			b = (t*w/a);		
  	    }
  	}
  	if(sgn==1) return -b; else return b;

 -- 
 Steve
 http://troutmask.apl.washington.edu/~kargl/
_______________________________________________
freebsd-bugs@freebsd.org mailing list
http://lists.freebsd.org/mailman/listinfo/freebsd-bugs
To unsubscribe, send any mail to "freebsd-bugs-unsubscribe@freebsd.org"