Ensure SingleFloat is 32-bit, DoubleFloat 64-bit.

* interp/types.boot: Tidy. * interp/sys-constants.boot: Define limits for builtin numeric types. * interp/sys-os.boot: Import sys-constants. * interp/foam_l.lisp: Import sys-constants. Align builtin numeric types with OpenAxiom runtime's. * interp/fortcall.boot: Tidy. * interp/i-output.boot: Likewise. * interp/Makefile.pamphlet (foam_l.$(FASLEXT)): Now require sys-constants.$(FASLEXT). * algebra/si.spad.pamphlet: Use $ShortMinimum and $ShortMaximum for SingleInteger limits. * algebra/sf.spad.pamphlet: Use $DoubleFloatMaximum, $DoubleFloatMinimum for DoubleFloat limits. * algebra/strap: Update cached Lisp translation.
author: dos-reis <gdr@axiomatics.org> 2008-09-07 05:41:35 +0000
committer: dos-reis <gdr@axiomatics.org> 2008-09-07 05:41:35 +0000
commit: 18b7b4b5d3511211092796bd8d593d5c9debec02 (patch)
tree: a90da73966bbb0311bf51c62e9ab8356465be694
parent: 52642d936f26a18f7d3818fcabe3fad2bff157ea (diff)
download: open-axiom-18b7b4b5d3511211092796bd8d593d5c9debec02.tar.gz
17 files changed, 160 insertions, 97 deletions
diff --git a/TODO b/TODO
index 0552b2cb..b0ae4949 100644
--- a/TODO
+++ b/TODO
@@ -17,9 +17,6 @@
 * Improve Boot, Spad
 * Integrate to development environments
 
-* Revisit the underlying definition of the types DFlo in foam
-  interface, and DoubleFloat in OpenAxiom.
-
 * Audit all codes that manipulate sockets.
 
 * Automate the process of "optimized" OpenAxiom: this requires a two-pass
@@ -49,6 +46,10 @@
 === DONE ===
 ============
 
+* Revisit the underlying definition of the types DFlo in foam
+  interface, and DoubleFloat in OpenAxiom.  Now, we ensure that
+  SingleFloat is 32-bit, and DoubleFloat is 64-bit.
+
 * Support more Lisp systems (GCL, SBCL, CLisp)
 
 * cleanup the Makefile pamphlets
diff --git a/src/ChangeLog b/src/ChangeLog
index 28fbb4b9..aaac1a7b 100644
--- a/src/ChangeLog
+++ b/src/ChangeLog
@@ -1,3 +1,21 @@
+2008-09-07  Gabriel Dos Reis  <gdr@cs.tamu.edu>
+
+	Ensure SingleFloat is 32-bit, DoubleFloat 64-bit.
+	* interp/types.boot: Tidy.
+	* interp/sys-constants.boot: Define limits for builtin numeric types.
+	* interp/sys-os.boot: Import sys-constants.
+	* interp/foam_l.lisp: Import sys-constants.  Align builtin numeric
+	types with OpenAxiom runtime's.
+	* interp/fortcall.boot: Tidy.
+	* interp/i-output.boot: Likewise.
+	* interp/Makefile.pamphlet (foam_l.$(FASLEXT)): Now require
+	sys-constants.$(FASLEXT).
+	* algebra/si.spad.pamphlet: Use $ShortMinimum and $ShortMaximum
+	for SingleInteger limits.
+	* algebra/sf.spad.pamphlet: Use $DoubleFloatMaximum,
+	$DoubleFloatMinimum for DoubleFloat limits.
+	* algebra/strap: Update cached Lisp translation.
+
 2008-09-06  Gabriel Dos Reis  <gdr@cs.tamu.edu>
 
 	* include/sockio.h (openaxiom_filedesc): New.  Abstract over
