¿Por qué GHCi da una respuesta incorrecta a continuación?
GHCi
λ> ((-20.24373193905347)^12)^2 - ((-20.24373193905347)^24)
4.503599627370496e15
Python3
>>> ((-20.24373193905347)**12)**2 - ((-20.24373193905347)**24)
0.0
ACTUALIZACIÓN Implementaría la función de Haskell (^) de la siguiente manera.
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
Aunque mi versión no parece más correcta que la proporcionada a continuación por @WillemVanOnsem, extrañamente da la respuesta correcta para este caso en particular al menos.
Python es similar.
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
que tiene un pequeño error de redondeo que se amplifica. Parece un problema de cancelación.
a^24
es aproximadamente2.2437e31
, y por lo tanto hay un error de redondeo que produce esto.