GHCiが以下の不正解を返すのはなぜですか?
GHCi
λ> ((-20.24373193905347)^12)^2 - ((-20.24373193905347)^24)
4.503599627370496e15
Python3
>>> ((-20.24373193905347)**12)**2 - ((-20.24373193905347)**24)
0.0
更新 Haskellの(^)関数を次のように実装します。
powerXY :: Double -> Int -> Double
powerXY x 0 = 1
powerXY x y
| y < 0 = powerXY (1/x) (-y)
| otherwise =
let z = powerXY x (y `div` 2)
in if odd y then z*z*x else z*z
main = do
let x = -20.24373193905347
print $ powerXY (powerXY x 12) 2 - powerXY x 24 -- 0
print $ ((x^12)^2) - (x ^ 24) -- 4.503599627370496e15
私のバージョンは@WillemVanOnsemによって以下に提供されているバージョンよりも正確に見えませんが、少なくともこの特定のケースでは奇妙な正解を示しています。
Pythonも同様です。
def pw(x, y):
if y < 0:
return pw(1/x, -y)
if y == 0:
return 1
z = pw(x, y//2)
if y % 2 == 1:
return z*z*x
else:
return z*z
# prints 0.0
print(pw(pw(-20.24373193905347, 12), 2) - pw(-20.24373193905347, 24))
2.243746917640863e31 - 2.2437469176408626e31
これは、増幅される小さな丸め誤差を計算します。キャンセルの問題のようです。
a^24
は約2.2437e31
であるため、これを生成する丸め誤差があります。