From fa3ed5eaf67bd54e8f1fa4a1bc2c7ffbd1898d6f Mon Sep 17 00:00:00 2001 From: "Eugeniy E. Mikhailov" Date: Wed, 20 Nov 2024 17:47:16 -0500 Subject: remove rounding of significant figures from takagi and bloch_messiah code --- takagi.m | 8 ++------ 1 file changed, 2 insertions(+), 6 deletions(-) (limited to 'takagi.m') 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); -- cgit v1.2.3