并行 STXXL 向量初始化
Parallel STXXL Vector Initialization
以下最小示例说明了 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 容器?
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;
}
以下最小示例说明了 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 容器?
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;
}