matlab mex并行.docxVIP

利用MEX文件在Matlab中创建并行应用# include mex.h# include math.h# include windows.h# define _MAX_THREAD 100/*判断线程是否执行完毕*/static volatile intWaitForThread[_MAX_THREAD];/*判断各线程执行情况*/static volatile double percentage[_MAX_THREAD];/*保存程序执行信息*/static intnThreadNum=0;static intisOpen=0;HANDLE * ThreadLIST=NULL;double * pOut=NULL;static int MM=0;static int NN=0;/*线程参数结构体*/typedefstruct _ThreadArg{int N;//输入阵列元素个数 double * pr1;//输入阵列1 double * pr2;//输入阵列2 double * pOut;//计算结果intthreadld;//线程编号}threadArg;unsigned __stdcallthreadFunc(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+ + ) {pOut[kk]=pr1[kk]+pr2[kk];//完成计算 percentage[id]=(kk-nb)*/(ne-nb); }WaitForThread[id]=0; return 0 }void mexFunction(intnlhs,mxArray * plhs[]){ if(nrhs= =1)//关闭线程、查询状态和返回结果 {int flag; flag=mxGetScalar(prhs[0]); if(flag==1)//关闭 {intkk;intisOK;for(kk=0;kknThreadNum;kk+ +) {isOK=TerminateThread(ThreadList[kk],0);CloseHandle(Thread[kk]); if(isOK) {mexPrintf(线程 %d 成功结束\n,kk); } else {mexPrint(线程 %d 成功失败\n,kk); } } } if(pOut) free(pOut);pOut=NULL;isOpen=0;}else if(flag= =2)//查询状态{intkk=0;for(kk=0;kknThreadNum;kk+ +) {mexPrintf(线程 %d，是否等候：%s，状态%d%%\n,WaitForThread[kk]?(是)：(否), (int)(precentage[kk]*100)); }} else//获取结果 {int flag=1;intkk;for(kk=0;kknThreadNum;kk+ +) { flag=flagWaitForThread[kk]; } if(!flag) {plhs[0]=mxCreateDoubleMatrix(MM,NN,mxREAL);memcpy(mxGetPr(plhs[0]),pOut,MM*NN*sizeof(double)); } else {plhs[0]=mxCreateDoubleMatrix(0,0,mxREAL); } }}else//启动{threadArg * pArg;intkk;if(nrhs= =2) {nThreadNum=4 }nThreadNum=mxGetScalar(prhs[nrhs-1]); if(nT