использование parfor в matlab для кода решетки boltszmann

Question

Я работаю над методом решетчатой ​​болтцман, и я написали код matlab. Я хотел бы распараллелить некоторые части кода, но я новичок в этом, поэтому я был бы признателен за вашу помощь. Я хотел бы знать, если это возможно, использовать parfor для этой части (оператор столкновений):

for i=1:lx 
    for j=1:ly 
     fork=1:9 
      f(k,i,j)=f(k,i,j) .* (1 - omega) + omega .* feq(k,i,j); 

     end 
    end 
end  

Я пытался заменить удаленные от центра цикл с parfor, но код, кажется медленнее.

любые предложения?

заранее спасибо

Answer 1

Вы должны быть в состоянии сделать всю эту операцию с одной строки кода без петель:

f = f.*(1 - omega) + omega .* feq; 

На моем компьютере с 2 ядрами и, начиная с:

f = rand(9,400,400); 
feq = rand(9,400,400); 

[lx,ly,lz] = size(f); 

omega = rand(1); 

Ваша петля занимает 0.087933 секунды, цикл parfor занимает 1.166662 секунды, и этот метод принимает 0.009388 секунд. Если можете, всегда вектурируйте свой код.