Mentions légales du service

Skip to content
Snippets Groups Projects
Commit 795ff00a authored by hhakim's avatar hhakim
Browse files

Fix erroneous use of spdiag() coming from smallbox (replacing it by sparse(diag()).

parent e56ab7c0
Branches
Tags
No related merge requests found
...@@ -414,7 +414,7 @@ classdef FaustFactory ...@@ -414,7 +414,7 @@ classdef FaustFactory
else else
[lambda, core_obj, Ddiag] = mexHierarchical_factCplx(U, mex_params, Lap); [lambda, core_obj, Ddiag] = mexHierarchical_factCplx(U, mex_params, Lap);
end end
D = spdiag(Ddiag); D = sparse(diag(Ddiag));
F = Faust(core_obj, isreal(U)); F = Faust(core_obj, isreal(U));
varargout = {F, D, lambda, p}; varargout = {F, D, lambda, p};
end end
...@@ -488,7 +488,7 @@ classdef FaustFactory ...@@ -488,7 +488,7 @@ classdef FaustFactory
end end
end end
[core_obj, D] = mexfgftgivensReal(Lap, J, t, verbosity); [core_obj, D] = mexfgftgivensReal(Lap, J, t, verbosity);
D = spdiag(D); D = sparse(diag(D));
FGFT = Faust(core_obj, true); FGFT = Faust(core_obj, true);
end end
...@@ -570,8 +570,8 @@ classdef FaustFactory ...@@ -570,8 +570,8 @@ classdef FaustFactory
[W1,D1] = matfaust.FaustFactory.eigtj(M*M', J, varargin{:}); [W1,D1] = matfaust.FaustFactory.eigtj(M*M', J, varargin{:});
[W2,D2] = matfaust.FaustFactory.eigtj(M'*M, J, varargin{:}); [W2,D2] = matfaust.FaustFactory.eigtj(M'*M, J, varargin{:});
S = diag(W1'*M*W2); S = diag(W1'*M*W2);
[~,I] = sort(abs(S), 'desc') [~,I] = sort(abs(S), 'desc');
S = spdiag(S(I)) S = sparse(diag(S(I)));
sign_S = sign(S); sign_S = sign(S);
S = S*sign_S; S = S*sign_S;
Id = eye(size(S)); Id = eye(size(S));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment