Daniel@0: #include "mexutil.h"
Daniel@0: 
Daniel@0: /* Functions to create uninitialized arrays. */
Daniel@0: 
Daniel@0: mxArray *mxCreateNumericArrayE(int ndim, const int *dims, 
Daniel@0:          mxClassID class, mxComplexity ComplexFlag)
Daniel@0: {
Daniel@0:   mxArray *a;
Daniel@0:   int i, *dims1 = mxMalloc(ndim*sizeof(int));
Daniel@0:   size_t sz = 1;
Daniel@0:   for(i=0;i<ndim;i++) {
Daniel@0:     sz *= dims[i];
Daniel@0:     dims1[i] = 1;
Daniel@0:   }
Daniel@0:   a = mxCreateNumericArray(ndim,dims1,class,ComplexFlag);
Daniel@0:   sz *= mxGetElementSize(a);
Daniel@0:   mxSetDimensions(a, dims, ndim);
Daniel@0:   mxFree(dims1);
Daniel@0:   mxSetData(a, mxRealloc(mxGetData(a), sz));
Daniel@0:   if(ComplexFlag == mxCOMPLEX) {
Daniel@0:     mxSetPi(a, mxRealloc(mxGetPi(a),sz));
Daniel@0:   }
Daniel@0:   return a;
Daniel@0: }
Daniel@0: mxArray *mxCreateNumericMatrixE(int m, int n, mxClassID class, 
Daniel@0: 				mxComplexity ComplexFlag)
Daniel@0: {
Daniel@0:   size_t sz = m*n*sizeof(double);
Daniel@0:   mxArray *a = mxCreateNumericMatrix(1, 1, class, ComplexFlag);
Daniel@0:   mxSetM(a,m);
Daniel@0:   mxSetN(a,n);
Daniel@0:   mxSetPr(a, mxRealloc(mxGetPr(a),sz));
Daniel@0:   if(ComplexFlag == mxCOMPLEX) {
Daniel@0:     mxSetPi(a, mxRealloc(mxGetPi(a),sz));
Daniel@0:   }
Daniel@0:   return a;
Daniel@0: }
Daniel@0: mxArray *mxCreateDoubleMatrixE(int m, int n, 
Daniel@0: 			       mxComplexity ComplexFlag)
Daniel@0: {
Daniel@0:   return mxCreateNumericMatrixE(m,n,mxDOUBLE_CLASS,ComplexFlag);
Daniel@0: }
Daniel@0: