[Openmcl-devel] representing infinity
Gary King
gwking at cs.umass.edu
Wed Mar 24 11:41:20 PST 2004
Hi Gary,
To follow up a bit on this, please look at the follow OpenMCL session.
I can create infinities one step at a time but not "all at once" in an
unwind protect. I'm afraid this goes beyond my debugging skills but I'm
hoping to learn more from your answer. FWIW, the unwind protect does
work in MCL 5.0
Welcome to OpenMCL Version (Beta: Darwin) 0.14.1-p1!
? (ccl:get-fpu-mode)
(:ROUNDING-MODE :NEAREST :OVERFLOW T :UNDERFLOW NIL :DIVISION-BY-ZERO T
:INVALID T :INEXACT NIL)
? (ccl:set-fpu-mode :division-by-zero nil)
(:ROUNDING-MODE :NEAREST :OVERFLOW T :UNDERFLOW NIL :DIVISION-BY-ZERO
NIL :INVALID T :INEXACT NIL)
? (setf positive-infinity (/ 0d0))
#<INFINITY 1.797693134862316D+308>
? (ccl:set-fpu-mode :division-by-zero t)
(:ROUNDING-MODE :NEAREST :OVERFLOW T :UNDERFLOW NIL :DIVISION-BY-ZERO T
:INVALID T :INEXACT NIL)
? positive-infinity
#<INFINITY 1.797693134862316D+308>
? (setf another-positive-infinity
(unwind-protect
(progn
(ccl:set-fpu-mode :division-by-zero nil)
(/ 0d0))
(ccl:set-fpu-mode :division-by-zero t))))
Unhandled exception 11 at 0x00004774, context->regs at #xf0134ad8
Read operation to unmapped address 0x004ba000
? for help
[5351] OpenMCL kernel debugger: ?
(S) Set specified GPR to new value
(A) Advance the program counter by one instruction (use with caution!)
(D) Describe the current exception in greater detail
(R) Show raw GPR/SPR register values
(L) Show Lisp values of tagged registers
(F) Show FPU registers
(B) Show backtrace
(X) Exit from this debugger, asserting that any exception was handled
(K) Kill OpenMCL process
(?) Show this help
[5351] OpenMCL kernel debugger: D
Read operation to unmapped address 0x004ba000
[5351] OpenMCL kernel debugger: R
r00 = 0x00000000 r08 = 0x00002015 r16 = 0x00002015 r24 = 0x05142A4E
r01 = 0xF0134FC0 r09 = 0x05396580 r17 = 0x050124EE r25 = 0x0000006C
r02 = 0x00101F50 r10 = 0x05390000 r18 = 0x050124DE r26 = 0x051429C6
r03 = 0x004B9FFC r11 = 0x051429C6 r19 = 0x00000033 r27 = 0x00000088
r04 = 0x00427930 r12 = 0x004B8730 r20 = 0x00000000 r28 = 0x00002015
r05 = 0x05141102 r13 = 0x00426070 r21 = 0x05399736 r29 = 0x7FFFFFFC
r06 = 0x004B87BE r14 = 0x0110C624 r22 = 0x000000E4 r30 = 0x0539A2FE
r07 = 0x00000000 r15 = 0x0510BD9E r23 = 0x00000340 r31 = 0x053983DE
PC = 0x00004774 LR = 0x0110C624 CTR = 0x01014B56 CCR = 0x42842482
XER = 0x00000000 MSR = 0x0000F930 DAR = 0x004BA000 DSISR = 0x40000000
[5351] OpenMCL kernel debugger: L
rcontext = 0x00101F50
nargs = 21301873
r15 (fn) = #<Function PRINT-A-NAN #x0510bd9e>
r23 (arg_z) = 208
r22 (arg_y) = 57
r21 (arg_x) = "performing CCL::UNKNOWN on (#<NaN
-2.696539702543562D+308>
r20 (temp0) = 0
r19 (temp1/next_method_context) = #<Unbound>
r18 (temp2/nfn) = #<Function %SET-FPSCR-CONTROL #x050124de>
r17 (temp3/fname) = %SET-FPSCR-CONTROL
r16 (temp4) = ()
r31 (save0) = #<3-element vector subtag = 72 @#x053983de>
r30 (save1) =
r29 (save2) = 536870911
r28 (save3) = ()
r27 (save4) = 34
r26 (save5) = #<2-element vector subtag = 9A @#x051429c6>
r25 (save6) = 27
r24 (save7) = FORMAT-ERROR
[5351] OpenMCL kernel debugger: F
f00 : 0xFFF80000000000D0 (nan)
f01 : 0xBFF0000000000000 (-1.000000)
f02 : 0x41D6403067400000 (1493221789.000000)
f03 : 0x412E848000000000 (1000000.000000)
f04 : 0x4173EA3800000000 (20882304.000000)
f05 : 0x3FE40C7A3A33C850 (0.626523)
f06 : 0xFFF8000000098F5B (nan)
f07 : 0xC24BC19587859393 (-238423838475.152924)
f08 : 0xC24BC19587859393 (-238423838475.152924)
f09 : 0xC24BC19587859393 (-238423838475.152924)
f10 : 0xC24BC19587859393 (-238423838475.152924)
f11 : 0xC24BC19587859393 (-238423838475.152924)
f12 : 0xC24BC19587859393 (-238423838475.152924)
f13 : 0xC24BC19587859393 (-238423838475.152924)
f14 : 0xC24BC19587859393 (-238423838475.152924)
f15 : 0xC24BC19587859393 (-238423838475.152924)
f16 : 0xC24BC19587859393 (-238423838475.152924)
f17 : 0xC24BC19587859393 (-238423838475.152924)
f18 : 0xC24BC19587859393 (-238423838475.152924)
f19 : 0xC24BC19587859393 (-238423838475.152924)
f20 : 0xC24BC19587859393 (-238423838475.152924)
f21 : 0xC24BC19587859393 (-238423838475.152924)
f22 : 0xC24BC19587859393 (-238423838475.152924)
f23 : 0xC24BC19587859393 (-238423838475.152924)
f24 : 0xC24BC19587859393 (-238423838475.152924)
f25 : 0xC24BC19587859393 (-238423838475.152924)
f26 : 0xC24BC19587859393 (-238423838475.152924)
f27 : 0xC24BC19587859393 (-238423838475.152924)
f28 : 0xC24BC19587859393 (-238423838475.152924)
f29 : 0xC24BC19587859393 (-238423838475.152924)
f30 : 0x4330000080000000 (4503601774854144.000000)
f31 : 0x0000000000000000 (0.000000)
FPSCR = 000000D0
[5351] OpenMCL kernel debugger: K
[Polanyi:~/CVS/gbbopen/gbbopen] gwking%
On Mar 19, 2004, at 8:28 PM, Gary Byers wrote:
>
>
> On Fri, 19 Mar 2004, Gary King wrote:
>
>> Is there a way in OpenMCL to represent IEEE positive and negative
>> infinity? An apropos on NAN gives me (among other things):
>>
>> CCL::IEEE-DOUBLE-FLOAT-SIGNALLING-NAN-BIT, Value: 51
>> CCL::IEEE-SINGLE-FLOAT-SIGNALLING-NAN-BIT, Value: 22
>> CCL::NAN-OR-INFINITY-P, Def: FUNCTION
>> CCL::NANNING
>> CCL::PRINT-A-NAN, Def: FUNCTION
>>
>> and one on "INFINITY" gives:
>>
>> CCL::INFINITY-P, Def: FUNCTION
>>
>> This is promising but even after looking at the source code it's not
>> clear to me how to produce something such that (infinity-p
>> <something>)
>> returns t.
>>
>> Thanks,
>> --
>> Gary Warren King, Lab Manager
>> EKSL East, University of Massachusetts * 413 577 0176
>>
>
> It's easy to "represent" an infinity or a NaN, since the hardwre does
> all of the hard work ...
>
> It's a little harder to produce one: FP operations whose results can't
> be represented (either because they're mathematically undefined or
> because their magnitude is too great) will generally either:
>
> a) raise an exception if that exception is enabled in the FPU in the
> current thread
> b) produce a NaN or an infinity otherwise.
>
> When a lisp thread starts up, IEEE OVERFLOW, DIVISION-BY-ZERO, and
> INVALID-OPERATION exceptions are enabled:
>
> ? (ccl:get-fpu-mode)
> (:ROUNDING-MODE :NEAREST :OVERFLOW T :UNDERFLOW NIL :DIVISION-BY-ZERO
> T :INVALID T :INEXACT NIL)
>
> so trying to do something that'd cause an overflow will generate a
> hardware
> exception (which is mapped to a lisp condition):
>
> ? (* 2 most-positive-single-float)
>> Error in process listener(1): FLOATING-POINT-OVERFLOW detected
>> performing * on (3.4028235E+38
>> 1.0842021E-19)
>> While executing: %SHORT-FLOAT*-2!
>> Type :POP to abort.
> Type :? for other options.
> 1 >
>
> [There's a bug there, in that the second operand is some more-or-less
> random
> number instead of 2.0.]
>
> If we disable the overflow exception, the result will be an infinity:
> ? (ccl:set-fpu-mode :overflow nil)
>
> then:
>
> ? (* 2 most-positive-single-float)
> #<INFINITY >
>
>
--
Gary Warren King, Lab Manager
EKSL East, University of Massachusetts * 413 577 0176
I believe that we can effectively defend ourselves abroad and at home
without dimming our principles. Indeed, I believe that our success in
defending ourselves depends precisely on not giving up what we stand
for.
-- Al Gore
More information about the Openmcl-devel
mailing list