SIGnal Floating Point Error – Pascal – Fórum – Programujte.com

Zdravím,

robím domácu úlohu na programovanie. Máme vykresliť graf lineárnej lomenej funkcie bod za bodom do grafickej plochy v Lazaruse. Keď však počítam funkčnú hodnotu ako y:= (a*x+b)/(c*x+d); tak mi vyskočí External: SIGFPE. Všetky hodnoty sú doubly, ktorým priradím hodnoty: počítacej premennej cyklu do x a prečítané inegery z editačných polí. SIGFPE nevyskočí vždy, ale odpozoroval som, že určite vyskočí, keď c=d alebo sa aspoň jedno z a;b;c;d rovná 0. Lenže všetko sú to doubly, takže by sa nikdy nemali rovnať, som z toho zmätený. Tuší niekto ako to vyriešiť? (okrem hintov: nepoužívaj lazarus, rád by som použil niečo iné alebo aspoň vykreslenie vyvolal v pyplote ale nesmiem použiť nič, okrem základného pascalu). Ďakujem

Nerobím v pascale, ale myslím, že floating point error by mal vyskakovať vtedy, keď sa robí niečo, čo je matematicky nemožné. Pre y:= (a*x+b)/(c*x+d); mi napadá len to, že ak c=d=0, delí sa tam nulou. Skús to ošetriť a ak stále bude nejaký vstup, pre ktorý to padá, napíš ten vstup.

To som ošetril podmienkou if cx+d=0 then y=10^6, čo pre moje účely stačí, ale aj v takom prípade vyhodí SIGFPE. Nejde teda len o bod x=1. Ako klinec do rakvy: ak by aj išlo o bod x=1 teda div zero, nemal by to byť problém, lebo pascal to zvláda a výsledok si pamätá ako "inf" resp. "-inf" a okrem zaokrúhľovania s tým vie robiť prakticky každú operáciu. Ak výpočet y=ax+b/cx+d rozdelím na w=ax+b, z=cx+d, y=w/z tak SIGFPE referuje na riadok y=w/z.

#3 SamuelKo
Mrkl sem na net a mělo by to být jak říká lukáš. Tj. dělíš nulou. Takže to ošetři a je to vyřešený.

btw. Nevím jestli jsem tě správně pochopil ale zdá se mi, že říkáš, že pascal zvládá dělit nulou. To je samozřejmě hloupost ale asi sem tě jen špatně pochopil.

Pointa je, že delenie nulou som ošetril a robí to stále. K tomu deleniu nulou... 

FloatToStr(1/0) = + Inf

Možná by bylo lepší, kdybys ukázal ten kód. Tohle je vaření z vody. Ale to co si sem hodil není ošetření jen ukázka jak shodit program.

To nemalo byť ošetrenie :D toto je moje ošetrenie (je to objektové, preto tie rozmery plochy, ale chcem vám to ukázať v pôvodnej podobe): 
 

procedure hodnota(a,b,c,d,x: double);
begin
  if ((x-(Image1.Width div 2))*c+d)=0 then y:=1000000
  else
  w:= ((x-(Image1.Width div 2))*a+b);
  z:= ((x-(Image1.Width div 2))*c+d);
  y:= w/z;
  y:= y+(Image1.Height div 2);

end;

SIGFPE odkazuje na riadok v bolde.

ten lazarus je priserna vec.. ani krokovat to evidentne neumi... a zobrazit co je v promennych v aktualnim kroku uz vubec ne... </sarkazmus>

Co je vubec ve w kdyz projde ta podminka ifu?

Co je v x? Je to vetsi nez polovina Width?

Co z? Kde to je definovane? (ostatne to same plati pro y)

Ospravedlňujem sa, že som to dostatočne nevysvetlil, ale všetky premenné (a;b;c;d;x;y;w;z) sú globálne doubly, ak prejde podmienka tak w je evidentne nenulové (vzhľadom na to, že sú to doubly, teda majú obmedzenú presnosť tak je tá podmienka prakticky zbytočná, lebo presnú nulu fakt nikdy nijak netrafíme). Celú túto procedúru vyvolávam v cykle pre i:=0 do Image1.Width, pričom i je longint a jeho hodnotu pred výpočtom priradím x. a;b;c;d sú parametre, ktoré načítam z editačného poľa a uložím do príslušných premenných. 

PS.: ja som za python, najradšej by som použil pyplot a hotovo. Bohužiaľ, som bol profesorkou uistený, že to je riešenie za 0b., takže sa môžem ďalej trápiť s asi najmenej user-friendly objektovým prostredím všetkých čias.

prvni zkus to krokovani...

jinak chapu ze clovek, ktery se jako prvni ucil treba python, bere cokoliv, co ma nejake striktnejsi typy jako "unfriendly"

kazdopadne to nehaze errory pro nic za nic...

+ globalni promenne jsou vetsinou prasarna :)

Floating point error nastane, když překročíš rozsah reálného čísla. Nemusí to být přímo dělení nulou, stačí i něco dostatečně blízko nuly. Takže místo "if jmenovatel=0" dej třeba "if abs(jmenovatel)<1e-10" (tu konstantu jsem jen tak plácnul, skutečnou hodnotu si nastav tak, aby to fungovalo).

Veľmi dobrý nápad. Ale za predpokladu, že použijem napr. koeficienty 1,2,3,3 a obmedzujúcu konštantu 1e-1 a moja polovica šírky plochy (časť osi, na ktorej vykresľujem) je 300, tak maximálna hodnota čitateľa je 302 a minimálna hodnota menovateľa je 0,1. Najvyššia funkčná hodnota nie je viac ako 3020. Aj tak vyskočí SIGFPE. 

PS.: nemal by double jednoducho pretiecť?

Není víc než 3020? Viděl jsi ty skutečné hodnoty 302 a 0.1 vypsané na obrazovce, nebo si jenom myslíš, že tam byly? To je, oč tu běží...

P.S.: Jestli se nemýlím, reálná čísla právě místo přetékání házejí errory.

302 je maximum, ktoré možno dostať zo vzťahu (x-(I1.W / 2))*a+b ak a=1, b=2, I1.W=600 a x je od 1 do 600. Menovaťeľ má ohraničenie 0,1, pretože to mám v podmienke ako "if abs(menovateľ)<0,1 then menovateľ=0,1". Double má horné ohraničenie 1,7e308. 

Vidím dvě možnosti. Buď máš porouchaný procesor (krajně nepravděpodobné), nebo se ti do toho výrazu nějak dostala taková čísla, že to těch teoretických 302 a 0.1 nedá. Nech si skutečné hodnoty vypsat na obrazovku nebo použij krokování a watch, jinak chybu nikdy nenajdeš.

#7 SamuelKo
Jé, teď mi to došlo. Všechny ty příkazy za else musíš zabalit do begin - end, jinak do toho else bude patřit jenom ten první a všechny ostatní pojedou vždycky, včetně toho dělení.

Ja som idiot :D fakt to môže byť také jednoduché? :D zlé návyky z pythonu... Dik moc :) teraz už ošetrenie funguje a plotí to pekne :)

#16 Mircosoft
jeste kdybych se ho na ten if neptal hned v mym prvnim prispevku :)