1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
/* some tests for non-numeric floats
* these are only meaningful for Lisps which support such values
* GCL, ECL, SBCL, maybe others;
*
* SBCL: allows both inf/-inf and NaN; enable via:
* :lisp (sb-vm::set-floating-point-modes :traps nil)
*
* GCL: allows both inf/-inf and NaN;
* requires *PRINT-NAN* = T but that's already enabled by Maxima
*
* ECL: allows only inf/-inf, not NaN
* (although maybe in a future version)
* no special action needed to enable, as far as I know
*/
kill (all);
done;
(is_sbcl : is (ssearch ("SBCL", build_info()@lisp_name) # false),
is_gcl : is (ssearch ("GCL", build_info()@lisp_name) # false),
is_ecl : is (ssearch ("ECL", build_info()@lisp_name) # false),
is (is_sbcl or is_gcl or is_ecl));
true;
(x:1e300, i:x*x, if not is_ecl then n:i-i, minusi:-i, 0);
0;
/* Maxima "=" calls ALIKE1 which compares atoms with EQUAL.
* As long as that's the case, we expect is(n=n) => true.
*/
[is(x=x), is(i=i), if is_ecl then [] else is(n=n), is(minusi=minusi)];
[true, true, true, true];
/* Lisp "=" yields false for NaN. */
[?\=(x, x), ?\=(i, i), if is_ecl then [] else ?\=(n, n), ?\=(minusi, minusi)];
[true, true, false, true];
/* bfloat and rat should trigger an error on non-numeric floats */
errcatch(bfloat(i));
[];
if is_ecl then false else errcatch(bfloat(n));
[];
errcatch(bfloat(minusi));
[];
/* GCL: rat(nonnumeric) => hangs; SF bug report # !! */
if not is_gcl then errcatch(rat(i));
[];
if is_ecl then false else if not is_gcl then errcatch(rat(n));
[];
if not is_gcl then errcatch(rat(minusi));
[];
/* Fortran output should produce INF/-INF or NAN */
map (fortran, ' ['x = x, 'i = i, 'n = n, 'minusi = minusi]);
[done, done, done, done];
/* At present, Maxima will output implementation-dependent
* representations of non-numeric floats.
*/
string('x = x);
''(if is_sbcl then "x = 1.0e+300"
elseif is_gcl then "x = 1.0E+300"
elseif is_ecl then "x = 1.0e+300");
string('i = i);
''(if is_sbcl then "i = #.SB-EXT:DOUBLE-FLOAT-POSITIVE-INFINITY"
elseif is_gcl then "i = #<inf>"
elseif is_ecl then "i = #.ext::double-float-positive-infinity");
string('n = n);
''(if is_sbcl then "n = #<DOUBLE-FLOAT quiet NaN>"
elseif is_gcl then "n = #<nan>");
string('minusi = minusi);
''(if is_sbcl then "minusi = #.SB-EXT:DOUBLE-FLOAT-NEGATIVE-INFINITY"
elseif is_gcl then "minusi = #<-inf>"
elseif is_ecl then "minusi = #.ext::double-float-negative-infinity");
/* SF bug report #2680: "1000!^0.01 gives i.nfE+6166368" */
errcatch (string (1000!^0.01));
''(if is_sbcl then ["#.SB-EXT:DOUBLE-FLOAT-POSITIVE-INFINITY"]
elseif is_gcl then ["#<inf>"]
elseif is_ecl then ["#.ext::double-float-positive-infinity"]);
/* SF bug report #2749: "Float evaluation of sinh causes Lisp error instead of overflowing to bigfloat" */
errcatch(string(sinh(1e3)));
''(if is_sbcl then ["#.SB-EXT:DOUBLE-FLOAT-POSITIVE-INFINITY"]
elseif is_gcl then ["#<inf>"]
elseif is_ecl then ["#.ext::double-float-positive-infinity"]);
errcatch(string(cosh(1e3)));
''(if is_sbcl then ["#.SB-EXT:DOUBLE-FLOAT-POSITIVE-INFINITY"]
elseif is_gcl then ["#<inf>"]
elseif is_ecl then ["#.ext::double-float-positive-infinity"]);
/* SF bug report #3081: abs(...) -> "bfloat: attempted conversion of floating-point infinity." */
[abs(%e^(10000/3)+1), abs(%e^10000+1), abs(%e^(10000/3))];
[abs(%e^(10000/3)+1), abs(%e^10000+1), abs(%e^(10000/3))];
errcatch(string(ev(exp(10000/3),numer)));
''(if is_sbcl then ["#.SB-EXT:DOUBLE-FLOAT-POSITIVE-INFINITY"]
elseif is_gcl then ["#<inf>"]
elseif is_ecl then ["#.ext::double-float-positive-infinity"]);
errcatch(string(ev(exp(10000/3),numer,%enumer)));
''(if is_sbcl then ["#.SB-EXT:DOUBLE-FLOAT-POSITIVE-INFINITY"]
elseif is_gcl then ["#<inf>"]
elseif is_ecl then ["#.ext::double-float-positive-infinity"]);
float_infinity_p (i);
true;
float_infinity_p (minusi);
true;
sign (i);
pos;
is (i > 0.0);
true;
sign (minusi);
neg;
is (minusi < 0.0);
true;
float_nan_p (n);
true;
map (float_infinity_p, [n, 0, 1.0, %pi, 'foo('blurf)]);
[false, false, false, false, false];
map (float_nan_p, [i, minusi, 0, 1.0, %pi, 'foo('blurf)]);
[false, false, false, false, false, false];
|
