刪除雙指針（矩陣）

Question

我正在解決一個量子機器問題，它需要我通過操縱一些矩陣來找到一些特徵值。這個問題的細節是不相關的，我只需要幫助C++的問題，我是這個語言的新手，幾個小時後，我發現任何嘗試自己解決這個問題都是徒勞的，所以我轉向你尋求幫助。

我有這個問題，其中glibc在我的程序結束時檢測到錯誤，我無法正確釋放，但在這裏複製粘貼太大，所以我只會複製實際出現錯誤的部分。

void hamiltonian(int, double **&); 

int i,j; 

int main() 
{ 
int N = 1000; double **A; 

hamiltonian(N, A); 

//Physics here 
. 
. 
. 
. 
. 
//Delete 
for(i=0; i<N; i++){delete []A[i];} 
delete []A; 

return 0; 
} 

void hamiltonian(int N, double **&A) 
{ 
A = new double *[N]; 
for(i=0; i<N; i++) 
{ 
A[i] = new double[N]; 
for(j=0; j<N; j++) 
{ 
if(i==j)A[i][j] = 2; 
if(i==j+1 || i==j-1){A[i][j] = 1;} 
} 
} 
} 

根據我的教授我有同樣的功能要解除我分配，但我壓根沒想到關於釋放被幾乎與我的項目完成後，所以我必須重新編寫大量的代碼，該問題是我不能在哈密爾頓函數內部釋放A，因爲我需要它在其他函數中（在物理內部）。

一定有辦法解決這個問題嗎？可能聽起來有點無知，但這聽起來像是一個效率較低的設計，如果我不得不按照我分配的相同函數進行釋放。

Answer 1

根據我的教授我有同樣的功能要解除我分配

這是純粹的愚蠢。有時（幾乎總是）您需要在函數外部使用分配的結構。對於對象來說肯定是錯誤的，因爲構造函數和析構函數是不同的函數。

任何方式，你可以逃脫不使用類，如果你犯了一個矩陣struct和相關newMatrix和deleteMatrix功能:)

#include <cstddef> 
#include <iostream> 

using namespace std; 

struct Matrix 
{ 
    int n; 
    int m; 
    double** v; 
}; 

Matrix newMatrix (int n, int m) 
{ 
    Matrix A; 
    A.n = n; 
    A.m = m; 
    A.v = new double*[n]; 
    for(int i = 0; i < n; i++){ 
     A.v[i] = new double[m]; 
    } 
    return A; 
} 

Matrix newHamiltonianMatrix (int n, int m) 
{ 
    Matrix A = newMatrix(n, m); 
    for(int i = 0; i < A.n; i++){ 
     for(int j = 0; j < A.m; j++){ 
      A.v[i][j] = 0.0; 
      if(i == j){ 
       A.v[i][j] = 2.0; 
      } 
      if(i == j + 1 or i == j - 1){ 
       A.v[i][j] = 1.0; 
      } 
     } 
    } 
    return A; 
} 

void deleteMatrix (Matrix A) 
{ 
    for(int i = 0; i < A.n; i++){ 
     delete [] A.v[i]; 
    } 
    delete [] A.v; 
    A.v = NULL; 
} 

int main() 
{ 
    Matrix A = newHamiltonianMatrix(10, 20); 
    for(int i = 0; i < A.n; i++){ 
     for(int j = 0; j < A.m; j++){ 
      cout << A.v[i][j] << " "; 
     } 
     cout << endl; 
    } 
    deleteMatrix(A); 
} 

Answer 2

delete A; 

需要是

delete[] A; 

如果你new[]，您必須delete[]它。此外，請使用vector - 他們會自己照顧自己。

vector<vector<double>> matrix; 

Answer 3

您的代碼有幾個問題。

（1）您不分配存儲器到的A指針成員。即A[i]未分配new[]。因此訪問它們是未定義的行爲。

（2）如果分配給new[]，您必須對指針做delete[]。在你的其他功能delete A;是錯誤的。使用delete[] A;

（3）使用new/new[]不是唯一的分配方式。事實上，你應該在沒有選擇的情況下使用這種動態分配。從你的代碼看來，你是硬編碼N=1000。所以最好使用2D數組。

const int N = 1000; // globally visible 
int main() 
{ 
    double A[N][N]; 
    ... 
} 
void hamiltonian (double (&A)[N][N]) 
{ 
    ... 
}