d/freebsd-dev
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Submitted by: J.T. Conklin <jtc@wimsey.com>

Browse Source 
First part of update to fdlibm 5.2: fix jn(n, x) and jnf(n, x).
jn(-1, x) was too large by a factor of 3.
This commit is contained in:
Bruce Evans 1995-04-07 23:13:43 +00:00
parent f9fc2cff0c
commit b29986e0f1
Notes: svn2git 2020-12-20 02:59:44 +00:00 
svn path=/head/; revision=7658
2 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
lib/msun/src/e_jn.c
View File

@ -11,7 +11,7 @@
*/
#ifndef lint
static char rcsid[] = "$Id: e_jn.c,v 1.6 1994/08/18 23:05:37 jtc Exp $";
static char rcsid[] = "$Id: e_jn.c,v 1.1.1.1 1994/08/19 09:39:44 jkh Exp $";
#endif
/*
@ -238,7 +238,7 @@ static double zero = 0.00000000000000000000e+00;
sign = 1;
if(n<0){
n = -n;
sign = 1 - ((n&1)<<2);
sign = 1 - ((n&1)<<1);
}
if(n==0) return(__ieee754_y0(x));
if(n==1) return(sign*__ieee754_y1(x));

4
lib/msun/src/e_jnf.c
View File

@ -14,7 +14,7 @@
*/
#ifndef lint
static char rcsid[] = "$Id: e_jnf.c,v 1.2 1994/08/18 23:05:39 jtc Exp $";
static char rcsid[] = "$Id: e_jnf.c,v 1.1.1.1 1994/08/19 09:39:55 jkh Exp $";
#endif
#include "math.h"
@ -192,7 +192,7 @@ static float zero = 0.0000000000e+00;
sign = 1;
if(n<0){
n = -n;
sign = 1 - ((n&1)<<2);
sign = 1 - ((n&1)<<1);
}
if(n==0) return(__ieee754_y0f(x));
if(n==1) return(sign*__ieee754_y1f(x));