Use the two-prod algorithm from crlibm documentation.
This handles overflows better and (c::two-prod 1.7976931214684583d308 (1+ (scale-float 1d0 -28))) doesn't signal an overflow like the old algorithm.
Loading
Please register or sign in to comment