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