Skip to content
Snippets Groups Projects
Commit c1be34eb authored by Achilles Kappis's avatar Achilles Kappis
Browse files

Revert to using serial for-loops in obsFiltEst.m

parent bbaeee86
No related branches found
No related tags found
No related merge requests found
...@@ -3,7 +3,7 @@ ...@@ -3,7 +3,7 @@
% Author: Achilles Kappis % Author: Achilles Kappis
% e-mail: axilleaz@protonmail.com % e-mail: axilleaz@protonmail.com
% %
% Date: 19/12/2024 (DD/MM/YYYY) % Date: 26/02/2025 (DD/MM/YYYY)
% %
% Copyright: MIT % Copyright: MIT
% -------------------------------------------------- % --------------------------------------------------
...@@ -109,7 +109,7 @@ function [est, err, errSqr, normErrSqr, See] = obsFiltEst(Pm, O, Pe, Svv) ...@@ -109,7 +109,7 @@ function [est, err, errSqr, normErrSqr, See] = obsFiltEst(Pm, O, Pe, Svv)
% Sum of squared estimation errors % Sum of squared estimation errors
if nargout > 2 if nargout > 2
parfor eIdx = 1:size(Pe, 1) for eIdx = size(Pe, 1):-1:1
errSqr(eIdx, 1) = err(eIdx, :) * Svv * err(eIdx, :)'; errSqr(eIdx, 1) = err(eIdx, :) * Svv * err(eIdx, :)';
end end
end end
...@@ -117,7 +117,7 @@ function [est, err, errSqr, normErrSqr, See] = obsFiltEst(Pm, O, Pe, Svv) ...@@ -117,7 +117,7 @@ function [est, err, errSqr, normErrSqr, See] = obsFiltEst(Pm, O, Pe, Svv)
% Normalised squared errors % Normalised squared errors
if nargout > 3 if nargout > 3
% Calculate power true spectral density at virtual microphones % Calculate power true spectral density at virtual microphones
parfor eIdx = 1:length(errSqr) for eIdx = length(errSqr):-1:1
See(eIdx, 1) = Pe(eIdx, :) * Svv * Pe(eIdx, :)'; See(eIdx, 1) = Pe(eIdx, :) * Svv * Pe(eIdx, :)';
end end
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment