Correct error in parfor loop

Question

xiaohuo 2022-3-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1668684-correct-error-in-parfor-loop

评论： xiaohuo 2022-3-10

I am working on a parfor loop in MatLab R2021b, shown below. But there are errors: Valid indices for 'Time_h1' and 'HD_Top' are restricted in PARFOR loops. I would like to how to solve this issue? Thanks.

        
parfor j = 1:nk2
    if abs(ID_sat) < 1e-5 
        T_mesh(nt,j) = 11;          
        Time_h1 (1:11,j) = 0:t_final/10:t_final; % error: restricted in PARFOR loops    
        HD_Top (1:11,j) = 0;  % error: restricted in PARFOR loops                    
    else                                  
        ...        
    end
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2022-3-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1668684-correct-error-in-parfor-loop#answer_914694

在 MATLAB Online 中打开

        
parfor j = 1:nk2
    if abs(ID_sat) < 1e-5 
        T_mesh(nt,j) = 11;
        h1 = Time_h1(:,j);
        HT = HD_top(:,j);
        h1(1:11) = 0:t_final/10:t_final;
        Time_h1(:,j) = h1;
        HT(1:11) = 0;
        HD_Top(:,j) = HT;                    
    else                                  
        ...        
    end
end

4 个评论
显示 2更早的评论隐藏 2更早的评论

Walter Roberson 2022-3-10

在 MATLAB Online 中打开

When you assign to an array that is indexed by a parfor variable, each index has to be either scalar or the : operator; when you read from an array that is indexed by a parfor variable, each index has to be either scalar or the : operator.

t1 = Time_h1(:,j);      %okay because it uses :
t1(1:T_mesh00) = Surface_flux(1:T_mesh00,1);  %okay because it is not indexed by j
Time_h1(:,j) = t1;     %okay because it uses :