% macierz rzadka
load bcsstk38.mat
S = Problem.A;
p = symamd(S);
tic; L = chol(S(p,p)); toc
nnz(L)/NN
spy(L'); pause
spy(S(p,p)); pause
p = colperm(S);
tic; L = chol(S(p,p)); toc
nnz(L)/NN
spy(L'); pause
spy(S(p,p)); pause
S
A = full(S);
A
A(1:10,1:10)
tic; L = chol(A); toc
