Zadanie 1. Porównaj efektywnoæ dwóch metod rozwi¹zywania uk³adów równañ liniowych dla uk³adów ró¿nej wielkoci (ró¿nej liczby równañ).
Mo¿esz naladowaæ powy¿szy eksperyment wykorzystuj¹c obie metody w jednej pêtli na losowo wybranych macierzy i wektorze. Wyniki zapisz w dwóch wektorach np. T1 i T2.
Mój kod:
Kod: Zaznacz cały
> T1 := Vector();
T2 := Vector();
for i from 20 by 20 to 100 do
b := RandomVector(i);
A1 := RandomMatrix(i);
s := time(); GaussJordan(A1, b);
T1 := <T1, time()-s>:
G := RandomVector(i);
CH := Cholesky(A1);
CH1 := SubMatrix(CH, 1 .. i, 1 .. i);
CH2 := SubMatrix(CH, 1 .. i, 1+i .. 2*i);
G := RozwiazTrojkatny(CH2, G);
WYNIK := RozwiazTrojkatny(CH2, G);
s := time(); #inicjalizuje stoper
T2 := <T2,i> : #zapisuje rozmiar macierzy w wektorze T2
end do;
> print(T2, T1);