% Primerjamo LU brez pivotiranja, LU z delnim pivotiranjem in LU s
% kompletnim pivotiranjem na razliènih matrikah

% 1) nakljuèna matrika 250 x 250

disp('Prvi test : nakljuèna matrika 400 x 400')

randn('state',3); 
A = randn(400);
x0 = randn(400,1);
b = A*x0;

format short e

[x1,x2,x3,res,pg]=lucompare(A,b,x0);

res
pg
cnd = cond(A)

pause

% 2) matrika, kjer imamo maksimalno pivotno rast pri delnem pivotiranju

disp('Oblika matrike z maksimalno pivotno rastjo pri delnem pivotiranju')

A5 = maxpivot(5)

pause

disp('Test z matriko 50 x 50 z maksimalno pivotno rastjo pri delnem pivotiranju')

A = maxpivot(50);
x0 = randn(50,1);
b = A*x0;

[x1,x2,x3,res,pg]=lucompare(A,b,x0);

res
pg
cnd = cond(A)

pause

% 3) matrika, kje imamo pivotno rast n pri kompletnem pivotiranju

disp('Hadarmardova matrika 8 x 8');

A8 = hadamard(8)

pause

disp('Test z Hadamardovo matriko 40 x 40');

A = hadamard(40);
x0 = randn(40,1);
b = A*x0;

[x1,x2,x3,res,pg]=lucompare(A,b,x0);

res
pg
cnd = cond(A)

pause

% 4) Golubova matrika

disp('Golubova matrika')

A6 = golub(8,10)

pause

disp('Test z matriko 8 x 8 z visoko obcutljivostjo')

A = A6;
x0 = randn(8,1);
b = A*x0;

[x1,x2,x3,res,pg]=lucompare(A,b,x0);

res
pg
cnd = cond(A)


