matlab mex并行.docxVIP

利用MEX文件在Matlab中创建并行应用# include mex.h# include math.h# include # define _MAX_THREAD 100/*判断线程是否执行完毕*/static volatile int WaitForThread_MAX_THREAD;/*判断各线程执行情况*/static volatile double percentage_MAX_THREAD;/*保存程序执行信息*/static int nThreadNum=0;static int isOpen=0;HANDLE * ThreadLIST=NULL;double * pOut=NULL;static int MM=0;static int NN=0;/*线程参数结构体*/typedef struct _ThreadArg int N;/输入阵列元素个数 double * pr1;/输入阵列1 double * pr2;/输入阵列2 double * pOut;/计算结果 int threadld;/线程编号threadArg;unsigned _stdcall threadFunc(threadArg * pArg) int id; int kk-0; int N=0; double * pr1; double * pr2; double * pOut; /分配数据 N=(int)(pArg-N/nThreadNum)+; nb=id*N; ne=(id+1)*N-; if(ne=pArg-N) ne=pArg-N-; for(kk=nb;kk=ne;kk+ + ) pOutkk=pr1kk+pr2kk;/完成计算 percentageid=(kk-nb)*/(ne-nb); WaitForThreadid=0; return 0 void mexFunction(int nlhs,mxArray * plhs) if(nrhs= =1)/关闭线程、查询状态和返回结果 int flag; flag=mxGetScalar(prhs0); if(flag=1)/关闭 int kk; int isOK; for(kk=0;kknThreadNum;kk+ +) isOK=TerminateThread(ThreadListkk,0); CloseHandle(Threadkk); if(isOK) mexPrintf(线程 %d 成功结束n,kk); else mexPrint(线程 %d 成功失败n,kk); if(pOut) free(pOut); pOut=NULL; isOpen=0;else if(flag= =2)/查询状态 int kk=0; for(kk=0;kknThreadNum;kk+ +) mexPrintf(线程 %d，是否等候：%s，状态%d%n, WaitForThreadkk?(是)：(否), (int)(precentagekk*100); else/获取结果 int flag=1; int kk; for(kk=0;kknThreadNum;kk+ +) flag=flag&WaitForThreadkk; if(!flag) plhs0=mxCreateDoubleMatrix(MM,NN,mxREAL); memcpy(mxGetPr(plhs0),pOut,MM*NN*sizeof(double); else plhs0=mxCreateDoubleMatrix(0,0,mxREAL); else/启动 threadArg * pArg; int kk; if(nrhs= =2) nThreadNum=4 nThreadNum=mxGetScalar(prhsnrhs-1); if(nThreadNum=0) return; ThreadList=(HANDLE* )malloc(nThreadNum* sizeof(HANDLE); pARG=(ThreadArg*)malloc(nTreadNum*sizeof(threadArg); if(pOut! =NULL) Free(pOut); pOut=NULL; /*初始化参数*/ pOut=(double*)malloc( sizeof(double)*mxGetM(prh0)*mxGetN(prhs0); for(kk=0;kknThreadNum;kk+ +) pArgkk.N=mxGetM(prhs0*mxGetN(prhs0); pArgkk.pr1=mxGetPr(prhs0); pArgkk.pr2=mxGetPr(prhs1); pArgkk.pOut=pOut; pArgkk.threadld=kk; WaitForThreadkk=1; percentagekk=0; /*保存输入数值列阵大小。假定输入前两个参数维数和各维大小相同*/ MM=mxGetM(prhs0); NN=mxGetN(prhs0); /*启动多线程*/ for(kk=0;kknThr