イテレータopenMPのループ

Question

OpenMPを使用してコードを並列化しています。私はループを2つのイテレータで並列化しようとします。

#include <iostream> 
#include <vector> 
#include <omp.h> 

using namespace std; 

int main(int argc, char *argv[]) 
{ 
    vector<float> v = {1, 2, 3, 4}; 
    vector<float> d = {0, 0, 0, 0}; 
    vector<float>::iterator iterV, iterD; 

    for(iterV = v.begin(), iterD = d.begin(); iterV < v.end(); ++iterV, ++iterD) 
    { 
     *iterD = *iterV; 
    } 


    for(iterD = d.begin(); iterD < d.end(); ++iterD) 
    cout << *iterD << endl; 

    return 0; 
} 

このコードの私の並列バージョン：

#include <iostream> 
#include <vector> 
#include <omp.h> 

using namespace std; 

int main(int argc, char *argv[]) 
{ 
    vector<float> v = {1, 2, 3, 4}; 
    vector<float> d = {0, 0, 0, 0}; 
    vector<float>::iterator iterV, iterD; 

    iterV = v.begin(); 
    iterD = d.begin(); 

#pragma omp parallel for 
    for(int i = 0; i < v.size(); ++i) 
    { 
     *(iterD + i) = *(iterV + i) ; 
    } 


    for(iterD = d.begin(); iterD < d.end(); ++iterD) 
    cout << *iterD << endl; 

    return 0; 
}` 

Answer 1

あなたの例では、ほとんどすべてのパフォーマンスの最適化を必要としないこれは非常に簡単です、私は私の実装では、この順次コードを並列化するための最良の方法であるかどうかを知りたいです。メモリをコピーするだけです（代わりにstd::copyを使用して最適化できます）。

あなたが書いたコードは正しいものであり、パフォーマンスを得るために他の方法で書くことはほとんどありません。しかし、より洗練されたコードを維持するために、私はループイテレータを各スレッド専用に維持しようとします。これにより、コードはきれいになります（個人的な好み）。

vector<float> v = {1, 2, 3, 4}; 
    vector<float> d = {0, 0, 0, 0}; 
    #pragma omp parallel for 
    for(auto it_v = v.begin(),it_d = d.begin(); it_v!=v.end();++it_v,++it_d) 
    { 
     *it_d = *it_v; 
    } 

EDIT

のOpenMP 3.1は、この場合に、複数のループ初期化を許可しません。彼らの仕様には準拠していません。ですから、1つの方法は次のとおりです。

#pragma omp parallel 
    { 
    auto it_v = v.begin(),it_d = d.begin(); 
    #pragma openmp for 
    for(; it_v!=v.end();++it_v) 
    { 
     *it_d = *it_v; 
    } 
    }