比cell2mat快
我现在正在优化一些代码,并且无法想出一种比MATLAB的cell2mat
更快的cell2mat
。 目前,我的代码中cell2mat
的多次使用代表了超过15%的处理时间 。
我认为它可以比这更快,因为我知道我将传递给函数的单元阵列的结构。
基本上,单元阵列是NxN
,其中:
左上(N-1)x(N-1)
块在每个单元中包含6x6
双矩阵
右下(N,N)
单元是一个MxM
双矩阵。
其他单元格具有用于连接的正确尺寸,即:
单元(1:(N-1),N)
是6 6xM
双矩阵,单元(N,1:(N-1))
是Mx6
双矩阵。 (为了清楚起见添加了图像,其中N=207
和M=300
)
由于单元格总是会被双倍填充,并且始终为2维,因此我已经只使用了一小部分cell2mat
代码,即:
function m = myCell2Mat(c)
rows = size(c,1);
m = cell(rows,1);
% Concatenate one dim first
for n=1:rows
m{n} = cat(2,c{n,:});% 73% of the time spent here
end
% Now concatenate the single column of cells into a matrix
m = cat(1,m{:});% 25.2% of the time spent there
end
这不会改变花费的时间(就像人们可以想象的那样),大部分时间花在这些线上。
我的问题是:有没有人有关于如何去除那里的循环的想法? 我尝试了以下几点:
N=207;
M=300;
C=cell(N,N);
for ii=1:N-1
for jj=1:N-1
C{ii,jj}=rand(6);
end
end
for kk=1:(N-1)
C{N,kk}=rand(M,6);
C{kk,N}=rand(6,M);
end
C{end}=rand(M,M);
tmp1=cat(1,C{:,1:(end-1)});
LeftPart=reshape(tmp1,[],6*(size(C,2)-1));
RightPart=cat(1,C{:,end});
Res=[LeftPart RightPart];
但它没有及时显示任何改进......(并且在reshape
操作时产生错误结果)
我曾考虑过使用递归函数,但它似乎没有去任何地方。
提前致谢!
评论的结论是,没有办法显着提高cell2mat
的性能。 相反,我会建议一个更好的数据结构。
你写道,数据实际上代表了一个二维矩阵。 不要将它分成块,而是将其重新组合成4维矩阵。 假设您的Cell被称为C
, M
是相应的2D矩阵
%Old Code
q=C{1,1}
%Faster way to index the same
%1) Convert 2D matrix to 4D
blocksize=6
M=reshape(M,blocksize,size(M,1)/blocksize,blocksize,size(M,2)/blocksize);
%2) Index a block
q=squeeze(M(:,1,:,1))
重塑基本上需要零时间,您可以切换回2D视图,例如处理最后一行/列。 如果你想避免转换回来,可以一次索引多个bock。 对于最后一列使用:对应于C{206,207}
reshape(M(:,206,:,207:end),6,[])
下一篇: VS Team Services and Azure Continuous delivery subdirectory