Spatny vypocet v Cecku - datove typy?

[masterboy]

Zdravim Vas,

zacinam si hrat s programovani uP ARM Cortex M3, ale narazil jsem zrejme na muj zacatecnicky nedostatek / neznamost programovani v C. Pro jednoduchost sem vkladam jen samotny kod bez hlavickovych souboru:

Kód: Vybrat vše

int main(void) {

uint32_t x[2] = {0,237154};
uint32_t y[2] = {0,225296};
int16_t K,Q;

K = 1000 * (y[0] - y[1]) / (x[0] - x[1]);
Q = y[0] - (K / 1000) * x[0];

printf("%i\t%i\n",K,Q);

while(1);

}

Problem je v tom, ze mi to nepocita spravne a zrejme tusim, ze to bude asi zalezitost ohledne datovych typu, ale nejak tomu nerozumim ... Pomuze nekdo?

Jinak spravny vysledek: K = 949 a Q = 0
Je to vlastne pocitani parametru smernice primky (y = K.x + Q), zadane dvema body.

Kód: Vybrat vše

[Prochy]

Nejsem si jistý, nemůže být problém v tom zlomku? Záleží, co udělá procesor jako první, jestli to vydělí či vynásobí. Ani si nenapsal, kolik ti to vyhodilo za hodnoty.

Když první vydělí, tak to je: 225296/237154 a to odřízne a je z toho nula.

Když vynásobí, tak by to mělo vyjít, tak jak píšeš.

Zkusil bych to dát ještě do jedný závorky: (1000 * (y[0] - y[1])) / (x[0] - x[1]);

Tohle jediný mě tedy napadlo, ikdyž nemám až takové zkušenosti s programováním procesorů.

Jinak nebylo by lepší použít float? Směrnice přímky přece nemusí být vždy celé číslo.

[masterboy]

To co mi to vypocitalo bylo K = 0 a Q = 0, takze toto neni urcite spravny vysledek

Vyzkousel jsem take to co jsi navrhoval, ale vysledek je uplne stejny.
Co se tyka vyuziti toho FLOATu, tak to me napadlo jako prvni, ale je potreba si uvedomit ze uP nejsou zase tak vykonne a ze je pro vypocet rychlejsi pouzit cele cisla (proto to nasobeni 1000). Na vypocet ve FLOATu to sezere skoro 6x vice casu nez v celych cislech, i kdyz se dopoustim chyby zaokrouhlovani.

[Panda38]

Změň všechny typy na int32_t

a:

Q = y[0] - K*x[0] / 1000;

[masterboy]

To Panda38: diky za tip, s tim to opravdu funguje, ale proc?

PS.: Jaky je rozdil v mem a tvem zapisu vypocti hodnoty Q?

[Panda38]

uint32_t je číslo bez znaménka. y[0] - y[1] pak není -225296, ale 4294742000, to násobíš 1000 (přeteče mezivýsledek) a dělíš 4294730142, což při dělení dá jiný výsledek, který se navíc nevleze do int16_t.

U Tvého výpočtu Q když dělíš 949 / 1000, je výsledkem 0 (protože je to celočíselné dělení). Proto je třeba změnit pořadí, nejdříve K*x[0] a až potom / 1000.

[Andrea]

Problém je ten unsigned int.
K = 1000 * ((int)y[0] - (int)y[1]) / ((int)x[0] - (int)x[1]);
A tady je potřeba závorka.
Q = y[0] - (K * x[0]) / 1000;

Jsem to tu pro jistotu zkoušela, než to odešlu a tak pozdě. No máš to tu dvakrát.

[masterboy]

Diky moc za vysvetleni, ted je mi to konecne jasne proc to nefungovalo Ale mel bych jeste jeden dotaz, kdyz bych chtel deklaraci nechat s uint32_t, je mozne nejakym zpusobem provest vypocet aby to vyslo tak jak ma? Nebo to musim udelat vzdy v deklaraci a pak se s tim uz neda nic udelat? Napadlo me to zkusit pretypovat, ale nevim jak na to a zda vy to vubec slo.

[Andrea]

Ode mě to tu máš s přetypováním.

[Panda38]

Přetypování Ti píše Andrea, jen raději asi na int32_t než int, protože int je typ závislý na procesoru, u některých může mít jiný rozměr ... ale u těch M3 by to bylo v pořádku 32 bitů:

K = 1000 * ((int32_t)y[0] - (int32_t)y[1]) / ((int32_t)x[0] - (int32_t)x[1]);

(...pozdě)

[Andrea]

Já to tu zkoušela v obecném Cčku, kde uint32_t není. Šlo ostatně o znaménkovost ne o počet bitů. Správně by asi i to Q mělo být int32_t, protože by se to do int16_t nemuselo vejít.

[masterboy]

A co je lepsi? Pouzivat primo typ v deklaraci nebo to radeji nadeklarovt jako unsigned a pak ve vypoctu to pretypovat? Jde mi hlavne o ticky

Jinak dekuji Vam obema mockrat za vysvetleni

PS.: Existuje nejake pravidlo na co si dat pozor kdyz pisu nejaky vzorec aby to nedopadlo jako v mem pripade? Preci jen u uP se to spatne zjistuje

[ZdenekHQ]

Mám podobnou zkušenost s výpočtama v PLC, když jsem kámošovi dal tři vzorečky pro výpočet teploty a výsledek byl katastrofa, protože to jeho PLC umí jen klasickej integer a všechno zaokrouhlí.

[mtajovsky]

PS.: Existuje nejake pravidlo na co si dat pozor kdyz pisu nejaky vzorec aby to nedopadlo jako v mem pripade?

Používat unsigned jen tam, kde to má opodstatnění, například tam, kde na nějaké čítače nebo délky je 2^31 málo a nebudou se dělat aritmetické operace, nemíchat signed a unsigned. Jinak si zaděláváte na problémy. Nadeklarovat si něco jako unsigned a pak to při použití přetypovávat na signed nemá moc smysl. Pokud váš překladač nemá dostatečnou diagnostiku a nedá varování, pak existují nástroje, které tyhle věci dokáží odhalit - http://www.coverity.com/

Jo, a přetypovat se dá celá závorka:
(int)(y[0] - y[1])/(int)(x[0] - x[1])

[Panda38]

S unsigned si člověk nadělá někdy i další problémy, když zapomene na použitý typ a napíše např.:

for (i = 100; i >= 0; i--)

a pak se diví proč mu program zamrzl.