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

[Git][cmucl/cmucl][issue-459-more-accurate-dd-complex-div] Clean up a comment
by Raymond Toy (＠rtoy) 19 Jan '26

19 Jan '26

Raymond Toy pushed to branch issue-459-more-accurate-dd-complex-div at cmucl / cmucl Commits: 466dee14 by Raymond Toy at 2026-01-19T11:26:01-08:00 Clean up a comment - - - - - 1 changed file: - src/code/numbers.lisp Changes: ===================================== src/code/numbers.lisp ===================================== @@ -751,12 +751,12 @@ (d (imagpart y)) (ab (max (abs a) (abs b))) (cd (max (abs c) (abs d))) - ;; S is always multipled by a power of 2. It's either - ;; 2, 0.5, or +be+, which is also a power of two. - ;; (Either 2^105 or 2^209). Hence, we can just use a - ;; double-float instead of a double-double-float - ;; because the double-double always has 0d0 for the lo - ;; part. + ;; S is always multipled by a power of 2. It's + ;; multiplied by either 2, 0.5, or +be+, which is + ;; also a power of two. (Either 2^105 or 2^209). + ;; Hence, we can just use a double-float instead + ;; of a double-double-float because the + ;; double-double always has 0d0 for the lo part. (s 1d0)) (declare (double-float s)) ;; If a or b is big, scale down a and b. View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/466dee1428a1b6ed3a439e7… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/466dee1428a1b6ed3a439e7… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-459-more-accurate-dd-complex-div] 2 commits: Let bind the constants as needed in each routine
by Raymond Toy (＠rtoy) 19 Jan '26

19 Jan '26

