リターンダブルポインタの2次元マトリックス

Question

私は、二重のポインタで定義された2次元配列で使用e.g、

double** array; 
array = (double**) calloc(numRows, sizeof(double*)); 
for (int i = 0; i < numRows; i++) 
{ 
    array[i] = (double*) calloc(numCols, sizeof(double)); 
    /* then array[i][j] = values */ 
} 

// code to return matlab array 
plhs[0] = mxCreateDoubleMatrix(numRows, numCols, mxREAL); 

// memory copy 
// ????????? 

for (i = 0; i < numRows; i++){ 
     free(array[i]); 
    }  
free(array); 

ています私はMathWorks社のMATLABでarrayを返したいです。作業されていない、我々はこのmxSetPr(OUT, *stacked2D);ようmem-copy機能でそれを行うことができる場合、私は疑問に思っています

stacked1D = mxGetPr(plhs[0]); 
int n = 0; 
for (int r = 0; r < max_degree; r++) 
    for (int c = 0; c < n_vars; c++) 
     stacked1D[n++] = stacked2D[r][c]; 

：私は// memory copy部分が今まで持っていると私はそれは大丈夫だと思うの実行は、私はあるを修正してくださいこの構文では

ヒント説明や回答がありますか？

Answer 1

コード内で行と列の反復を逆にする必要があります。PaulMcKenzieが提案したことは、基本的には良い考えですが、Mexマトリックスでは機能しません（列ごとにレイアウトされています。不自然で混乱するM[column][row]であなたのマトリックスにアクセスしなければならないでしょう）。

また、次のような単純なラッパー使用することができます

template <class T, class I = unsigned> 
struct MatrixWrapper 
{ 
    T *data; 
    I nrows, ncols; 

    MatrixWrapper() 
     { clear(); } 
    MatrixWrapper(T *data_, I nrows_, I ncols_) 
     { set(data_,nrows_,ncols_); } 

    inline void clear() 
     { data = NULL; nrows = ncols = 0; } 

    inline void set(T *data_, I nrows_, I ncols_) 
     { data = data_; nrows = nrows_; ncols = ncols_; } 

    inline T& operator() (I r, I c) const 
     { return data[ r + nrows*c ]; } 
}; 

をして、MEX関数は次のようになります。

// allocate a temporary matrix 
double *M_data = new double[ nrows*ncols ]; 
MatrixWrapper<double> M(M_data, nrows, ncols); 

// access values with M(i,j) as a "normal" matrix 

// allocate a Mex output 
plhs[0] = mxCreateDoubleMatrix(nrows, ncols, mxREAL); 
MatrixWrapper<double> out0(mxGetPr(plhs[0]), nrows, ncols); 

// copy to the Mex output 
for (unsigned c = 0; c < ncols; c++) 
for (unsigned r = 0; r < nrows; r++) 
    out0(r,c) = M(r,c); 

// free temporary allocation 
M.clear(); 
delete[] M_data;