0

我试图证明代码可以并行运行，代码运行良好，除非我将数据分配给切片变量，我相信我尊重 Matlab 的 parfor 限制，实际上它对于一些 parfor 迭代运行顺利，但是某些点 Matlab 崩溃要求退出。

将值分配给“功能”时，问题出在最后一行（我也尝试过使用单元格但同样的问题）：

if (matlabpool('size')==0) 
    matlabpool();
end


%Prepare iterations
numTrain = size(Xnorm(1000:1250,:),1);
numFeatures = size(Xnorm,2);
classData=Y(1000:1250,:);

%Linear Space for Grid
zz = linspace(0.1,20,10);
zsize=length(zz);

%Definitions to use inside parfor
features=zeros(numFeatures,1);
fmax=@(x) max(x);

parfor i = 1:numFeatures

    %Internal variables definition
    CsCross=zeros(zsize,1);
    SigmaCsMax=zeros(zsize,1);

    fprintf('\nFeature:%i...',i);

    for j=1:zsize

    %Precompute kernel

       K =  [ (1:numTrain)' , gpuKernel(Xnorm(1000:1250,i),zz(j)) ];

        for iter = 1:zsize

            options = sprintf('-c %d -t 4 -v 10 -q',zz(iter));
            model = svmtrain(classData, K, options);
            CsCross(iter)=model;

        end

        SigmaCsMax(j)=fmax(CsCross);
    end
   temp=fmax(SigmaCsMax);

   %%% HERE THE ERROR IS CAUSED. When removed code runs smoothly
   features(i)=temp;

end

有什么建议吗？我究竟做错了什么？（只是说代码可以在没有 parfor 的情况下工作）

提前致谢，

塞尔吉

4

1 回答 1

0

我找到了让它工作的方法。我将“特征”预定义为 features=[] （没有任何预定义的大小），然后让 parfor 构建整个数组。令人惊讶的是，结果尊重了顺序。这里的代码（我省略了初始部分）：[.....]

%HERE I DEFINE AND EMPTY ARRAY
features=[];


parfor i = 1:numFeatures

    %Internal variables definition
    CsCross=zeros(zsize,1);
    SigmaCsMax=zeros(zsize,1);
    temp=0;

    fprintf('\nFeature:%i...',i);

    for j=1:zsize

    %Precompute kernel

       K =  [ (1:numTrain)' , gpuKernel(Xnorm(1000:1250,i),zz(j)) ];

        for iter = 1:zsize

            options = sprintf('-c %d -t 4 -v 10 -q',zz(iter));
            model = svmtrain(classData, K, options);
            CsCross(iter)=model;

        end

        SigmaCsMax(j)=fmax(CsCross);
    end
   temp=fmax(SigmaCsMax);

   %%% EVERY ITERATION ADD NEW RESULTS

   features=[features,temp];

end

在此处找到有关其工作原理的更详细说明：http: //www.mathworks.es/es/help/distcomp/getting-started-with-parfor.html#brdqn6p-1

希望它可能会有所帮助。

塞尔吉

于 2014-07-02T16:07:08.270 回答