diff --git a/src/algebra/sf.spad.pamphlet b/src/algebra/sf.spad.pamphlet
index 66714b07..ddbedb5b 100644
--- a/src/algebra/sf.spad.pamphlet
+++ b/src/algebra/sf.spad.pamphlet
@@ -481,13 +481,13 @@ DoubleFloat(): Join(FloatingPointSystem, DifferentialRing, OpenMath,
      base() = 2 => precision()
      base() = 16 => 4*precision()
      wholePart(precision()*log2(base()::%))::PositiveInteger
-   max()            == MOST_-POSITIVE_-LONG_-FLOAT$Lisp
-   min()            == MOST_-NEGATIVE_-LONG_-FLOAT$Lisp
+   max()            == _$DoubleFloatMaximum$Lisp
+   min()            == _$DoubleFloatMinimum$Lisp
    order(a) == precision() + exponent a - 1
-   0                == FLOAT(0$Lisp,MOST_-POSITIVE_-LONG_-FLOAT$Lisp)$Lisp
-   1                == FLOAT(1$Lisp,MOST_-POSITIVE_-LONG_-FLOAT$Lisp)$Lisp
+   0                == FLOAT(0$Lisp,_$DoubleFloatMaximum$Lisp)$Lisp
+   1                == FLOAT(1$Lisp,_$DoubleFloatMaximum$Lisp)$Lisp
    -- rational approximation to e accurate to 23 digits
-   exp1()           == FLOAT(534625820200,MOST_-POSITIVE_-LONG_-FLOAT$Lisp)$Lisp / FLOAT(196677847971,MOST_-POSITIVE_-LONG_-FLOAT$Lisp)$Lisp
+   exp1()           == FLOAT(534625820200,_$DoubleFloatMaximum$Lisp)$Lisp / FLOAT(196677847971,_$DoubleFloatMaximum$Lisp)$Lisp
    pi()             == PI$Lisp
    coerce(x:%):OutputForm == 
      outputForm(FORMAT(NIL$Lisp,format,x)$Lisp pretend DoubleFloat)
@@ -506,7 +506,7 @@ DoubleFloat(): Join(FloatingPointSystem, DifferentialRing, OpenMath,
    log10 x          == checkComplex log(x)$Lisp
    x:% ** i:Integer == EXPT(x,i)$Lisp
    x:% ** y:%       == checkComplex EXPT(x,y)$Lisp
-   coerce(i:Integer):% == FLOAT(i,MOST_-POSITIVE_-LONG_-FLOAT$Lisp)$Lisp
+   coerce(i:Integer):% == FLOAT(i,_$DoubleFloatMaximum$Lisp)$Lisp
    exp x            == EXP(x)$Lisp
    log x            == checkComplex LN(x)$Lisp
    log2 x           == checkComplex LOG2(x)$Lisp
diff --git a/src/algebra/si.spad.pamphlet b/src/algebra/si.spad.pamphlet
index 65253830..48bf6ac4 100644
--- a/src/algebra/si.spad.pamphlet
+++ b/src/algebra/si.spad.pamphlet
@@ -236,8 +236,8 @@ SingleInteger(): Join(IntegerNumberSystem,Logic,OpenMath) with
  == add
 
    seed : % := 1$Lisp               -- for random()
-   MAXINT ==> MOST_-POSITIVE_-FIXNUM$Lisp
-   MININT ==> MOST_-NEGATIVE_-FIXNUM$Lisp
+   MAXINT ==> _$ShortMaximum$Lisp
+   MININT ==> _$ShortMinimum$Lisp
    BASE ==> 67108864$Lisp           -- 2**26
    MULTIPLIER ==> 314159269$Lisp    -- from Knuth's table
    MODULUS ==> 2147483647$Lisp      -- 2**31-1
diff --git a/src/algebra/strap/DFLOAT.lsp b/src/algebra/strap/DFLOAT.lsp
index 43c27b6f..88f168b4 100644
--- a/src/algebra/strap/DFLOAT.lsp
+++ b/src/algebra/strap/DFLOAT.lsp
@@ -88,31 +88,31 @@
            (|check-subtype| (> #0# 0) '(|PositiveInteger|) #0#))))))) 
 
 (PUT '|DFLOAT;max;$;12| '|SPADreplace|
-     '(XLAM NIL MOST-POSITIVE-LONG-FLOAT)) 
+     '(XLAM NIL |$DoubleFloatMaximum|)) 
 
-(DEFUN |DFLOAT;max;$;12| ($) MOST-POSITIVE-LONG-FLOAT) 
+(DEFUN |DFLOAT;max;$;12| ($) |$DoubleFloatMaximum|) 
 
 (PUT '|DFLOAT;min;$;13| '|SPADreplace|
-     '(XLAM NIL MOST-NEGATIVE-LONG-FLOAT)) 
+     '(XLAM NIL |$DoubleFloatMinimum|)) 
 
-(DEFUN |DFLOAT;min;$;13| ($) MOST-NEGATIVE-LONG-FLOAT) 
+(DEFUN |DFLOAT;min;$;13| ($) |$DoubleFloatMinimum|) 
 
 (DEFUN |DFLOAT;order;$I;14| (|a| $)
   (- (+ (FLOAT-DIGITS 0.0) (SPADCALL |a| (|getShellEntry| $ 28))) 1)) 
 
 (PUT '|DFLOAT;Zero;$;15| '|SPADreplace|
-     '(XLAM NIL (FLOAT 0 MOST-POSITIVE-LONG-FLOAT))) 
+     '(XLAM NIL (FLOAT 0 |$DoubleFloatMaximum|))) 
 
-(DEFUN |DFLOAT;Zero;$;15| ($) (FLOAT 0 MOST-POSITIVE-LONG-FLOAT)) 
+(DEFUN |DFLOAT;Zero;$;15| ($) (FLOAT 0 |$DoubleFloatMaximum|)) 
 
 (PUT '|DFLOAT;One;$;16| '|SPADreplace|
-     '(XLAM NIL (FLOAT 1 MOST-POSITIVE-LONG-FLOAT))) 
+     '(XLAM NIL (FLOAT 1 |$DoubleFloatMaximum|))) 
 
-(DEFUN |DFLOAT;One;$;16| ($) (FLOAT 1 MOST-POSITIVE-LONG-FLOAT)) 
+(DEFUN |DFLOAT;One;$;16| ($) (FLOAT 1 |$DoubleFloatMaximum|)) 
 
 (DEFUN |DFLOAT;exp1;$;17| ($)
-  (/ (FLOAT 534625820200 MOST-POSITIVE-LONG-FLOAT)
-     (FLOAT 196677847971 MOST-POSITIVE-LONG-FLOAT))) 
+  (/ (FLOAT 534625820200 |$DoubleFloatMaximum|)
+     (FLOAT 196677847971 |$DoubleFloatMaximum|))) 
 
 (PUT '|DFLOAT;pi;$;18| '|SPADreplace| '(XLAM NIL PI)) 
 
@@ -179,10 +179,10 @@
   (|DFLOAT;checkComplex| (EXPT |x| |y|) $)) 
 
 (PUT '|DFLOAT;coerce;I$;35| '|SPADreplace|
-     '(XLAM (|i|) (FLOAT |i| MOST-POSITIVE-LONG-FLOAT))) 
+     '(XLAM (|i|) (FLOAT |i| |$DoubleFloatMaximum|))) 
 
 (DEFUN |DFLOAT;coerce;I$;35| (|i| $)
-  (FLOAT |i| MOST-POSITIVE-LONG-FLOAT)) 
+  (FLOAT |i| |$DoubleFloatMaximum|)) 
 
 (PUT '|DFLOAT;exp;2$;36| '|SPADreplace| 'EXP) 
 
@@ -315,7 +315,7 @@
 (DEFUN |DFLOAT;wholePart;$I;71| (|x| $) (FIX |x|)) 
 
 (DEFUN |DFLOAT;float;2IPi$;72| (|ma| |ex| |b| $)
-  (* |ma| (EXPT (FLOAT |b| MOST-POSITIVE-LONG-FLOAT) |ex|))) 
+  (* |ma| (EXPT (FLOAT |b| |$DoubleFloatMaximum|) |ex|))) 
 
 (PUT '|DFLOAT;convert;2$;73| '|SPADreplace| '(XLAM (|x|) |x|)) 
 
@@ -372,7 +372,7 @@
     (RETURN
       (SEQ (LETT |n| (FIX |x|) |DFLOAT;retract;$I;79|)
            (EXIT (COND
-                   ((= |x| (FLOAT |n| MOST-POSITIVE-LONG-FLOAT)) |n|)
+                   ((= |x| (FLOAT |n| |$DoubleFloatMaximum|)) |n|)
                    ('T (|error| "Not an integer")))))))) 
 
 (DEFUN |DFLOAT;retractIfCan;$U;80| (|x| $)
@@ -380,7 +380,7 @@
     (RETURN
       (SEQ (LETT |n| (FIX |x|) |DFLOAT;retractIfCan;$U;80|)
            (EXIT (COND
-                   ((= |x| (FLOAT |n| MOST-POSITIVE-LONG-FLOAT))
+                   ((= |x| (FLOAT |n| |$DoubleFloatMaximum|))
                     (CONS 0 |n|))
                    ('T (CONS 1 "failed")))))))) 
 
@@ -404,17 +404,16 @@
                          (LETT |x| (FLOAT-SIGN 1.0 |x|)
                                |DFLOAT;manexp|)
                          (COND
-                           ((< MOST-POSITIVE-LONG-FLOAT |x|)
+                           ((< |$DoubleFloatMaximum| |x|)
                             (PROGN
                               (LETT #0#
                                     (CONS
                                      (+
                                       (* |s|
-                                       (SPADCALL
-                                        MOST-POSITIVE-LONG-FLOAT
+                                       (SPADCALL |$DoubleFloatMaximum|
                                         (|getShellEntry| $ 27)))
                                       1)
-                                     (SPADCALL MOST-POSITIVE-LONG-FLOAT
+                                     (SPADCALL |$DoubleFloatMaximum|
                                       (|getShellEntry| $ 28)))
                                     |DFLOAT;manexp|)
                               (GO #0#))))
@@ -610,9 +609,9 @@
                                      (SPADCALL |x|
                                       (/
                                        (FLOAT |n|
-                                        MOST-POSITIVE-LONG-FLOAT)
+                                        |$DoubleFloatMaximum|)
                                        (FLOAT |d|
-                                        MOST-POSITIVE-LONG-FLOAT))
+                                        |$DoubleFloatMaximum|))
                                       (|getShellEntry| $ 59)))))))))))
            #0# (EXIT #0#))))) 
 
diff --git a/src/algebra/strap/INS-.lsp b/src/algebra/strap/INS-.lsp
index a342958f..1d3cdcca 100644
--- a/src/algebra/strap/INS-.lsp
+++ b/src/algebra/strap/INS-.lsp
@@ -52,8 +52,7 @@
       (|getShellEntry| $ 30))) 
 
 (DEFUN |INS-;convert;SDf;11| (|x| $)
-  (FLOAT (SPADCALL |x| (|getShellEntry| $ 27))
-         MOST-POSITIVE-LONG-FLOAT)) 
+  (FLOAT (SPADCALL |x| (|getShellEntry| $ 27)) |$DoubleFloatMaximum|)) 
 
 (DEFUN |INS-;convert;SIf;12| (|x| $)
   (SPADCALL (SPADCALL |x| (|getShellEntry| $ 27))
diff --git a/src/algebra/strap/INT.lsp b/src/algebra/strap/INT.lsp
index 61a4dfb8..bee38736 100644
--- a/src/algebra/strap/INT.lsp
+++ b/src/algebra/strap/INT.lsp
@@ -135,10 +135,9 @@
   (SPADCALL |x| (|getShellEntry| $ 45))) 
 
 (PUT '|INT;convert;$Df;24| '|SPADreplace|
-     '(XLAM (|x|) (FLOAT |x| MOST-POSITIVE-LONG-FLOAT))) 
+     '(XLAM (|x|) (FLOAT |x| |$DoubleFloatMaximum|))) 
 
-(DEFUN |INT;convert;$Df;24| (|x| $)
-  (FLOAT |x| MOST-POSITIVE-LONG-FLOAT)) 
+(DEFUN |INT;convert;$Df;24| (|x| $) (FLOAT |x| |$DoubleFloatMaximum|)) 
 
 (DEFUN |INT;convert;$If;25| (|x| $)
   (SPADCALL |x| (|getShellEntry| $ 50))) 
diff --git a/src/algebra/strap/SINT.lsp b/src/algebra/strap/SINT.lsp
index 3dc9ea4f..c89f6600 100644
--- a/src/algebra/strap/SINT.lsp
+++ b/src/algebra/strap/SINT.lsp
@@ -80,13 +80,13 @@
 
 (DEFUN |SINT;base;$;12| ($) 2) 
 
-(PUT '|SINT;max;$;13| '|SPADreplace| '(XLAM NIL MOST-POSITIVE-FIXNUM)) 
+(PUT '|SINT;max;$;13| '|SPADreplace| '(XLAM NIL |$ShortMaximum|)) 
 
-(DEFUN |SINT;max;$;13| ($) MOST-POSITIVE-FIXNUM) 
+(DEFUN |SINT;max;$;13| ($) |$ShortMaximum|) 
 
-(PUT '|SINT;min;$;14| '|SPADreplace| '(XLAM NIL MOST-NEGATIVE-FIXNUM)) 
+(PUT '|SINT;min;$;14| '|SPADreplace| '(XLAM NIL |$ShortMinimum|)) 
 
-(DEFUN |SINT;min;$;14| ($) MOST-NEGATIVE-FIXNUM) 
+(DEFUN |SINT;min;$;14| ($) |$ShortMinimum|) 
 
 (PUT '|SINT;=;2$B;15| '|SPADreplace| 'EQL) 
 
@@ -240,8 +240,8 @@
 
 (DEFUN |SINT;coerce;I$;51| (|x| $)
   (SEQ (COND
-         ((NULL (< MOST-POSITIVE-FIXNUM |x|))
-          (COND ((NULL (< |x| MOST-NEGATIVE-FIXNUM)) (EXIT |x|)))))
+         ((NULL (< |$ShortMaximum| |x|))
+          (COND ((NULL (< |x| |$ShortMinimum|)) (EXIT |x|)))))
        (EXIT (|error| "integer too large to represent in a machine word")))) 
 
 (DEFUN |SINT;random;$;52| ($)
diff --git a/src/interp/Makefile.in b/src/interp/Makefile.in
index 3563e592..4c3441a0 100644
--- a/src/interp/Makefile.in
+++ b/src/interp/Makefile.in
@@ -425,15 +425,15 @@ buildom.$(FASLEXT): buildom.boot sys-macros.$(FASLEXT)
 diagnostics.$(FASLEXT): sys-globals.$(FASLEXT) vmlisp.$(FASLEXT)
 sys-driver.$(FASLEXT): sys-driver.boot types.$(FASLEXT)
 sys-globals.$(FASLEXT): sys-constants.$(FASLEXT) hash.$(FASLEXT)
+sys-os.$(FASLEXT): sys-constants.$(FASLEXT) cfuns.$(FASLEXT) 
 sys-constants.$(FASLEXT): types.$(FASLEXT)
 hash.$(FASLEXT): types.$(FASLEXT)
 union.$(FASLEXT): vmlisp.$(FASLEXT)
 ggreater.$(FASLEXT): vmlisp.$(FASLEXT)
 axext_l.$(FASLEXT): foam_l.$(FASLEXT)
-foam_l.$(FASLEXT): vmlisp.$(FASLEXT)
+foam_l.$(FASLEXT): vmlisp.$(FASLEXT) sys-constants.$(FASLEXT)
 sys-utility.$(FASLEXT): vmlisp.$(FASLEXT) sys-os.$(FASLEXT)
 vmlisp.$(FASLEXT): types.$(FASLEXT) sys-globals.$(FASLEXT)
-sys-os.$(FASLEXT): types.$(FASLEXT) cfuns.$(FASLEXT) 
 types.$(FASLEXT): boot-pkg.$(FASLEXT)
 boot-pkg.$(FASLEXT): boot-pkg.lisp
 
diff --git a/src/interp/Makefile.pamphlet b/src/interp/Makefile.pamphlet
index 7f159117..f46eab6c 100644
--- a/src/interp/Makefile.pamphlet
+++ b/src/interp/Makefile.pamphlet
@@ -674,15 +674,15 @@ buildom.$(FASLEXT): buildom.boot sys-macros.$(FASLEXT)
 diagnostics.$(FASLEXT): sys-globals.$(FASLEXT) vmlisp.$(FASLEXT)
 sys-driver.$(FASLEXT): sys-driver.boot types.$(FASLEXT)
 sys-globals.$(FASLEXT): sys-constants.$(FASLEXT) hash.$(FASLEXT)
+sys-os.$(FASLEXT): sys-constants.$(FASLEXT) cfuns.$(FASLEXT) 
 sys-constants.$(FASLEXT): types.$(FASLEXT)
 hash.$(FASLEXT): types.$(FASLEXT)
 union.$(FASLEXT): vmlisp.$(FASLEXT)
 ggreater.$(FASLEXT): vmlisp.$(FASLEXT)
 axext_l.$(FASLEXT): foam_l.$(FASLEXT)
-foam_l.$(FASLEXT): vmlisp.$(FASLEXT)
+foam_l.$(FASLEXT): vmlisp.$(FASLEXT) sys-constants.$(FASLEXT)
 sys-utility.$(FASLEXT): vmlisp.$(FASLEXT) sys-os.$(FASLEXT)
 vmlisp.$(FASLEXT): types.$(FASLEXT) sys-globals.$(FASLEXT)
-sys-os.$(FASLEXT): types.$(FASLEXT) cfuns.$(FASLEXT) 
 types.$(FASLEXT): boot-pkg.$(FASLEXT)
 boot-pkg.$(FASLEXT): boot-pkg.lisp
 
diff --git a/src/interp/foam_l.lisp b/src/interp/foam_l.lisp
index c948b274..8863d2d6 100644
--- a/src/interp/foam_l.lisp
+++ b/src/interp/foam_l.lisp
@@ -75,6 +75,7 @@
 #+:gcl (in-package "BOOT")
 (in-package "AxiomCore")
 (import-module "vmlisp")
+(import-module "sys-constants")
 
 (in-package "FOAM")
 
@@ -184,24 +185,18 @@
 
 
 ;; type defs for Foam types
-(deftype |Char| () 'character)
+(deftype |Char| () 'BOOT::|%Char|)
 (deftype |Clos| () 'list)
 (deftype |Bool| () '(member t nil))
-(deftype |Byte| () 'unsigned-byte)
+(deftype |Byte| () 'BOOT::|%Byte|)
 (deftype |HInt| () '(integer #.(- (expt 2 15)) #.(1- (expt 2 15))))
-(deftype |SInt| () 'fixnum)
+(deftype |SInt| () 'BOOT::|%Short|)
 
-#+:AKCL
-(deftype |BInt| () t)
-#-:AKCL
-(deftype |BInt| () 'integer)
+(deftype |BInt| () 'BOOT::|%Integer|)
 
-(deftype |SFlo| () 'short-float)
+(deftype |SFlo| () 'BOOT::|%SingleFloat|)
 
-#+:AKCL
-(deftype |DFlo| () t)
-#-:AKCL
-(deftype |DFlo| () 'long-float)
+(deftype |DFlo| () 'BOOT::|%DoubleFloat|)
 
 (deftype |Level| () t) ;; structure??
 
@@ -221,9 +216,8 @@
 (defconstant |HIntInit| (the |HInt| 0))
 (defconstant |SIntInit| (the |SInt| 0))
 (defconstant |BIntInit| (the |BInt|  0))
-(defconstant |SFloInit| (the |SFlo| 0.0s0))
-;; FIXME: Revisit the definition of DFlo as long-double.
-(defconstant |DFloInit| (the |DFlo| 0.0l0))
+(defconstant |SFloInit| (coerce 0 '|SFlo|))
+(defconstant |DFloInit| (coerce 0 '|DFlo|))
 (defconstant |PtrInit|  (the |Ptr|  nil))
 (defconstant |ArrInit|  (the |Arr|  nil))
 (defconstant |RecordInit|  (the |Record|  nil))
@@ -261,11 +255,11 @@
 (defmacro |CharOrd|     (x) `(CHAR-INT (the |Char| ,x)))
 (defmacro |CharNum|     (x) `(INT-CHAR (the |SInt| ,x)))
 
-(defmacro |SFlo0|        () 0.0s0)
-(defmacro |SFlo1|        () 1.0s0)
-(defmacro |SFloMin|      () most-negative-short-float)
-(defmacro |SFloMax|      () most-positive-short-float)
-(defmacro |SFloEpsilon|  () short-float-epsilon)
+(defmacro |SFlo0|        () (coerce 0 '|SFlo|))
+(defmacro |SFlo1|        () (coerce 1 '|SFlo|))
+(defmacro |SFloMin|      () BOOT::|$SingleFloatMinimum|)
+(defmacro |SFloMax|      () BOOT::|$SingleFloatMaximum|)
+(defmacro |SFloEpsilon|  () BOOT::|$SingleFloatEpsilon|)
 (defmacro |SFloIsZero|  (x) `(zerop (the |SFlo| ,x)))
 (defmacro |SFloIsNeg|   (x) `(minusp (the |SFlo| ,x)))
 (defmacro |SFloIsPos|   (x) `(plusp (the |SFlo| ,x)))
@@ -274,8 +268,8 @@
 (defmacro |SFloEQ|    (x y) `(= (the |SFlo| ,x) (the |SFlo| ,y)))
 (defmacro |SFloNE|    (x y) `(/= (the |SFlo| ,x) (the |SFlo| ,y)))
 (defmacro |SFloNegate|  (x) `(the |SFlo| (- (the |SFlo| ,x))))
-(defmacro |SFloNext|    (x) `(the |SFlo| (+ (the |SFlo| ,x) 1.0s0)))
-(defmacro |SFloPrev|    (x) `(the |SFlo| (- (the |SFlo| ,x) 1.0s0)))
+(defmacro |SFloNext|    (x) `(the |SFlo| (+ (the |SFlo| ,x) |SFlo1|)))
+(defmacro |SFloPrev|    (x) `(the |SFlo| (- (the |SFlo| ,x) |SFlo1|)))
 (defmacro |SFloMinus| (x y) `(the |SFlo| (- (the |SFlo| ,x) (the |SFlo| ,y))))
 (defmacro |SFloTimes| (x y) `(the |SFlo| (* (the |SFlo| ,x) (the |SFlo| ,y))))
 (defmacro |SFloTimesPlus| (x y z)
@@ -295,11 +289,11 @@
 ;;(defmacro |SFloFloor|    (x) `(the |BInt| (floor (the |SFlo| ,x))))
 ;;(defmacro |SFloCeiling|  (x) `(the |BInt| (ceiling (the |SFlo| ,x))))
 
-(defmacro |DFlo0|         () 0.0d0)
-(defmacro |DFlo1|         () 1.0d0)
-(defmacro |DFloMin|       () most-negative-long-float)
-(defmacro |DFloMax|       () most-positive-long-float)
-(defmacro |DFloEpsilon|   () long-float-epsilon)
+(defmacro |DFlo0|         () (coerce 0 '|DFlo|))
+(defmacro |DFlo1|         () (coerce 1 '|DFlo|))
+(defmacro |DFloMin|       () BOOT::|$DoubleFloatMinimum|)
+(defmacro |DFloMax|       () BOOT::|$DoubleFloatMaximum|)
+(defmacro |DFloEpsilon|   () BOOT::|$DoubleFloatEpsilon|)
 (defmacro |DFloIsZero|   (x) `(zerop (the |DFlo| ,x)))
 (defmacro |DFloIsNeg|    (x) `(minusp (the |DFlo| ,x)))
 (defmacro |DFloIsPos|    (x) `(plusp (the |DFlo| ,x)))
@@ -308,8 +302,8 @@
 (defmacro |DFloLT|     (x y) `(< (the |DFlo| ,x) (the |DFlo| ,y)))
 (defmacro |DFloNE|     (x y) `(/= (the |DFlo| ,x) (the |DFlo| ,y)))
 (defmacro |DFloNegate|   (x) `(the |DFlo| (- (the |DFlo| ,x))))
-(defmacro |DFloNext|     (x) `(the |DFlo| (+ (the |DFlo| ,x) 1.0d0)))
-(defmacro |DFloPrev|     (x) `(the |DFlo| (- (the |DFlo| ,x) 1.0d0)))
+(defmacro |DFloNext|     (x) `(the |DFlo| (+ (the |DFlo| ,x) |DFlo1|)))
+(defmacro |DFloPrev|     (x) `(the |DFlo| (- (the |DFlo| ,x) |DFlo1|)))
 (defmacro |DFloPlus|   (x y) `(the |DFlo| (+ (the |DFlo| ,x) (the |DFlo| ,y))))
 (defmacro |DFloMinus|  (x y) `(the |DFlo| (- (the |DFlo| ,x) (the |DFlo| ,y))))
 (defmacro |DFloTimes|  (x y) `(the |DFlo| (* (the |DFlo| ,x) (the |DFlo| ,y))))
@@ -344,8 +338,8 @@
 
 (defmacro |SInt0|         () 0)
 (defmacro |SInt1|         () 1)
-(defmacro |SIntMin|       () `(the |SInt| most-negative-fixnum))
-(defmacro |SIntMax|       () `(the |SInt| most-positive-fixnum))
+(defmacro |SIntMin|       () `(the |SInt| BOOT::|$ShortMinimum|))
+(defmacro |SIntMax|       () `(the |SInt| BOOT::|$ShortMaximum|))
 (defmacro |SIntIsZero|   (x) `(zerop (the |SInt| ,x)))
 (defmacro |SIntIsNeg|    (x) `(minusp (the |SInt| ,x)))
 (defmacro |SIntIsPos|    (x) `(plusp (the |SInt| ,x)))
diff --git a/src/interp/fortcall.boot b/src/interp/fortcall.boot
index 9d7d0d50..a8322d8f 100644
--- a/src/interp/fortcall.boot
+++ b/src/interp/fortcall.boot
@@ -399,11 +399,11 @@ vec2Lists u == [vec2Lists1 ELT(u,i) for i in 0..#u-1]
 spad2lisp(u) ==
   -- Turn complexes into arrays of floats
   first first(u)="Complex" =>
-    makeVector([makeVector([CADR u,CDDR u],'DOUBLE_-FLOAT)],NIL)
+    makeVector([makeVector([CADR u,CDDR u],"%DoubleFloat")],NIL)
   -- Turn arrays of complexes into arrays of floats so that tarnsposing
   -- them puts them in the correct fortran order
   first first(u)="Matrix" and first SECOND first(u) = "Complex" =>
-    makeVector([makeVector(complexRows vec2Lists rest u,'DOUBLE_-FLOAT)],NIL)
+    makeVector([makeVector(complexRows vec2Lists rest u,"%DoubleFloat")],NIL)
   rest(u)
 
 invokeFortran(objFile,args,dummies,decls,results,actual) ==
@@ -526,14 +526,14 @@ spadify(l,results,decls,names,actual) ==
 lispType u ==
   -- Return the lisp type equivalent to the given Fortran type.
   LISTP u => lispType first u
-  u = "real" => "SHORT-FLOAT"
-  u = "double" => "DOUBLE-FLOAT"
-  u = "double precision" => "DOUBLE-FLOAT"
+  u = "real" => "%SingleFloat"
+  u = "double" => "%DoubleFloat"
+  u = "double precision" => "DoubleFloat"
   u = "integer" => "FIXNUM"
   u = "logical" => "BOOLEAN"
   u = "character" => "CHARACTER"
-  u = "complex" => "SHORT-FLOAT"
-  u = "double complex" => "DOUBLE-FLOAT"
+  u = "complex" => "%SingleFloat"
+  u = "double complex" => "%DoubleFloat"
   error ['"Unrecognised Fortran type: ",u]
 
 getVal(u,names,values) ==
@@ -552,11 +552,11 @@ prepareData(args,dummies,values,decls) ==
 
 
 checkForBoolean u ==
-  u = "BOOLEAN" => "FIXNUM"
+  u = "BOOLEAN" => "%Short"
   u
 
-shortZero == COERCE(0.0,'SHORT_-FLOAT)
-longZero == COERCE(0.0,'DOUBLE_-FLOAT)
+shortZero == COERCE(0,"%SingleFloat")
+longZero == COERCE(0,"%DoubleFloat")
 
 prepareResults(results,args,dummies,values,decls) ==
   -- Create the floating point zeros (boot doesn't like 0.0d0, 0.0D0 etc)
@@ -570,7 +570,7 @@ prepareResults(results,args,dummies,values,decls) ==
           makeVector(  makeList(
             2*APPLY('_*,[getVal(tt,argNames,actual) for tt in rest(type)]),_
             if first(type)="complex" then shortZero else longZero),_
-          if first(type)="complex" then "SHORT-FLOAT" else "DOUBLE-FLOAT" )
+          if first(type)="complex" then "%SingleFloat" else "%DoubleFloat" )
         LISTP type => makeVector(_
           makeList(
             APPLY('_*,[getVal(tt,argNames,actual) for tt in rest(type)]),_
@@ -582,8 +582,8 @@ prepareResults(results,args,dummies,values,decls) ==
         type = "double precision" => longZero
         type = "logical" => 0
         type = "character" => MAKE_-STRING(1)
-        type = "complex" => makeVector([shortZero,shortZero],'SHORT_-FLOAT)
-        type = "double complex" => makeVector([longZero,longZero],'LONG_-FLOAT)
+        type = "complex" => makeVector([shortZero,shortZero],"%SingleFloat")
+        type = "double complex" => makeVector([longZero,longZero],"%DoubleFloat")
         error ['"Unrecognised Fortran type: ",type]
   NREVERSE data
 
diff --git a/src/interp/i-output.boot b/src/interp/i-output.boot
index a3151587..a9a26e81 100644
--- a/src/interp/i-output.boot
+++ b/src/interp/i-output.boot
@@ -876,7 +876,7 @@ WIDTH u ==
       negative := 0
     -- Try and be fairly exact for smallish integers:
     u = 0 => 1
-    u < MOST_-POSITIVE_-LONG_-FLOAT => 1+negative+FLOOR ((log10 u) + 0.0000001)
+    u < $DoubleFloatMaximum => 1+negative+FLOOR ((log10 u) + 0.0000001)
     -- Rough guess: integer-length returns log2 rounded up, so divide it by
     -- roughly log2(10). This should return an over-estimate, but for objects
     -- this big does it matter?
diff --git a/src/interp/sfsfun.boot b/src/interp/sfsfun.boot
index 8118fb4b..a740f3d5 100644
--- a/src/interp/sfsfun.boot
+++ b/src/interp/sfsfun.boot
@@ -370,7 +370,7 @@ rPsiW(n,x) ==
         xmin := float(FLOOR(alpha + beta*n) + 1)
         if n>0
         then
-                a := MIN(0,1.0/float(n)*LOG(DOUBLE_-FLOAT_-EPSILON/MIN(1.0,x)))
+                a := MIN(0,1.0/float(n)*LOG($DoubleFloatPrecision/MIN(1.0,x)))
                 c := EXP(a)
                 if ABS(a) >= 0.001
                 then
diff --git a/src/interp/sys-constants.boot b/src/interp/sys-constants.boot
index 12b11104..4bf36c9e 100644
--- a/src/interp/sys-constants.boot
+++ b/src/interp/sys-constants.boot
@@ -40,6 +40,55 @@
 import types
 namespace BOOT
 
+--%
+--% Numeric limits
+--%
+
+++ Minimum for %Short values.
+$ShortMinimum ==
+  MOST_-NEGATIVE_-FIXNUM
+
+++ Maximum for %Short values.
+$ShortMaximum ==
+  MOST_-POSITIVE_-FIXNUM
+
+++ Minimum for %SingleFloat values.
+$SingleFloatMinimum ==
+)if %hasFeature KEYWORD::GCL
+  MOST_-NEGATIVE_-SHORT_-FLOAT
+)else
+  MOST_-NEGATIVE_-SINGLE_-FLOAT
+)endif
+
+++ Maximum for %SingleFloat values.
+$SingleFloatMaximum ==
+)if %hasFeature KEYWORD::GCL
+  MOST_-POSITIVE_-SHORT_-FLOAT
+)else
+  MOST_-POSITIVE_-SINGLE_-FLOAT
+)endif
+
+++ Machine precision for %SingleFloat
+$SingleFloatEpsilon ==
+)if %hasFeature KEYWORD::GCL
+  SHORT_-FLOAT_-EPSILON
+)else
+  SINGLE_-FLOAT_-EPSILON
+)endif
+
+++ Maximum for %DoubleFloat values
+$DoubleFloatMinimum ==
+  MOST_-NEGATIVE_-DOUBLE_-FLOAT
+
+++ Maximum for %DoubleFloat values
+$DoubleFloatMaximum ==
+  MOST_-POSITIVE_-DOUBLE_-FLOAT
+
+++ Machine precision for %DoubleFloat
+$DoubleFloatEpsilon ==
+  DOUBLE_-FLOAT_-EPSILON
+
+--% 
 
 ++ Clock time unit per second.
 $timerTicksPerSecond == 
diff --git a/src/interp/sys-os.boot b/src/interp/sys-os.boot
index 9f22adf0..95342e27 100644
--- a/src/interp/sys-os.boot
+++ b/src/interp/sys-os.boot
@@ -38,7 +38,7 @@
 -- supporting C runtime libopen-axiom-core.
 --
 
-import types
+import sys_-constants
 import cfuns
 namespace BOOT
 module sys_-os
@@ -181,7 +181,7 @@ $minusInfinity == SB_-EXT::DOUBLE_-FLOAT_-NEGATIVE_-INFINITY
 )else
 -- In general Common Lisp does not provide support for infinities
 -- and the like.
-$plusInfinity == MOST_-POSITIVE_-DOUBLE_-FLOAT
+$plusInfinity == $DoubleFloatMaximum
 
 $minusInfinity == -$plusInfinity
 )endif
diff --git a/src/interp/types.boot b/src/interp/types.boot
index 748e5070..6bb4a555 100644
--- a/src/interp/types.boot
+++ b/src/interp/types.boot
@@ -64,12 +64,16 @@ namespace BOOT
   INTEGER
 
 ++ Type of single precision floating point numbers.  Most of the
-++ time, this is a 32-bit datatype.
+++ time, this is a 32-bit datatype on IEEE-754 host.
 %SingleFloat <=>
+)if %hasFeature KEYWORD::GCL
+  SHORT_-FLOAT
+)else
   SINGLE_-FLOAT
+)endif
 
 ++ Type of double precision floating point numbers.  Most of the time,
-++ this is a 64-bit sized datatype.
+++ this is a 64-bit sized datatype on IEEE-756 host.
 %DoubleFloat <=>
   DOUBLE_-FLOAT
author	dos-reis <gdr@axiomatics.org>	2008-09-07 05:41:35 +0000
committer	dos-reis <gdr@axiomatics.org>	2008-09-07 05:41:35 +0000
commit	18b7b4b5d3511211092796bd8d593d5c9debec02 (patch)
tree	a90da73966bbb0311bf51c62e9ab8356465be694
parent	52642d936f26a18f7d3818fcabe3fad2bff157ea (diff)
download	open-axiom-18b7b4b5d3511211092796bd8d593d5c9debec02.tar.gz