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

[Git][cmucl/cmucl][issue-425-correctly-rounded-math-functions] 3 commits: Define CORE64_SRCS instead of CORE64_OBJS
by Raymond Toy (＠rtoy) 28 Jan '26

28 Jan '26

Raymond Toy pushed to branch issue-425-correctly-rounded-math-functions at cmucl / cmucl Commits: 6f100c06 by Raymond Toy at 2026-01-28T08:21:07-08:00 Define CORE64_SRCS instead of CORE64_OBJS Define the list of sources for core-math routines. The rules in GNUMakefile handle converting the list of sources into a list of object files. Then we can clean up the makefile a bit to reduce unneeded complexity. - - - - - 08f23a38 by Raymond Toy at 2026-01-28T08:22:47-08:00 Use fdlibm_scalbn even when core-math is used. Core-math doesn't have a scalbn function, but fdlibm does so use that. We could use libm's scalbn, but that doesn't signal errors on overflow like we want. - - - - - 34837462 by Raymond Toy at 2026-01-28T08:28:42-08:00 Disable some tests for core-math This is so we can pass CI. We'll need to come back and look at these again after we decide what we want to do with these failures for core-math which pass when using fdlibm. - - - - - 4 changed files: - src/lisp/Config.x86_common - src/lisp/GNUmakefile - src/lisp/irrat.c - tests/fdlibm.lisp Changes: ===================================== src/lisp/Config.x86_common ===================================== @@ -42,29 +42,29 @@ endif ifdef FEATURE_CORE_MATH CORE_MATH_64=core-math/src/binary64 -CORE64_OBJS=\ - $(CORE_MATH_64)/sin/sin.o \ - $(CORE_MATH_64)/cos/cos.o \ - $(CORE_MATH_64)/tan/tan.o \ - $(CORE_MATH_64)/atan2/atan2.o \ - $(CORE_MATH_64)/asin/asin.o \ - $(CORE_MATH_64)/acos/acos.o \ - $(CORE_MATH_64)/atan/atan.o \ - $(CORE_MATH_64)/sinh/sinh.o \ - $(CORE_MATH_64)/cosh/cosh.o \ - $(CORE_MATH_64)/tanh/tanh.o \ - $(CORE_MATH_64)/asinh/asinh.o \ - $(CORE_MATH_64)/acosh/acosh.o \ - $(CORE_MATH_64)/atanh/atanh.o \ - $(CORE_MATH_64)/exp/exp.o \ - $(CORE_MATH_64)/log/log.o \ - $(CORE_MATH_64)/log10/log10.o \ - $(CORE_MATH_64)/log2/log2.o \ - $(CORE_MATH_64)/pow/pow.o \ - $(CORE_MATH_64)/hypot/hypot.o \ - $(CORE_MATH_64)/log1p/log1p.o \ - $(CORE_MATH_64)/expm1/expm1.o \ - $(CORE_MATH_64)/sincos/sincos.o +CORE64_SRCS=\ + $(CORE_MATH_64)/sin/sin.c \ + $(CORE_MATH_64)/cos/cos.c \ + $(CORE_MATH_64)/tan/tan.c \ + $(CORE_MATH_64)/atan2/atan2.c \ + $(CORE_MATH_64)/asin/asin.c \ + $(CORE_MATH_64)/acos/acos.c \ + $(CORE_MATH_64)/atan/atan.c \ + $(CORE_MATH_64)/sinh/sinh.c \ + $(CORE_MATH_64)/cosh/cosh.c \ + $(CORE_MATH_64)/tanh/tanh.c \ + $(CORE_MATH_64)/asinh/asinh.c \ + $(CORE_MATH_64)/acosh/acosh.c \ + $(CORE_MATH_64)/atanh/atanh.c \ + $(CORE_MATH_64)/exp/exp.c \ + $(CORE_MATH_64)/log/log.c \ + $(CORE_MATH_64)/log10/log10.c \ + $(CORE_MATH_64)/log2/log2.c \ + $(CORE_MATH_64)/pow/pow.c \ + $(CORE_MATH_64)/hypot/hypot.c \ + $(CORE_MATH_64)/log1p/log1p.c \ + $(CORE_MATH_64)/expm1/expm1.c \ + $(CORE_MATH_64)/sincos/sincos.c endif ifeq ($(filter 2% 3%, $(shell $(CC) -dumpversion)),) ===================================== src/lisp/GNUmakefile ===================================== @@ -36,7 +36,7 @@ SRCS = lisp.c coreparse.c alloc.c monitor.c print.c interr.c \ dynbind.c breakpoint.c regnames.c backtrace.c save.c purify.c \ runprog.c time.c case-mapping.c exec-init.c \ irrat.c \ - ${CORE64_OBJS:.o=.c} \ + ${CORE64_SRCS} \ ${FDLIBM} ${ARCH_SRC} ${ASSEM_SRC} ${OS_SRC} ${GC_SRC} @@ -61,8 +61,8 @@ lisp: ${OBJS} ${CORE64_OBJS} version.o # Create a library out of all the object files so we can build an # executable. However, we need to remove exec-init.o from the library -lisp.a: version.o ${OBJS} ${CORE64_OBJS} ${EXEC_FINAL_OBJ} - $(AR) crs lisp.a ${OBJS} ${CORE64_OBJS} version.o +lisp.a: version.o ${OBJS} ${EXEC_FINAL_OBJ} + $(AR) crs lisp.a ${OBJS} version.o ifneq (${EXEC_FINAL_OBJ},) $(AR) d lisp.a exec-init.o endif ===================================== src/lisp/irrat.c ===================================== @@ -7,6 +7,8 @@ #include "lisp.h" #include "internals.h" +#include "fdlibm.h" + #ifdef FEATURE_CORE_MATH #include <math.h> extern double cr_sin(double); @@ -29,8 +31,6 @@ extern double cr_pow(double, double); extern double cr_hypot(double, double); extern double cr_log1p(double); extern double cr_expm1(double); -#else -#include "fdlibm.h" #endif @@ -241,9 +241,5 @@ lisp_expm1(double x) double lisp_scalbn(double x, int n) { -#ifdef FEATURE_CORE_MATH - return scalbn(x, n); -#else return fdlibm_scalbn(x, n); -#endif } ===================================== tests/fdlibm.lisp ===================================== @@ -171,6 +171,7 @@ (ext:with-float-traps-masked (:invalid) (assert-true (ext::float-nan-p (kernel:%expm1 *snan*)))) ;; expm1(x) = -1 for x < -56*log(2), signaling inexact + #-core-math (let ((x (* -57 (log 2d0)))) (ext:with-float-traps-enabled (:inexact) (assert-error 'floating-point-inexact @@ -217,6 +218,7 @@ (ext:with-float-traps-masked (:overflow) (assert-equal ext:double-float-positive-infinity (kernel:%exp 710d0))) + #-core-math (let ((modes (ext:get-floating-point-modes))) (unwind-protect (progn View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/66d306ffceda908c6d56b8… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/66d306ffceda908c6d56b8… 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] 6 commits: Fix #466: Add C wrapper for the special functions
by Raymond Toy (＠rtoy) 28 Jan '26

