这两天抽点时间研究matlab调用cuda程序的问题，，发现网上这方面的资料比较少，而且nvidia提供的白皮书描述的不够详尽，因此对于开发流程进行总结，希望大家有用。
总的来说，matlab中调用cuda程序有两种方法，第一是利用开勇的cuda wizard创建cuda程序的dll，然后结合matlab的mex编译器在matlab中进行调用；第二种是根据白皮书中所写，将cuda程写成标准的mex文件格式，然后利用mex编译器生成二进制代码，在matlab中进行调用。两种方法生成的都是以.mexw32或.mexw64结尾的二进制程序，直接在m文件中使用即可。
1. Matlab白皮书方法
A. 仅仅使用cufft，将cuda代码和C代码函数写在mex标准文件中，则利用下面的编译指令可以完成二进制程序的生成。
mex fft2_cuda.c -IC:/CUDA/include -LC:/CUDA/lib -lcudart -lcufft
B. 代码里面包含核函数，readme中说明使用如下指令可以完成，但是总是报错，说nvmex.pl有错误，可能是vs或者matlab版本不对？（我用的vs2005 sp1 and Matlab R2010）
nvmex -f nvmexopts.bat Szeta.cu -IC:/cuda/include -LC:/cuda/lib -lcufft -lcudart
C. 从网上下载了一个代码，非常不错，里面包含matlab编译脚本文件my_compile.m, cuda文件test.cu,mex标准文件test.cpp.这是一个非常好的例子，参照即可。首先运行my_compile test.cpp, 然后就会生成二进制程序文件以供matlab调用。
my_compile.m
///////////////////////
function my_compile(varargin)!”%VS80COMNTOOLS%vsvars32.bat” & nvcc -c -arch compute_13 test.cun=getenv(‘CUDA_LIB_PATH’); if n(1)=='”‘, n=n(2:end); end, if n(end)=='”‘, n=n(1:end-1);endmex([‘-L’ n],’-lcudart’,’test.obj’,varargin{:});
test.cu
//////////////////
extern “C” void gpuAdd(double* A, double* B, double *C);
__global__ void vecAdd(double* A, double* B, double* C){ int i = threadIdx.x; C[i]=A[i]+B[i];}
void gpuAdd(double* A, double* B, double *C){ double *Ad,*Bd,*Cd; cudaMalloc((void**)&Ad,5*sizeof(double)); cudaMalloc((void**)&Bd,5*sizeof(double)); cudaMalloc((void**)&Cd,5*sizeof(double)); cudaMemcpy(Ad,A,5*sizeof(double),cudaMemcpyHostToDevice); cudaMemcpy(Bd,B,5*sizeof(double),cudaMemcpyHostToDevice); vecAdd<<<1,5>>>(Ad,Bd,Cd); cudaMemcpy(C,Cd,5*sizeof(double),cudaMemcpyDeviceToHost); cudaFree(Ad);cudaFree(Bd);cudaFree(Cd);}
test.cpp
/////////////
#include “mex.h”
extern “C” void gpuAdd(double* A, double* B, double *C);
void mexFunction( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]){ double *C; if (mxGetNumberOfElements(prhs[0])!=5) mexErrMsgTxt(“Wrong number of elements in A!”); if (mxGetNumberOfElements(prhs[1])!=5) mexErrMsgTxt(“Wrong number of elements in B!”); C=(double*)mxGetData(plhs[0]=mxCreateDoubleMatrix(1,5,mxREAL)); gpuAdd((double*)mxGetData(prhs[0]),(double*)mxGetData(prhs[1]),C);}
2. 动态库的方法
首先按照开勇的dll开发方法，完成cuda程序的dll封装。然后创建一个包含mexFunction的标准mex文件，将调用的cuda函数写入其中即可。编译的方法类似：
mex -L. -ltest test.c
3. 直接把cuda程序编译成exe文件，在matlab中通过system函数来调用也可，（懒人办法）
真凉爽啊！青山绿水映入我的眼中，景色怡人啊！
