🚀 go-pugleaf

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
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;

	mov ebx, tab				; zaladowanie adresu tablicy
	fld QWORD PTR [ebx+8]	; zaladowanie wartosci f(x)
	fstp t1					; zachowanie wartosci do zmiennej lokalnej, zwolnienie
stosu
	mov ebx, beta			; zaladowanie adresu zmiennej beta (liczonej przez
funkcje "Beta(j, n)"
	fld QWORD PTR [ebx]		; zaladowanie wartosci spod wskaznika beta
	fcos						; cosinus (beta)
	fadd st(0), st(0)			; 2 * cos (beta)
	fstp t2					; zachowanie powyzszej wartosci do lokalnej zmiennej

	mov eax, k				; zachowanie zmiennej k
	mov temp, eax
	inc temp					; temp = k+1
	mov eax, tab				; zachowanie adresu tablicy
	mov tempptr, eax
	add tempptr, 32			; przesuniecie sie do nastepnego elementu tablicy

	INVOKE wspU, n, temp, beta, tempptr	; wywolanie funkcji wspU (n, k+1,
beta, tab[k+1])

	fld QWORD PTR [eax]		; odczytanie wyniku zwroconego przez wspU
	fstp t3					; zapisanie wyniku

	mov eax, k				; zachowanie zmiennej k
	mov temp, eax
	inc temp					; temp = k+1
	inc temp					; temp = k+2
	mov eax, tab				; zachowanie adresu tablicy
	mov tempptr, eax
	add tempptr, 64			; przesuniecie sie w tablicy o 2 elementy

	INVOKE wspU, n, temp, beta, tempptr	; wywolanie funkcji wspU (n, k+2,
beta, tab[k+2])

	fld QWORD PTR [eax]		; odczytanie wyniku zwroconego przez wspU
	fstp t4					; zapisanie wyniku

	fld t4					; st2 = U(k+2)
	fld t3					; st1 = U(k+1)
	fld t2					; st0 = 2 * cos(beta)

	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

	jmp Koniec

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

Koniec:
	ret ;koniec funkcji



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...

wspU ENDP