common-lisp/sbcl.sbcl
1
0
Fork 0
mirror of git://git.code.sf.net/p/sbcl/sbcl synced 2026-09-10 07:26:40 -04:00
Code Issues Actions 8 Packages Projects Releases Wiki Activity

Make the cross float sqrt work on -0.0

Browse source
This commit is contained in:
Stas Boukarev 2025-09-11 02:54:02 +03:00
parent c76b206937
commit d526dd1968
1 changed files with 8 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

11
src/code/cross-float.lisp
View file

@ -336,9 +336,14 @@
(%%sqrt rational (cl:/ (isqrt (numerator rational)) (isqrt (denominator rational))))))
(defun sb-xc:sqrt (arg)
(let ((format (if (rationalp arg) 'single-float (type-of arg))))
(with-memoized-math-op (sqrt arg)
(flonum-from-rational (%sqrt (rational arg)) format))))
(cond ((eql arg 0)
0.0)
((sb-xc:= arg 0)
arg)
(t
(let ((format (if (rationalp arg) 'single-float (type-of arg))))
(with-memoized-math-op (sqrt arg)
(flonum-from-rational (%sqrt (rational arg)) format))))))
;;; There seems to be no portable way to mask float traps, so right
;;; now we ignore them and hardcode special cases.