OpenMP threadprivate 指令詳解：線程私有變量與數據隔離機制

代碼示例

#include 
int alpha[10], beta[10], i;

#pragma omp threadprivate(alpha)  // alpha 為線程私有

int main() {
    // 第一個並行區域：初始化
    #pragma omp parallel private(i, beta)
    {
        int id = omp_get_thread_num();
        for (i = 0; i < 10; i++)
            alpha[i] = beta[i] = id * i;  // 每線程獨立副本
    }

    // 第二個並行區域：驗證持久性
    #pragma omp parallel
    printf("I am thread %d : alpha[3] = %d, beta[3] = %d\n",
           omp_get_thread_num(), alpha[3], beta[3]);
    return 0;
}

輸出結果分析

I am thread 0 : alpha[3] = 0, beta[3] = 0
I am thread 3 : alpha[3] = 9, beta[3] = 0
I am thread 2 : alpha[3] = 6, beta[3] = 0
I am thread 1 : alpha[3] = 3, beta[3] = 0

核心機制解析

關鍵規則

1. 聲明位置：#pragma omp threadprivate 必須放在文件作用域，變量定義之後、首次使用之前。
2. 初始化：首次進入並行區時，各線程副本按主線程值複製初始化（若有初始化表達式）。
3. 持久性：同一線程在後續並行區見到上次修改的值。
4. 動態線程：若線程數變化，新線程副本未定義；建議配合 omp_set_dynamic(0) 固定線程數。
5. 不支持：類類型構造/析構函數、引用類型、不完全類型。

典型應用場景

• 線程局部緩衝區：避免鎖競爭，如各線程獨立 rand() 種子、臨時數組。
• 狀態保持：跨並行區累積統計（計數器、求和）。
• 第三方庫上下文：非線程安全庫的上下文隔離。

對比：threadprivate vs private vs firstprivate

threadprivate: 全局/靜態變量，跨並行區持久，線程隔離
private:       任意變量，僅當前並行區，每次進入未初始化
firstprivate:  任意變量，僅當前並行區，按主線程值初始化

結語

threadprivate 是 OpenMP 實現線程局部存儲（TLS）的標準方式。理解其「全局聲明、線程隔離、跨區持久」三大特性，配合 private/firstprivate 靈活選擇，能高效解決並行區間的數據隔離與狀態傳遞問題。