NUMA対応Cppコンテナ

Question

NUMA対応の方法でcpp標準コンテナAPIを使用する便利で効率的な方法はありますか？

私は、cpp環境でOpenMP並列スパース行列ベクトル乗算を行いたいと考えています。 NUMAドメインに関して、ベクトルと行列の値を割り当て、初期化する 、Cのコードが何らかの形で次のようになり：CPPで

size_t N = 1000000; 
double* vecVal = malloc (N*sizeof(double)); 

#pragma OMP parallel for 
for (size_i=0; i<N; ++i) 
{ 
    vecVal[i] = 0.; 
} 
/* do spMV */ 
delete vecVal; 

を私はスタンダードを使用したい::ベクトル（STD ::の配列固定サイズもOKです）。 std :: vector :: reserve（）はそのトリックを行いますか？ それはこのような何かをするために有効です。

私はその後のstd ::ベクトルに正しいサイズを設定するにはどうすればよい

std::vector<double> vec; 
vec.reserve(N); 
double *vecVal = vec.data(); 

#pragma OMP parallel for 
for (size_i=0; i<N; ++i) 
{ 
    vecVal[i] = 0.; 
} 
/* do spMV */ 

？

誰かがもっとエレガントな方法を知っていますか？

Answer 1

ここでは、特別なnuma対応のアロケータを使用する必要があります。 https://github.com/STEllAR-GROUP/hpx/blob/master/hpx/parallel/util/numa_allocator.hpp

基本的な考え方は、アロケータの割り当て関数の内部に最初のタッチを置くことです。 HPXエグゼキュータのものを#omp parallel for schedule(static)ループに置き換えれば問題ありません。