OpenMP разреженного Jacobi

Question

Я пытаюсь определить, есть ли способ распараллеливания методы Якобите с использованием форматов разреженных матриц (в частности, сжатый формат строка)

У меня есть рабочая разреженная матрица Якоби. Я не знаю, могу ли я разместить

! $ OMP PARALLEL DO

директив по середине делают петлю, потому что х время как запись и чтение с. Я предполагаю, что внутренний цикл do может иметь его, но тот же t перезаписывается, поэтому я не знаю, возможно ли это там. Я что-то забыл? Благодарю.

x(:) = 0 
do p = 1, numIterations 
    do i=1, n 
     t=b(i) 
     do j = IA(i), IA(i+1) - 1 
      if j=i 
       d=A(j) 
      else 
       t = t - A(j) * x(jA(j)) 
      end if 
     end do 
     x(i) = t/d 
    end do 
end do 

Answer 1

Это правда, у вас есть зависимость от t во внутреннем цикле, так как он используется в качестве аккумулятора. Однако это также означает, что у вас может быть личная копия t в каждом из потоков (так как массивы A и x не записываются в цикле, значение t зависит только от значения j, что также является потоком частного).

Следующие должны работать:

x(:) = 0 
do p = 1, numIterations 
    do i=1, n 
     t=0 
     !$OMP PARALLEL DO 
     !$OMP REDUCTION(+:t) 
     do j = IA(i), IA(i+1) - 1 
      if j=i 
       d=A(j) 
      else 
       t = A(j) * x(jA(j)) 
      end if 
     end do 
     x(i) = (b(i)-t)/d 
    end do 
end do 

Обратите внимание, что d можно записать только один из потоков, так что переменная может совместно использоваться betewen нитей, без каких-либо зависимостей петель переносимых на d.