Параллельный STXXL вектор инициализации

Question

Ниже минимальный пример иллюстрирует поведение stxxl при инициализации контейнеров параллельно (с использованием OpenMP):

#include <omp.h> 
#include <stdio.h> 
#include <stxxl.h> 

typedef stxxl::VECTOR_GENERATOR<float>::result VEC_T; 

int main(int argc, char* argv[]) { 
    const unsigned long NUM = 8; 
    #pragma omp parallel num_threads(NUM) 
    {  
      VEC_T v; 
      printf("%d\t%p\n", omp_get_thread_num(), &v); 
    } 
    return EXIT_SUCCESS; 
} 

, работающие в любой

[STXXL-ERROR] File too large 

или

[SYSTEM-ERROR]Segmentation fault 

Как я могу выделить stxxl-контейнеры в нескольких потоках?

Answer 1

Инициализация контейнеров stxxl не является потокобезопасной, поэтому требуется взаимное исключение для потока, инициализирующего контейнер. Использование OpenMP, это будет следующим образом:

#include <omp.h> 
#include <stdio.h> 
#include <stxxl.h> 

typedef stxxl::VECTOR_GENERATOR<float>::result VEC_T; 

int main(int argc, char* argv[]) { 
    const unsigned long NUM = 8; 
    #pragma omp parallel num_threads(NUM) 
    { 
      VEC_T* v; 
      #pragma omp critical 
      { 
        v = new VEC_T(); 
      } 
      printf("%d\t%p\n", omp_get_thread_num(), &v); 
      delete v; 
    } 
    return EXIT_SUCCESS; 
}