CUDAを使用して、ベクトル(私の場合はint2)で構成される配列を使用してコンパイル時の最適化を試みていますが、これをクリーンな方法で実現することはできません。具体的には、2つの定数配列cとwを使用する問題に取り組んでいます。配列wはfloatで構成され、配列cはint2で構成されます。これらの配列は一定であるため、コンパイラにコンパイル時の最適化を実行させ、それによって配列アクセスを効果的に最適化してもらいたいと思います。たとえば、次の2つのデバイス関数の場合、コンパイラはループを展開し、配列アクセスをcとwの値に直接置き換えることで配列アクセスを最適化します。
__forceinline__ __device__ float someFunction1() {
const int2 c[9] = {make_int2(0, 0), make_int2(1, 0), make_int2(0, 1), make_int2(-1, 0), make_int2(0, -1),
make_int2(1, 1), make_int2(-1, 1), make_int2(-1, -1), make_int2(1, -1)};
const float w[9] = {4.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/36.0f, 1.0f/36.0f, 1.0f/36.0f, 1.0f/36.0f};
#pragma unroll
for (int i = 0; i < 9; ++i) {
//Do something here, accessing c[i] and w[i]
}
}
__forceinline__ __device__ float someFunction2() {
const int2 c[9] = {make_int2(0, 0), make_int2(1, 0), make_int2(0, 1), make_int2(-1, 0), make_int2(0, -1),
make_int2(1, 1), make_int2(-1, 1), make_int2(-1, -1), make_int2(1, -1)};
const float w[9] = {4.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/36.0f, 1.0f/36.0f, 1.0f/36.0f, 1.0f/36.0f};
#pragma unroll
for (int i = 0; i < 9; ++i) {
//Do something here, accessing c[i] and w[i]
}
}
さて、問題は、cとwを使用する各デバイス関数でcとwを継続的に宣言したくないということです。wをグローバルに宣言することはできますが、cをグローバルに宣言することはできません。CUDAではグローバル変数でmake_int2コンストラクターを呼び出すことができないためです。つまり、以下のプログラムは、「デバイス上で空でないコンストラクタまたはデストラクタのコードを生成できません」というエラーを出します。
//Declaring array c like this is not allowed
__device__ const int2 c[9] = {make_int2(0, 0), make_int2(1, 0), make_int2(0, 1), make_int2(-1, 0), make_int2(0, -1),
make_int2(1, 1), make_int2(-1, 1), make_int2(-1, -1), make_int2(1, -1)};
__device__ const float w[9] = {4.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/9.0f, 1.0f/36.0f, 1.0f/36.0f, 1.0f/36.0f, 1.0f/36.0f};
__forceinline__ __device__ float someFunction() {
#pragma unroll
for (int i = 0; i < 9; ++i) {
//Do something here, accessing c[i] and w[i]
}
}
私の質問は、これらの変数にアクセスする各関数でcとwの宣言を防ぎ、必要なコンパイル時の最適化を維持するにはどうすればよいですか?または、別の言い方をすれば、ベクトルの配列をグローバルに宣言するための回避策はありますか?
注意:cとwをグローバルメモリまたは__constant__メモリに格納できることは承知していますが、コンパイル時の最適化はできません。__constant__メモリも、不定期にアクセスすると問題になる可能性があります。