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

[Git][cmucl/cmucl] Pushed new branch issue-432-fix-man-option-for-make-extra-dist
by Raymond Toy (＠rtoy) 18 Aug '25

18 Aug '25

Raymond Toy pushed new branch issue-432-fix-man-option-for-make-extra-dist at cmucl / cmucl -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/tree/issue-432-fix-man-option-… 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] Update release notes
by Raymond Toy (＠rtoy) 17 Aug '25

17 Aug '25

Raymond Toy pushed to branch issue-431-fix-set-floating-point-modes at cmucl / cmucl Commits: 02a3fe80 by Raymond Toy at 2025-08-17T07:18:07-07:00 Update release notes - - - - - 1 changed file: - src/general-info/release-21f.md Changes: ===================================== src/general-info/release-21f.md ===================================== @@ -139,6 +139,7 @@ public domain. 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))` + * $431 Fix setting of x87 FP modes in `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/02a3fe80b25f6d508d59363… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/02a3fe80b25f6d508d59363… 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] More cleanups
by Raymond Toy (＠rtoy) 17 Aug '25

17 Aug '25

Raymond Toy pushed to branch issue-431-fix-set-floating-point-modes at cmucl / cmucl Commits: ef367c52 by Raymond Toy at 2025-08-17T07:14:14-07:00 More cleanups Remove old comment and debugging prints. Wrap comments neatly. Clean up printing of the mode word for x87. - - - - - 2 changed files: - src/code/float-trap.lisp - src/compiler/x86/float-sse2.lisp Changes: ===================================== src/code/float-trap.lisp ===================================== @@ -109,37 +109,17 @@ (ash x87-enables 7)))) final-mode)) - #+nil - (defun (setf floating-point-modes) (new-mode) - (declare (type (unsigned-byte 32) new-mode)) - ;; Set the floating point modes for both X87 and SSE2. This - ;; include the rounding control bits. - (let* ((rc (ldb float-rounding-mode new-mode)) - (new-exceptions (logand #x3f new-mode)) - (new-enables (logand #x3f (ash new-mode -7))) - (x87-modes - (logior new-exceptions - (ash rc 10) - (ash new-enables 16) - ;; Set precision control to be 64-bit, always. We - ;; don't use the x87 registers with sse2, so this - ;; is ok and would be the correct setting if we - ;; ever support long-floats. - (ash 3 (+ 8 16))))) - (setf (vm::sse2-floating-point-modes) (ldb (byte 24 0) new-mode)) - (setf (vm::x87-floating-point-modes) (ldb (byte 24 0) x87-modes))) - new-mode) + (defun (setf floating-point-modes) (new-mode) (declare (type (unsigned-byte 32) new-mode)) ;; Set the floating point modes for both X87 and SSE2. This ;; include the rounding control bits. - #+nil - (format t "new-mode = ~8,'0x~%" new-mode) (let* ((rc (ldb float-rounding-mode new-mode)) (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. + ;; 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)) @@ -152,8 +132,7 @@ (setf (vm::sse2-floating-point-modes) (ldb (byte 16 0) new-mode)) (setf (vm::x87-floating-point-modes) (ldb (byte 30 0) x87-modes))) - new-mode) - ) + new-mode)) #+(and sse2 darwin) (progn ===================================== src/compiler/x86/float-sse2.lisp ===================================== @@ -1422,7 +1422,7 @@ ;; 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 float-rounding-mode sse-mode))) - (format t "Rounding control: ~30T~A ~A~%" + (format t "Rounding control: ~30T#b~2,'0b ~S~%" rc (car (rassoc rc rounding-mode-alist)))) (print-fp-exceptions-enabled (ldb float-traps-byte sse-mode)) @@ -1432,26 +1432,28 @@ "Print X87 floating modes word in a readable fashion." ;; 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 "SW: FPU busy: ~30T~A~%" (ldb (byte 1 15) x87-mode)) - (format t "SW: Condition code C3: ~30T~A~%" (ldb (byte 1 14) x87-mode)) - (format t "SW: Top of stack: ~30T~D~%" (ldb (byte 3 11) x87-mode)) - (format t "SW: Condition code C2: ~30T~A~%" (ldb (byte 1 10) x87-mode)) - (format t "SW: Condition code C1: ~30T~A~%" (ldb (byte 1 9) x87-mode)) - (format t "SW: Condition code C0: ~30T~A~%" (ldb (byte 1 8) x87-mode)) - (format t "SW: Error summary: ~30T~A~%" (ldb (byte 1 7) x87-mode)) - (format t "SW: Stack fault: ~30T~A~%" (ldb (byte 1 6) x87-mode)) + (format t "Status word:~%") + (format t "FPU busy: ~30T~A~%" (ldb (byte 1 15) x87-mode)) + (format t "Condition code C3: ~30T~A~%" (ldb (byte 1 14) x87-mode)) + (format t "Top of stack: ~30T~D~%" (ldb (byte 3 11) x87-mode)) + (format t "Condition code C2: ~30T~A~%" (ldb (byte 1 10) x87-mode)) + (format t "Condition code C1: ~30T~A~%" (ldb (byte 1 9) x87-mode)) + (format t "Condition code C0: ~30T~A~%" (ldb (byte 1 8) x87-mode)) + (format t "Error summary: ~30T~A~%" (ldb (byte 1 7) x87-mode)) + (format t "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 x87-float-infinity-control-byte x87-mode)) + (format t "~%Control word:~%") + (format t "Reserved: ~30T#b~2,'0b~%" (ldb (byte 2 (+ 13 16)) x87-mode)) + (format t "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~%" + (format t "Rounding control: ~30T#b~2,'0b ~S~%" rc (car (rassoc rc rounding-mode-alist)))) (let ((pc (ldb x87-float-precision-control-byte x87-mode))) - (format t "CW: Precision control: ~30T~A ~A~%" + (format t "Precision control: ~30T#b~2,'0b ~S~%" pc (car (rassoc pc x87-float-precision-control-alist)))) - (format t "CW: Reserved: ~30T~A~%" (ldb (byte 2 (+ 6 16)) x87-mode)) + (format t "Reserved: ~30T#b~2,'0b~%" (ldb (byte 2 (+ 6 16)) x87-mode)) (print-fp-exceptions-enabled (ldb x87-float-traps-byte x87-mode))) (defknown x87-floating-point-modes () float-modes (flushable)) View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/ef367c524ca6e359970c5fe… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/ef367c524ca6e359970c5fe… 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] 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