Primena GCC-a u numerici

Priloženi program ne daje iste rezultate kada se kompajlira sa i bez optimizacije. Ispravan rezultat je sum=0 i on se mora dobiti shewchuk-ovim algoritmom. Međutim, "optimizacije" ga upropaste. Pokušao sam da kompajliram sa



ali neće. I dalje izbacuje sumu različitu od nule. Kada kompajliram sa



sve radi fino. Sa druge strane VC++ optimizuje kod, ali kako treba - dobija se rezultat 0. Kako ovo kompajlirati GCC-om sa optimizacijom da da se dobije ispravan rezultat?

Potpis

Nije bitno koji su zaključci izvučeni, već kako se do njih došlo.

Kod mene sa GCC-om 4.4.5 (tačnije: Debijanov paket 4.4.5-10, amd64), obe komandne linije koje si naveo rezultuju nultom sumom. Takođe izlazi nula i kada kompilujem čisto kao

Nenultu vrednost dobijem tek ako pokušam

Isto se ponaša i GCC 4.3.5 (Deb. 4.3.5-4).

Hvala. Ja imam verziju 4.5.0 (openSUSE 11.3) i radi kako sam naveo. Probaću noviju verziju GCC-a. Koliko vidim, izašla je verzija 4.5.2.

Potpis

Nije bitno koji su zaključci izvučeni, već kako se do njih došlo.

Upravo isprobah najnoviju verziju MinGW-a (4.5.2) pod Windows-om i ima opisanu grešku.

Potpis

Nije bitno koji su zaključci izvučeni, već kako se do njih došlo.

Onda će biti defekt u grani 4.5, koji bi mogao da prijaviš razvojnom timu. Ovde su potanko opisali kako se podnosi izveštaj o defektu: http://gcc.gnu.org/bugs/ . (Jedna stavka je i da se izgradi trenutna verzija GCC-a iz izvorne riznice i proveri na njoj, pošto ne prihvataju izveštaje na osnovu paketa treće strane; malo oštro, ali bih rekao da ima smisla za nešto poput kompilatora.)

Prijavljeno.

Potpis

Nije bitno koji su zaključci izvučeni, već kako se do njih došlo.

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.“