🚀 go-pugleaf

MrokU pisze:
[...]
> Generalnie wyliczać mam taki wzorek:
>
> U(k) = f(x) + 2 * U(k+1) * cos(B) - U(k+2)
[...]

> Na razie, po pierwszych próbach funkcja (oczywiście nie działa dobrze)
> wygląda tak:
> 
> ;******************************************************************
> ;****** Procedura obliczajaca pierwsze 2 wyrazy ciagu Un
> ;****** Wg wzoru:
> ;****** double U (int k, double beta, Element *tab)  //obliczanie
> wspólczynnika U - wzór 28
> ;****** {
> ;******		if ( k > n )
> ;******			return (0.0);
> ;******		else
> ;******			return (tab[k].f_x + (2.0) * ( U(k+1, beta, tab) * cosf
> (beta) ) - U(k+2, beta, tab));
> ;****** }
> 
> ;****** definicja elementu tablicy
> 
> ;****** struct Element
> ;****** {
> ;****** 	double x;
> ;******	double f_x;
> ;****** 	double g_x;
> ;****** 	double blad;
> ;****** };
> 
> wspU PROC n:DWORD, k:DWORD, beta:PTR, tab:PTR

 Wiesz, że ten prototyp wygląda inaczej niż ten z C i musisz tę
funkcję (znacznie) inaczej uruchamiać, prawda?

> LOCAL t1:QWORD, t2:QWORD, t3:QWORD, t4:QWORD, temp:DWORD, tempptr:PTR
> 	pushad
> 	finit
> 
> 	mov eax, k
> 	cmp eax, n				; if (k > n)
> 
> 	ja Zwroc0				; return 0;

 Załóżmy najpierw, że w tym punkcie warunek jest spełniony i skok JA
wykonany. W końcu kiedyś będzie się to musiało zdarzyć. Spójrz teraz,
co się dzieje pod "Zwroc0".

[...]
> 	fmul st(0), st(1)			; st0 = 2 * cos(beta) * U(k+1)
> 	fld t1					; st0 = f(x)
> 	fadd st(0), st(1)			; st0 = f(x) + 2 * cos(beta) * U(k+1)
> 	fsub st(0), st(3)			; st0 = f(x) + 2 * cos(beta) * U(k+1) - U(k+2)
> 
> 	popad
> 
> 	fstp QWORD PTR [eax]	; zapisanie wyniku

 Ta sama uwaga, co poniżej dla zwracania zera.

> 	jmp Koniec
> 
> Zwroc0:
> 	fldz						; ladujemy 0 do st0
> 	popad
> 	fstp QWORD PTR [eax]; zapisanie zera pod adres wynikowy

 Wielkie pytanie, co jest w EAX w tym miejscu, gdy "ja Zwroc0"
zostanie wykonane. A z resztą, nawet jak nie będzie wykonane, to co
jest w EAX? Wykonując POPAD, przywracasz wartość z momentu
uruchomienia funkcji, a rejestr EAX jest używany dowolnie, więc jego
wartość jest nieokreślona. Piszesz więc pod nieznany adres (przy
pierwszym uruchomieniu funkcji). Inicjalizuj EAX na tyle, ile trzeba
(czyli pewnie tab+8).

> No i teraz nie wiem dlaczego, ale modyfikowane są wartości tab[].x
> czyli te, pod adresami tab + x * 32 (bo tyle zajmują 4 double).
> 
> na szczęście błąd ten jest mało ważny. Jednak dla identycznych
> parametrów, funkcja w C i w ASMie zwraca zupełnie inne wyniki i tego
> już nie kumam... Patrze w te instrukcje i nie widzę gdzie może być
> błąd...

 Pewnie wszystko wina niezainicjalizowanych wskaźników. Jeśli nie, to
pokaż kod wywołujący.

-- 
Pozdrawiam/Regards - Bogdan                     (GNU/Linux & FreeDOS)
Kurs asemblera x86 (DOS, GNU/Linux):http://rudy.mif.pg.gda.pl/~bogdro
Grupy dyskusyjne o asm:  pl.comp.lang.asm alt.pl.asm alt.pl.asm.win32
www.JabberPL.org www.TorProject.org Soft (EN): miniurl.pl/bogdro-soft