August 2025 - cmucl-cvs - mailman3.common-lisp.net

[Git][cmucl/cmucl][issue-431-fix-set-floating-point-modes] 5 commits: Add some more tests for set-floating-point-modes
by Raymond Toy (＠rtoy) 16 Aug '25

16 Aug '25

Raymond Toy pushed to branch issue-431-fix-set-floating-point-modes at cmucl / cmucl Commits: e80d2fe4 by Raymond Toy at 2025-08-15T15:15:00-07:00 Add some more tests for set-floating-point-modes - - - - - 1dfdaba2 by Raymond Toy at 2025-08-15T19:17:20-07:00 Rename things and update corresponding uses. We don't need x87-exceptions-byte anymore because it's the same as float-exceptions-byte. Rename the defconstants to include the word "float". Also make the x87 constant names more like the corresponding names in parms.lisp, which are basically for sse2. - - - - - 3f9a71d1 by Raymond Toy at 2025-08-16T07:52:41-07:00 Renamed some x87 FP constants x87-exceptions-mask-byte -> x87-float-traps-byte x87-rounding-control-byte -> x87-float-rounding-mode These new names match the existing float-traps-byte and float-rounding-mode constants that exist in parms.lisp. - - - - - a1adbe86 by Raymond Toy at 2025-08-16T07:54:15-07:00 rename x87-traps-byte to x87-float-traps-byte - - - - - 0ac12cb6 by Raymond Toy at 2025-08-16T07:54:43-07:00 Update pot files for updated docstrings - - - - - 7 changed files: - src/bootfiles/21e/boot-2025-07.lisp - src/code/float-trap.lisp - src/compiler/x86/float-sse2.lisp - src/compiler/x86/parms.lisp - src/i18n/locale/cmucl-sse2.pot - src/i18n/locale/cmucl-x86-vm.pot - tests/float.lisp Changes: ===================================== src/bootfiles/21e/boot-2025-07.lisp ===================================== @@ -2,18 +2,24 @@ (in-package :x86) -(defconstant x87-exceptions-byte - (byte 6 0)) -(defconstant x87-infinity-control-byte - (byte 1 (+ 12 16))) -(defconstant x87-rounding-control-byte - (byte 2 (+ 10 16))) -(defconstant x87-precision-control-byte - (byte 2 (+ 8 16))) -(defconstant x87-exceptions-mask-byte - (byte 6 16)) -(defconstant x87-precision-control-alist +(defconstant x87-float-infinity-control-byte + (byte 1 (+ 12 16)) + "The bit in the x87 FPU control word that controls the infinity mode.") +(defconstant x87-float-rounding-mode + (byte 2 (+ 10 16)) + "The bits in the x87 FPU control word for the rounding mode.") +(defconstant x87-float-precision-control-byte + (byte 2 (+ 8 16)) + "The bits in the x87 FPU contol word for the FP operation precision.") +(defconstant x87-float-traps-byte + (byte 6 16) + "The bits in the x87 FPU control word indicating the exceptions that + are enabled.") +(defconstant x87-float-precision-control-alist `((:24-bits . 0) (:reserved . 1) (:53-bits . 2) - (:64-bits . 3))) + (:64-bits . 3)) + "Alist for the x87 precison control. The car is the symbolic + precision and the cdr is the value for the precision control field.") + ===================================== src/code/float-trap.lisp ===================================== @@ -139,16 +139,15 @@ (new-exceptions (ldb float-exceptions-byte new-mode)) (new-enables (ldb float-traps-byte new-mode)) (x87-modes 0)) - ;; From the new mode, update the x87 FPU mode with the same set - ;; of accrued exceptions, enabled exceptions, and rounding mode. - (setf x87-modes (dpb new-exceptions x87-exceptions-byte x87-modes)) - (setf x87-modes (dpb new-enables x87-exceptions-mask-byte x87-modes)) - (setf x87-modes (dpb rc x87-rounding-control-byte x87-modes)) + ;; From the new mode, update the x87 FPU mode with the same set of accrued exceptions, enabled exceptions, and rounding mode. + (setf x87-modes (dpb new-exceptions float-exceptions-byte x87-modes)) + (setf x87-modes (dpb new-enables x87-float-traps-byte x87-modes)) + (setf x87-modes (dpb rc x87-float-rounding-mode x87-modes)) ;; Set precision to 64-bits (double-extended). (setf x87-modes - (dpb (cdr (assoc :64-bits x87-precision-control-alist)) - x87-precision-control-byte + (dpb (cdr (assoc :64-bits x87-float-precision-control-alist)) + x87-float-precision-control-byte x87-modes)) (setf (vm::sse2-floating-point-modes) (ldb (byte 16 0) new-mode)) ===================================== src/compiler/x86/float-sse2.lisp ===================================== @@ -1379,21 +1379,26 @@ ;; 10 double precision (53 bits) ;; 11 double extended precision (64 bits) -(defconstant x87-exceptions-byte - (byte 6 0)) -(defconstant x87-infinity-control-byte - (byte 1 (+ 12 16))) -(defconstant x87-rounding-control-byte - (byte 2 (+ 10 16))) -(defconstant x87-precision-control-byte - (byte 2 (+ 8 16))) -(defconstant x87-exceptions-mask-byte - (byte 6 16)) -(defconstant x87-precision-control-alist +(defconstant x87-float-infinity-control-byte + (byte 1 (+ 12 16)) + "The bit in the x87 FPU control word that controls the infinity mode.") +(defconstant x87-float-rounding-mode + (byte 2 (+ 10 16)) + "The bits in the x87 FPU control word for the rounding mode.") +(defconstant x87-float-precision-control-byte + (byte 2 (+ 8 16)) + "The bits in the x87 FPU contol word for the FP operation precision.") +(defconstant x87-float-traps-byte + (byte 6 16) + "The bits in the x87 FPU control word indicating the exceptions that + are enabled.") +(defconstant x87-float-precision-control-alist `((:24-bits . 0) (:reserved . 1) (:53-bits . 2) - (:64-bits . 3))) + (:64-bits . 3)) + "Alist for the x87 precison control. The car is the symbolic + precision and the cdr is the value for the precision control field.") (defun print-fp-exceptions-enabled (enabled) (format t "Precision enable: ~30T~A~%" (ldb (byte 1 5) enabled)) @@ -1416,7 +1421,7 @@ ;; Note that Intel uses masks to disable the exception, but to match ;; the rest of cmucl, these bits are represented as enable bits. (format t "Flush-to-zero: ~30T~A~%" (ldb (byte 1 15) sse-mode)) - (let ((rc (ldb (byte 2 13) sse-mode))) + (let ((rc (ldb float-rounding-mode sse-mode))) (format t "Rounding control: ~30T~A ~A~%" rc (car (rassoc rc rounding-mode-alist)))) @@ -1437,17 +1442,17 @@ (format t "SW: Stack fault: ~30T~A~%" (ldb (byte 1 6) x87-mode)) (print-fp-current-exceptions (ldb float-exceptions-byte x87-mode)) (format t "CW: Reserved: ~30T~A~%" (ldb (byte 2 (+ 13 16)) x87-mode)) - (format t "CW: Infinity control: ~30T~A~%" (ldb (byte 1 (+ 12 16)) x87-mode)) - (let ((rc (ldb (byte 2 (+ 10 16)) x87-mode))) + (format t "CW: Infinity control: ~30T~A~%" (ldb x87-float-infinity-control-byte x87-mode)) + (let ((rc (ldb x87-float-rounding-mode x87-mode))) (format t "CW: Rounding control: ~30T~A ~A~%" rc (car (rassoc rc rounding-mode-alist)))) - (let ((pc (ldb (byte 2 (+ 8 16)) x87-mode))) + (let ((pc (ldb x87-float-precision-control-byte x87-mode))) (format t "CW: Precision control: ~30T~A ~A~%" pc - (car (rassoc pc x87-precision-control-alist)))) + (car (rassoc pc x87-float-precision-control-alist)))) (format t "CW: Reserved: ~30T~A~%" (ldb (byte 2 (+ 6 16)) x87-mode)) - (print-fp-exceptions-enabled (ldb x87-exceptions-mask-byte x87-mode))) + (print-fp-exceptions-enabled (ldb x87-float-traps-byte x87-mode))) (defknown x87-floating-point-modes () float-modes (flushable)) (defknown ((setf x87-floating-point-modes)) (float-modes) ===================================== src/compiler/x86/parms.lisp ===================================== @@ -205,9 +205,15 @@ ;; need to do it this way because the interface assumes the modes are ;; in the same order as the MXCSR register. (defconstant float-rounding-mode (byte 2 13)) -(defconstant float-sticky-bits (byte 6 0)) -(defconstant float-traps-byte (byte 6 7)) -(defconstant float-exceptions-byte (byte 6 0)) +(defconstant float-sticky-bits (byte 6 0) + "The bits in the FP mode that hold the accrued exceptions that have + occurred since the bits were reset.") +(defconstant float-traps-byte (byte 6 7) + "The bits in the FP mode that hold FP exceptions that are enabled.") +(defconstant float-exceptions-byte (byte 6 0) + "The bits in the FP mode that hold the current exception. However + for x86, there aren't separate bits for this, so we use the stick + bits from the accrued exceptions") (progn ;; SSE2 has a flush-to-zero flag, which we use as the fast bit. Some ===================================== src/i18n/locale/cmucl-sse2.pot ===================================== @@ -86,6 +86,30 @@ msgstr "" msgid "inline ftruncate" msgstr "" +#: src/compiler/x86/float-sse2.lisp +msgid "The bit in the x87 FPU control word that controls the infinity mode." +msgstr "" + +#: src/compiler/x86/float-sse2.lisp +msgid "The bits in the x87 FPU control word for the rounding mode." +msgstr "" + +#: src/compiler/x86/float-sse2.lisp +msgid "The bits in the x87 FPU contol word for the FP operation precision." +msgstr "" + +#: src/compiler/x86/float-sse2.lisp +msgid "" +"The bits in the x87 FPU control word indicating the exceptions that\n" +" are enabled." +msgstr "" + +#: src/compiler/x86/float-sse2.lisp +msgid "" +"Alist for the x87 precison control. The car is the symbolic\n" +" precision and the cdr is the value for the precision control field." +msgstr "" + #: src/compiler/x86/float-sse2.lisp msgid "Print SSE2 modes word in a readable fashion." msgstr "" ===================================== src/i18n/locale/cmucl-x86-vm.pot ===================================== @@ -136,6 +136,23 @@ msgstr "" msgid "Maximum number of bits in a positive fixnum" msgstr "" +#: src/compiler/x86/parms.lisp +msgid "" +"The bits in the FP mode that hold the accrued exceptions that have\n" +" occurred since the bits were reset." +msgstr "" + +#: src/compiler/x86/parms.lisp +msgid "The bits in the FP mode that hold FP exceptions that are enabled." +msgstr "" + +#: src/compiler/x86/parms.lisp +msgid "" +"The bits in the FP mode that hold the current exception. However\n" +" for x86, there aren't separate bits for this, so we use the stick\n" +" bits from the accrued exceptions" +msgstr "" + #: src/compiler/x86/vm.lisp msgid "Redefining SC number ~D from ~S to ~S." msgstr "" ===================================== tests/float.lisp ===================================== @@ -343,3 +343,48 @@ (assert-true (typep new-mode 'x86::float-modes)) (assert-equal new-mode (setf (x86::x87-floating-point-modes) new-mode)))) +(define-test set-floating-point-modes + (let ((old-x87-modes (x86::x87-floating-point-modes)) + (old-sse2-modes (x86::sse2-floating-point-modes))) + (unwind-protect + (progn + ;; Set some new traps and rounding mode + (ext:set-floating-point-modes :traps '(:underflow + :overflow + :invalid + :divide-by-zero) + :rounding-mode :zero) + (let* ((x87-modes (x86::x87-floating-point-modes)) + (x87-exceptions-enabled (ldb x86::x87-float-traps-byte x87-modes)) + (x87-rc (ldb x86::x87-float-rounding-mode x87-modes)) + (sse2-modes (x86::sse2-floating-point-modes)) + (sse2-exceptions-enabled (ldb x86::float-traps-byte sse2-modes)) + (sse2-rc (ldb x86::float-rounding-mode sse2-modes))) + (x86::print-x87-fp-modes x87-modes) + ;; Verify that we set the enabled exceptions + ;; correctly. First for sse2, then for x87. + (assert-false (logbitp 5 sse2-exceptions-enabled)) ; precision + (assert-true (logbitp 4 sse2-exceptions-enabled)) ; underflow + (assert-true (logbitp 3 sse2-exceptions-enabled)) ; overflow + (assert-true (logbitp 2 sse2-exceptions-enabled)) ; divide-by-zero + (assert-false (logbitp 1 sse2-exceptions-enabled)) ; denormal + (assert-true (logbitp 0 sse2-exceptions-enabled)) ; invalid + + (assert-false (logbitp 5 x87-exceptions-enabled)) ; precision + (assert-true (logbitp 4 x87-exceptions-enabled)) ; underflow + (assert-true (logbitp 3 x87-exceptions-enabled)) ; overflow + (assert-true (logbitp 2 x87-exceptions-enabled)) ; divide-by-zero + (assert-false (logbitp 1 x87-exceptions-enabled)) ; denormal + (assert-true (logbitp 0 x87-exceptions-enabled)) ; invalid + + ;; Verify the rounding mode is set to zero + (assert-eql :zero (car (rassoc sse2-rc x86::rounding-mode-alist))) + (assert-eql :zero (car (rassoc x87-rc x86::rounding-mode-alist))) + + ;; Verify precision for x87 + (assert-eql :64-bits + (car (rassoc (ldb x86::x87-float-precision-control-byte x87-modes) + x86::x87-float-precision-control-alist))))) + + (setf (x86::x87-floating-point-modes) old-x87-modes) + (setf (x86::sse2-floating-point-modes) old-sse2-modes)))) View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/c370c74e74004ba8d25a09… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/c370c74e74004ba8d25a09… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-431-fix-set-floating-point-modes] Forgot to add bootstrap for CI
by Raymond Toy (＠rtoy) 15 Aug '25

15 Aug '25

Raymond Toy pushed to branch issue-431-fix-set-floating-point-modes at cmucl / cmucl Commits: c370c74e by Raymond Toy at 2025-08-15T14:16:12-07:00 Forgot to add bootstrap for CI Needed to build the code so that the new x87 defconstants are available when compile `set-floating-point-modes`. - - - - - 1 changed file: - .gitlab-ci.yml Changes: ===================================== .gitlab-ci.yml ===================================== @@ -4,7 +4,7 @@ variables: download_url: "https://common-lisp.net/project/cmucl/downloads/snapshots/$year/$month" version: "$year-$month-x86" tar_ext: "xz" - bootstrap: "" + bootstrap: "-B boot-2025-07" # Default install configuration to download the cmucl tarballs to use # for building. View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/c370c74e74004ba8d25a099… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/c370c74e74004ba8d25a099… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl] Pushed new branch issue-431-fix-set-floating-point-modes
by Raymond Toy (＠rtoy) 15 Aug '25

15 Aug '25

Raymond Toy pushed new branch issue-431-fix-set-floating-point-modes at cmucl / cmucl -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/tree/issue-431-fix-set-floatin… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][master] Update release notes for #426
by Raymond Toy (＠rtoy) 14 Aug '25

14 Aug '25

Raymond Toy pushed to branch master at cmucl / cmucl Commits: eb98570f by Raymond Toy at 2025-08-14T07:55:12-07:00 Update release notes for #426 [skip-ci] - - - - - 1 changed file: - src/general-info/release-21f.md Changes: ===================================== src/general-info/release-21f.md ===================================== @@ -138,6 +138,7 @@ public domain. * #418 Update asdf to version 3.3.7.4 (for package-local-nicknames) * #424 Use fdlibm `hypot` to fix bug in snapshot 2025-07 + * #426 Define float-modes type correctly for `(setf (x87-set-floating-point-modes))` * Other changes: * Improvements to the PCL implementation of CLOS: * Changes to building procedure: View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/eb98570f9cc44474f64dae6… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/eb98570f9cc44474f64dae6… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][master] 2 commits: Fix #426: Define float-modes type correctly
by Raymond Toy (＠rtoy) 14 Aug '25

14 Aug '25

Raymond Toy pushed to branch master at cmucl / cmucl Commits: d2c96f1b by Raymond Toy at 2025-08-14T07:45:18-07:00 Fix #426: Define float-modes type correctly - - - - - 2db294b1 by Raymond Toy at 2025-08-14T07:45:18-07:00 Merge branch 'issue-426-fix-float-modes-type-def' into 'master' Fix #426: Define float-modes type correctly Closes #426 See merge request cmucl/cmucl!313 - - - - - 2 changed files: - src/compiler/x86/float-sse2.lisp - tests/float.lisp Changes: ===================================== src/compiler/x86/float-sse2.lisp ===================================== @@ -1265,7 +1265,7 @@ ;;;; Float mode hackery: -(deftype float-modes () '(unsigned-byte 24)) +(deftype float-modes () '(unsigned-byte 30)) ;; For the record, here is the format of the MXCSR register. ;; ===================================== tests/float.lisp ===================================== @@ -333,3 +333,13 @@ (declare (ignore c)) (values (invoke-restart 'lisp::largest-float))))) (read-from-string string)))))) + +#+x86 +(define-test x87-set-floating-point-modes + (:tag :issues) + (let ((new-mode (setf (x86::x87-floating-point-modes) + (dpb 2 (byte 2 24) + (x86::x87-floating-point-modes))))) + (assert-true (typep new-mode 'x86::float-modes)) + (assert-equal new-mode (setf (x86::x87-floating-point-modes) new-mode)))) + View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/cb874f7fa2cd9d0dd38435… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/cb874f7fa2cd9d0dd38435… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][master] Update release notes for #424
by Raymond Toy (＠rtoy) 14 Aug '25

14 Aug '25

Raymond Toy pushed to branch master at cmucl / cmucl Commits: cb874f7f by Raymond Toy at 2025-08-14T07:42:48-07:00 Update release notes for #424 - - - - - 1 changed file: - src/general-info/release-21f.md Changes: ===================================== src/general-info/release-21f.md ===================================== @@ -137,6 +137,7 @@ public domain. * #417 PCL complains about repeated aux variables in defmethod * #418 Update asdf to version 3.3.7.4 (for package-local-nicknames) + * #424 Use fdlibm `hypot` to fix bug in snapshot 2025-07 * Other changes: * Improvements to the PCL implementation of CLOS: * Changes to building procedure: View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/cb874f7fa2cd9d0dd38435c… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/cb874f7fa2cd9d0dd38435c… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][master] 2 commits: Fix #424: Use fdlibm hypot
by Raymond Toy (＠rtoy) 14 Aug '25

14 Aug '25

Raymond Toy pushed to branch master at cmucl / cmucl Commits: 4660c011 by Raymond Toy at 2025-08-14T07:40:17-07:00 Fix #424: Use fdlibm hypot - - - - - 6a2028c5 by Raymond Toy at 2025-08-14T07:40:17-07:00 Merge branch 'issue-424-fix-hypot' into 'master' Fix #424: Use fdlibm hypot Closes #424 See merge request cmucl/cmucl!311 - - - - - 4 changed files: - src/code/irrat.lisp - src/lisp/GNUmakefile - + src/lisp/e_hypot.c - tests/irrat.lisp Changes: ===================================== src/code/irrat.lisp ===================================== @@ -87,7 +87,7 @@ (def-math-rtn ("__ieee754_pow" %pow) 2) #-(or x86 sparc-v7 sparc-v8 sparc-v9) (def-math-rtn "sqrt" 1) -(def-math-rtn "hypot" 2) +(def-math-rtn ("__ieee754_hypot" %hypot) 2) (def-math-rtn ("fdlibm_log1p" %log1p) 1) (def-math-rtn ("fdlibm_expm1" %expm1) 1) ===================================== src/lisp/GNUmakefile ===================================== @@ -28,6 +28,7 @@ FDLIBM = k_sin.c k_cos.c k_tan.c s_sin.c s_cos.c s_tan.c sincos.c \ e_rem_pio2.c k_rem_pio2.c \ e_log10.c s_scalbn.c \ setexception.c \ + e_hypot.c \ log2.c SRCS = lisp.c coreparse.c alloc.c monitor.c print.c interr.c \ ===================================== src/lisp/e_hypot.c ===================================== @@ -0,0 +1,144 @@ + +/* @(#)e_hypot.c 1.3 95/01/18 */ +/* + * ==================================================== + * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. + * + * Developed at SunSoft, a Sun Microsystems, Inc. business. + * Permission to use, copy, modify, and distribute this + * software is freely granted, provided that this notice + * is preserved. + * ==================================================== + */ + +/* __ieee754_hypot(x,y) + * + * Method : + * If (assume round-to-nearest) z=x*x+y*y + * has error less than sqrt(2)/2 ulp, than + * sqrt(z) has error less than 1 ulp (exercise). + * + * So, compute sqrt(x*x+y*y) with some care as + * follows to get the error below 1 ulp: + * + * Assume x>y>0; + * (if possible, set rounding to round-to-nearest) + * 1. if x > 2y use + * x1*x1+(y*y+(x2*(x+x1))) for x*x+y*y + * where x1 = x with lower 32 bits cleared, x2 = x-x1; else + * 2. if x <= 2y use + * t1*y1+((x-y)*(x-y)+(t1*y2+t2*y)) + * where t1 = 2x with lower 32 bits cleared, t2 = 2x-t1, + * y1= y with lower 32 bits chopped, y2 = y-y1. + * + * NOTE: scaling may be necessary if some argument is too + * large or too tiny + * + * Special cases: + * hypot(x,y) is INF if x or y is +INF or -INF; else + * hypot(x,y) is NAN if x or y is NAN. + * + * Accuracy: + * hypot(x,y) returns sqrt(x^2+y^2) with error less + * than 1 ulps (units in the last place) + */ + +#include "fdlibm.h" + +static inline void +set_hiword(double* d, int hi) +{ + union { int i[2]; double d; } ux; + + ux.d = *d; + ux.i[HIWORD] = hi; + *d = ux.d; +} + +static inline int +get_loword(double d) +{ + union { int i[2]; double d; } ux; + ux.d = d; + return ux.i[LOWORD]; +} + +static inline int +get_hiword(double d) +{ + union { int i[2]; double d; } ux; + ux.d = d; + + return ux.i[HIWORD]; +} + + +#ifdef __STDC__ + double __ieee754_hypot(double x, double y) +#else + double __ieee754_hypot(x,y) + double x, y; +#endif +{ + double a=x,b=y,t1,t2,y1,y2,w; + int j,k,ha,hb; + + ha = get_hiword(x)&0x7fffffff; /* high word of x */ + hb = get_hiword(y)&0x7fffffff; /* high word of y */ + if(hb > ha) {a=y;b=x;j=ha; ha=hb;hb=j;} else {a=x;b=y;} + set_hiword(&a,ha); /* a <- |a| */ + set_hiword(&b, hb); /* b <- |b| */ + + if((ha-hb)>0x3c00000) {return a+b;} /* x/y > 2**60 */ + k=0; + if(ha > 0x5f300000) { /* a>2**500 */ + if(ha >= 0x7ff00000) { /* Inf or NaN */ + w = a+b; /* for sNaN */ + if(((ha&0xfffff)|get_loword(a))==0) w = a; + if(((hb^0x7ff00000)|get_loword(b))==0) w = b; + return w; + } + /* scale a and b by 2**-600 */ + ha -= 0x25800000; hb -= 0x25800000; k += 600; + set_hiword(&a, ha); + set_hiword(&b, hb); + } + if(hb < 0x20b00000) { /* b < 2**-500 */ + if(hb <= 0x000fffff) { /* subnormal b or 0 */ + if((hb|(get_loword(b)))==0) return a; + t1=0; + set_hiword(&t1, 0x7fd00000); /* t1=2^1022 */ + b *= t1; + a *= t1; + k -= 1022; + } else { /* scale a and b by 2^600 */ + ha += 0x25800000; /* a *= 2^600 */ + hb += 0x25800000; /* b *= 2^600 */ + k -= 600; + set_hiword(&a, ha); + set_hiword(&b, hb); + } + } + /* medium size a and b */ + w = a-b; + if (w>b) { + t1 = 0; + set_hiword(&t1, ha); + t2 = a-t1; + w = sqrt(t1*t1-(b*(-b)-t2*(a+t1))); + } else { + a = a+a; + y1 = 0; + set_hiword(&y1, hb); + y2 = b - y1; + t1 = 0; + set_hiword(&t1, ha+0x00100000); + t2 = a - t1; + w = sqrt(t1*y1-(w*(-w)-(t1*y2+t2*b))); + } + if(k!=0) { + t1 = 1.0; + set_hiword(&t1, get_hiword(t1) + (k<<20)); + return t1*w; + } else return w; +} ===================================== tests/irrat.lisp ===================================== @@ -245,3 +245,13 @@ (tanh #c(200w0 200w0)))) +;; See bug #424 +(define-test hypot + (:tag :issues) + (assert-eql 3.8950612975366328d0 + (kernel:%hypot 2.302585092994046d0 3.141592653589793d0)) + (let ((result (expt #C(2.302585092994046d0 3.141592653589793d0) 4))) + (assert-eql -188.4466069439329d0 + (realpart result)) + (assert-eql -132.16721026205448d0 + (imagpart result)))) View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/9a8fd3c477f3bbfa1dec82… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/9a8fd3c477f3bbfa1dec82… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-425-correctly-rounded-math-functions] Set 64-bit precision (80-bit float) in call-out
by Raymond Toy (＠rtoy) 13 Aug '25

13 Aug '25

Raymond Toy pushed to branch issue-425-correctly-rounded-math-functions at cmucl / cmucl Commits: 5ff177cb by Raymond Toy at 2025-08-13T15:06:30-07:00 Set 64-bit precision (80-bit float) in call-out Add some tests for the special functions where we were computing the incorrect value when the precision was not set to 64-bit (80-bit float). - - - - - 2 changed files: - src/compiler/x86/sse2-c-call.lisp - tests/irrat.lisp Changes: ===================================== src/compiler/x86/sse2-c-call.lisp ===================================== @@ -41,7 +41,7 @@ #+core-math (:temporary (:sc unsigned-stack) fpu-cw) #+core-math - (:temporary (:sc unsigned-reg) temp-cw) + (:temporary (:sc unsigned-reg :offset esi-offset) temp-cw) (:node-var node) (:vop-var vop) (:save-p t) @@ -51,11 +51,14 @@ #+core-math (progn ;; Save the x87 FPU control word. Then modify it to set the - ;; precision bits to double (2). + ;; precision bits to 3 for 64-bit mantissas for 80-bit + ;; arithmetic. If we don't some of some special functions + ;; return incorrect values because the x87 precision was set to + ;; single. (inst fnstcw save-fpu-cw) (move temp-cw save-fpu-cw) (inst and temp-cw (dpb 0 (byte 2 8) #xffff)) ; Zap the precision control bits - (inst or temp-cw (dpb 2 (byte 2 8) 0)) ; Set precision control to double + (inst or temp-cw (dpb 3 (byte 3 8) 0)) ; Set precision control bits (move fpu-cw temp-cw) (inst fldcw fpu-cw) ; New CW ) ===================================== tests/irrat.lisp ===================================== @@ -245,3 +245,18 @@ (tanh #c(200w0 200w0)))) +(define-test cos-tiny + (:tag issues) + ;; This test comes from the Maxima testsuite where core-math was not + ;; computing cos(8.881784197001252d-16) correctly. It should be + ;; exactly 1. We were getting one bit less. + (assert-eql 1d0 + (cos (kernel:make-double-float 1020264448 0)))) + +(define-test log.special-value + (:tag issues) + ;; This test comes from the Maxima testsuite where core-math was not + ;; computing log(0.6899991035461426d0) => -0.37106498060016496d0. + ;; The computed result looked it had like a single-precision accuracy. + (assert-eql -0.37106498060016496d0 + (kernel:%log (kernel:make-double-float (+ 1071644672 398457) 0)))) View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/5ff177cb9f69fcd1042dbdf… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/5ff177cb9f69fcd1042dbdf… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-425-correctly-rounded-math-functions] Linux cross-compile needs x86_linux_clang
by Raymond Toy (＠rtoy) 12 Aug '25

12 Aug '25

Raymond Toy pushed to branch issue-425-correctly-rounded-math-functions at cmucl / cmucl Commits: 0cca89a6 by Raymond Toy at 2025-08-12T15:52:01-07:00 Linux cross-compile needs x86_linux_clang Forgot to set the build config when running bin/build.sh. - - - - - 1 changed file: - .gitlab-ci.yml Changes: ===================================== .gitlab-ci.yml ===================================== @@ -182,10 +182,10 @@ linux:cross-build: script: # Use x86_linux_clang to build because x86_linux with gcc gives an # internal compiler error building core-math. - - bin/create-target.sh xtarget x86_linux_clang - - bin/create-target.sh xcross x86_linux_clang + - bin/create-target.sh xtarget "$build_config" + - bin/create-target.sh xcross "$build_config" - bin/cross-build-world.sh -crl xtarget xcross src/tools/cross-scripts/cross-x86-x86.lisp dist/bin/lisp - - bin/build.sh -b xlinux $bootstrap -R -C "" -o "xtarget/lisp/lisp -lib xtarget/lisp" + - bin/build.sh -b xlinux $bootstrap -R -C "$build_config" -o "xtarget/lisp/lisp -lib xtarget/lisp" - bin/make-dist.sh -I xdist xlinux-4 linux:test: View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/0cca89a6ac031562f71b294… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/0cca89a6ac031562f71b294… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-425-correctly-rounded-math-functions] Linux cross-compile needs x86_linux_clang
by Raymond Toy (＠rtoy) 12 Aug '25

12 Aug '25

Raymond Toy pushed to branch issue-425-correctly-rounded-math-functions at cmucl / cmucl Commits: eed071ef by Raymond Toy at 2025-08-12T07:19:06-07:00 Linux cross-compile needs x86_linux_clang As we did for the linux and opensuse builders, we need to set the build_config to Config.x86_linux_clang to build the core-math routines. - - - - - 1 changed file: - .gitlab-ci.yml Changes: ===================================== .gitlab-ci.yml ===================================== @@ -174,6 +174,11 @@ linux:cross-build: - job: linux:build #- job: linux:install artifacts: true + variables: + # The CI on linux needs to use clang to build the core-math + # routines. Gcc gets an internal compiler error for the OS + # running the linux builder. + build_config: "x86_linux_clang" script: # Use x86_linux_clang to build because x86_linux with gcc gives an # internal compiler error building core-math. View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/eed071ef3e30185ea44794e… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/eed071ef3e30185ea44794e… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0