mxshowsparsematrix.c

这是机械工业张德丰《MATLAB与外部程序接口编程》随书源码

/* #include "mex.h"
#include "matrix.h"
void showSparseMatrix(mxArray *pa)
{
	if(pa==NULL)
	{
		return;
	}
	if(mxlsSparse(pa))
	{
		mwIndex *ir, *jc;
		double *pr, *pi;
		size_t M,N;
		int i,j;
		mwSize curColumNZnum=0;    //*每列的非零元素个数
		mwSize nznum=0;     //已知的非零元素个数
		mwSize nzmax=mxGetNzmat(pa);
		N=mxGetN(pa);
		M=mxGetM(pa);
		ir=mxGetir(pa);
		jc=mxGetjc(pa);
		pr=mxGetpr(pa);
		pi=mxGetpi(pa);
		for(i=0;i<N;i++)
		{
			curColumNZnum=jc[j+1]-jc[i];
			while(curColumNZnum>0)
			{
				if(pi==NULL)
				{
					mexPrintf("\t<%d,%d>=%6.2f\n",ir[nznum]+1,i+1,pr[nznum]);
				}
				else
				{
					mexPrintf("\t<%d,%d>=%6.2f+%6.2fi\n",ir[nznum]+1,i+1,pr[nznum],pi[nznum]);
				}
				curColumNZnum--;
			nznum++;
			}		
		}
	}
	else
	{
		mexPrintf("\t输入变量非稀疏矩阵!\n");
	}
}

void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[])
{
	int i=0;
	for(i=0;i<nrhs;i++)
	{
		mexPrintf("第%d个输入参数：\n",i++);
		showSparseMatrix(prhs[i]);
	}
}
文件内容提要*/
#include "mex.h"
#include "matrix.h"
void mexFunction(int nlhs,mxArray * plhs90,int nrhs,const mxArray * prhs[])
{
	int i,j,k;
	mwIndex index;
	double* pr=UNLL;
	double* pr=UNLL;
	size_t M,N;
	size_t ndim;
	mwSize dims[2];
	for(i=0,i<nrhs;i++)
	{
		if((mxIsDouble(prhs[i])) && (mxGetNumberOfDimensions(prhs[i]==2))
		{
			pr=mxGetPr(prhs[i]);
			pi=mxGetPi(prhs[i]);
			M=mxGetM(prhs[i]);
			N=mxGetN(prhs[i]);
			ndim=mxGetNumberOfDimensions(prhs[i]);

			mexPrintf("变量%d:\n",i)；
				for(j=0;j<M;j++)
				{
					dims[0]=j;
					for(k=0,k<n;k++)
					{
						dims[1]=k;
						index=mxCalcSingleSubscript(prhs[i],ndim,dims);
						if(pi==NULL)
						{
							mexPrintf("%6.2f,",pr[index]);
						}
						else
						{
							mexPrintf("%6.2f+%6.2fj,",pr[index],pi[index]);
						}
					}
					mexPrintf("\n");
				}
			}
			else
			{
			mexPrintf("输入%d个变量不是二维double数值阵列!\n",i);
			}
	}
}