Vidim da je tvoj izveštaj preusmeren na drugi, i to
vrlo zanimljiv. Ja, dakle, problem nisam video jer sam isprobavao na amd64 procesoru, kod kojeg GCC automatski koristi SSE za operacije u pokretnom zarezu. Sad sam probao na ne-SSE ia32 procesoru, gde GCC spada na iks87, i tamo vidim problem; reši se, kao što ti reče R. Ginter, tek kad nametnem kopiranje u memoriju (
-ffloat-store) ili snizim tačnost (
-mpc32; ali ne i
-mpc64).
Kad bolje pogledam, i sam Šučak je u
svom radu upozorio da proširena tačnost iks87 FPU-ova dovodi do kraha algoritma: „
[...] the presence of extended precision internal floating-point registers, such as those on the Intel 80486 and Pentium processors [...] cause the methods described herein for determining the roundoff of an operation to fail.“