Dotaz: Jak si overit, jestli vysledne cislo neni -nan || nan

[Bundas]

U teto matematicke operace stranaB = sqrt(pow(stranaC,2) - pow(stranaA,2)); obcas vyjde, ze vysledek je -nan nebo nan. To znamena, ze je to matematicky nemozny, nebo tk neco.


Nevite, jak si overit, jestli vysledek je nebo neni -nan || nan? Nejak takhle?

if((stranaB = sqrt(pow(stranaC,2) - pow(stranaA,2)) == -nan || (stranaB = sqrt(pow(stranaC,2) - pow(stranaA,2)) == nan){}

 

nebo jak?

 

predem diky za pomoc

Odpovědi

[vdusek]

Začal bych zde http://www.mathworks.com/help/matlab/ref/nan.html a ošetřil proěnné podle tyu operace.

 

Např. kontroloval zda nedělím nulou, nedělám odmocninu ze záporného čísla apod. Poté výsledek +-NaN nemohu získat

[chrono]

Ak nevadí, že sa použije C99, tak sa dá použiť funkcia isnan.

[Sněhulák]

Spíš si zařiď, abys nedělal odmocninu ze záporné ho čísla, pak ti nebude vycházet NaN. Stačí neodečítat od menšího čísla větší a žádné ověřování na NaN nebude potřeba. Taky bys to mohl trochu upravit, aby se to všechno nepočítalo zbytečně dvakrát...

[potato]

Pokud to záporné čislo vzniká kvůli konečné přesnosti, je to obecně problém, pokud nechceš striktní IEEE sémantiku (tedy pomalou matematiku): splnění nějaké podmínky dvěma floating point čísly nezaručuje splnění téže podmínky v budoucnu, když se třeba mezitím přesnesou z registrů do paměti a zpět...

Vím-li, že se potýkám jen s konečnou přesností, tak normálně stačí

sqrt(fmax(stranaC*stranaC - stranaA*stranaA, 0.0))

(superpomalou funkci pow() na vynásobení dvou čísel používá jen prase).

[lertimir]

Předpokádám, že jsou to zase ty trojúhelníky a to, co počítáš je Pythagorova věta. NaN ti vyjde, když sqrt má záporný argument, což je když A>C. A to je správně i z pohledu trojúhelníků, protože pravoúhlý trojúhelník neuděláš, když jedna odvěsna bude delší než přepona. Takže posoudit ifem na C>A a jen v tom případě vzorec počítat.

[mluno]

Ještě pitomnější je, když se NaN zacyklí. Mějjme třeba algoritmus filtru prvního řádu y_nové = k1*x + k2*y_předchozí je jasné, že pokud se tam jednou jedinkrát dostane NaN tak už nikdy nic jiného nedostanete (NaN + cokoliv je zase NaN). A to samý platí pro inf, které se při dělení nulou obdrží (a to spíš měl být náš případ, ale mám pocit, že to bylo NaN). To se mi to dobře mluví, ale dva dny jsme pitvali program a už jsme chtěli do mathworksu volat o pomoc a chtít vrátit peníze a já nevim co, než jsme toto objevili. V našem případě to vznikalo dělením nulou, tak jsme zajistili, že pokud je číslo menší než něco, nahradí se tím něco a problém zmizel (vzniklo sice veliké číslo, ale to se potom o kus dál násobí tim malym a výsledkem je skoro nula). Jinak samozřejmně jsou tu funkce isnan, isfinite, isfloat a další, kterými by se to také vyřešilo, pro náš případ bylo výhodnější eliminovat to dělení nulou tak, jak jsem psal. Ve vašem případě by se dalo přidat pod to odmocnítko abs() což by sice nezpůsobilo NaN ale dávalo by to chybné výsledky v případě chybného zadání. Takže spíš to zachitit před tím odmocněním např. xxx = pow(stranaC,2) - pow(stranaA,2);
if xxx>0
...stranaB = sqrt(xxx);
else
...///něco je špatně, když vychází pod odmocnítkem záporné
číslo
...// případně stranaB = j*sqrt(-xxx) pokud to má smysl - u délek stran asi - imaginární délka neexistuje
end