diff options
Diffstat (limited to 'takagi.m')
-rw-r--r-- | takagi.m | 8 |
1 files changed, 2 insertions, 6 deletions
@@ -1,10 +1,9 @@ -function [rl, U] = takagi(N, tol, rounding) +function [rl, U] = takagi(N, tol) % Autonne-Takagi decomposition of a complex symmetric (not Hermitian!) matrix. % % Args: % N (complex matrix): square, symmetric matrix N % tol (double): tolerance for symmetry check (default: 1e-13) - % rounding (integer): decimal places for rounding singular values (default: 13) % % Returns: % rl (vector): rounded singular values @@ -13,9 +12,6 @@ function [rl, U] = takagi(N, tol, rounding) if nargin < 2 tol = 1e-13; end - if nargin < 3 - rounding = 13; - end [n, m] = size(N); @@ -55,7 +51,7 @@ function [rl, U] = takagi(N, tol, rounding) [v, l, ws] = svd(N); w = ws'; - rl = round(diag(l), rounding); + rl = diag(l); % Group degenerate singular values [~, ~, ic] = unique(rl); |