2015-05-05 69 views
0

不幸的是,我不得不运行一个简单但巨大的3次嵌套for循环。循环填充表示三维坐标中函数强度的三维矩阵。 (需要填充整个矩阵,因为我需要它来绘制3D强度云图...)为什么matlab不能并行运行我的简单嵌套for循环?

我想并行运行它,因为使用我的笔记本电脑需要数周时间才能顺序运行:

% Initialize vars 
c0values = 400:1:600; 
n0values = 1:1:10000; 
phivalues = 0:0.1:2*pi; 
mtxHeight = zeros(numel(c0values),numel(n0values),numel(phivalues)); 
% Run over c0 
parfor c0ctr=1:numel(c0values) 
    % Run over n0 
    for n0ctr=1:numel(n0values) 
     % Run over phi 
     for phictr=1:numel(phivalues) 
      % Sum over ell 
      dHeight = 0; 
      for ell=1:10 
        dHeight = dHeight + fnToMaximize(mtxObs(ell,:), ... 
         c0values(c0ctr), ... 
         n0values(n0ctr), ... 
         phivalues(phictr)); 
      end  
      mtxHeight(c0ctr,n0ctr,phictr) = dHeight; 
     end 
    end 
end 

由于某些原因,Matlab抱怨mtxHeight无法分类并拒绝运行并行化的代码。

我在做什么是不允许的?

+2

并行化并不是一件容易的事。特别是对于一台机器。 – knedlsepp

回答

0

MATLAB不能处理3个嵌套循环与parfor我认为。我建议你在内部for循环之外更新mtxHeight
类似这样的:

parfor c0ctr=1:numel(c0values) 
    % Run over n0 
    for n0ctr=1:numel(n0values) 
     % Run over phi 
     dHeight = zeros(1,numel(phivalues)); 
     for phictr=1:numel(phivalues) 
      % Sum over ell 
      for ell=1:10 
        dHeight(phictr) = dHeight(phictr) + fnToMaximize(mtxObs(ell,:), ... 
         c0values(c0ctr), ... 
         n0values(n0ctr), ... 
         phivalues(phictr)); 
      end  
     end 
     mtxHeight(c0ctr,n0ctr,:) = dHeight;  
    end 
end