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