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

[git] CMU Common Lisp branch master updated. snapshot-2014-02-5-gf04a802
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via f04a80287e33b51b08242d53e7ad7c1c2ec9c59e (commit) from 8970ad11b03f863b40a9ee352d3aced52b783cce (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit f04a80287e33b51b08242d53e7ad7c1c2ec9c59e Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Sun Feb 9 09:59:01 2014 -0800 Update from commit logs. diff --git a/src/general-info/release-20f.txt b/src/general-info/release-20f.txt index 628236b..fdc8837 100644 --- a/src/general-info/release-20f.txt +++ b/src/general-info/release-20f.txt @@ -68,6 +68,7 @@ New in this release: * Trac Tickets: * Ticket #90 fixed. + * Ticket #92 fixed. * Other changes: ----------------------------------------------------------------------- Summary of changes: src/general-info/release-20f.txt | 1 + 1 file changed, 1 insertion(+) hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp branch master updated. snapshot-2014-06-13-gc7bf8f8
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via c7bf8f8fe16c0bc1a9c9c65060033046855eca82 (commit) from 203b0b0144dec1a5f481ada67792040e4a485f21 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit c7bf8f8fe16c0bc1a9c9c65060033046855eca82 Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Tue Jul 22 10:28:30 2014 -0700 ppc doesn't have sincos in libm so use our own version. diff --git a/src/lisp/Config.ppc_darwin b/src/lisp/Config.ppc_darwin index 0902b89..a4516a8 100644 --- a/src/lisp/Config.ppc_darwin +++ b/src/lisp/Config.ppc_darwin @@ -72,4 +72,8 @@ OS_LIBS = -lSystem -lc -lm # According to Config.x86_common, this no longer has aliasing # problems, so we don't need any additional compilation options. e_rem_pio2.o : e_rem_pio2.c - $(CC) -c $(CFLAGS) $< + $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< + +OS_SRC += sincos.c +sincos.o : sincos.c + $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< ----------------------------------------------------------------------- Summary of changes: src/lisp/Config.ppc_darwin | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp branch master updated. snapshot-2014-06-12-g203b0b0
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via 203b0b0144dec1a5f481ada67792040e4a485f21 (commit) from 90530eaf3540c66bbcd445535b1cd9941bccbc3e (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 203b0b0144dec1a5f481ada67792040e4a485f21 Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Tue Jul 22 10:25:07 2014 -0700 For consistency with sincos, we need to include the fdlibm trig functions on Darwin/x86 even though Darwin's libm is accurate. diff --git a/src/lisp/Config.x86_darwin b/src/lisp/Config.x86_darwin index 2af8d77..3dcdc61 100644 --- a/src/lisp/Config.x86_darwin +++ b/src/lisp/Config.x86_darwin @@ -18,13 +18,22 @@ OS_LIBS = EXEC_FINAL_OBJ = exec-final.o -OS_SRC += sincos.c k_sin.o k_cos.o +OS_SRC += k_sin.c k_cos.c k_tan.c s_sin.c s_cos.c s_tan.c sincos.c -sincos.o : sincos.c +k_sin.o : k_sin.c + $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< +k_cos.o : k_cos.c + $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< +k_tan.o : k_tan.c $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_sin.o : k_sin.c +s_sin.o : s_sin.c + $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< +s_cos.o : s_cos.c + $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< +s_tan.o : s_tan.c $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_cos.o : k_cos.c +sincos.o : sincos.c $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< + ----------------------------------------------------------------------- Summary of changes: src/lisp/Config.x86_darwin | 17 +++++++++++++---- 1 file changed, 13 insertions(+), 4 deletions(-) hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp annotated tag snapshot-2014-05 created. snapshot-2014-05
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The annotated tag, snapshot-2014-05 has been created at 4d75cc008ac57b0861309b40ad779b0c9ed9938d (tag) tagging fc6ec6d77cdde9c73bb3b6c8a2c5d1d6ced9094f (commit) replaces begin-x87-removal tagged by Raymond Toy on Sat May 3 10:19:25 2014 -0700 - Log ----------------------------------------------------------------- Snapshot 2014-05 -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.13 (Darwin) iEYEABECAAYFAlNlJSQACgkQJ5IjUmgZO7LSXQCgi31+oAc0VszJfYmt4dW+r6nZ sKoAn1NsRHCovupvqZBbJkvwqYHf3eXG =nuX9 -----END PGP SIGNATURE----- Raymond Toy (23): Remove support for x87 cores from build scripts. Remove lisp -fpu command-line switch. Remove -fpu command-line switch. Simplify call-out for sse2. Remove commented-out stuff and clean up comments. Simplify call_into_c more, putting more of the functionality into More x87 removal from sse2-c-call.lisp Comment out x87 stuff from call_into_lisp. Remove commented code and some redundant jmps. Remove support for compiling with x87 versions. Print more digits for % passed since we have 75K tests now. Make sse2 the default and remove the conditioals. Make sse2 the default and remove the conditionals. Make sse2 the default and always call alloc_overflow_sse2. Add comment that 16-byte alignment is required on Darwin. Remove alloc_overflow_x87. Regenerated. Rmove more x87 fpu mode stuff. Remove more sse2 stuff. Don't merge the x87 FP modes with the SSE2 modes when getting and Oops. Only compile float-sse2 on x86 machines! Oops. Other archs still need float.lisp. Ignore fasls in this directory. ----------------------------------------------------------------------- hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp branch master updated. snapshot-2013-12-a-26-gb71b7c8
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via b71b7c854469191e462162c9d11d0c35222a284a (commit) from 712df0bc4e655226bc5c9ed91aa9c875b4a5eb0d (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit b71b7c854469191e462162c9d11d0c35222a284a Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Fri Dec 20 18:12:03 2013 -0800 Increase accuracy of tan near multiples of pi/4. code/irrat-dd.lisp o Make REDUCE-ARG return an extra result. Since the reduction returns 3 double's, return a double-double result and the third double result for extra accuracy o Update dd-%tan and dd-%%tan to take the extra arg. o Add new constant dd-pi/4-lo. o Increase accuracy of tan by using the relationship tan(pi/4-y)=(1-tan(y))/(1+tan(y)). tests/trig.lisp: o Update the allowed error threshold for two tests to reflect the increased accuracy. diff --git a/src/code/irrat-dd.lisp b/src/code/irrat-dd.lisp index 381d678..2eabb58 100644 --- a/src/code/irrat-dd.lisp +++ b/src/code/irrat-dd.lisp @@ -56,6 +56,12 @@ 0.7853981633974483096156608458198757210492923w0 _N"Pi/4") +;; dd-pi/4-lo is such that dd-pi/4-lo + dd-pi/4 is equal to pi/4 to +;; twice the precision of a double-double-float. +(defconstant dd-pi/4-lo + -7.486924524295848886603985669688687133352026408988280860093283232w-34 + ) + ;; log2-c1 and log-c2 are log(2) arranged in such a way that log2-c1 + ;; log2-c2 is log(2) to an accuracy greater than double-double-float. (defconstant log2-c1 @@ -1072,9 +1078,24 @@ pi/4 11001001000011111101101010100010001000010110100011000 010001101001100010 (/ y) y)))) -(defun dd-%%tan (x) - (declare (type double-double-float x)) - (dd-tancot x nil)) +(defun dd-%%tan (x extra) + (declare (type double-double-float x) + (double-float extra)) + (cond ((>= (abs x) 0.6744) + ;; For 0.6744 <= |x| <= pi/4, we want to use the relationship + ;; + ;; tan(x) = tan(pi/4-y) = (1 - tan(y))/(1 + tan(y)) + ;; = 1 - 2*(tan(y) - tan(y)^2)/(1+tan(y)) + (if (minusp x) + (- (dd-%%tan (- x) (- extra))) + (let* ((z (- dd-pi/4 x)) + (w (- dd-pi/4-lo extra)) + (tan (dd-tancot (+ z w) nil))) + (- 1 + (/ (* 2 (- tan (* tan tan))) + (+ 1 tan)))))) + (t + (dd-tancot x nil)))) (declaim (inline %kernel-rem-pi/2)) (alien:def-alien-routine ("__kernel_rem_pio2" %kernel-rem-pi/2) c-call:int @@ -1139,10 +1160,9 @@ pi/4 11001001000011111101101010100010001000010110100011000 010001101001100010 (length parts) 3 (vector-sap two-over-pi)))) - (sum (+ (coerce (aref y 2) 'double-double-float) - (coerce (aref y 1) 'double-double-float) + (sum (+ (coerce (aref y 1) 'double-double-float) (coerce (aref y 0) 'double-double-float)))) - (values n sum)))) + (values n sum (aref y 2))))) (declaim (ftype (function (double-double-float) double-double-float) @@ -1181,15 +1201,15 @@ pi/4 11001001000011111101101010100010001000010110100011000 010001101001100010 dd-%tan)) (defun dd-%tan (x) (declare (double-double-float x)) - (cond ((< (abs x) (/ pi 4)) - (dd-%%tan x)) + (cond ((<= (abs x) (/ pi 4)) + (dd-%%tan x 0d0)) (t ;; Argument reduction needed - (multiple-value-bind (n reduced) + (multiple-value-bind (n reduced extra) (reduce-arg x) (if (evenp n) - (dd-%%tan reduced) - (- (/ (dd-%%tan reduced)))))))) + (dd-%%tan reduced extra) + (- (/ (dd-%%tan reduced extra)))))))) (defun dd-%sincos (x) (declare (double-double-float x)) diff --git a/src/tests/trig.lisp b/src/tests/trig.lisp index b46c904..58d5440 100644 --- a/src/tests/trig.lisp +++ b/src/tests/trig.lisp @@ -376,12 +376,12 @@ (assert-eq t (rel-or-abs-error (tan (* 7/4 kernel:dd-pi)) -1.000000000000000000000000000000001844257310064121018312678894979w0 - 6.467w-33)) + 3.422w-49)) ;; Test for argument reduction with n odd (assert-eq t (rel-or-abs-error (tan (* 9/4 kernel:dd-pi)) 1.000000000000000000000000000000025802415787810837455445433037983w0 - 5.773w-33)) + 0w0)) ;; Test for argument reduction, big value (assert-eq t (rel-or-abs-error (tan (scale-float 1w0 120)) ----------------------------------------------------------------------- Summary of changes: src/code/irrat-dd.lisp | 42 +++++++++++++++++++++++++++++++----------- src/tests/trig.lisp | 4 ++-- 2 files changed, 33 insertions(+), 13 deletions(-) hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp branch master updated. snapshot-2014-06-23-g06ca7d3
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via 06ca7d326f688d40ad8730bbd2faa8ca7813d2f0 (commit) from 3533a13a8cc53e43a7c4d7d27b3e6b8bcc004672 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 06ca7d326f688d40ad8730bbd2faa8ca7813d2f0 Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Sat Jul 26 08:29:22 2014 -0700 Add include guard. diff --git a/src/lisp/fdlibm.h b/src/lisp/fdlibm.h index 61794b3..467daa5 100644 --- a/src/lisp/fdlibm.h +++ b/src/lisp/fdlibm.h @@ -10,6 +10,8 @@ * ==================================================== */ +#ifndef FDLIBM_H_ +#define FDLIBM_H_ /* Sometimes it's necessary to define __LITTLE_ENDIAN explicitly but these catch some common cases. */ @@ -48,3 +50,5 @@ extern double __kernel_tan(double x, double y, int iy); extern double sin(double x); extern double cos(double x); extern double tan(double x); + +#endif ----------------------------------------------------------------------- Summary of changes: src/lisp/fdlibm.h | 4 ++++ 1 file changed, 4 insertions(+) hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp branch master updated. snapshot-2014-08-1-g00d977a
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via 00d977a1b07c525948d1fa0e33caefe1536cdec3 (commit) from 98319835e05b93ea62e17f6bc3da4144f96ac2c5 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 00d977a1b07c525948d1fa0e33caefe1536cdec3 Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Fri Aug 8 19:08:37 2014 -0700 Forgot a few items that were done in the snapshot. Update it now (post snapshot). The wiki entry has the right info. diff --git a/src/general-info/release-20f.txt b/src/general-info/release-20f.txt index b49ecc6..b6ce0d5 100644 --- a/src/general-info/release-20f.txt +++ b/src/general-info/release-20f.txt @@ -73,7 +73,7 @@ New in this release: for atanh in the CLHS. * CLEAR-OUTPUT was not actually doing anything. Now, CLEAR-OUTPUT will discard any buffered data that has not been written out - yet. + yet. * Bugfixes: * Fix error in pi reduction on x87. It was not noticed previously @@ -94,7 +94,9 @@ New in this release: * For Gray streams, CLEAR-OUTPUT was calling STREAM-FORCE-OUTPUT instead of STREAM-CLEAR-OUTPUT, so the output wasn't actually cleared. - + * In some situations KERNEL:DOUBLE-FLOAT-BITS on x86 would cause a + segfault. This has been fixed. + * For linux, motifd is no longer a 64-bit app. * Trac Tickets: * Ticket #90 fixed. @@ -105,6 +107,8 @@ New in this release: * Ticket #98 fixed. * Ticket #104 fixed. * Ticket #101, item 1 fixed. + * Ticket #105, fixed. + * Ticket #84 fixed on x86. * Other changes: ----------------------------------------------------------------------- Summary of changes: src/general-info/release-20f.txt | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp branch master updated. snapshot-2014-06-78-g07a07a7
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via 07a07a799922a073fde0093867fbb4bef756b219 (commit) via b19192f76254973d133297de39c5594cf7438514 (commit) from e5bc2c586e644fc7fe7a1e48d495600c8d27a19f (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 07a07a799922a073fde0093867fbb4bef756b219 Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Sat Aug 2 19:44:32 2014 -0700 Remove the list of fdlibm files from here because they're in GNUMakefile now. diff --git a/src/lisp/Config.ppc_darwin b/src/lisp/Config.ppc_darwin index 4edf0eb..b0768e0 100644 --- a/src/lisp/Config.ppc_darwin +++ b/src/lisp/Config.ppc_darwin @@ -69,36 +69,3 @@ OS_LIBS = -lSystem -lc -lm #darwin-lispadjuster: darwin-lispadjuster.c -# According to Config.x86_common, this no longer has aliasing -# problems, so we don't need any additional compilation options. -e_rem_pio2.o : e_rem_pio2.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -OS_SRC += k_sin.c k_cos.c k_tan.c s_sin.c s_cos.c s_tan.c sincos.c s_log1p.c s_expm1.c e_pow.c - -k_sin.o : k_sin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_cos.o : k_cos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -k_tan.o : k_tan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -s_sin.o : s_sin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_cos.o : s_cos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_tan.o : s_tan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -sincos.o : sincos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -s_log1p.o : s_log1p.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_exmp1.o : s_expm1.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_pow.o : e_pow.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - diff --git a/src/lisp/Config.sparc_sunc b/src/lisp/Config.sparc_sunc index 1578b56..2e0b4e0 100644 --- a/src/lisp/Config.sparc_sunc +++ b/src/lisp/Config.sparc_sunc @@ -28,31 +28,3 @@ CFLAGS += -xlibmieee -O DEPEND_FLAGS = -xM ASFLAGS = $(AS_V8PLUS) OS_LINK_FLAGS = -M /usr/lib/ld/map.noexstk - -OS_SRC += k_sin.c k_cos.c k_tan.c s_sin.c s_cos.c s_tan.c sincos.c s_log1p.c s_expm1.c e_pow.c - -k_sin.o : k_sin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_cos.o : k_cos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_tan.o : k_tan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -s_sin.o : s_sin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_cos.o : s_cos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_tan.o : s_tan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -sincos.o : sincos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -s_log1p.o : s_log1p.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_exmp1.o : s_expm1.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_pow.o : e_pow.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - commit b19192f76254973d133297de39c5594cf7438514 Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Sat Aug 2 19:40:28 2014 -0700 Clean up config files. * GNUMakefile: * List all of the fdlibm files here instead of the individual configs because we use fdlibm everywhere now. * Config.x86_common: * Remove e_rem_pio2 and k_rem_pio2 which is in GNUMakefile now. * Config.x86_darwin: * Remove all of the fdlibm files since they're in GNUMakefile now. * Config.x86_linux: * Remove all of the fdlibm files since they're in GNUMakefile now. diff --git a/src/lisp/Config.x86_common b/src/lisp/Config.x86_common index f263d40..1873f7d 100644 --- a/src/lisp/Config.x86_common +++ b/src/lisp/Config.x86_common @@ -59,19 +59,7 @@ ASFLAGS = -g ASSEM_SRC = x86-assem.S ARCH_SRC = x86-arch.c -OS_SRC = os-common.c e_rem_pio2.c k_rem_pio2.c +OS_SRC = os-common.c NM = nm -gp DEPEND_FLAGS = -MM - -# This no longer has aliasing problems, so no need to use -# -fno-strict-aliasing anymore. However, if we're building with x87, -# we MUST use -ffloat-store to get proper double-float rounding. -e_rem_pio2.o : e_rem_pio2.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -k_rem_pio2.o : k_rem_pio2.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - - - diff --git a/src/lisp/Config.x86_darwin b/src/lisp/Config.x86_darwin index a76280d..39c09bd 100644 --- a/src/lisp/Config.x86_darwin +++ b/src/lisp/Config.x86_darwin @@ -17,62 +17,3 @@ OS_LINK_FLAGS = -m32 $(MIN_VER) OS_LIBS = EXEC_FINAL_OBJ = exec-final.o - -OS_SRC += k_sin.c k_cos.c k_tan.c s_sin.c s_cos.c s_tan.c sincos.c \ - s_log1p.c s_expm1.c e_pow.c e_exp.c e_log.c \ - e_acos.c e_asin.c s_atan.c \ - e_cosh.c e_sinh.c s_tanh.c \ - e_acosh.c s_asinh.c e_atanh.c \ - e_atan2.c - -k_sin.o : k_sin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_cos.o : k_cos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_tan.o : k_tan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -s_sin.o : s_sin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_cos.o : s_cos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_tan.o : s_tan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -sincos.o : sincos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -s_log1p.o : s_log1p.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_exmp1.o : s_expm1.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_pow.o : e_pow.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_exp.o : e_exp.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_log.o : e_log.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_acos.o : e_acos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_asin.o : e_asin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_atan.o : s_atan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_cosh.o : e_cosh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_sinh.o : e_sinh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_tanh.o : s_tanh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_acosh.o : e_acosh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_asinh.o : s_asinh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_atanh.o : e_atanh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_atan2.o : e_atan2.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< diff --git a/src/lisp/Config.x86_linux b/src/lisp/Config.x86_linux index 8c06bae..213e42c 100644 --- a/src/lisp/Config.x86_linux +++ b/src/lisp/Config.x86_linux @@ -11,62 +11,3 @@ OS_LINK_FLAGS = -m32 -rdynamic -Xlinker --export-dynamic -Xlinker -Map -Xlinker OS_LINK_FLAGS += -Wl,-z,noexecstack EXEC_FINAL_OBJ = exec-final.o - -OS_SRC += k_sin.c k_cos.c k_tan.c s_sin.c s_cos.c s_tan.c sincos.c \ - s_log1p.c s_expm1.c e_pow.c e_exp.c e_log.c \ - e_acos.c e_asin.c s_atan.c \ - e_cosh.c e_sinh.c s_tanh.c \ - e_acosh.c s_asinh.c e_atanh.c \ - e_atan2.c - -k_sin.o : k_sin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_cos.o : k_cos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -k_tan.o : k_tan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -s_sin.o : s_sin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_cos.o : s_cos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_tan.o : s_tan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -sincos.o : sincos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -s_log1p.o : s_log1p.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_exmp1.o : s_expm1.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_pow.o : e_pow.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_exp.o : e_exp.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_log.o : e_log.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_acos.o : e_acos.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_asin.o : e_asin.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_atan.o : s_atan.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_cosh.o : e_cosh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_sinh.o : e_sinh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -s_tanh.o : s_tanh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_acosh.o : e_acosh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_asinh.o : s_asinh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< -e_atanh.o : e_atanh.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< - -e_atan2.o : e_atan2.c - $(CC) -c $(CFLAGS) $(CPPFLAGS) $(CC_REM_PIO2) $< diff --git a/src/lisp/GNUmakefile b/src/lisp/GNUmakefile index a9ab587..55a0b5b 100644 --- a/src/lisp/GNUmakefile +++ b/src/lisp/GNUmakefile @@ -5,11 +5,19 @@ all: lisp.nm -include internals.inc include Config +FDLIBM = k_sin.c k_cos.c k_tan.c s_sin.c s_cos.c s_tan.c sincos.c \ + s_log1p.c s_expm1.c e_pow.c e_exp.c e_log.c \ + e_acos.c e_asin.c s_atan.c \ + e_cosh.c e_sinh.c s_tanh.c \ + e_acosh.c s_asinh.c e_atanh.c \ + e_atan2.c \ + e_rem_pio2.c k_rem_pio2.c + 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 exec-init.c \ - ${ARCH_SRC} ${ASSEM_SRC} ${OS_SRC} ${GC_SRC} + ${FDLIBM} ${ARCH_SRC} ${ASSEM_SRC} ${OS_SRC} ${GC_SRC} OBJS = $(patsubst %.c,%.o,$(patsubst %.S,%.o,$(patsubst %.s,%.o,$(SRCS)))) ----------------------------------------------------------------------- Summary of changes: src/lisp/Config.ppc_darwin | 33 ------------------------- src/lisp/Config.sparc_sunc | 28 --------------------- src/lisp/Config.x86_common | 14 +---------- src/lisp/Config.x86_darwin | 59 -------------------------------------------- src/lisp/Config.x86_linux | 59 -------------------------------------------- src/lisp/GNUmakefile | 10 +++++++- 6 files changed, 10 insertions(+), 193 deletions(-) hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp branch master updated. snapshot-2014-08-39-g074f4d0
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via 074f4d0bac83c1f2f8aff49ec9dd6e23950106fc (commit) from 284043eee5e68291c9f66be8f81cedcf44f7adef (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 074f4d0bac83c1f2f8aff49ec9dd6e23950106fc Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Wed Aug 27 20:23:01 2014 -0700 Move the test cases from compiler/srctran.lisp to tests/srctran.lisp. diff --git a/src/compiler/srctran.lisp b/src/compiler/srctran.lisp index ec37fd4..c2e1a83 100644 --- a/src/compiler/srctran.lisp +++ b/src/compiler/srctran.lisp @@ -2033,41 +2033,6 @@ (list (- limit)) limit) :high (list limit)))))) -#| Test cases -(floor-quotient-bound (make-interval :low 0.3 :high 10.3)) -=> #S(INTERVAL :LOW 0 :HIGH 10) -(floor-quotient-bound (make-interval :low 0.3 :high '(10.3))) -=> #S(INTERVAL :LOW 0 :HIGH 10) -(floor-quotient-bound (make-interval :low 0.3 :high 10)) -=> #S(INTERVAL :LOW 0 :HIGH 10) -(floor-quotient-bound (make-interval :low 0.3 :high '(10))) -=> #S(INTERVAL :LOW 0 :HIGH 9) -(floor-quotient-bound (make-interval :low '(0.3) :high 10.3)) -=> #S(INTERVAL :LOW 0 :HIGH 10) -(floor-quotient-bound (make-interval :low '(0.0) :high 10.3)) -=> #S(INTERVAL :LOW 0 :HIGH 10) -(floor-quotient-bound (make-interval :low '(-1.3) :high 10.3)) -=> #S(INTERVAL :LOW -2 :HIGH 10) -(floor-quotient-bound (make-interval :low '(-1.0) :high 10.3)) -=> #S(INTERVAL :LOW -1 :HIGH 10) -(floor-quotient-bound (make-interval :low -1.0 :high 10.3)) -=> #S(INTERVAL :LOW -1 :HIGH 10) - - -(floor-rem-bound (make-interval :low 0.3 :high 10.3)) -=> #S(INTERVAL :LOW 0 :HIGH '(10.3)) -(floor-rem-bound (make-interval :low 0.3 :high '(10.3))) -=> #S(INTERVAL :LOW 0 :HIGH '(10.3)) -(floor-rem-bound (make-interval :low -10 :high -2.3)) -#S(INTERVAL :LOW (-10) :HIGH 0) -(floor-rem-bound (make-interval :low 0.3 :high 10)) -=> #S(INTERVAL :LOW 0 :HIGH '(10)) -(floor-rem-bound (make-interval :low '(-1.3) :high 10.3)) -=> #S(INTERVAL :LOW '(-10.3) :HIGH '(10.3)) -(floor-rem-bound (make-interval :low '(-20.3) :high 10.3)) -=> #S(INTERVAL :LOW (-20.3) :HIGH (20.3)) -|# - ;;; Same functions for CEILING (defun ceiling-quotient-bound (quot) @@ -2133,45 +2098,7 @@ limit) :high (list limit)))))) -#| Test cases -(ceiling-quotient-bound (make-interval :low 0.3 :high 10.3)) -=> #S(INTERVAL :LOW 1 :HIGH 11) -(ceiling-quotient-bound (make-interval :low 0.3 :high '(10.3))) -=> #S(INTERVAL :LOW 1 :HIGH 11) -(ceiling-quotient-bound (make-interval :low 0.3 :high 10)) -=> #S(INTERVAL :LOW 1 :HIGH 10) -(ceiling-quotient-bound (make-interval :low 0.3 :high '(10))) -=> #S(INTERVAL :LOW 1 :HIGH 10) -(ceiling-quotient-bound (make-interval :low '(0.3) :high 10.3)) -=> #S(INTERVAL :LOW 1 :HIGH 11) -(ceiling-quotient-bound (make-interval :low '(0.0) :high 10.3)) -=> #S(INTERVAL :LOW 1 :HIGH 11) -(ceiling-quotient-bound (make-interval :low '(-1.3) :high 10.3)) -=> #S(INTERVAL :LOW -1 :HIGH 11) -(ceiling-quotient-bound (make-interval :low '(-1.0) :high 10.3)) -=> #S(INTERVAL :LOW 0 :HIGH 11) -(ceiling-quotient-bound (make-interval :low -1.0 :high 10.3)) -=> #S(INTERVAL :LOW -1 :HIGH 11) - - -(ceiling-rem-bound (make-interval :low 0.3 :high 10.3)) -=> #S(INTERVAL :LOW (-10.3) :HIGH 0) -(ceiling-rem-bound (make-interval :low 0.3 :high '(10.3))) -=> #S(INTERVAL :LOW 0 :HIGH '(10.3)) -(ceiling-rem-bound (make-interval :low -10 :high -2.3)) -=> #S(INTERVAL :LOW 0 :HIGH (10)) -(ceiling-rem-bound (make-interval :low 0.3 :high 10)) -=> #S(INTERVAL :LOW (-10) :HIGH 0) -(ceiling-rem-bound (make-interval :low '(-1.3) :high 10.3)) -=> #S(INTERVAL :LOW (-10.3) :HIGH (10.3)) -(ceiling-rem-bound (make-interval :low '(-20.3) :high 10.3)) -=> #S(INTERVAL :LOW (-20.3) :HIGH (20.3)) -|# - - - - (defun truncate-quotient-bound (quot) ;; For positive quotients, truncate is exactly like floor. For ;; negative quotients, truncate is exactly like ceiling. Otherwise, diff --git a/tests/srctran.lisp b/tests/srctran.lisp new file mode 100644 index 0000000..f88e977 --- /dev/null +++ b/tests/srctran.lisp @@ -0,0 +1,82 @@ +;; Tests for various source transformations. + +(defpackage :srctran-tests + (:use :cl :lisp-unit)) + +(in-package "SRCTRAN-TESTS") + +(define-test floor-quotient-bound + "Test the first value of FLOOR returns the correct interval" + (assert-equalp (c::make-interval :low 0 :high 10) + (c::floor-quotient-bound + (c::make-interval :low 0.3 :high 10.3))) + (assert-equalp (c::make-interval :LOW 0 :HIGH 10) + (c::floor-quotient-bound + (c::make-interval :low 0.3 :high '(10.3)))) + (assert-equalp (c::make-interval :LOW 0 :HIGH 10) + (c::floor-quotient-bound (c::make-interval :low 0.3 :high 10))) + (assert-equalp (c::make-interval :LOW 0 :HIGH 9) + (c::floor-quotient-bound (c::make-interval :low 0.3 :high '(10)))) + (assert-equalp (c::make-interval :LOW 0 :HIGH 9) + (c::floor-quotient-bound (c::make-interval :low 0.3 :high '(10)))) + (assert-equalp (c::make-interval :LOW 0 :HIGH 10) + (c::floor-quotient-bound (c::make-interval :low '(0.3) :high 10.3))) + (assert-equalp (c::make-interval :LOW 0 :HIGH 10) + (c::floor-quotient-bound (c::make-interval :low '(0.0) :high 10.3))) + (assert-equalp (c::make-interval :LOW -2 :HIGH 10) + (c::floor-quotient-bound (c::make-interval :low '(-1.3) :high 10.3))) + (assert-equalp (c::make-interval :LOW -1 :HIGH 10) + (c::floor-quotient-bound (c::make-interval :low '(-1.0) :high 10.3))) + (assert-equalp (c::make-interval :LOW -1 :HIGH 10) + (c::floor-quotient-bound (c::make-interval :low -1.0 :high 10.3)))) + +(define-test floor-rem-bound + "Test the second value of FLOOR returns the correct interval" + (assert-equalp (c::make-interval :low 0 :high '(10.3)) + (c::floor-rem-bound (c::make-interval :low 0.3 :high 10.3))) + (assert-equalp (c::make-interval :low 0 :high '(10.3)) + (c::floor-rem-bound (c::make-interval :low 0.3 :high '(10.3)))) + (assert-equalp (c::make-interval :low '(-10) :high 0) + (c::floor-rem-bound (c::make-interval :low -10 :high -2.3))) + (assert-equalp (c::make-interval :low 0 :high '(10)) + (c::floor-rem-bound (c::make-interval :low 0.3 :high 10))) + (assert-equalp (c::make-interval :low '(-10.3) :high '(10.3)) + (c::floor-rem-bound (c::make-interval :low '(-1.3) :high 10.3))) + (assert-equalp (c::make-interval :low '(-20.3) :high '(20.3)) + (c::floor-rem-bound (c::make-interval :low '(-20.3) :high 10.3)))) + +(define-test ceiling-quotient-bound + "Test the first value of CEILING returns the correct interval" + (assert-equalp (c::make-interval :low 1 :high 11) + (c::ceiling-quotient-bound (c::make-interval :low 0.3 :high 10.3))) + (assert-equalp (c::make-interval :low 1 :high 11) + (c::ceiling-quotient-bound (c::make-interval :low 0.3 :high '(10.3)))) + (assert-equalp (c::make-interval :low 1 :high 10) + (c::ceiling-quotient-bound (c::make-interval :low 0.3 :high 10))) + (assert-equalp (c::make-interval :low 1 :high 10) + (c::ceiling-quotient-bound (c::make-interval :low 0.3 :high '(10)))) + (assert-equalp (c::make-interval :low 1 :high 11) + (c::ceiling-quotient-bound (c::make-interval :low '(0.3) :high 10.3))) + (assert-equalp (c::make-interval :low 1 :high 11) + (c::ceiling-quotient-bound (c::make-interval :low '(0.0) :high 10.3))) + (assert-equalp (c::make-interval :low -1 :high 11) + (c::ceiling-quotient-bound (c::make-interval :low '(-1.3) :high 10.3))) + (assert-equalp (c::make-interval :low 0 :high 11) + (c::ceiling-quotient-bound (c::make-interval :low '(-1.0) :high 10.3))) + (assert-equalp (c::make-interval :low -1 :high 11) + (c::ceiling-quotient-bound (c::make-interval :low -1.0 :high 10.3)))) + +(define-test ceiling-rem-bound + "Test the second value of CEILING returns the correct interval" + (assert-equalp (c::make-interval :low '(-10.3) :high 0) + (c::ceiling-rem-bound (c::make-interval :low 0.3 :high 10.3))) + (assert-equalp (c::make-interval :low '(-10.3) :high 0) + (c::ceiling-rem-bound (c::make-interval :low 0.3 :high '(10.3)))) + (assert-equalp (c::make-interval :low 0 :high '(10)) + (c::ceiling-rem-bound (c::make-interval :low -10 :high -2.3))) + (assert-equalp (c::make-interval :low '(-10) :high 0) + (c::ceiling-rem-bound (c::make-interval :low 0.3 :high 10))) + (assert-equalp (c::make-interval :low '(-10.3) :high '(10.3)) + (c::ceiling-rem-bound (c::make-interval :low '(-1.3) :high 10.3))) + (assert-equalp (c::make-interval :low '(-20.3) :high '(20.3)) + (c::ceiling-rem-bound (c::make-interval :low '(-20.3) :high 10.3)))) ----------------------------------------------------------------------- Summary of changes: src/compiler/srctran.lisp | 73 ---------------------------------------- tests/srctran.lisp | 82 +++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 82 insertions(+), 73 deletions(-) create mode 100644 tests/srctran.lisp hooks/post-receive -- CMU Common Lisp

