cuda, pycuda - 복소수 작성법 - 오류 : 클래스 "cuComplex"에는 "i"회원이 없습니다

Question

나는 cuda, pycuda에서 복소수를 사용하는 데 어려움이 있습니다.

나는이 C에 있습니다

또한

#include <complex> 
typedef std::complex<double> cmplx; 
.... 
cmplx j(0.,1.); 

이 같은 코드 :

#include <boost/python.hpp> 
#include <boost/array.hpp> 
... 
typedef std::vector< boost::array<std::complex<double>,3 > > ComplexFieldType; 
typedef std::vector< boost::array<double,3> > RealFieldType; 
... 
__global__ void compute(RealFieldType const & Rs,ComplexFieldType const & M,..) 
... 

어떻게 내가이 pycuda와 함께 사용하기 위해 변환 할 수 있습니까? 난 ('예에 의해 쿠다'책에 따라)이 추천 STH 시도 :

struct cuComplex { 
    float real; 
    float imag; 
    cuComplex(float a,float b): real(a),imag(b){} 
    cuComplex operator *(const cuComplex& a) { 
    return cuComplex(real*a.real -imag*a.imag ,imag*a.real +real*a.imag); 
    } 
cuComplex operator +(const cuComplex& a) { 
    return cuComplex(real+a.real ,imag+a.imag); 
    }; 

cuComplex j(0.,1.); //instead of cmplx j(0.,1.); 

__global__ void compute(float *Rs,cuComplex * M,..) //instead of RealFieldType const & Rs,ComplexFieldType const & M 
.... 

제가 가지고 일부 오류

은 :

데이터 부재 초기화는 허용되지

이 선언에는 저장 클래스 또는 형식 지정자가 없습니다.

고마워요!

-------------------- -EDIT- --------------------- ------------------------

I했던 (위 상대적) #include <pycuda-complex.hpp>를 사용하여 다음

pycuda::complex<float> cmplx; 

cmplx j(0.,1.); 

등을 typedef std::vector< boost::array<std::complex<double>,3 > > ComplexFieldType;

및 ComplexFieldType const & M에 대해 글로벌 기능 내에서 나는 "float * M"또는 "cmplx * M"을 시도했습니다. 지금까지

, 내가 점점 오전 오류 :

변수 "cmplx"는 유형 이름

내가 pycuda :: 복잡한 cmplx 사용하는 경우 아니다; 다음, 내가 얻을 :

식별자 "cmplx은"다음에 정의되지 않은

이름입니다 "::"클래스 또는 네임 스페이스 이름 또한

해야합니다 :

표현식에 포인터 대 객체 유형이 있어야합니다 (코드의 다른 부분에서 온 것일 수도 있음)

Answer 1

실제로 무엇을하려고하는지 명확하게 알지 못하지만 (실제로 자신을 아는 사람이라면) 편집 및 댓글이 계속 올라감에 따라 점차 혼란스러워지고 있습니다.

#include <pycuda-complex.hpp> 

template<typename T> 
__global__ void kernel(const T * x, const T *y, T *z) 
{ 
    int tid = threadIdx.x + blockDim.x * blockIdx.x; 

    z[tid] = x[tid] + y[tid]; 
} 


typedef pycuda::complex<float> scmplx; 
typedef pycuda::complex<double> dcmplx; 

template void kernel<float>(const float *, const float *, float *); 
template void kernel<double>(const double *, const double *, double *); 
template void kernel<scmplx>(const scmplx *, const scmplx *, scmplx *); 
template void kernel<dcmplx>(const dcmplx *, const dcmplx *, dcmplx *); 

이 당신에게 사소한 커널 컴파일의 단일 및 이중 실수와 복소수 버전을 제공 :하지만, 여기에 제대로 pycuda 기본 복합 형식을 사용 CUDA 코드의 간단한, 컴파일 가능한 부분이며, 안드레아스의 대답은 조금 확장 이 같은 NVCC의 뭔가 :

$ nvcc -arch=sm_20 -Xptxas="-v" -I$HOME/pycuda-2011.1.2/src/cuda -c scmplx.cu 
ptxas info : Compiling entry function '_Z6kernelIN6pycuda7complexIdEEEvPKT_S5_PS3_' for 'sm_20' 
ptxas info : Used 12 registers, 44 bytes cmem[0], 168 bytes cmem[2], 4 bytes cmem[16] 
ptxas info : Compiling entry function '_Z6kernelIN6pycuda7complexIfEEEvPKT_S5_PS3_' for 'sm_20' 
ptxas info : Used 8 registers, 44 bytes cmem[0], 168 bytes cmem[2] 
ptxas info : Compiling entry function '_Z6kernelIdEvPKT_S2_PS0_' for 'sm_20' 
ptxas info : Used 8 registers, 44 bytes cmem[0], 168 bytes cmem[2] 
ptxas info : Compiling entry function '_Z6kernelIfEvPKT_S2_PS0_' for 'sm_20' 
ptxas info : Used 4 registers, 44 bytes cmem[0], 168 bytes cmem[2] 

아마도이 구조 정의에 따르면 귀하의 질문에 ....

Answer 2

사용하는의 된 STLport 버전에서 파생 된 사실 std::complex<>로

#include <pycuda-complex.hpp> 

{ 
    pycuda::complex<float> x(5, 17); 
} 

같은 인터페이스.