mirrors/ecl
1
Fork 0
mirror of https://gitlab.com/embeddable-common-lisp/ecl.git synced 2026-01-14 05:12:38 -08:00
Code Issues Projects Releases Packages Wiki Activity

Use si:trap-fpe also when computing the small and large constants

Browse source
This commit is contained in:
Juan Jose Garcia Ripoll 2009-06-09 11:49:04 +02:00
parent d6e98053f6
commit c9b2bfbf82
1 changed files with 3 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

4
src/lsp/numlib.lsp
View file

@ -34,6 +34,7 @@
(epsilon- (x)
(/= (float 1 x) (- (float 1 x) x))))
`(eval-when (compile load eval)
(let ((bits (si::trap-fpe 'last nil)))
(defconstant short-float-epsilon
,(binary-search #'epsilon+ (coerce 0 'short-float) (coerce 1 'short-float))
"The smallest postive short-float E that satisfies
@ -65,7 +66,8 @@
(defconstant long-float-negative-epsilon
,(binary-search #'epsilon- (coerce 0 'long-float) (coerce 1 'long-float))
"The smallest positive long-float E that satisfies
(not (= (float 1 E) (- (float 1 E) E)))")))
(not (= (float 1 E) (- (float 1 E) E)))")
(si::trap-fpe 'last nil))))
#+IEEE-FLOATING-POINT
(locally