28 Jan '26

Raymond Toy pushed to branch issue-425-correctly-rounded-math-functions at cmucl / cmucl Commits: 7f708431 by Raymond Toy at 2026-01-27T08:45:32-08:00 Fix #466: Add C wrapper for the special functions - - - - - 318d0560 by Raymond Toy at 2026-01-27T08:45:32-08:00 Merge branch 'issue-466-c-wrapper-specfun' into 'master' Fix #466: Add C wrapper for the special functions Closes #466 See merge request cmucl/cmucl!345 - - - - - 8ceef426 by Raymond Toy at 2026-01-27T18:42:27-08:00 Merge branch 'master' into issue-425-correctly-rounded-math-functions - - - - - 5a960d14 by Raymond Toy at 2026-01-27T19:44:00-08:00 Fix up incorrect merge We need to get rid of most of the #+core-math stuff. - - - - - f364498b by Raymond Toy at 2026-01-27T20:01:10-08:00 Add core-math routines to SRCS Just for completeness. - - - - - 66d306ff by Raymond Toy at 2026-01-27T20:02:09-08:00 Call the core-math routines when FEATURE_CORE_MATH is set Otherwise, we use the old fdlibm routines. - - - - - 5 changed files: - src/code/irrat.lisp - src/lisp/GNUmakefile - src/lisp/fdlibm.h - + src/lisp/irrat.c - src/lisp/log2.c Changes: ===================================== src/code/irrat.lisp ===================================== @@ -75,73 +75,38 @@ ;;; Please refer to the Unix man pages for details about these routines. -#+core-math -(progn -;;; Trigonometric. -(def-math-rtn ("cr_sin" %sin) 1) -(def-math-rtn ("cr_cos" %cos) 1) -(def-math-rtn ("cr_tan" %tan) 1) -(def-math-rtn ("cr_atan" %atan) 1) -(def-math-rtn ("cr_atan2" %atan2) 2) -(def-math-rtn ("cr_asin" %asin) 1) -(def-math-rtn ("cr_acos" %acos) 1) -(def-math-rtn ("cr_sinh" %sinh) 1) -(def-math-rtn ("cr_cosh" %cosh) 1) -(def-math-rtn ("cr_tanh" %tanh) 1) -(def-math-rtn ("cr_asinh" %asinh) 1) -(def-math-rtn ("cr_acosh" %acosh) 1) -(def-math-rtn ("cr_atanh" %atanh) 1) - -;;; Exponential and Logarithmic. -(def-math-rtn ("cr_exp" %exp) 1) -(def-math-rtn ("cr_log" %log) 1) -(def-math-rtn ("cr_log10" %log10) 1) -(def-math-rtn ("cr_log2" %log2) 1) - -(def-math-rtn ("__ieee754_pow" %pow) 2) -#-(or x86 sparc-v7 sparc-v8 sparc-v9) -(def-math-rtn "sqrt" 1) -(def-math-rtn ("cr_hypot" %hypot) 2) - -(def-math-rtn ("cr_log1p" %log1p) 1) -(def-math-rtn ("cr_expm1" %expm1) 1) -) - -#-core-math -(progn ;;; Trigonometric. -(def-math-rtn ("fdlibm_sin" %sin) 1) -(def-math-rtn ("fdlibm_cos" %cos) 1) -(def-math-rtn ("fdlibm_tan" %tan) 1) -(def-math-rtn ("fdlibm_atan" %atan) 1) -(def-math-rtn ("__ieee754_atan2" %atan2) 2) -(def-math-rtn ("__ieee754_asin" %asin) 1) -(def-math-rtn ("__ieee754_acos" %acos) 1) -(def-math-rtn ("__ieee754_sinh" %sinh) 1) -(def-math-rtn ("__ieee754_cosh" %cosh) 1) -(def-math-rtn ("fdlibm_tanh" %tanh) 1) -(def-math-rtn ("fdlibm_asinh" %asinh) 1) -(def-math-rtn ("__ieee754_acosh" %acosh) 1) -(def-math-rtn ("__ieee754_atanh" %atanh) 1) +(def-math-rtn ("lisp_sin" %sin) 1) +(def-math-rtn ("lisp_cos" %cos) 1) +(def-math-rtn ("lisp_tan" %tan) 1) +(def-math-rtn ("lisp_atan" %atan) 1) +(def-math-rtn ("lisp_atan2" %atan2) 2) +(def-math-rtn ("lisp_asin" %asin) 1) +(def-math-rtn ("lisp_acos" %acos) 1) +(def-math-rtn ("lisp_sinh" %sinh) 1) +(def-math-rtn ("lisp_cosh" %cosh) 1) +(def-math-rtn ("lisp_tanh" %tanh) 1) +(def-math-rtn ("lisp_asinh" %asinh) 1) +(def-math-rtn ("lisp_acosh" %acosh) 1) +(def-math-rtn ("lisp_atanh" %atanh) 1) ;;; Exponential and Logarithmic. -(def-math-rtn ("__ieee754_exp" %exp) 1) -(def-math-rtn ("__ieee754_log" %log) 1) -(def-math-rtn ("__ieee754_log10" %log10) 1) -(def-math-rtn ("cmucl_log2" %log2) 1) +(def-math-rtn ("lisp_exp" %exp) 1) +(def-math-rtn ("lisp_log" %log) 1) +(def-math-rtn ("lisp_log10" %log10) 1) +(def-math-rtn ("lisp_log2" %log2) 1) -(def-math-rtn ("__ieee754_pow" %pow) 2) +(def-math-rtn ("lisp_pow" %pow) 2) #-(or x86 sparc-v7 sparc-v8 sparc-v9) (def-math-rtn "sqrt" 1) -(def-math-rtn ("__ieee754_hypot" %hypot) 2) +(def-math-rtn ("lisp_hypot" %hypot) 2) -(def-math-rtn ("fdlibm_log1p" %log1p) 1) -(def-math-rtn ("fdlibm_expm1" %expm1) 1) -) +(def-math-rtn ("lisp_log1p" %log1p) 1) +(def-math-rtn ("lisp_expm1" %expm1) 1) (declaim (inline %scalbn)) (export '%scalbn) -(alien:def-alien-routine ("fdlibm_scalbn" %scalbn) c-call:double +(alien:def-alien-routine ("lisp_scalbn" %scalbn) c-call:double (x double-float) (n c-call:int)) ===================================== src/lisp/GNUmakefile ===================================== @@ -35,6 +35,8 @@ SRCS = lisp.c coreparse.c alloc.c monitor.c print.c interr.c \ vars.c parse.c interrupt.c search.c validate.c globals.c \ dynbind.c breakpoint.c regnames.c backtrace.c save.c purify.c \ runprog.c time.c case-mapping.c exec-init.c \ + irrat.c \ + ${CORE64_OBJS:.o=.c} \ ${FDLIBM} ${ARCH_SRC} ${ASSEM_SRC} ${OS_SRC} ${GC_SRC} ===================================== src/lisp/fdlibm.h ===================================== @@ -54,8 +54,21 @@ extern double fdlibm_tan(double x); extern double fdlibm_expm1(double x); extern double fdlibm_log1p(double x); extern double fdlibm_atan(double x); +extern double fdlibm_tanh(double x); +extern double fdlibm_asinh(double x); extern double __ieee754_exp(double x); extern double __ieee754_log(double x); +extern double __ieee754_atan2(double y, double x); +extern double __ieee754_asin(double x); +extern double __ieee754_acos(double x); +extern double __ieee754_sinh(double x); +extern double __ieee754_cosh(double x); +extern double __ieee754_atanh(double x); +extern double __ieee754_acosh(double x); +extern double __ieee754_log10(double x); +extern double __ieee754_pow(double x, double y); +extern double __ieee754_hypot(double x, double y); +extern double fdlibm_scalbn(double x, int n); enum FDLIBM_EXCEPTION { FDLIBM_DIVIDE_BY_ZERO, ===================================== src/lisp/irrat.c ===================================== @@ -0,0 +1,249 @@ +/* + + This code was written as part of the CMU Common Lisp project at + Carnegie Mellon University, and has been placed in the public domain. + +*/ + +#include "lisp.h" +#include "internals.h" +#ifdef FEATURE_CORE_MATH +#include <math.h> +extern double cr_sin(double); +extern double cr_cos(double); +extern double cr_tan(double); +extern double cr_atan(double); +extern double cr_atan2(double, double); +extern double cr_asin(double); +extern double cr_acos(double); +extern double cr_sinh(double); +extern double cr_cosh(double); +extern double cr_tanh(double); +extern double cr_asinh(double); +extern double cr_acosh(double); +extern double cr_atanh(double); +extern double cr_exp(double); +extern double cr_log(double); +extern double cr_log10(double); +extern double cr_pow(double, double); +extern double cr_hypot(double, double); +extern double cr_log1p(double); +extern double cr_expm1(double); +#else +#include "fdlibm.h" +#endif + + +/* + * Wrappers for the special functions + */ + +double +lisp_sin(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_sin(x); +#else + return fdlibm_sin(x); +#endif +} + +double +lisp_cos(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_cos(x); +#else + return fdlibm_cos(x); +#endif +} + +double +lisp_tan(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_tan(x); +#else + return fdlibm_tan(x); +#endif +} + +double +lisp_atan(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_atan(x); +#else + return fdlibm_atan(x); +#endif +} + +double +lisp_atan2(double y, double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_atan2(y, x); +#else + return __ieee754_atan2(y, x); +#endif +} + +double +lisp_asin(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_asin(x); +#else + return __ieee754_asin(x); +#endif +} + +double +lisp_acos(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_acos(x); +#else + return __ieee754_acos(x); +#endif +} + +double +lisp_sinh(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_sinh(x); +#else + return __ieee754_sinh(x); +#endif +} + +double +lisp_cosh(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_cosh(x); +#else + return __ieee754_cosh(x); +#endif +} + +double +lisp_tanh(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_tanh(x); +#else + return fdlibm_tanh(x); +#endif +} + +double +lisp_asinh(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_asinh(x); +#else + return fdlibm_asinh(x); +#endif +} + +double +lisp_acosh(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_acosh(x); +#else + return __ieee754_acosh(x); +#endif +} + +double +lisp_atanh(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_atanh(x); +#else + return __ieee754_atanh(x); +#endif +} + +double +lisp_exp(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_exp(x); +#else + return __ieee754_exp(x); +#endif +} + +double +lisp_log(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_log(x); +#else + return __ieee754_log(x); +#endif +} + +double +lisp_log10(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_log10(x); +#else + return __ieee754_log10(x); +#endif +} + +double +lisp_pow(double x, double y) +{ +#ifdef FEATURE_CORE_MATH + return cr_pow(x, y); +#else + return __ieee754_pow(x, y); +#endif +} + +double +lisp_hypot(double x, double y) +{ +#ifdef FEATURE_CORE_MATH + return cr_hypot(x, y); +#else + return __ieee754_hypot(x, y); +#endif +} + +double +lisp_log1p(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_log1p(x); +#else + return fdlibm_log1p(x); +#endif +} + +double +lisp_expm1(double x) +{ +#ifdef FEATURE_CORE_MATH + return cr_expm1(x); +#else + return fdlibm_expm1(x); +#endif +} + +double +lisp_scalbn(double x, int n) +{ +#ifdef FEATURE_CORE_MATH + return scalbn(x, n); +#else + return fdlibm_scalbn(x, n); +#endif +} ===================================== src/lisp/log2.c ===================================== @@ -39,7 +39,7 @@ cp = 9.61796693925975554329e-01, /* 0x3FEEC709, 0xDC3A03FD =2/(3ln2) */ cp_h = 9.61796700954437255859e-01, /* 0x3FEEC709, 0xE0000000 =(float)cp */ cp_l = -7.02846165095275826516e-09; /* 0xBE3E2FE0, 0x145B01F5 =tail of cp_h*/ -double cmucl_log2(double x) +double lisp_log2(double x) { double ax; int k, hx, lx, ix; View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/a5440432136bbcc0de718c… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/a5440432136bbcc0de718c… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][master] 2 commits: Fix #466: Add C wrapper for the special functions
by Raymond Toy (＠rtoy) 27 Jan '26

