bin/144306: Nasty bug in jn(3)

[Steven G. Kargl]

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

From: "Steven G. Kargl" <kargl at troutmask.apl.washington.edu>
To: FreeBSD-gnats-submit at FreeBSD.org, freebsd-bugs at 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/