cuComplex.hとexp（）

Question

Q0：

cuComplex.hはexp（）をサポートしていますか？

Q1：

A = Bの*のEXPを作成する方法（私はCを*）、Aは、B、Cは、実数の同じサイズの配列でありますか？これは正しいですか？メイン

：

cuComplex A; 
float B; 
cuComplex c; 

カーネル：

c[idx] = (0, C[idx]); 

A[idx] = B[idx] * exp(c[idx]); 

Q2：

cuComplex私は元の行列、次いで2回より多くのメモリを割り当てる必要が手段2つのフロート番号を含んでいます。純粋な虚数を作る方法はありますか？

Answer 1

cuComplex.hは、cuComplex（主にCUBLASおよびCUFFTライブラリ内で使用されるもの）で基本操作をいくつか提供していますが、指数関数はサポートされていません。

コンポーネント指向の算術を使用して指数関数を実装することができます。 cuComplexは、複素数の実数部をx成分に格納し、虚数部をy成分に格納します。所与の複素数z = X + I yを*指数は次のように計算することができる。

EXP（Z）= EXP（X）*（COS（Y）+ iが（Y）罪*）

これは、（未テスト）以下のCUDAコードにつながる：何

cuComplex my_complex_exp (cuComplex arg) 
{ 
    cuComplex res; 
    float s, c; 
    float e = expf(arg.x); 
    sincosf(arg.y, &s, &c); 
    res.x = c * e; 
    res.y = s * e; 
    return res; 
}