編程問題使用OpenMP

Question

我有如下使用OpenMP問題描述：

我有這樣

subroutine ... 
    ... 
    do i=1,N 
    .... 
    end do 
end subroutine ... 

和OpenMP代碼的串行代碼是

subroutine ... 
    use omp_lib 
    ... 
    call omp_set_num_threads(omp_get_num_procs()) 
    !$omp parallel do 
    do i=1,N 
    .... 
    end do 
    !$omp end parallel do 
end subroutine ... 

沒有問題與編譯，但是當我運行該程序時，有兩個主要問題與串行代碼的結果相比：

1. 該程序的運行速度比串行代碼還要慢（假設在do循環中執行矩陣乘法（matmul））
2. 與串行代碼相比，數字精度似乎有所下降（我對它進行了檢查）

任何想法可能會發生什麼？
謝謝，
小玉

Answer 1

在使用OpenMP並行的情況下，你需要指定線程程序是使用數量。您可以通過使用環境變量OMP_NUM_THREADS（例如，通過撥打您的程序

OMP_NUM_THREADS=5 ./myprogram 

使用5個線程執行它。

或者，您可以設置運行時的線程數omp_set_num_threads（documentation）。

旁註

1. 不要忘記設置私有變量，如果在循環中有什麼！
   實施例：

   $ OMP平行做私人（prelimRes）
   做I = 1，N
   prelimRes = myFunction的（ⅰ）
   RES（ⅰ）= prelimRes + someValue中
   端做
   ！$ OMP到底並行做

   注變量如何prelimRes聲明private讓每個線程都有自己的workspac即

2. 根據循環內的實際情況（即使用OpenBLAS），由於不同的相似處理，您的結果可能確實有所不同（變化應該小於1e-8關於雙精度變量）。

3. 如果您不確定發生了什麼，則應在程序運行時使用htop或類似程序檢查CPU負載。

附錄：設置線程數，以自動匹配的CPU 的數量如果您想使用有用的最大線程數，例如使用盡可能多的線程，因爲有CPU的，你可以這樣做，通過使用（就像你在你的問題說）：

subroutine ... 
    use omp_lib 
    ... 
    call omp_set_num_threads(omp_get_num_procs()) 
    !$omp parallel do 
    do i=1,N 
     .... 
    end do 
    !$omp end do 
    !$omp end parallel 
end subroutine ...