[U,X] = qr(randn(80));
[V,X] = qr(randn(80));
S = diag(2.^(-1:-1:-80));
A = U*S*V;
[QC,RC] = clgs(A)
semilogy(diag(RC));
