Update ENTERI() macro

Steve Kargl sgk at troutmask.apl.washington.edu
Fri Mar 8 03:04:20 UTC 2019


On Thu, Mar 07, 2019 at 04:22:45PM -0800, Steve Kargl wrote:
> 
> Make that 3 bugs.  
> 
> Index: ld80/k_expl.h
> ===================================================================
> --- ld80/k_expl.h	(revision 344600)
> +++ ld80/k_expl.h	(working copy)
> @@ -285,7 +285,7 @@
>  	y = cimagl(z);
>  	__k_expl(x, &hi, &lo, &k);
>  
> -	exp_x = (lo + hi) * 0x1p16382;
> +	exp_x = (lo + hi) * 0x1p16382L;
>  	expt += k - 16382;
>  
>  	scale1 = 1;
> @@ -292,9 +292,9 @@
>  	half_expt = expt / 2;
>  	SET_LDBL_EXPSIGN(scale1, BIAS + half_expt);
>  	scale2 = 1;
> -	SET_LDBL_EXPSIGN(scale1, BIAS + expt - half_expt);
> +	SET_LDBL_EXPSIGN(scale2, BIAS + expt - half_expt);

This bug appears in ld128/k_expl.h.  The other two bugs
are not present or were fixed at some point.

>  
> -	return (CMPLXL(cos(y) * exp_x * scale1 * scale2,
> +	return (CMPLXL(cosl(y) * exp_x * scale1 * scale2,
>  	    sinl(y) * exp_x * scale1 * scale2));
>  }
>  #endif /* _COMPLEX_H */
> 
> -- 
> Steve
> _______________________________________________
> freebsd-numerics at freebsd.org mailing list
> https://lists.freebsd.org/mailman/listinfo/freebsd-numerics
> To unsubscribe, send any mail to "freebsd-numerics-unsubscribe at freebsd.org"

-- 
Steve
20170425 https://www.youtube.com/watch?v=VWUpyCsUKR4
20161221 https://www.youtube.com/watch?v=IbCHE-hONow


More information about the freebsd-numerics mailing list