1 #include "cs_mex.h"
2 /* cs_print: print the contents of a sparse matrix. */
mexFunction(int nargout,mxArray * pargout[],int nargin,const mxArray * pargin[])3 void mexFunction
4 (
5 int nargout,
6 mxArray *pargout [ ],
7 int nargin,
8 const mxArray *pargin [ ]
9 )
10 {
11 CS_INT brief ;
12 if (nargout > 0 || nargin < 1 || nargin > 2)
13 {
14 mexErrMsgTxt ("Usage: cs_print(A,brief)") ;
15 }
16 brief = (nargin < 2) ? 0 : mxGetScalar (pargin [1]) ; /* get brief */
17 if (mxIsComplex (pargin [0]))
18 {
19 #ifndef NCOMPLEX
20 cs_cl Amatrix, *A ;
21 A = cs_cl_mex_get_sparse (&Amatrix, 0, pargin [0]) ; /* get A */
22 cs_cl_print (A, brief) ; /* print A */
23 cs_free (A->x) ;
24 #else
25 mexErrMsgTxt ("complex matrices not supported") ;
26 #endif
27 }
28 else
29 {
30 cs_dl Amatrix, *A ;
31 A = cs_dl_mex_get_sparse (&Amatrix, 0, 1, pargin [0]) ;/* get A */
32 cs_print (A, brief) ; /* print A */
33 }
34 }
35