How to resolve indexing the sliced variable error in the case of for nested with parfor?

Question

Rahul Bhadani 2021-4-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/817330-how-to-resolve-indexing-the-sliced-variable-error-in-the-case-of-for-nested-with-parfor

评论： Qitr 2023-11-11

I have following use case:

Nth = 102;
p = 0.0000001:0.1:0.99999;
NthArray = (Nth - 100 > 0)*Nth + (Nth - 100 <=0)*1:1:Nth+100;
III = zeros( length(p), length(NthArray));
parfor ip = 1:numel(p)
    for inth = 1:numel(NthArray)
        III(ip, inth) = g( p(ip), ntharray(inth));
    end
end

However, upon running it I get following error:

Error: When indexing the sliced variable 'III', the range of the for-loop variable 'inth' must be a row vector of positive constant numbers or variables. For more information, see Parallel
for Loops in MATLAB, "Nested for-Loops with Sliced Variables".

How should I resolve this error?

Thanks.

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

请先登录，再进行评论。

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

Answer 1

Edric Ellis 2021-4-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/817330-how-to-resolve-indexing-the-sliced-variable-error-in-the-case-of-for-nested-with-parfor#answer_688440

在 MATLAB Online 中打开

The problem here is that the inner loop bounds must be a constant - basically this means that you must compute it outside the parfor loop, like so:

Nth = 102;
p = 0.0000001:0.1:0.99999;
NthArray = (Nth - 100 > 0)*Nth + (Nth - 100 <=0)*1:1:Nth+100;
III = zeros( length(p), length(NthArray));
% Must calculate inner loop bounds outside the PARFOR loop:
n_Nth = numel(NthArray);
g = @plus; % Dummy definition of "g"
parfor ip = 1:numel(p)
    for inth = 1:n_Nth
        III(ip, inth) = g( p(ip), NthArray(inth));
    end
end
disp('Success!')
Success!

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

Qitr 2023-11-10

移动：Edric Ellis 2023-11-10

在 MATLAB Online 中打开

I have the same problem in the following code. Even though, both indexes are constant I obtaine error:

When indexing the sliced variable 'RintE_mat', the range of the for-loop variable 'itm2' must

be a row vector of positive constant numbers or variables. For more information, see Parallel

for Loops in MATLAB, "Nested for-Loops with Sliced Variables".

What am I doing wrong?

parfor itm=1:2*NmaxAMM
        gz_real_coeff = squeeze(gzvec_real(itm,:));
        gz_real_interp(itm,:) = sqrt(interp1(grid_vec,gz_real_coeff.^2,fvec,'spline','extrap'));
        for itm2=1:2*NmaxAMM
            RintE_coeff = squeeze(RintE_mat0(itm,itm2,:));
            RintH_coeff = squeeze(RintH_mat0(itm,itm2,:));
            RintE_interp_real = interp1(grid_vec,real(RintE_coeff),fvec,'spline','extrap');
            RintE_interp_imag = interp1(grid_vec,imag(RintE_coeff),fvec,'spline','extrap');
            RintE_mat(itm,itm2,:) = RintE_interp_real + 1i*RintE_interp_imag;
            RintH_interp_real = interp1(grid_vec,real(RintH_coeff),fvec,'spline','extrap');
            RintH_interp_imag = interp1(grid_vec,imag(RintH_coeff),fvec,'spline','extrap');
            RintH_mat(itm,itm2,:) = RintH_interp_real+1i*RintH_interp_imag;
        end
end

Edric Ellis 2023-11-10

在 MATLAB Online 中打开

Here's an executable version of my previous comment. Hopefully this includes the salient points from your code, and this does work.

limit = 2*3;
out = zeros(limit, limit, 2);
parfor i = 1:limit
    for j = 1:limit
        out(i,j,:) = [i j];
    end
end
Starting parallel pool (parpool) using the 'Processes' profile ...
Parallel pool using the 'Processes' profile is shutting down.
disp(out)
(:,:,1) =

     1     1     1     1     1     1
     2     2     2     2     2     2
     3     3     3     3     3     3
     4     4     4     4     4     4
     5     5     5     5     5     5
     6     6     6     6     6     6


(:,:,2) =

     1     2     3     4     5     6
     1     2     3     4     5     6
     1     2     3     4     5     6
     1     2     3     4     5     6
     1     2     3     4     5     6
     1     2     3     4     5     6