在内核中调用cublas

我想并行使用Zgemv。在内核中调用cublas

__global__ void S_Cphir(cuDoubleComplex *S,cuDoubleComplex *A,cuDoubleComplex *B, int n,int l) 
{ 
    .... 
cublasZgemv(handle,CUBLAS_OP_N,n,n,&alpha,S+i*n*n,n,A+n*i,1,&beta,B+i*n,1);} 

void S_Cphir_(cuDoubleComplex *S,cuDoubleComplex *A,cuDoubleComplex *B, int n,int l){ 
dim3 grid = dim3(1,1,1); 
dim3 block = dim3(32,1,1); 
S_Cphir<<<grid,block>>>(S,A,B,n,l);}

我的编译命令是

nvcc -c -arch=compute_30 -code=sm_35 time_propagation_cublas.cu --relocatable-device-code true 
nvcc -o ./main.v2 time_propagation_cublas.o -lcublas

第一行是工作。但第二行是错误的！

In function`__sti____cudaRegisterAll_58_tmpxft_000032b7_00000000_6_time_propagation_cublas_cpp1_ii_0d699356()';tmpxft_000032b7_00000000-3_time_propagation_cublas.cudafe1.cpp:(.text+0x17a4): 
undefined reference to `__cudaRegisterLinkedBinary_58_tmpxft_000032b7_00000000_6_time_propagation_cublas_cpp1_ii_0d699356' 
collect2: ld returned 1 exit status

我搜索“cudaRegisterLinkedBinary”，但我什么都没有！

我知道nvcc支持在内核中调用cublas。

来源

2013-10-19 周繼暉

使用CUBLAS Device Library sample code作为参考。在一个标准的CUDA 5.5安装，你会在找到它：

/usr/local/cuda/samples/7_CUDALibraries/simpleDevLibCUBLAS

参考该目录中的Makefile，你的编译命令应该是这样的：

nvcc -arch=sm_35 -rdc=true -o main.v2 time_propagation_cublas.cu -lcublas -lcublas_device -lcudadevrt

来源

2013-10-19 06:02:33

感谢。我可以编译。但我仍然不能在内核中并行Zgemv ...我想说的是我有32个矩阵向量乘法，我想并行执行。 –

我也发现这个[链接]（https://solarianprogrammer.com/2012/05/31/matrix-multiplication-cuda-cublas-curand-thrust/）有帮助的例子。 –

在内核中调用cublas

回答

相关问题