Dotaz: FPU - odečítání

[Jardík]

Mám takový problém s odečítáním v assembleru (pomocí FPU).

Mám takovýto stav st registrů:

st(0) = 1/5! = 8,3333e-3
st(1) = 0,25/7! = 4,96031e-5
st(2) = 0,25
st(3) = 1

Teď potřebuju provést odčítání st(0) = st(0) - st(1) a dostat

st(0) = 1/5! - 0,25/7! = 8,2837e-3
st(1) = 0,25
st(2) = 1

Provedl jsem tedy

fsubrp %st(0), %st(1)  ; fsubrp st(1), st(0) v intel syntax

Což by mělo udělat to, co potřebuji:

Subtract ST(i) from ST(0), store result in ST(i), and pop register stack

Jenže nějakým zázrakem dostanu v st(0) hodnotu -8,2837e-3, takže to udělalo pravý opak. Když instrukci změním na fsubp, tak dostanu správný výsledek, ale podle manuálu je to blbost (?). Je to bug v gasu nebo ve mně?

---

---

Odpovědi

[Andrej Herceg]

Pravdepodobne ide o tento problém 9.13.13 AT&T Syntax bugs (a teda gas má kvôli gcc vymenený zdrojový a cieľový register).

[Jardík]

Super ... a co teď s tím, přece nebudu kvůli chybnému gasu odevzdávat chybný úkol v assembleru, který má být zkompilovatelný gasem. Takže buď:

udělám správný kód a budu cvičícímu vysvětlovat, že gas je chybný

udělám špatný úkol, který gas zkompiluje jako správný, ale cvičící mě pošle někam, protože je úkol špatně

vykašlu se na gas, použiju yasm a budu přesvědčovat cvičícího, že gas je zabugovaný, není schopen správně přeložit kód

zadání

Napište podprogram pro výpočet funkce sinus. sin(x) = x/1! - x^3/3! + x^5/5! - x^7/7!

není jednoznačné a proto použiju instrukci fsin (v čemž mi zadání nebrání), která s gasem vrátí správný výsledek ... a hádat se s cvičícím, že měl zakázat použití instrukce fsin, ten mi to stejně neuzná, jak je na ČVUT zvykem.

Co byste udělal vy?

[Ladicek]

Použij druhou možnost, nad tu instrukci napiš komentář, a rozeber to v dokumentaci s odkazem na tu bugu. Ohavná realita sometimes happens.

[Jardík]

Asi jo, ale stejně mi to příde dost úchylný psát kvůli chybnému kompilátoru úmyslně chybný kód, aby byl výsledek nakonec správný.

[Petr Šobáň]

Asi by u té instrukce mělo být makro které bude testovat verzi překladače a podle toho tam vloží správnou instrukci aby to počítalo správně.

[Petr Šobáň]

Teď mě ještě napadlo není to test od učitele jak si poradíte s kompilátorem?

A s tím makrem by to mělo být čisté. (Teoreticky)

[l4m4]

Doufám, že si uvědomuješ, že je to celé 32bitová šmejďárna. Normálně se už dávno pro floating point aritmetiku používají SSE instrukce na XMM registrech, protože je to podstatně efektivnější než 387 instrukce.

[Jardík]

Samozřejmě bych radši použil SSE2, ale zadání je zadání.