Параллельная независимая итерация элементов в массиве

Question

Я хочу, чтобы итератировать все элементы независимо (одинаковое условие для всех элементов). Я хочу остановить итерацию, когда значение перестает меняться. Я отправлю часть своего кода. Согласно некоторым исследованиям, я подумал, что это также можно сделать с помощью цикла parfor, но я не знаю, как его реализовать. Может кто-нибудь исправить мой код? Заранее спасибо.

probability = (ones(1,2048) .* 1/2048); 

Tij = sum(StateTransitionfwd); %gives us an array of 2048 elements. 

probability = ((probability * StateTransitionbwd) - (Tij .* probability)); 

threshold = (ones(1,2048) .* 0.05); 

old = zeros(1,2048); 

new = zeros(1,2048); 

while(1) 

    probability = ((probability * StateTransitionbwd) - (Tij .* probability)); 

    new = probability; 

    if old-new <= threshold 

     break 
    end 
    old = probability; 

end 

Поэтому в основном я хочу, устойчивое состояние вероятности (где она не меняющаяся больше)

Answer 1

Для одного вы не можете parallellise в то время как петля.

Для реализации цикла parfor убедитесь в следующем: - Все переменные должны быть определены до цикла парфор. Таким образом, определите выходные переменные перед циклом. - Все итерации должны быть абсолютно независимыми друг от друга.

Второе условие не выполняется в вашем случае, так как вы постоянно обновляете переменную new, делая ее не независимой.

Я не хотел даже parallellise этого кода, так как он предназначен, чтобы разбить, прежде чем он прошел через все элементы, но ради аргумента Постараюсь это:

parfor ii = 1:length(Tij) 
    probability{ii} = ((probability * StateTransitionbwd) - (Tij .* probability)); 
end 

Либо использовать клетку или матрица, в зависимости от того, что более полезно для вашей работы. Теперь вы можете пойти дальше и использовать простой if найти свой порог:

for ii = 1:length(probability)-1 
    if probability(ii+1)-probability(ii)<threshold 
     break 
    end 
end 

Как для остановки в то время как петля условия: это опасно, чтобы бежать цикл неконтролируемого некоторое условия, как и вы. Достаточно легко застревать в бесконечном цикле. Используйте дополнительное условие MaxIt для максимального количества итераций, которые необходимо выполнить перед завершением, независимо от другого условия.

while probability<threshold && MaxIt<someNumber 
    your code 
end 

Также вы проверяете достоверность полного массива.

A = logical([1;0]); 
while A 
    B=1; 
end 
C = logical([1;1]); 
while C 
    D=1; 
end 

Первый цикл while не будет работать, поскольку одна из его записей неверна; второй цикл представляет собой бесконечный цикл. Если вы хотите прервать каждую строку, как только условие будет выполнено, поместите цикл while внутри цикла for через каждую строку. Этот forloop может быть параллелизирован, если вы поместите цикл while в функцию, см. Документацию parfor для этого.