openMP v. 2.0을 사용하여 필터 반복자를 병렬화하십시오.

Question

큰 벡터가 있습니다. 긍정적 인 회원만을 평행하게 늘리고 싶습니다.

class is_positive_number { 
public: 
    bool operator()(int x) { return 0 < x; } 
}; 

int main() 
{ 
    int numbers_[] = { 0, -1, 4, -3, 5, 8, -2, 3, 4, -1 }; 
    const int N = sizeof(numbers_)/sizeof(int); 
    std::vector<int> test_vec(numbers_, numbers_ + N); 

    std::cout << "max: " << omp_get_max_threads() << std::endl; 
    int number_of_threads = 4; 
    int SubVecLen = N/number_of_threads, 
     LeftOvers = N % number_of_threads; 
    int k = 0, i = 0, n = 0; 
    std::vector<int> start_vec, end_vec; 
    while (n < number_of_threads) 
    { 
     k += i; 
     i = (SubVecLen + (LeftOvers-- > 0 ? 1 : 0)); 
     start_vec.push_back(k); 
     end_vec.push_back((k + i - 1)); 

     ++n; 
    } 
    for (size_t n = 0; n < number_of_threads; n++) 
    { 
     std::cout << start_vec[n] << "\t" << end_vec[n] << "\n"; 
    } 

    is_positive_number predicate; 
    typedef boost::filter_iterator<is_positive_number, std::vector<int>::iterator> 
     FilterIter; 

#pragma omp parallel 
    { 
#pragma omp for 
     for (int s = 0; s < 4; s++) 
     { 
      FilterIter filter_iter_first(predicate, test_vec.begin(), (test_vec.begin() + start_vec[omp_get_thread_num()])); 
      FilterIter filter_iter_last(predicate, test_vec.begin(), (test_vec.begin() + end_vec[omp_get_thread_num()] + 1)); 

      for (auto iter = filter_iter_first; 
       iter != filter_iter_last; iter++) 
      { 
       std::cout << "num: " << *iter << std::endl; 
       (*iter) = (*iter) + 1; 
      } 

     } 
    } 

    for (size_t n = 0; n < test_vec.size(); n++) 
    { 
     std::cout << test_vec[n] << "\n"; 
    } 
    return 0; 
} 

결과가 잘못되었습니다.

그러나 내가 해결하려고하는 실제 문제는이 것보다 복잡합니다. 나는 더 쉽게 이해할 수 있도록 내 질문을 단순화하려고 노력했다. openMP v.2.0을 사용하도록 제한되었습니다. 그 이유는 내 벡터를 분할하고 벡터의 각 부분을 특정 스레드 전용으로 사용하려고했기 때문입니다.

어떤 아이디어라도 감사 할 것입니다.

Answer 1

필터 반복기를 만들 때 잘못된 시작 및 중지 인덱스를 사용합니다. 작업 예제를 첨부했습니다. 나는 병렬 for 문 안에서 스레드를 선택하기 위해 s를 사용했다. 버전이 제대로 작동하는지 확신 할 수 없으므로 이전 버전으로 변경했습니다.

#pragma omp parallel 
{ 
    #pragma omp for 
    for (int s = 0; s < 4; s++) 
    { 
     /* The following lines are important */ 
     FilterIter filter_iter_first(predicate, test_vec.begin() + start_vec[s], (test_vec.begin() + end_vec[s] + 1)); 
     FilterIter filter_iter_last(predicate, test_vec.begin() + end_vec[s] + 1, (test_vec.begin() + end_vec[s] + 1)); 

     for (auto iter = filter_iter_first; iter != filter_iter_last; iter++) 
     { 
      std::cout << "num: " << *iter << std::endl; 
      (*iter) = (*iter) + 1; 
     } 
    } 
}