27 Jan '26

Raymond Toy pushed to branch master at cmucl / cmucl Commits: 7f708431 by Raymond Toy at 2026-01-27T08:45:32-08:00 Fix #466: Add C wrapper for the special functions - - - - - 318d0560 by Raymond Toy at 2026-01-27T08:45:32-08:00 Merge branch 'issue-466-c-wrapper-specfun' into 'master' Fix #466: Add C wrapper for the special functions Closes #466 See merge request cmucl/cmucl!345 - - - - - 5 changed files: - src/code/irrat.lisp - src/lisp/GNUmakefile - src/lisp/fdlibm.h - + src/lisp/irrat.c - src/lisp/log2.c Changes: ===================================== src/code/irrat.lisp ===================================== @@ -76,37 +76,37 @@ ;;; Please refer to the Unix man pages for details about these routines. ;;; Trigonometric. -(def-math-rtn ("fdlibm_sin" %sin) 1) -(def-math-rtn ("fdlibm_cos" %cos) 1) -(def-math-rtn ("fdlibm_tan" %tan) 1) -(def-math-rtn ("fdlibm_atan" %atan) 1) -(def-math-rtn ("__ieee754_atan2" %atan2) 2) -(def-math-rtn ("__ieee754_asin" %asin) 1) -(def-math-rtn ("__ieee754_acos" %acos) 1) -(def-math-rtn ("__ieee754_sinh" %sinh) 1) -(def-math-rtn ("__ieee754_cosh" %cosh) 1) -(def-math-rtn ("fdlibm_tanh" %tanh) 1) -(def-math-rtn ("fdlibm_asinh" %asinh) 1) -(def-math-rtn ("__ieee754_acosh" %acosh) 1) -(def-math-rtn ("__ieee754_atanh" %atanh) 1) +(def-math-rtn ("lisp_sin" %sin) 1) +(def-math-rtn ("lisp_cos" %cos) 1) +(def-math-rtn ("lisp_tan" %tan) 1) +(def-math-rtn ("lisp_atan" %atan) 1) +(def-math-rtn ("lisp_atan2" %atan2) 2) +(def-math-rtn ("lisp_asin" %asin) 1) +(def-math-rtn ("lisp_acos" %acos) 1) +(def-math-rtn ("lisp_sinh" %sinh) 1) +(def-math-rtn ("lisp_cosh" %cosh) 1) +(def-math-rtn ("lisp_tanh" %tanh) 1) +(def-math-rtn ("lisp_asinh" %asinh) 1) +(def-math-rtn ("lisp_acosh" %acosh) 1) +(def-math-rtn ("lisp_atanh" %atanh) 1) ;;; Exponential and Logarithmic. -(def-math-rtn ("__ieee754_exp" %exp) 1) -(def-math-rtn ("__ieee754_log" %log) 1) -(def-math-rtn ("__ieee754_log10" %log10) 1) -(def-math-rtn ("cmucl_log2" %log2) 1) +(def-math-rtn ("lisp_exp" %exp) 1) +(def-math-rtn ("lisp_log" %log) 1) +(def-math-rtn ("lisp_log10" %log10) 1) +(def-math-rtn ("lisp_log2" %log2) 1) -(def-math-rtn ("__ieee754_pow" %pow) 2) +(def-math-rtn ("lisp_pow" %pow) 2) #-(or x86 sparc-v7 sparc-v8 sparc-v9) (def-math-rtn "sqrt" 1) -(def-math-rtn ("__ieee754_hypot" %hypot) 2) +(def-math-rtn ("lisp_hypot" %hypot) 2) -(def-math-rtn ("fdlibm_log1p" %log1p) 1) -(def-math-rtn ("fdlibm_expm1" %expm1) 1) +(def-math-rtn ("lisp_log1p" %log1p) 1) +(def-math-rtn ("lisp_expm1" %expm1) 1) (declaim (inline %scalbn)) (export '%scalbn) -(alien:def-alien-routine ("fdlibm_scalbn" %scalbn) c-call:double +(alien:def-alien-routine ("lisp_scalbn" %scalbn) c-call:double (x double-float) (n c-call:int)) ===================================== src/lisp/GNUmakefile ===================================== @@ -35,6 +35,7 @@ SRCS = lisp.c coreparse.c alloc.c monitor.c print.c interr.c \ vars.c parse.c interrupt.c search.c validate.c globals.c \ dynbind.c breakpoint.c regnames.c backtrace.c save.c purify.c \ runprog.c time.c case-mapping.c exec-init.c \ + irrat.c \ ${FDLIBM} ${ARCH_SRC} ${ASSEM_SRC} ${OS_SRC} ${GC_SRC} OBJS = $(patsubst %.c,%.o,$(patsubst %.S,%.o,$(patsubst %.s,%.o,$(SRCS)))) ===================================== src/lisp/fdlibm.h ===================================== @@ -54,8 +54,21 @@ extern double fdlibm_tan(double x); extern double fdlibm_expm1(double x); extern double fdlibm_log1p(double x); extern double fdlibm_atan(double x); +extern double fdlibm_tanh(double x); +extern double fdlibm_asinh(double x); extern double __ieee754_exp(double x); extern double __ieee754_log(double x); +extern double __ieee754_atan2(double y, double x); +extern double __ieee754_asin(double x); +extern double __ieee754_acos(double x); +extern double __ieee754_sinh(double x); +extern double __ieee754_cosh(double x); +extern double __ieee754_atanh(double x); +extern double __ieee754_acosh(double x); +extern double __ieee754_log10(double x); +extern double __ieee754_pow(double x, double y); +extern double __ieee754_hypot(double x, double y); +extern double fdlibm_scalbn(double x, int n); enum FDLIBM_EXCEPTION { FDLIBM_DIVIDE_BY_ZERO, ===================================== src/lisp/irrat.c ===================================== @@ -0,0 +1,138 @@ +/* + + This code was written as part of the CMU Common Lisp project at + Carnegie Mellon University, and has been placed in the public domain. + +*/ + +#include "fdlibm.h" + +/* + * Wrappers for the special functions + */ + +double +lisp_sin(double x) +{ + return fdlibm_sin(x); +} + +double +lisp_cos(double x) +{ + return fdlibm_cos(x); +} + +double +lisp_tan(double x) +{ + return fdlibm_tan(x); +} + +double +lisp_atan(double x) +{ + return fdlibm_atan(x); +} + +double +lisp_atan2(double y, double x) +{ + return __ieee754_atan2(y, x); +} + +double +lisp_asin(double x) +{ + return __ieee754_asin(x); +} + +double +lisp_acos(double x) +{ + return __ieee754_acos(x); +} + +double +lisp_sinh(double x) +{ + return __ieee754_sinh(x); +} + +double +lisp_cosh(double x) +{ + return __ieee754_cosh(x); +} + +double +lisp_tanh(double x) +{ + return fdlibm_tanh(x); +} + +double +lisp_asinh(double x) +{ + return fdlibm_asinh(x); +} + +double +lisp_acosh(double x) +{ + return __ieee754_acosh(x); +} + +double +lisp_atanh(double x) +{ + return __ieee754_atanh(x); +} + +double +lisp_exp(double x) +{ + return __ieee754_exp(x); +} + +double +lisp_log(double x) +{ + return __ieee754_log(x); +} + +double +lisp_log10(double x) +{ + return __ieee754_log10(x); +} + +double +lisp_pow(double x, double y) +{ + return __ieee754_pow(x, y); +} + +double +lisp_hypot(double x, double y) +{ + return __ieee754_hypot(x, y); +} + +double +lisp_log1p(double x) +{ + return fdlibm_log1p(x); +} + +double +lisp_expm1(double x) +{ + return fdlibm_expm1(x); +} + +double +lisp_scalbn(double x, int n) +{ + return fdlibm_scalbn(x, n); +} ===================================== src/lisp/log2.c ===================================== @@ -39,7 +39,7 @@ cp = 9.61796693925975554329e-01, /* 0x3FEEC709, 0xDC3A03FD =2/(3ln2) */ cp_h = 9.61796700954437255859e-01, /* 0x3FEEC709, 0xE0000000 =(float)cp */ cp_l = -7.02846165095275826516e-09; /* 0xBE3E2FE0, 0x145B01F5 =tail of cp_h*/ -double cmucl_log2(double x) +double lisp_log2(double x) { double ax; int k, hx, lx, ix; View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/ea927341913119e29bf891… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/ea927341913119e29bf891… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-466-c-wrapper-specfun] Fix typo
by Raymond Toy (＠rtoy) 27 Jan '26

27 Jan '26

Raymond Toy pushed to branch issue-466-c-wrapper-specfun at cmucl / cmucl Commits: e2cb39cc by Raymond Toy at 2026-01-27T07:47:05-08:00 Fix typo - - - - - 1 changed file: - src/code/irrat.lisp Changes: ===================================== src/code/irrat.lisp ===================================== @@ -49,7 +49,7 @@ (string-upcase name))))) `(progn (declaim (inline ,lisp-name)) - (export ,lisp-name) + (export ',lisp-name) (alien:def-alien-routine (,c-name ,lisp-name) double-float ,@(let ((results nil)) (dotimes (i num-args (nreverse results)) View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/e2cb39cc083c3156b760522… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/e2cb39cc083c3156b760522… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-466-c-wrapper-specfun] Finish reversion of single-float changes
by Raymond Toy (＠rtoy) 27 Jan '26

27 Jan '26

Raymond Toy pushed to branch issue-466-c-wrapper-specfun at cmucl / cmucl Commits: c55685ee by Raymond Toy at 2026-01-27T07:42:55-08:00 Finish reversion of single-float changes Remove the parens in the export in def-math-rtn. This returns us to what was there originally instead of being gratuitously different. - - - - - 1 changed file: - src/code/irrat.lisp Changes: ===================================== src/code/irrat.lisp ===================================== @@ -49,7 +49,7 @@ (string-upcase name))))) `(progn (declaim (inline ,lisp-name)) - (export '(,lisp-name)) + (export ,lisp-name) (alien:def-alien-routine (,c-name ,lisp-name) double-float ,@(let ((results nil)) (dotimes (i num-args (nreverse results)) View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/c55685eeac58664c350a026… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/c55685eeac58664c350a026… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-466-c-wrapper-specfun] Address review comments
by Raymond Toy (＠rtoy) 27 Jan '26

27 Jan '26

Raymond Toy pushed to branch issue-466-c-wrapper-specfun at cmucl / cmucl Commits: e2d1c242 by Raymond Toy at 2026-01-27T07:29:36-08:00 Address review comments Remove the changes for single-float which will be done in a different MR. - - - - - 4 changed files: - src/code/exports.lisp - src/code/irrat.lisp - src/compiler/float-tran.lisp - src/lisp/irrat.c Changes: ===================================== src/code/exports.lisp ===================================== @@ -2184,28 +2184,7 @@ "STANDARD-PPRINT-DISPATCH-TABLE-MODIFIED-ERROR" "%IEEE754-REM-PI/2" - "%SINCOS" - "%SINF" - "%COSF" - "%TANF" - "%ATANF" - "%ATAN2F" - "%ASINF" - "%ACOSF" - "%SINHF" - "%COSHF" - "%TANHF" - "%ASINHF" - "%ACOSHF" - "%ATANHF" - "%EXPF" - "%LOGF" - "%LOG10F" - "%LOG2F" - "%POWF" - "%HYPOTF" - "%LOG1PF" - "%EXPM1F") + "%SINCOS") #+heap-overflow-check (:export "DYNAMIC-SPACE-OVERFLOW-WARNING-HIT" "DYNAMIC-SPACE-OVERFLOW-ERROR-HIT" ===================================== src/code/irrat.lisp ===================================== @@ -47,29 +47,21 @@ (intern (concatenate 'simple-string "%" (string-upcase name))))) - (let ((c-name-f (concatenate 'simple-string c-name "f")) - (lisp-name-f (symbolicate lisp-name "F"))) `(progn - (declaim (inline ,lisp-name ,lisp-name-f)) - (export '(,lisp-name ,lisp-name-f)) + (declaim (inline ,lisp-name)) + (export '(,lisp-name)) (alien:def-alien-routine (,c-name ,lisp-name) double-float ,@(let ((results nil)) (dotimes (i num-args (nreverse results)) (push (list (intern (format nil "ARG-~D" i)) 'double-float) - results)))) - (alien:def-alien-routine (,c-name-f ,lisp-name-f) single-float - ,@(let ((results nil)) - (dotimes (i num-args (nreverse results)) - (push (list (intern (format nil "ARG-~D" i)) - 'single-float) - results)))))))) + results))))))) (eval-when (compile load eval) (defun handle-reals (function var) `((((foreach fixnum single-float bignum ratio)) - (,(symbolicate function "F") (coerce ,var 'single-float))) + (coerce (,function (coerce ,var 'double-float)) 'single-float)) ((double-float) (,function ,var)) #+double-double ===================================== src/compiler/float-tran.lisp ===================================== @@ -620,84 +620,42 @@ (double-float) double-float (movable foldable flushable)) -(defknown (%tanf %sinhf %asinhf %atanhf %logf %log10f) - (single-float) single-float - (movable foldable flushable)) - (defknown (%sin %cos %tanh #+x87 %sin-quick #+x87 %cos-quick) (double-float) (double-float -1.0d0 1.0d0) (movable foldable flushable)) -(defknown (%sinf %cosf %tanhf) - (single-float) (single-float -1.0f0 1.0f0) - (movable foldable flushable)) - (defknown (%asin %atan) (double-float) (double-float #.(- (/ pi 2)) #.(/ pi 2)) (movable foldable flushable)) -(defknown (%asinf %atanf) - (single-float) (single-float #.(coerce (- (/ pi 2)) 'single-float) - #.(coerce (/ pi 2) 'single-float)) - (movable foldable flushable)) - (defknown (%acos) (double-float) (double-float 0.0d0 #.pi) (movable foldable flushable)) -(defknown (%acosf) - (single-float) (single-float 0.0f0 #.(coerce pi 'single-float)) - (movable foldable flushable)) - (defknown (%cosh) (double-float) (double-float 1.0d0) (movable foldable flushable)) -(defknown (%coshf) - (single-float) (single-float 1.0f0) - (movable foldable flushable)) - (defknown (%acosh %exp %sqrt) (double-float) (double-float 0.0d0) (movable foldable flushable)) -(defknown (%acoshf %expf) - (single-float) (single-float 0.0f0) - (movable foldable flushable)) - (defknown %expm1 (double-float) (double-float -1d0) (movable foldable flushable)) -(defknown %expm1f - (single-float) (single-float -1f0) - (movable foldable flushable)) - (defknown (%hypot) (double-float double-float) (double-float 0d0) (movable foldable flushable)) -(defknown (%hypotf) - (single-float single-float) (single-float 0f0) - (movable foldable flushable)) - (defknown (%pow) (double-float double-float) double-float (movable foldable flushable)) -(defknown (%powf) - (single-float single-float) single-float - (movable foldable flushable)) - (defknown (%atan2) (double-float double-float) (double-float #.(- pi) #.pi) (movable foldable flushable)) -(defknown (%atan2f) - (single-float single-float) (single-float #.(coerce (- pi) 'single-float) - #.(coerce pi 'single-float)) - (movable foldable flushable)) - (defknown (%scalb) (double-float double-float) double-float (movable foldable flushable)) @@ -710,12 +668,9 @@ (double-float) double-float (movable foldable flushable)) -(defknown (%log1pf) - (single-float) single-float - (movable foldable flushable)) - (dolist (stuff '((exp %exp *) (log %log float) + (sqrt %sqrt float) (sin %sin float) (cos %cos float) (tan %tan float) @@ -729,19 +684,10 @@ (acosh %acosh float) (atanh %atanh float))) (destructuring-bind (name prim rtype) stuff - (let ((primf (symbolicate prim "F"))) - (deftransform name ((x) '(single-float) rtype :eval-name t) - `(,primf x)) - (deftransform name ((x) '(double-float) rtype :eval-name t :when :both) - `(,prim x))))) - -(deftransform sqrt ((x) (double-float) double-float :when :both) - `(%sqrt x)) - -;; We don't currently have sqrt specialized for single-floats, so use -;; the double-float version. -(deftransform sqrt ((x) (single-float) single-float) - `(coerce (%sqrt (coerce x 'double-float)) 'single-float)) + (deftransform name ((x) '(single-float) rtype :eval-name t) + `(coerce (,prim (coerce x 'double-float)) 'single-float)) + (deftransform name ((x) '(double-float) rtype :eval-name t :when :both) + `(,prim x)))) (defknown (%sincos) (double-float) (values double-float double-float) ===================================== src/lisp/irrat.c ===================================== @@ -136,126 +136,3 @@ lisp_scalbn(double x, int n) { return fdlibm_scalbn(x, n); } - -/* - * Wrappers for the single-float versions - */ -float -lisp_sinf(float x) -{ - return (float) fdlibm_sin((double) x); -} - -float -lisp_cosf(float x) -{ - return (float) fdlibm_cos((double) x); -} - -float -lisp_tanf(float x) -{ - return (float) fdlibm_tan((double) x); -} - -float -lisp_atanf(float x) -{ - return (float) fdlibm_atan((double) x); -} - -float -lisp_atan2f(float y, float x) -{ - return (float) __ieee754_atan2((double) y, (double) x); -} - -float -lisp_asinf(float x) -{ - return (float) __ieee754_asin((double) x); -} - -float -lisp_acosf(float x) -{ - return (float) __ieee754_acos((double) x); -} - -float -lisp_sinhf(float x) -{ - return (float) __ieee754_sinh((double) x); -} - -float -lisp_coshf(float x) -{ - return (float) __ieee754_cosh((double) x); -} - -float -lisp_tanhf(float x) -{ - return (float) fdlibm_tanh((double) x); -} - -float -lisp_asinhf(float x) -{ - return (float) fdlibm_asinh((double) x); -} - -float -lisp_acoshf(float x) -{ - return (float) __ieee754_acosh((double) x); -} - -float -lisp_atanhf(float x) -{ - return (float) __ieee754_atanh((double) x); -} - -float -lisp_expf(float x) -{ - return (float) __ieee754_exp((double) x); -} - -float -lisp_logf(float x) -{ - return (float) __ieee754_log((double) x); -} - -float -lisp_log10f(float x) -{ - return (float) __ieee754_log10((double) x); -} - -float -lisp_powf(float x, float y) -{ - return (float) __ieee754_pow((double) x, (double) y); -} - -float -lisp_hypotf(float x, float y) -{ - return (float) __ieee754_hypot((double) x, (double) y); -} - -float -lisp_log1pf(float x) -{ - return (float) fdlibm_log1p((double) x); -} - -float -lisp_expm1f(float x) -{ - return (float) fdlibm_expm1((double) x); -} View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/e2d1c2429cf7b0de17fbc44… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/e2d1c2429cf7b0de17fbc44… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-466-c-wrapper-specfun] Handle sqrt specially since we don't have a sqrtf function
by Raymond Toy (＠rtoy) 27 Jan '26

27 Jan '26

Raymond Toy pushed to branch issue-466-c-wrapper-specfun at cmucl / cmucl Commits: 5eda6726 by Raymond Toy at 2026-01-26T17:30:13-08:00 Handle sqrt specially since we don't have a sqrtf function For now, sqrt of a single-float is converted to %sqrt that handles double-floats. Obvious conversions done handle single-float. - - - - - 1 changed file: - src/compiler/float-tran.lisp Changes: ===================================== src/compiler/float-tran.lisp ===================================== @@ -716,7 +716,6 @@ (dolist (stuff '((exp %exp *) (log %log float) - (sqrt %sqrt float) (sin %sin float) (cos %cos float) (tan %tan float) @@ -736,6 +735,14 @@ (deftransform name ((x) '(double-float) rtype :eval-name t :when :both) `(,prim x))))) +(deftransform sqrt ((x) (double-float) double-float :when :both) + `(%sqrt x)) + +;; We don't currently have sqrt specialized for single-floats, so use +;; the double-float version. +(deftransform sqrt ((x) (single-float) single-float) + `(coerce (%sqrt (coerce x 'double-float)) 'single-float)) + (defknown (%sincos) (double-float) (values double-float double-float) (movable foldable flushable)) View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/5eda6726b04abf3a826c726… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/5eda6726b04abf3a826c726… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][issue-466-c-wrapper-specfun] Support single-float versions of the special functions
by Raymond Toy (＠rtoy) 26 Jan '26

26 Jan '26

Raymond Toy pushed to branch issue-466-c-wrapper-specfun at cmucl / cmucl Commits: c970737c by Raymond Toy at 2026-01-26T12:58:33-08:00 Support single-float versions of the special functions Define C wrapper for the single-float versions of the special functions. This just converts the arg to a double and then converts the result to a single. Update `def-math-rtn` to generate the single-float interface by appending an "f" to the the lisp name as well as the C name. Export all the new single-float functions. Finally, update the deftransfrom to call the single-float version and add defknowns for the new single-float functions. - - - - - 4 changed files: - src/code/exports.lisp - src/code/irrat.lisp - src/compiler/float-tran.lisp - src/lisp/irrat.c Changes: ===================================== src/code/exports.lisp ===================================== @@ -2184,7 +2184,28 @@ "STANDARD-PPRINT-DISPATCH-TABLE-MODIFIED-ERROR" "%IEEE754-REM-PI/2" - "%SINCOS") + "%SINCOS" + "%SINF" + "%COSF" + "%TANF" + "%ATANF" + "%ATAN2F" + "%ASINF" + "%ACOSF" + "%SINHF" + "%COSHF" + "%TANHF" + "%ASINHF" + "%ACOSHF" + "%ATANHF" + "%EXPF" + "%LOGF" + "%LOG10F" + "%LOG2F" + "%POWF" + "%HYPOTF" + "%LOG1PF" + "%EXPM1F") #+heap-overflow-check (:export "DYNAMIC-SPACE-OVERFLOW-WARNING-HIT" "DYNAMIC-SPACE-OVERFLOW-ERROR-HIT" ===================================== src/code/irrat.lisp ===================================== @@ -47,21 +47,29 @@ (intern (concatenate 'simple-string "%" (string-upcase name))))) + (let ((c-name-f (concatenate 'simple-string c-name "f")) + (lisp-name-f (symbolicate lisp-name "F"))) `(progn - (declaim (inline ,lisp-name)) - (export ',lisp-name) + (declaim (inline ,lisp-name ,lisp-name-f)) + (export '(,lisp-name ,lisp-name-f)) (alien:def-alien-routine (,c-name ,lisp-name) double-float ,@(let ((results nil)) (dotimes (i num-args (nreverse results)) (push (list (intern (format nil "ARG-~D" i)) 'double-float) - results))))))) + results)))) + (alien:def-alien-routine (,c-name-f ,lisp-name-f) single-float + ,@(let ((results nil)) + (dotimes (i num-args (nreverse results)) + (push (list (intern (format nil "ARG-~D" i)) + 'single-float) + results)))))))) (eval-when (compile load eval) (defun handle-reals (function var) `((((foreach fixnum single-float bignum ratio)) - (coerce (,function (coerce ,var 'double-float)) 'single-float)) + (,(symbolicate function "F") (coerce ,var 'single-float))) ((double-float) (,function ,var)) #+double-double ===================================== src/compiler/float-tran.lisp ===================================== @@ -620,42 +620,84 @@ (double-float) double-float (movable foldable flushable)) +(defknown (%tanf %sinhf %asinhf %atanhf %logf %log10f) + (single-float) single-float + (movable foldable flushable)) + (defknown (%sin %cos %tanh #+x87 %sin-quick #+x87 %cos-quick) (double-float) (double-float -1.0d0 1.0d0) (movable foldable flushable)) +(defknown (%sinf %cosf %tanhf) + (single-float) (single-float -1.0f0 1.0f0) + (movable foldable flushable)) + (defknown (%asin %atan) (double-float) (double-float #.(- (/ pi 2)) #.(/ pi 2)) (movable foldable flushable)) +(defknown (%asinf %atanf) + (single-float) (single-float #.(coerce (- (/ pi 2)) 'single-float) + #.(coerce (/ pi 2) 'single-float)) + (movable foldable flushable)) + (defknown (%acos) (double-float) (double-float 0.0d0 #.pi) (movable foldable flushable)) +(defknown (%acosf) + (single-float) (single-float 0.0f0 #.(coerce pi 'single-float)) + (movable foldable flushable)) + (defknown (%cosh) (double-float) (double-float 1.0d0) (movable foldable flushable)) +(defknown (%coshf) + (single-float) (single-float 1.0f0) + (movable foldable flushable)) + (defknown (%acosh %exp %sqrt) (double-float) (double-float 0.0d0) (movable foldable flushable)) +(defknown (%acoshf %expf) + (single-float) (single-float 0.0f0) + (movable foldable flushable)) + (defknown %expm1 (double-float) (double-float -1d0) (movable foldable flushable)) +(defknown %expm1f + (single-float) (single-float -1f0) + (movable foldable flushable)) + (defknown (%hypot) (double-float double-float) (double-float 0d0) (movable foldable flushable)) +(defknown (%hypotf) + (single-float single-float) (single-float 0f0) + (movable foldable flushable)) + (defknown (%pow) (double-float double-float) double-float (movable foldable flushable)) +(defknown (%powf) + (single-float single-float) single-float + (movable foldable flushable)) + (defknown (%atan2) (double-float double-float) (double-float #.(- pi) #.pi) (movable foldable flushable)) +(defknown (%atan2f) + (single-float single-float) (single-float #.(coerce (- pi) 'single-float) + #.(coerce pi 'single-float)) + (movable foldable flushable)) + (defknown (%scalb) (double-float double-float) double-float (movable foldable flushable)) @@ -668,6 +710,10 @@ (double-float) double-float (movable foldable flushable)) +(defknown (%log1pf) + (single-float) single-float + (movable foldable flushable)) + (dolist (stuff '((exp %exp *) (log %log float) (sqrt %sqrt float) @@ -684,10 +730,11 @@ (acosh %acosh float) (atanh %atanh float))) (destructuring-bind (name prim rtype) stuff - (deftransform name ((x) '(single-float) rtype :eval-name t) - `(coerce (,prim (coerce x 'double-float)) 'single-float)) - (deftransform name ((x) '(double-float) rtype :eval-name t :when :both) - `(,prim x)))) + (let ((primf (symbolicate prim "F"))) + (deftransform name ((x) '(single-float) rtype :eval-name t) + `(,primf x)) + (deftransform name ((x) '(double-float) rtype :eval-name t :when :both) + `(,prim x))))) (defknown (%sincos) (double-float) (values double-float double-float) ===================================== src/lisp/irrat.c ===================================== @@ -136,3 +136,126 @@ lisp_scalbn(double x, int n) { return fdlibm_scalbn(x, n); } + +/* + * Wrappers for the single-float versions + */ +float +lisp_sinf(float x) +{ + return (float) fdlibm_sin((double) x); +} + +float +lisp_cosf(float x) +{ + return (float) fdlibm_cos((double) x); +} + +float +lisp_tanf(float x) +{ + return (float) fdlibm_tan((double) x); +} + +float +lisp_atanf(float x) +{ + return (float) fdlibm_atan((double) x); +} + +float +lisp_atan2f(float y, float x) +{ + return (float) __ieee754_atan2((double) y, (double) x); +} + +float +lisp_asinf(float x) +{ + return (float) __ieee754_asin((double) x); +} + +float +lisp_acosf(float x) +{ + return (float) __ieee754_acos((double) x); +} + +float +lisp_sinhf(float x) +{ + return (float) __ieee754_sinh((double) x); +} + +float +lisp_coshf(float x) +{ + return (float) __ieee754_cosh((double) x); +} + +float +lisp_tanhf(float x) +{ + return (float) fdlibm_tanh((double) x); +} + +float +lisp_asinhf(float x) +{ + return (float) fdlibm_asinh((double) x); +} + +float +lisp_acoshf(float x) +{ + return (float) __ieee754_acosh((double) x); +} + +float +lisp_atanhf(float x) +{ + return (float) __ieee754_atanh((double) x); +} + +float +lisp_expf(float x) +{ + return (float) __ieee754_exp((double) x); +} + +float +lisp_logf(float x) +{ + return (float) __ieee754_log((double) x); +} + +float +lisp_log10f(float x) +{ + return (float) __ieee754_log10((double) x); +} + +float +lisp_powf(float x, float y) +{ + return (float) __ieee754_pow((double) x, (double) y); +} + +float +lisp_hypotf(float x, float y) +{ + return (float) __ieee754_hypot((double) x, (double) y); +} + +float +lisp_log1pf(float x) +{ + return (float) fdlibm_log1p((double) x); +} + +float +lisp_expm1f(float x) +{ + return (float) fdlibm_expm1((double) x); +} View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/c970737c75c9273bf89cfb8… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/commit/c970737c75c9273bf89cfb8… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl] Pushed new branch issue-466-c-wrapper-specfun
by Raymond Toy (＠rtoy) 26 Jan '26

26 Jan '26

Raymond Toy pushed new branch issue-466-c-wrapper-specfun at cmucl / cmucl -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/tree/issue-466-c-wrapper-specf… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0

[Git][cmucl/cmucl][master] 2 commits: Fix #428: Add some constants for common values in irrat.lisp
by Raymond Toy (＠rtoy) 20 Jan '26

20 Jan '26

Raymond Toy pushed to branch master at cmucl / cmucl Commits: 6d264493 by Raymond Toy at 2026-01-20T07:02:39-08:00 Fix #428: Add some constants for common values in irrat.lisp - - - - - ea927341 by Raymond Toy at 2026-01-20T07:02:40-08:00 Merge branch 'issue-428-add-constants-for-irrat' into 'master' Fix #428: Add some constants for common values in irrat.lisp Closes #428 See merge request cmucl/cmucl!344 - - - - - 3 changed files: - src/code/irrat.lisp - src/general-info/release-22a.md - src/i18n/locale/cmucl.pot Changes: ===================================== src/code/irrat.lisp ===================================== @@ -24,6 +24,18 @@ (defconstant pi 3.14159265358979323846264338327950288419716939937511L0) ;(defconstant e 2.71828182845904523536028747135266249775724709369996L0) +(defconstant +log-2d0+ + (log 2d0) + "Natural log of 2d0") + +(defconstant +log-2w0+ + (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. ;;; @@ -488,14 +500,14 @@ complex) (if (and (zerop base) (plusp (realpart power))) (* base power) - (exp (* power (* (log2 base 1w0) (log 2w0)))))) + (exp (* power (* (log2 base 1w0) +log-2w0+))))) (((foreach fixnum (or bignum ratio) single-float double-float) (foreach (complex double-float))) ;; Result should have double-float accuracy. Use log2 in ;; case the base won't fit in a double-float. (if (and (zerop base) (plusp (realpart power))) (* base power) - (exp (* power (* (log2 base) (log 2d0)))))) + (exp (* power (* (log2 base) +log-2d0+))))) ((double-float (foreach (complex rational) (complex single-float))) (if (and (zerop base) (plusp (realpart power))) @@ -508,7 +520,7 @@ ;; in case the base won't fit in a double-float. (if (and (zerop base) (plusp (realpart power))) (* base power) - (exp (* power (* (log2 base 1w0) (log 2w0)))))) + (exp (* power (* (log2 base 1w0) +log-2w0+))))) (((foreach fixnum (or bignum ratio) single-float) (foreach (complex rational) (complex single-float))) (if (and (zerop base) (plusp (realpart power))) @@ -566,18 +578,18 @@ ;; log(2), with the precision specified by the type of F (number-dispatch ((f real)) ((double-float) - #.(log 2d0)) + +log-2d0+) #+double-double ((double-double-float) - #.(log 2w0)))) + +log-2w0+))) (log-2-pi (f) ;; log(pi), with the precision specified by the type of F (number-dispatch ((f real)) ((double-float) - #.(/ pi (log 2d0))) + #.(/ pi +log-2d0+)) #+double-double ((double-double-float) - #.(/ dd-pi (log 2w0))))) + #.(/ dd-pi +log-2w0+)))) (log1p (x) ;; log(1+x), with the precision specified by the type of ;; X @@ -827,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) @@ -852,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)? @@ -1393,7 +1405,7 @@ This is for use with J /= 0 only when |z| is huge." (let ((t0 #.(/ 1 (sqrt 2.0d0))) (t1 1.2d0) (t2 3d0) - (ln2 #.(log 2d0)) + (ln2 +log-2d0+) (x (float (realpart z) 1.0d0)) (y (float (imagpart z) 1.0d0))) (multiple-value-bind (rho k) ===================================== 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 ===================================== @@ -5008,6 +5008,18 @@ 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 "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/3f733bae60b20dadee5f4f… -- View it on GitLab: https://gitlab.common-lisp.net/cmucl/cmucl/-/compare/3f733bae60b20dadee5f4f… You're receiving this email because of your account on gitlab.common-lisp.net.

1 0