Matlab中的平行迎风方案

Question

声明

假设我们在Matlab中考虑了一个简单的迎风方案

for k = 1:kFin % time-loop
    % space-loop
    un(2:nx) = u(2:nx) - cfl*(u(2:nx) - u(1:nx-1)); 
    % BC
    un(1) = un(nx);
    % update
    u = un; 
end

问题

如何通过 parfor 或 spmd 或其他方式并行化它？