diff --git a/svd.m b/svd.m index 2a60e00..86455c9 100644 --- a/svd.m +++ b/svd.m @@ -42,10 +42,8 @@ if p < tol q = q + 1; V{q} = Q_0' * d; - fprintf("here") else if q > 0 - fprintf("there.") k = size(S,1); Y = [S, cell2mat(V)]; [Qy, Sy, Ry] = svd(Y, 'econ'); @@ -53,7 +51,6 @@ R1 = Ry(1:k,:); R2 = Ry(k+1:end,:); R = [R*R1; R2]; d = Qy' * d; - %clear Y Qy Sy Ry end V = {}; q = 0;