Skip to content

Commit

Permalink
Update svd.m
Browse files Browse the repository at this point in the history
  • Loading branch information
bchaoss authored Nov 22, 2023
1 parent 299b745 commit f0f3d5f
Showing 1 changed file with 0 additions and 3 deletions.
3 changes: 0 additions & 3 deletions svd.m
Original file line number Diff line number Diff line change
Expand Up @@ -42,18 +42,15 @@
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');
Q = Q*(Q_0*Qy); S = Sy;
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;

Expand Down

0 comments on commit f0f3d5f

Please sign in to comment.