OpenMPを使用したポインタ

Question

私のプログラムでOpenMPを使用しようとしています（私はOpenMPを使用している初心者です）。プログラムは2つの場所でエラーを返します。ここで

はサンプルコードです：

#include <iostream> 
#include <cstdint> 
#include <vector> 
#include <boost/multi_array.hpp> 
#include <omp.h> 

class CNachbarn { 
public: 
    CNachbarn() { a = 0; } 
    uint32_t Get_Next_Neighbor() { return a++; } 

private: 
    uint32_t a; 
}; 

class CNetwork { 
public: 
    CNetwork (uint32_t num_elements_); 
    ~CNetwork(); 
    void Validity(); 
    void Clean(); 

private: 
    uint32_t num_elements; 
    uint32_t nachbar; 

    std::vector<uint32_t> remove_node_v; 
    CNachbarn *Nachbar; 
}; 

CNetwork::CNetwork(uint32_t num_elements_ ) { 
    num_elements = num_elements_; 
    Nachbar = new CNachbarn(); 

    remove_node_v.reserve(num_elements); 
} 

CNetwork::~CNetwork() { 
    delete Nachbar; 
} 

inline void CNetwork::Validity() { 
    #pragma omp parallel for 
    for (uint32_t i = 0 ; i < num_elements ; i++) { 
     #pragma omp critical 
     remove_node_v.push_back(i); 
    } 
} 

void CNetwork::Clean() { 
    #pragma omp parallel for 
    for (uint8_t j = 0 ; j < 2 ; j++) { 
     nachbar = Nachbar->Get_Next_Neighbor(); 
     std::cout << "i: " << i << ", neighbor: " << nachbar << std::endl; 
    } 

    remove_node_v.clear(); 
} 

int main() { 
    uint32_t num_elements = 1u << 3; 
    uint32_t i   = 0; 
    CNetwork Network(num_elements); 

    do { 
     Network.Validity(); 
     Network.Clean(); 
    } while (++i < 2); 

    return 0; 
} 

私はの#pragmaが重要OMP場合

1. を知っていただきたいと思いpush_back()に適したソリューションですか？ （この問題は解決しますか？）スレッドごとに独自のベクトルを定義し、それらを結合する方がよいでしょうか？または何らかの種類のlock？

2. 私の元のコードでは、実行中のエラー：nachbar = Nachbar->Get_Next_Neighbor(&remove_node_v[i]);が表示されますが、この例では表示されません。私はOpenMPがCNachbarnコアの数として使用することを望んでいます。CNachbarnは再帰的な計算であり、他のスレッドの影響を受けてはならないからです。問題はそれをいかにスマートにするかです。 （私は私が私のシミュレーションと時間で百万回以上この関数を呼び出すので、私は、forループを開始するたびに重要であるCNachbarnを定義するために、スマートであると思ういけない

Answer 1

あなたの最初の問題に関して：。 あなた機能の有効性は、並列ループ内でシリアルパフォーマンスの下に達成するための最適な方法ですが、あなたはすでに正しい答えを与えたあなたは、スレッドごとに独立したベクトルを記入し、その後それらをマージする必要があり

inline void CNetwork::Validity() { 
#pragma omp parallel for 
for (uint32_t i = 0 ; i < num_elements ; i++) { 
    #pragma omp critical 
    remove_node_v.push_back(i); 
} 
} 

EDIT：。。。 A可能な救済策はこのように見えるかもしれません（あなたがあなたの要素へのシリアルアクセスをequire、あなたはOMPの最大数は可能なスレッドの大きさ、および個別のアクセス権を持つCNachbarnの配列を定義することによって、あなたの第二の問題は解決できる

inline void CNetwork::Validity() { 
remove_node_v.reserve(num_elements); 
#pragma omp parallel 
{ 
    std::vector<uint32_t> remove_node_v_thread_local; 
    uint32_t thread_id=omp_get_thread_num(); 
    uint32_t n_threads=omp_get_num_threads(); 
    for (uint32_t i = thread_id ; i < num_elements ; i+=n_threads) 
    remove_node_v_thread_local.push_back(i); 

    #pragma omp critical 
    remove_node_v.insert(remove_node_v.end(), remove_node_v_thread_local.begin(), remove_node_v_thread_local.end()); 
} 
} 

）ループを少し変更する必要があります

CNachbarn* meine_nachbarn=alle_meine_nachbarn[omp_get_thread_num()]