Raymond Toy pushed to branch issue-459-more-accurate-dd-complex-div at cmucl / cmucl Commits: b29c585b by Raymond Toy at 2026-01-19T11:07:59-08:00 Let bind the constants as needed in each routine This also means the macro needs args to define the values. Address review comments. - - - - - 184720ee by Raymond Toy at 2026-01-19T11:18:04-08:00 Remove unneeded defconstants and update/clarify some comments - - - - - 1 changed file: - src/code/numbers.lisp Changes: ===================================== src/code/numbers.lisp ===================================== @@ -602,53 +602,11 @@ ;; In particular iteration 1 and 3 are added. Iteration 2 and 4 were ;; not added. The test examples from iteration 2 and 4 didn't change ;; with or without changes added. -(defconstant +cdiv-rmin+ least-positive-normalized-double-float) -(defconstant +cdiv-rbig+ (/ most-positive-double-float 2)) -(defconstant +cdiv-rmin2+ (scale-float 1d0 -53)) -(defconstant +cdiv-rminscal+ (scale-float 1d0 51)) -(defconstant +cdiv-rmax2+ (* +cdiv-rbig+ +cdiv-rmin2+)) -;; This is the value of %eps from Scilab -(defconstant +cdiv-eps+ (scale-float 1d0 -52)) -(defconstant +cdiv-be+ (/ 2 (* +cdiv-eps+ +cdiv-eps+))) -(defconstant +cdiv-2/eps+ (/ 2 +cdiv-eps+)) - -;; Same constants but for DOUBLE-DOUBLE-FLOATS. Some of these aren't -;; well-defined for DOUBLE-DOUBLE-FLOATS so we make our best guess at -;; what they might be. Since double-doubles have about twice as many -;; bits of precision as a DOUBLE-FLOAT, we generally just double the -;; exponent (for SCALE-FLOAT) of the corresponding DOUBLE-FLOAT values -;; above. ;; -;; Note also that both LEAST-POSITIVE-NORMALIZED-DOUBLE-DOUBLE-FLOAT -;; and MOST-POSITIVE-DOUBLE-DOUBLE-FLOAT have a low component of 0, so -;; there's no loss of precision if we use the corresponding -;; DOUBLE-FLOAT values. Likewise, all the other constants here are -;; powers of 2, so the DOUBLE-DOUBLE-FLOAT values can be represented -;; exactly as a DOUBLE-FLOAT. We can use DOUBLE-FLOAT values. -(defconstant +cdiv-dd-rmin+ - least-positive-normalized-double-float) -(defconstant +cdiv-dd-rbig+ - (/ most-positive-double-float 2)) -(defconstant +cdiv-dd-rmin2+ - (scale-float 1d0 -106)) -(defconstant +cdiv-dd-rminscal+ - (scale-float 1d0 102)) -(defconstant +cdiv-dd-rmax2+ - (* +cdiv-dd-rbig+ +cdiv-dd-rmin2+)) -;; Epsilon for double-doubles isn't really well-defined because things -;; like (+ 1w0 1w-200) is a valid double-double float. -(defconstant +cdiv-dd-eps+ - (scale-float 1d0 -104)) -(defconstant +cdiv-dd-be+ - (/ 2 (* +cdiv-dd-eps+ +cdiv-dd-eps+))) -(defconstant +cdiv-dd-2/eps+ - (/ 2 +cdiv-dd-eps+)) - - -;; Make these functions accessible. cdiv-double-float and -;; cdiv-single-float are used by deftransforms. Of course, two-arg-/ -;; is the interface to division. cdiv-generic isn't used anywhere -;; else. +;; Make these functions accessible outside the block compilation. +;; CDIV-DOUBLE-FLOAT and CDIV-SINGLE-FLOAT are used by deftransforms. +;; Of course, TWO-ARG-/ is the interface to division. CDIV-GENERIC +;; isn't used anywhere else. (declaim (ext:start-block cdiv-double-float cdiv-single-float cdiv-double-double-float two-arg-/)) @@ -657,27 +615,15 @@ ;; implementation is basically identical except for the constants. ;; use a macro to generate both versions. (eval-when (:compile-toplevel) -(defmacro define-cdiv (type) - (let* ((dd (if (eq type 'double-float) - "" - "DD-")) - (opt '((optimize (speed 3) (safety 0)))) +(defmacro define-cdiv (type &key rmin rbig rmin2 rminscal eps) + (let* ((opt '((optimize (speed 3) (safety 0)))) (name (symbolicate "CDIV-" type)) (compreal (symbolicate "CDIV-" type "-INTERNAL-COMPREAL")) (subinternal (symbolicate "CDIV-" type "-ROBUST-SUBINTERNAL")) (internal (symbolicate "CDIV-" type "-ROBUST-INTERNAL")) (docstring (let ((*print-case* :downcase)) (format nil "Accurate division of complex ~A numbers x and y." - type))) - ;; Create the correct symbols for the constants we need, - ;; as defined above. - (+rmin+ (symbolicate "+CDIV-" dd "RMIN+")) - (+rbig+ (symbolicate "+CDIV-" dd "RBIG+")) - (+rmin2+ (symbolicate "+CDIV-" dd "RMIN2+")) - (+rminscal+ (symbolicate "+CDIV-" dd "RMINSCAL+")) - (+rmax2+ (symbolicate "+CDIV-" dd "RMAX2+")) - (+be+ (symbolicate "+CDIV-" dd "BE+")) - (+2/eps+ (symbolicate "+CDIV-" dd "2/EPS+"))) + type)))) `(progn (defun ,compreal (a b c d r tt) (declare (,type a b c d r tt) @@ -699,7 +645,7 @@ ;; = (a + b*r)/(c + d*r). ;; ;; Thus tt = (c + d*r). - (cond ((>= (abs r) ,+rmin+) + (cond ((>= (abs r) ,rmin) (let ((br (* b r))) (if (/= br 0) (/ (+ a br) tt) @@ -719,115 +665,154 @@ ,@opt) (let* ((r (/ d c)) (tt (+ c (* d r)))) - ;; e is the real part and f is the imaginary part. We - ;; can use internal-compreal for the imaginary part by - ;; noticing that the imaginary part of (a+i*b)/(c+i*d) is - ;; the same as the real part of (b-i*a)/(c+i*d). + ;; e is the real part and f is the imaginary part. We can + ;; use COMPREAL for the imaginary part by noticing that the + ;; imaginary part of (a+i*b)/(c+i*d) is the same as the + ;; real part of (b-i*a)/(c+i*d). (let ((e (,compreal a b c d r tt)) (f (,compreal b (- a) c d r tt))) (values e f)))) (defun ,internal (x y) - (declare (type (complex ,type) x y) - ,@opt) - (let ((a (realpart x)) - (b (imagpart x)) - (c (realpart y)) - (d (imagpart y))) - (declare (,type a b c d)) - (flet ((maybe-scale (abs-tst a b c d) - (declare (,type a b c d)) - ;; This implements McGehearty's iteration 3 to - ;; handle the case when some values are too big - ;; and should be scaled down. Also if some - ;; values are too tiny, scale them up. - (let ((abs-a (abs a)) - (abs-b (abs b))) - (if (or (> abs-tst ,+rbig+) - (> abs-a ,+rbig+) - (> abs-b ,+rbig+)) - (setf a (* a 0.5d0) - b (* b 0.5d0) - c (* c 0.5d0) - d (* d 0.5d0)) - (if (< abs-tst ,+rmin2+) - (setf a (* a ,+rminscal+) - b (* b ,+rminscal+) - c (* c ,+rminscal+) - d (* d ,+rminscal+)) - (if (or (and (< abs-a ,+rmin+) - (< abs-b ,+rmax2+) - (< abs-tst ,+rmax2+)) - (and (< abs-b ,+rmin+) - (< abs-a ,+rmax2+) - (< abs-tst ,+rmax2+))) - (setf a (* a ,+rminscal+) - b (* b ,+rminscal+) - c (* c ,+rminscal+) - d (* d ,+rminscal+))))) - (values a b c d)))) - (cond - ((<= (abs d) (abs c)) - ;; |d| <= |c|, so we can use robust-subinternal to - ;; perform the division. - (multiple-value-bind (a b c d) - (maybe-scale (abs c) a b c d) - (,subinternal a b c d))) - (t - ;; |d| > |c|. So, instead compute - ;; - ;; (b + i*a)/(d + i*c) = ((b*d+a*c) + (a*d-b*c)*i)/(d^2+c^2) - ;; - ;; Compare this to (a+i*b)/(c+i*d) and we see that - ;; realpart of the former is the same, but the - ;; imagpart of the former is the negative of the - ;; desired division. - (multiple-value-bind (a b c d) - (maybe-scale (abs d) a b c d) - (multiple-value-bind (e f) - (,subinternal b a d c) - (values e (- f))))))))) - (defun ,name (x y) - ,docstring - (declare (type (complex ,type) x y) - ,@opt) - (let* ((a (realpart x)) - (b (imagpart x)) - (c (realpart y)) - (d (imagpart y)) - (ab (max (abs a) (abs b))) - (cd (max (abs c) (abs d))) - ;; S is always multipled by a power of 2. It's either - ;; 2, 0.5, or +be+, which is also a power of two. - ;; (Either 2^105 or 2^209). Hence, we can just use a - ;; double-float instead of a double-double-float - ;; because the double-double always has 0d0 for the lo - ;; part. - (s 1d0)) - (declare (double-float s)) - ;; If a or b is big, scale down a and b. - (when (>= ab ,+rbig+) - (setf x (* x 0.5d0) - s (* s 2))) - ;; If c or d is big, scale down c and d. - (when (>= cd ,+rbig+) - (setf y (* y 0.5d0) - s (* s 0.5d0))) - ;; If a or b is tiny, scale up a and b. - (when (<= ab (* ,+rmin+ ,+2/eps+)) - (setf x (* x ,+be+) - s (/ s ,+be+))) - ;; If c or d is tiny, scale up c and d. - (when (<= cd (* ,+rmin+ ,+2/eps+)) - (setf y (* y ,+be+) - s (* s ,+be+))) - (multiple-value-bind (e f) - (,internal x y) - (complex (* s e) - (* s f)))))))) + (declare (type (complex ,type) x y) + ,@opt) + (let ((a (realpart x)) + (b (imagpart x)) + (c (realpart y)) + (d (imagpart y))) + (declare (,type a b c d)) + (flet ((maybe-scale (abs-tst a b c d) + (declare (,type a b c d)) + ;; This implements McGehearty's iteration 3 to + ;; handle the case when some values are too big + ;; and should be scaled down. Also if some + ;; values are too tiny, scale them up. + (let ((+rbig+ ,rbig) + (+rmin2+ ,rmin2) + (+rminscal+ ,rminscal) + (+rmax2+ (* ,rbig ,rmin2)) + (+rmin+ ,rmin) + (abs-a (abs a)) + (abs-b (abs b))) + (if (or (> abs-tst +rbig+) + (> abs-a +rbig+) + (> abs-b +rbig+)) + (setf a (* a 0.5d0) + b (* b 0.5d0) + c (* c 0.5d0) + d (* d 0.5d0)) + (if (< abs-tst +rmin2+) + (setf a (* a +rminscal+) + b (* b +rminscal+) + c (* c +rminscal+) + d (* d +rminscal+)) + (if (or (and (< abs-a +rmin+) + (< abs-b +rmax2+) + (< abs-tst +rmax2+)) + (and (< abs-b +rmin+) + (< abs-a +rmax2+) + (< abs-tst +rmax2+))) + (setf a (* a +rminscal+) + b (* b +rminscal+) + c (* c +rminscal+) + d (* d +rminscal+))))) + (values a b c d)))) + (cond + ((<= (abs d) (abs c)) + ;; |d| <= |c|, so we can use robust-subinternal to + ;; perform the division. + (multiple-value-bind (a b c d) + (maybe-scale (abs c) a b c d) + (,subinternal a b c d))) + (t + ;; |d| > |c|. So, instead compute + ;; + ;; (b + i*a)/(d + i*c) = ((b*d+a*c) + (a*d-b*c)*i)/(d^2+c^2) + ;; + ;; Compare this to (a+i*b)/(c+i*d) and we see that + ;; realpart of the former is the same, but the + ;; imagpart of the former is the negative of the + ;; desired division. + (multiple-value-bind (a b c d) + (maybe-scale (abs d) a b c d) + (multiple-value-bind (e f) + (,subinternal b a d c) + (values e (- f))))))))) + (defun ,name (x y) + ,docstring + (declare (type (complex ,type) x y) + ,@opt) + (let ((+rmin+ ,rmin) + (+rbig+ ,rbig) + (+2/eps+ (/ 2 ,eps)) + (+be+ (/ 2 (* ,eps ,eps)))) + (let* ((a (realpart x)) + (b (imagpart x)) + (c (realpart y)) + (d (imagpart y)) + (ab (max (abs a) (abs b))) + (cd (max (abs c) (abs d))) + ;; S is always multipled by a power of 2. It's either + ;; 2, 0.5, or +be+, which is also a power of two. + ;; (Either 2^105 or 2^209). Hence, we can just use a + ;; double-float instead of a double-double-float + ;; because the double-double always has 0d0 for the lo + ;; part. + (s 1d0)) + (declare (double-float s)) + ;; If a or b is big, scale down a and b. + (when (>= ab +rbig+) + (setf x (* x 0.5d0) + s (* s 2))) + ;; If c or d is big, scale down c and d. + (when (>= cd +rbig+) + (setf y (* y 0.5d0) + s (* s 0.5d0))) + ;; If a or b is tiny, scale up a and b. + (when (<= ab (* +rmin+ +2/eps+)) + (setf x (* x +be+) + s (/ s +be+))) + ;; If c or d is tiny, scale up c and d. + (when (<= cd (* +rmin+ +2/eps+)) + (setf y (* y +be+) + s (* s +be+))) + (multiple-value-bind (e f) + (,internal x y) + (complex (* s e) + (* s f))))))))) ) -(define-cdiv double-float) -(define-cdiv double-double-float) +(define-cdiv double-float + :rmin least-positive-normalized-double-float + :rbig (/ most-positive-double-float 2) + :rmin2 (scale-float 1d0 -53) + :rminscal (scale-float 1d0 51) + ;; This is the value of %eps from Scilab + :eps (scale-float 1d0 -52)) + +;; Same constants but for DOUBLE-DOUBLE-FLOATS. Some of these aren't +;; well-defined for DOUBLE-DOUBLE-FLOATS (like eps) so we make our +;; best guess at what they might be. Since double-doubles have +;; twice as many bits of precision as a DOUBLE-FLOAT, we generally +;; just double the exponent (for SCALE-FLOAT) of the corresponding +;; DOUBLE-FLOAT values above. +;; +;; Note also that since both +;; LEAST-POSITIVE-NORMALIZED-DOUBLE-DOUBLE-FLOAT and +;; MOST-POSITIVE-DOUBLE-DOUBLE-FLOAT have a low component of 0, +;; there's no loss of precision if we use the corresponding +;; DOUBLE-FLOAT values. Likewise, all the other constants here are +;; powers of 2, so the DOUBLE-DOUBLE-FLOAT values can be represented +;; exactly as a DOUBLE-FLOAT. We can use DOUBLE-FLOAT values. This +;; also makes some of the operations a bit simpler since operations +;; between a DOUBLE-DOUBLE-FLOAT and a DOUBLE-FLOAT are simpler than +;; if both were DOUBLE-DOUBLE-FLOATs. +(define-cdiv double-double-float + :rmin least-positive-normalized-double-float + :rbig (/ most-positive-double-float 2) + :rmin2 (scale-float 1d0 -106) + :rminscal (scale-float 1d0 102) + :eps (scale-float 1d0 -104)) + ;; Smith's algorithm for complex division for (complex single-float). ;; We convert the parts to double-floats before computing the result. View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/7ef6184520ac5201dd878e… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/7ef6184520ac5201dd878e… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-459-more-accurate-dd-complex-div] 5 commits: Update cmucl.pot for new docstrings
by Raymond Toy (＠rtoy) 19 Jan '26

19 Jan '26

Raymond Toy pushed to branch issue-459-more-accurate-dd-complex-div at cmucl / cmucl Commits: 336b9322 by Raymond Toy at 2026-01-15T09:13:30-08:00 Update cmucl.pot for new docstrings Forgot to update this when accurate complex division was committed. [skip-ci] - - - - - 1b688607 by Raymond Toy at 2026-01-18T15:38:33-08:00 Fix #463: Handle double-double-float and double-float comparisons - - - - - 67dc5c9d by Raymond Toy at 2026-01-18T15:38:34-08:00 Merge branch 'issue-463-double-double-float-comparison' into 'master' Fix #463: Handle double-double-float and double-float comparisons Closes #463 See merge request cmucl/cmucl!342 - - - - - bde45741 by Raymond Toy at 2026-01-18T15:53:04-08:00 Merge branch 'master' into issue-459-more-accurate-dd-complex-div - - - - - 7ef61845 by Raymond Toy at 2026-01-19T06:41:38-08:00 Update cmucl.pot - - - - - 1 changed file: - src/compiler/float-tran-dd.lisp Changes: ===================================== src/compiler/float-tran-dd.lisp ===================================== @@ -678,6 +678,18 @@ (kernel:double-double-hi b) (kernel:double-double-lo b))) +(deftransform = ((a b) (vm::double-double-float double-float) *) + `(dd= (kernel:double-double-hi a) + (kernel:double-double-lo a) + b + 0d0)) + +(deftransform = ((a b) (double-float vm::double-double-float) *) + `(dd= a + 0d0 + (kernel:double-double-hi b) + (kernel:double-double-lo b))) + (deftransform < ((a b) (vm::double-double-float vm::double-double-float) *) `(dd< (kernel:double-double-hi a) @@ -685,10 +697,34 @@ (kernel:double-double-hi b) (kernel:double-double-lo b))) +(deftransform < ((a b) (vm::double-double-float double-float) *) + `(dd< (kernel:double-double-hi a) + (kernel:double-double-lo a) + b + 0d0)) + +(deftransform < ((a b) (double-float vm::double-double-float) *) + `(dd< a + 0d0 + (kernel:double-double-hi b) + (kernel:double-double-lo b))) + (deftransform > ((a b) (vm::double-double-float vm::double-double-float) *) `(dd> (kernel:double-double-hi a) (kernel:double-double-lo a) (kernel:double-double-hi b) (kernel:double-double-lo b))) + +(deftransform > ((a b) (vm::double-double-float double-float) *) + `(dd> (kernel:double-double-hi a) + (kernel:double-double-lo a) + b + 0d0)) + +(deftransform > ((a b) (double-float vm::double-double-float) *) + `(dd> a + 0d0 + (kernel:double-double-hi b) + (kernel:double-double-lo b))) ) ; end progn View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/3b75ab818edd40c13b9e1f… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/3b75ab818edd40c13b9e1f… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][master] 2 commits: Fix #463: Handle double-double-float and double-float comparisons
by Raymond Toy (＠rtoy) 18 Jan '26

18 Jan '26

Raymond Toy pushed to branch master at cmucl / cmucl Commits: 1b688607 by Raymond Toy at 2026-01-18T15:38:33-08:00 Fix #463: Handle double-double-float and double-float comparisons - - - - - 67dc5c9d by Raymond Toy at 2026-01-18T15:38:34-08:00 Merge branch 'issue-463-double-double-float-comparison' into 'master' Fix #463: Handle double-double-float and double-float comparisons Closes #463 See merge request cmucl/cmucl!342 - - - - - 1 changed file: - src/compiler/float-tran-dd.lisp Changes: ===================================== src/compiler/float-tran-dd.lisp ===================================== @@ -674,6 +674,18 @@ (kernel:double-double-hi b) (kernel:double-double-lo b))) +(deftransform = ((a b) (vm::double-double-float double-float) *) + `(dd= (kernel:double-double-hi a) + (kernel:double-double-lo a) + b + 0d0)) + +(deftransform = ((a b) (double-float vm::double-double-float) *) + `(dd= a + 0d0 + (kernel:double-double-hi b) + (kernel:double-double-lo b))) + (deftransform < ((a b) (vm::double-double-float vm::double-double-float) *) `(dd< (kernel:double-double-hi a) @@ -681,10 +693,34 @@ (kernel:double-double-hi b) (kernel:double-double-lo b))) +(deftransform < ((a b) (vm::double-double-float double-float) *) + `(dd< (kernel:double-double-hi a) + (kernel:double-double-lo a) + b + 0d0)) + +(deftransform < ((a b) (double-float vm::double-double-float) *) + `(dd< a + 0d0 + (kernel:double-double-hi b) + (kernel:double-double-lo b))) + (deftransform > ((a b) (vm::double-double-float vm::double-double-float) *) `(dd> (kernel:double-double-hi a) (kernel:double-double-lo a) (kernel:double-double-hi b) (kernel:double-double-lo b))) + +(deftransform > ((a b) (vm::double-double-float double-float) *) + `(dd> (kernel:double-double-hi a) + (kernel:double-double-lo a) + b + 0d0)) + +(deftransform > ((a b) (double-float vm::double-double-float) *) + `(dd> a + 0d0 + (kernel:double-double-hi b) + (kernel:double-double-lo b))) ) ; end progn View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/336b9322330eee9a4b44d9… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/336b9322330eee9a4b44d9… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-428-add-constants-for-irrat] 2 commits: Add constant for log base 2 of e
by Raymond Toy (＠rtoy) 18 Jan '26

18 Jan '26

Raymond Toy pushed to branch issue-428-add-constants-for-irrat at cmucl / cmucl Commits: 46e88b6a by Raymond Toy at 2026-01-18T07:51:21-08:00 Add constant for log base 2 of e - - - - - 688490f6 by Raymond Toy at 2026-01-18T07:57:48-08:00 Update release notes and cmucl.pot for docstrings - - - - - 3 changed files: - src/code/irrat.lisp - src/general-info/release-22a.md - src/i18n/locale/cmucl.pot Changes: ===================================== src/code/irrat.lisp ===================================== @@ -32,7 +32,10 @@ (log 2w0) "Natural log of 2w0") - +(defconstant +log2-e+ + (log (exp 1d0) 2d0) + "Base 2 logarithm of e") + ;;; Make these INLINE, since the call to C is at least as compact as a Lisp ;;; call, and saves number consing to boot. ;;; @@ -836,7 +839,7 @@ (if (minusp number) (complex (coerce (log (- number)) 'single-float) (coerce pi 'single-float)) - (coerce (/ (log2 number) #.(log (exp 1d0) 2d0)) 'single-float))) + (coerce (/ (log2 number) +log2-e+) 'single-float))) ((ratio) (if (minusp number) (complex (coerce (log (- number)) 'single-float) @@ -861,7 +864,7 @@ (coerce (%log1p (coerce (- number 1) 'double-float)) 'single-float) (coerce (/ (- (log2 top) (log2 bot)) - #.(log (exp 1d0) 2d0)) + +log2-e+) 'single-float))))) (((foreach single-float double-float)) ;; Is (log -0) -infinity (libm.a) or -infinity + i*pi (Kahan)? ===================================== src/general-info/release-22a.md ===================================== @@ -28,6 +28,7 @@ public domain. * #373: Add macros to handle temp files and directories * #387: Update CLX from upstream to version from 2024-09-11 * #408: Add C functions to get user name and home dir + * #428: Add some constants for common values in irrat.lisp * #437: Use `realpath` to find the lisp executable * #444: `unix-stat` and friends return wrong timestamp * #449: Touch 64-bit timestamp.txt with desired time for testing ===================================== src/i18n/locale/cmucl.pot ===================================== @@ -5012,6 +5012,10 @@ msgstr "" msgid "Natural log of 2w0" msgstr "" +#: src/code/irrat.lisp +msgid "Base 2 logarithm of e" +msgstr "" + #: src/code/irrat.lisp msgid "Return e raised to the power NUMBER." msgstr "" View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/3451a113880c01328881c9… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/3451a113880c01328881c9… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-428-add-constants-for-irrat] Update cmucl.pot for new docstrings
by Raymond Toy (＠rtoy) 18 Jan '26

18 Jan '26

Raymond Toy pushed to branch issue-428-add-constants-for-irrat at cmucl / cmucl Commits: 3451a113 by Raymond Toy at 2026-01-18T07:42:55-08:00 Update cmucl.pot for new docstrings - - - - - 1 changed file: - src/i18n/locale/cmucl.pot Changes: ===================================== src/i18n/locale/cmucl.pot ===================================== @@ -5004,6 +5004,14 @@ msgid "" " (= x (float (rationalize x) x))" msgstr "" +#: src/code/irrat.lisp +msgid "Natural log of 2d0" +msgstr "" + +#: src/code/irrat.lisp +msgid "Natural log of 2w0" +msgstr "" + #: src/code/irrat.lisp msgid "Return e raised to the power NUMBER." msgstr "" View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/3451a113880c01328881c90… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/3451a113880c01328881c90… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl] Pushed new branch issue-428-add-constants-for-irrat
by Raymond Toy (＠rtoy) 18 Jan '26

18 Jan '26

Raymond Toy pushed new branch issue-428-add-constants-for-irrat at cmucl / cmucl -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/tree/issue-428-add-constants-f… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-464-lisp-version-uses-branch-name] Update sed regexp for macos
by Raymond Toy (＠rtoy) 16 Jan '26

16 Jan '26

Raymond Toy pushed to branch issue-464-lisp-version-uses-branch-name at cmucl / cmucl Commits: 377248f2 by Raymond Toy at 2026-01-16T08:15:20-08:00 Update sed regexp for macos Use "*" instead of "+" in sed regex because macos sed doesn't recognize "+". - - - - - 1 changed file: - bin/git-version.sh Changes: ===================================== bin/git-version.sh ===================================== @@ -41,7 +41,7 @@ else # we have to remove everything before the first slash which # contains things like "tag/" or "head/". GIT_DESC="`git describe --all --dirty || git describe 2>/dev/null`" - GIT_HASH="`echo ${GIT_DESC} | sed 's;^[^/]\+/;;' 2>/dev/null`" + GIT_HASH="`echo ${GIT_DESC} | sed 's;^[^/]*/;;' 2>/dev/null`" BRANCH="`git rev-parse --abbrev-ref HEAD`*" case "$GIT_HASH" in View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/377248f29a553f46b5f64e8… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/377248f29a553f46b5f64e8… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-464-lisp-version-uses-branch-name] Replace backtick with single quote; add some double quotes
by Raymond Toy (＠rtoy) 16 Jan '26

16 Jan '26

Raymond Toy pushed to branch issue-464-lisp-version-uses-branch-name at cmucl / cmucl Commits: 6dac9de2 by Raymond Toy at 2026-01-15T18:18:23-08:00 Replace backtick with single quote; add some double quotes We aren't trying to run a command so use single quotes instead of backticks in the comment. Add some double quotes around shell var. - - - - - 1 changed file: - bin/git-version.sh Changes: ===================================== bin/git-version.sh ===================================== @@ -8,7 +8,7 @@ git-version.sh [-hfv] -f The version is printed as a C file #define expression. Otherwise, the version is just printed to stdout - -v Use this as the version instead of using `git describe`, which + -v Use this as the version instead of using 'git describe', which is the default Determine the version of cmucl. By @@ -40,7 +40,7 @@ else # name as appropriate. This is much more informative. However, # we have to remove everything before the first slash which # contains things like "tag/" or "head/". - GIT_DESC="`(git describe --all --dirty || git describe 2>/dev/null)`" + GIT_DESC="`git describe --all --dirty || git describe 2>/dev/null`" GIT_HASH="`echo ${GIT_DESC} | sed 's;^[^/]\+/;;' 2>/dev/null`" BRANCH="`git rev-parse --abbrev-ref HEAD`*" @@ -75,7 +75,7 @@ else fi if [ -z "$FILE" ]; then - echo $DEFAULT_VERSION + echo "$DEFAULT_VERSION" else cat <<EOF /* View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/6dac9de2d8325668589140d… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/6dac9de2d8325668589140d… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-464-lisp-version-uses-branch-name] 2 commits: Prefix the version with CI if it was built via CI
by Raymond Toy (＠rtoy) 16 Jan '26

16 Jan '26

Raymond Toy pushed to branch issue-464-lisp-version-uses-branch-name at cmucl / cmucl Commits: 47db53d7 by Raymond Toy at 2026-01-15T18:08:03-08:00 Prefix the version with CI if it was built via CI - - - - - f41a1277 by Raymond Toy at 2026-01-15T18:13:10-08:00 Use CASE statement instead of sequence of if/elif This makes it easier to read and lets us do some globbing to match the git hash which is easier than using regexp with `expr`. - - - - - 1 changed file: - bin/git-version.sh Changes: ===================================== bin/git-version.sh ===================================== @@ -44,22 +44,29 @@ else GIT_HASH="`echo ${GIT_DESC} | sed 's;^[^/]\+/;;' 2>/dev/null`" BRANCH="`git rev-parse --abbrev-ref HEAD`*" - if [ `expr "X$GIT_HASH" : 'Xsnapshot-[0-9][0-9][0-9][0-9]-[01][0-9]'` != 0 ]; then - # The git hash looks like snapshot-yyyy-mm-<stuff>. Remove the - # "snapshot-" part. - DEFAULT_VERSION=`expr "$GIT_HASH" : "snapshot-$.*$"` - elif [ `expr "X$GIT_HASH" : 'X[0-9][0-9][a-f]'` != 0 ]; then - # The git hash looks like a release which is 3 hex digits. Use it as is. - DEFAULT_VERSION="${GIT_HASH}" - elif [ `expr "${GIT_HASH}" : "${BRANCH}"` != 0 ]; then - # The hash looks like the current branch with possibly more - # stuff at the end. Use the hash as the version. - DEFAULT_VERSION="${GIT_HASH}" - elif [ `expr "${GIT_HASH}" : "[0-9]*"` != 0 ]; then - # Assuming this is CI which seems to produce a githash like - # "pipeline/<digits>". - DEFAULT_VERSION="${GIT_HASH}" - fi + case "$GIT_HASH" in + snapshot-[0-9][0-9][0-9][0-9]-[01][0-9]) + # The git hash looks like snapshot-yyyy-mm-<stuff>. Remove the + # "snapshot-" part. + DEFAULT_VERSION=`expr "$GIT_HASH" : "snapshot-$.*$"` + ;; + [0-9][0-9][a-f]) + # The git hash looks like a release which is 3 hex digits. + # Use it as is. + DEFAULT_VERSION="${GIT_HASH}" + ;; + ${BRANCH}*) + # The hash looks like the current branch with possibly more + # stuff at the end. Use the hash as the version. + DEFAULT_VERSION="${GIT_HASH}" + ;; + [0-9]*) + # Assuming this is CI which seems to produce a githash like + # "pipeline/<digits>". Make the version include "ci-" so we + # know this was done via CI. + DEFAULT_VERSION="ci-${GIT_HASH}" + ;; + esac if [ -z "$DEFAULT_VERSION" ]; then echo "Unable to determine a default version from git describe: $GIT_DESC" View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/800cf0dbf497725bbb432f… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/800cf0dbf497725bbb432f… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0