Fast aggregates over overlapping subsets of a vector

Question

dymitr ruta 2023-2-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1912065-fast-aggregates-over-overlapping-subsets-of-a-vector

回答： Swaraj 2023-3-7

I have a long vector of values say x=rand(1e5,1), and a long but shorter cellarray with indices to the elements of x. I just want to compute means or sums over the vector subsets using these indices. There is a variable number of indices in cells and the indices are exclusive within each cell but may not be exclusive among the cells i.e. different cells may contain the same indices. I am aware of the accumarray but I assumed it can work for mutually exclusive sets of indices i.e. each element is used only once in aggregation. Is there any faster way to do that than in a loop:

n=1e5;
%Vector
x=rand(n,1); 
%Random length indices
id=arrayfun(@(i) randperm(n, randi(n)), 1:round(n/10), 'UniformOutput',false)';
%And a loop test scenario
tic; y=zeros(size(id)); for i=1:numel(id) y(i)=mean(q(id{i})); end; t=toc
%t=2.788

Cheers

2 个评论
显示无隐藏无

Dyuman Joshi 2023-2-14

在 MATLAB Online 中打开

Is there a need to store those indices? If not then, you can call them in each iteration and store mean accordingly -

n=1e5;
%Vector
y=zeros(1,round(n/10));
for k=1:round(n/10)
    y(i)=mean(q(randperm(n,randi(n))));
end

However, You have not used x in the code and we don't know what q is, thus it's not clear what you exactly want to do.

dymitr ruta 2023-2-14

编辑：dymitr ruta 2023-2-14

Unfortunately yes, these indices come out from kdtree rangesearch() function and in a form of cellarray of variable length indices' vectors. Given the vector x and indices id as input I just need to get fastest computed means

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Swaraj 2023-3-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1912065-fast-aggregates-over-overlapping-subsets-of-a-vector#answer_1187365

You are correct in saying that “accumarray” only functions for sets of indices that are mutually exclusive. You can apply a function (like mean or sum) to each cell in your cell array that contains indices by using the “cellfun” function.

Please go through the below documentation for more details.

Apply function to each cell in cell array - MATLAB cellfun (mathworks.com)