1 0

[git] CMU Common Lisp branch master updated. snapshot-2013-12-a-35-g4e3589d
by rtoy＠common-lisp.net 08 Apr '15

08 Apr '15

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "CMU Common Lisp". The branch, master has been updated via 4e3589d9d175382fe85b6efe2e3e94ccf598a88e (commit) via 1ba45c8177e4a6ef66f7fac3a3870f477a3d03eb (commit) from a9598e606edc4490b8b58aedb4ecdca64b6ffed9 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 4e3589d9d175382fe85b6efe2e3e94ccf598a88e Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Sun Dec 22 22:41:03 2013 -0800 More documentation. o Derive the value of asin(2), lest we forget again why it is the way it is. o Use the formulas from the CLHS in the docstrings. diff --git a/src/code/irrat.lisp b/src/code/irrat.lisp index 078e56f..ff23a5e 100644 --- a/src/code/irrat.lisp +++ b/src/code/irrat.lisp @@ -1916,8 +1916,43 @@ Z may be any number, but the result is always a complex." (realpart sqrt-z+1)))))))) +;; What is the value of asin(2)? Here is a derivation. +;; +;; asin(2) = -i*log(i*2+sqrt(-3)) +;; = -i*log(2*i+sqrt(3)*i) +;; = -i*(log(2+sqrt(3)) + i*pi/2) +;; = pi/2 - i*log(2+sqrt(3)) +;; +;; Note that this differs from asin(2+0.0*i) because we support signed +;; zeroes. +;; +;; asin(2+0*i) = -i*log(i*(2+0*i) + sqrt(1-(4+0*i))) +;; = -i*log((2*i - 0) + sqrt(-3-0*i)) +;; = -i*log(-0 + 2*i - sqrt(3)*i) +;; = -i*log(-0 + i*(2-sqrt(3))) +;; = -i*(log(2-sqrt(3)) + i*pi/2) +;; = pi/2 - i*log(2-sqrt(3)) +;; = pi/2 + i*log(2+sqrt(3)) +;; +;; The last equation follows because (2-sqrt(3)) = 1/(2+sqrt(3)). +;; Hence asin(2) /= asin(2+0*i). +;; +;; Also +;; +;; asin(2-0*i) = -i*log(i*(2-0*i) + sqrt(1-(4-0*i))) +;; = -i*log((2*i + 0) + sqrt(-3+0*i)) +;; = -i*log(0 + 2*i + sqrt(3)*i) +;; = -i*log(0 + i*(2+sqrt(3))) +;; = -i*(log(2+sqrt(3)) + i*pi/2) +;; = pi/2 + i*log(2+sqrt(3)) +;; +;; Hence asin(2) = asin(2-0.0*i). +;; +;; Similar derivations will show that asin(-2) = asin(-2 + 0.0*i) and +;; asin(-2+0.0*i) is different from asin(-2-0.0*i) because of the +;; branch cut, of course. (defun complex-asin (z) - "Compute asin z = asinh(i*z)/i + "Compute asin z = -i*log(i*z + sqrt(1-z^2)) Z may be any number, but the result is always a complex." (declare (number z)) @@ -1950,7 +1985,7 @@ Z may be any number, but the result is always a complex." (- (realpart result))))) (defun complex-atan (z) - "Compute atan z = atanh (i*z) / i + "Compute atan z = (log(1+i*z) - log(1-i*z))/(2*i) Z may be any number, but the result is always a complex." (declare (number z)) commit 1ba45c8177e4a6ef66f7fac3a3870f477a3d03eb Author: Raymond Toy <toy.raymond(a)gmail.com> Date: Sun Dec 22 22:39:37 2013 -0800 Use the formulas from the CLHS in the docstrings. diff --git a/src/code/irrat-dd.lisp b/src/code/irrat-dd.lisp index 2eabb58..024ec93 100644 --- a/src/code/irrat-dd.lisp +++ b/src/code/irrat-dd.lisp @@ -1854,7 +1854,7 @@ Z may be any number, but the result is always a complex." (defun dd-complex-asin (z) - _N"Compute asin z = asinh(i*z)/i + _N"Compute asin z = -i*log(i*z + sqrt(1-z^2)) Z may be any number, but the result is always a complex." (declare (number z)) @@ -1894,7 +1894,7 @@ Z may be any number, but the result is always a complex." (- (realpart result))))) (defun dd-complex-atan (z) - _N"Compute atan z = atanh (i*z) / i + _N"Compute atan z = (log(1+i*z) - log(1-i*z))/(2*i) Z may be any number, but the result is always a complex." (declare (number z)) ----------------------------------------------------------------------- Summary of changes: src/code/irrat-dd.lisp | 4 ++-- src/code/irrat.lisp | 39 +++++++++++++++++++++++++++++++++++++-- 2 files changed, 39 insertions(+), 4 deletions(-) hooks/post-receive -- CMU Common Lisp

1 0