From: Eric M. <eri...@fr...> - 2011-10-21 21:03:09
|
Hi, ,---- | * (lisp-implementation-version) | "1.0.52.15-7cee246" | * (sqrt -4.9406564584124654d-320) | debugger invoked on a TYPE-ERROR in thread #<THREAD>: | The value 1061 is not of type (INTEGER -1074 1024). | 0: [ABORT] Exit debugger, returning to top level. | (SB-KERNEL::CSSQS #<unavailable argument>) `---- Running on Linux/AMD64. -- Eric Marsden |
From: Paul K. <pv...@pv...> - 2011-10-22 17:09:21
|
In article <87m...@fr...>, Eric Marsden <eri...@fr...> wrote: > Hi, > > ,---- > | * (lisp-implementation-version) > | "1.0.52.15-7cee246" > | * (sqrt -4.9406564584124654d-320) > | debugger invoked on a TYPE-ERROR in thread #<THREAD>: > | The value 1061 is not of type (INTEGER -1074 1024). > | 0: [ABORT] Exit debugger, returning to top level. > | (SB-KERNEL::CSSQS #<unavailable argument>) > `---- > > Running on Linux/AMD64. >From CSSQS ;; If we're here, neither x nor y are infinity and at ;; least one is non-zero.. Thus logb returns a nice ;; integer. (let ((k (- (logb-finite (max (abs x) (abs y)))))) (values (+ (square (scalb x k)) (square (scalb y k))) (- k)))) SCALB doesn't like denorms. Not sure what to do about this; I suppose we could clamp the scale factor? Paul Khuong |
From: Nikodemus S. <nik...@ra...> - 2011-10-26 11:02:16
|
On 22 October 2011 20:08, Paul Khuong <pv...@pv...> wrote: > SCALB doesn't like denorms. Not sure what to do about this; I suppose > we could clamp the scale factor? Can't we just loosen the exponent type in SCALB, or use SCALE-FLOAT directly? Cheers, -- nikodemus |