比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=207M=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被称为CM是相应的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,[])

    链接地址: http://www.djcxy.com/p/88563.html

    上一篇: Go faster than cell2mat

    下一篇: VS Team Services and Azure Continuous delivery subdirectory