如何解决MATLAB parfor循环中B_mat与C_mat的Broadcast变量警告

问题原因 #

当前代码通过线性索引n转换为(r,k)后再取B_mat{r}/C_mat{r}，MATLAB无法确定r的取值规律，会将整个B_mat/C_mat广播到所有并行工作进程，触发通信开销警告。

修改方案 #

核心思路是让parfor循环直接关联r的维度，让MATLAB明确每个工作进程仅需要B_mat/C_mat的对应元素，而非整个数组。

方案1：拆分循环（推荐，开销更低）

将原线性索引的parfor拆分为外层parfor遍历r，内层串行遍历k，这样每个工作进程仅获取对应r的B_mat{r}和C_mat{r}：

% 已定义变量
N_RIPETIZIONI=2;
K=201;
parfor r=1:N_RIPETIZIONI
    B = B_mat{r};
    C = C_mat{r};
    % 处理当前r对应的所有k
    for k=1:K
        % 写入原n循环中对应(r,k)的业务逻辑
    end
end

方案2：保留独立迭代（适合需每个(r,k)并行的场景）

预先生成所有(r,k)的组合列表，让parfor直接遍历该列表，MATLAB可分析出每个迭代仅需对应r的元素：

% 已定义变量
N_RIPETIZIONI=2;
K=201;
% 生成r和k的对应列表
r_list = repelem(1:N_RIPETIZIONI, K);
k_list = repmat(1:K, 1, N_RIPETIZIONI);

parfor idx=1:length(r_list)
    r = r_list(idx);
    k = k_list(idx);
    B = B_mat{r};
    C = C_mat{r};
    % 写入原n循环中的业务逻辑
end

效果说明 #

两种修改方式都能让MATLAB识别到B_mat/C_mat的访问是索引化的、每个进程仅需部分元素，从而避免广播整个数组，消除警告并减少不必要的通信开销。

内容的提问来源于stack exchange，提问作者Lusian