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