remove rounding of significant figures from takagi and bloch_messiah code

1 files changed, 2 insertions, 6 deletions

diff --git a/takagi.m b/takagi.m
index f9375a8..dafd632 100644
--- a/takagi.m
+++ b/takagi.m
@@ -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);