這是我的改寫，自己覺得很有道理拉 XD NR_END我覺得是為了記憶體宣告時不access overrun存在的 怎麼看我都覺得malloc時轉(size_t)很多餘，這樣做對於現在compiler應該很多餘吧 有錯請指指正，謝謝！ float **constraint_matrix(long nrl, long nrh, long ncl, long nch) /* allocate a float matrix with subscript range m[nrl..nrh][ncl..nch] */ { //利用記憶體是否初始的方式來執行range m[nrl..nrh][ncl..nch] //也就是說對於一個(nrh+1)*(nrh+1)的矩陣, 除了m[nrl..nrh][ncl..nch]部份,記憶 //體的存取為無效(access violation)，實際應用上只能以 m[rI][cI], rI<nrl //記憶體存取無效的方式來實作,至於cI部分卻只能自行判斷， //且為了保證range m[nrl..nrh][ncl..nch]在宣告記憶體時有效, //因此需要宣告m[nrl+nrh][ncl+nch]大小的連續記憶體空間,以避免 memory access //overrun const int rowStart = nrl; const int colStart = ncl; long nrow=nrh-nrl+1,ncol=nch-ncl+1; float **m; /* allocate pointers to rows */ m=(float **) malloc((size_t)((nrow+rowStart)*sizeof(float*))); if (!m) new Exception("allocation failure 1 in matrix()"); /* allocate rows and set pointers to them */ m[nrl]=(float *) malloc((size_t)((nrow*ncol+colStart)*sizeof(float))); if (!m[nrl]) new Exception("allocation failure 2 in matrix()"); //initial矩陣m中(row index>nrl)的記憶體起始位址 for(int rI=nrl+1;rI<=nrh;rI++) { m[rI]=m[rI-1]+ncol; memset(m[rI],-1, colStart*sizeof(float));//驗證使用 //理論上會將m[rI][0..colStart-1]設為0, //實際上是將m[rI][last..last-colStart] 設為0 } //承上, 所以若要將m[nrl..nrh][ncl..nch]設為0,應執行 for(int i=nrl; i<=nrh; i++) memset(m[i]+colStart, 0, ncol*sizeof(float)); /* return pointer to array of pointers to rows */ return m; } ※ 引述《realmeat (真肉)》之銘言： : ※ 引述《HeyScng ( )》之銘言： : : 剛剛在書中的程式看到的，不過目前還看不懂他在寫什麼， : : 大概跑了一下也是抓不到規律，有興趣的人可以看看! 謝謝！ : : : 跑了一下才知道它要幹嘛 : 難怪會取到軌異的位置去 : 簡單的說就是 : 取個矩陣只有指定的範圍才有宣告實體位置 : 其他則都沒宣告實體位置 : 寫法相當差, 很多無意義的程式碼 : 個人習慣, 全面性改寫 : float **matrix(long row_low, long row_high, long col_low, long col_high) : { : long i; : long row_sub = row_high-row_low+1; : long col_sub = col_high-col_low+1; : // 宣告放指標的空間 : float **m=(float **) malloc((size_t)((row_sub)*sizeof(float*))); : if (!m){ : //error ("allocation failure 1 in matrix()"); : return NULL; : } : // 宣告放float的空間 : float* m2=(float *) malloc((size_t)((row_sub*col_sub)*sizeof(float))); : if (!m2){ : //error ("allocation failure 2 in matrix()"); : return NULL; : } : m[0] = m2-col_low; // 做col 位移 : // 依序填入值 : for(i=1; i<row_sub; i++){ : m[i] = m[i-1]+col_sub; : } : return &m[-row_low]; //做row 位移 : } : 簡單的說就是宣告2次陣列的簡單程式 : 然後把位址移動一下欺騙程式就好了 : 寫那麼複雜是怎樣? --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 118.169.42.127 ※ 編輯: HeyScng 來自: 118.169.42.127 (11/05 21:13)