错误删除多维矩阵

Question

当用该析构函数删除多维矩阵：

matrix::~matrix(){ 
    int i; 
    for(i=0;i<n;i++){ 
     delete[] user_matrix[i];} 
    delete[] user_matrix;} 

我恢复这样的错误：

*** glibc detected *** ./a.out: free(): invalid pointer: 0x00007fdb33067778 *** 
======= Backtrace: ========= 
/lib/x86_64-linux-gnu/libc.so.6(+0x7eb96)[0x7fdb32d2db96] 
./a.out[0x40157c] 
./a.out[0x40172b] 
/lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed)[0x7fdb32cd076d] 
./a.out[0x400b09] 

随之而来的是“存储器映象”。当我通过Valgrind的运行它，它表明，有在操作过载为*错误：

matrix matrix::operator* (matrix param) { 
    //n is the size of the square matrix 
    if(n!=param.n){ 
    //returns an empty matrix if they are not of equal size 
    matrix blah; 
    return blah;} 
    //initiates a nxn matrix that is empty 
    matrix temp(n,0); 
    temp.user_matrix=matrix_mult(user_matrix,param.user_matrix); 
    return temp;} 

这符合析构函数被调用，也与析构函数注释掉运行，直到计算机内存用完或者计算足够小并结束。

#include"matrix.h" 
using namespace std; 

matrix::matrix(int n1,int initiate){ 
    srand(3534.34535); 
    n=n1; 
    //float** user_matrix 
    user_matrix=new float* [n]; 
    int i; 
    for(i=0;i<n;i++){ 
    user_matrix[i]=new float [n];} 

    if(initiate==1){ 

    int j; 
    for(i=0;i<n;i++){ 

    for(j=0;j<n;j++){ 

     cout<<"please ["<<i<<"]["<<j<<"]"<<endl; 

     cin>>user_matrix[j][i];} 
    } 
    }else if(initiate==2){ 
    user_matrix=random_matrix(n);} 

} 

float** matrix::inverse(){ 
    int i,k; 
    float sub_det,detin; 

    detin=det(); 

    if(detin==0){cout<<"uninvertable"<<endl;}; 

    float** inverse = new float* [n]; 
    for(i=0;i<n;i++){ 
    inverse[i]=new float [n];} 

    float invertdet=1.0/detin; 

    for(i=0;i<n;i++){ 
    for(k=0;k<n;k++){ 
     inverse[k][i]=invertdet*pow(-1,i+k)*determinant(sub_matrix(user_matrix,i,k,n),n-1); 
    } 
    } 

    return inverse;} 

void matrix::display(){ 
    //cout<<"lol"<<endl<<n<<endl; 
    int i,j; 
    cout.precision(5); 
    for(j=0;j<n;j++){ 
    cout<<"|"; 
    for(i=0;i<n;i++){ 

     cout<<user_matrix[i][j]<<" ";};cout<<"|"; 
    cout<<endl<<endl;} 
    cout<<endl<<endl; 
} 

float matrix::determinant(float** matrix,int n1){ 
    if(n1==1){return matrix[0][0];} 
    int i; 
    float det1=0; 
    i=0; 
    for(i=0;i<n1;i++){  
    float** temp_matrix=sub_matrix(matrix,i,0,n1); 
    det1 = det1 + pow(-1.0,i)*matrix[i][0]*determinant(temp_matrix,n1-1); 
    int j=0; 
    for(j=0;j<n1-1;j++){delete[] temp_matrix[j];} 
    delete[] temp_matrix; 
    } 
    return det1;} 

float matrix::det(){return determinant(user_matrix,n);} 




float** matrix::sub_matrix(float** matrix,int colum,int row,int n){ 
    float **sub_matrix=new float *[n-1]; 
    int iter; 

    for(iter=0;iter<(n-1);iter++){ 
    sub_matrix[iter]=new float [n-1];} 
    int iter2; 
    int placeholder1,placeholder2; 

    placeholder1=placeholder2=0; 

    for(iter=0;iter<n;iter++){ 

    if(iter==colum){continue;} 
    placeholder2=0; 
    for(iter2=0;iter2<n;iter2++){ 
     if(iter2==row){continue;} 
     sub_matrix[placeholder1][placeholder2]=matrix[iter][iter2]; 
     placeholder2++; 
    } 
    placeholder1++; 

    } 

    return sub_matrix;} 

float** matrix::random_matrix(int n){ 

    int i,j; 
    float** temp_mat=new float* [n]; 
    for(i=0;i<n;i++){ 
    temp_mat[i]=new float [n];} 

    for(i=0;i<n;i++){ 
    for(j=0;j<n;j++){ 
     temp_mat[i][j]=rand()%10 +1;} 
    } 
    return temp_mat;} 

float** matrix::matrix_mult(float** matrix1,float** matrix2){ 
    int i,j,k; 
    float subresult; 
    float** ret_mat; 
    ret_mat=new float* [n]; 
    for(i=0;i<n;i++){ 
    ret_mat[i]=new float [n];} 
    for(i=0;i<n;i++){ 
    for(j=0;j<n;j++){ 

     for(k=0;k<n;k++){ 
    subresult=subresult + matrix1[k][i]*matrix2[j][k]; 
     } 
     ret_mat[i][j]=subresult;} 
    } 
    return ret_mat;} 
matrix::~matrix(){ 
    int i;for(i=0;i<n;i++){delete[] user_matrix[i];};delete[] user_matrix;} 

matrix matrix::operator* (matrix param) { 
    if(n!=param.n){ 
    matrix blah; 
    return blah;} 

    matrix temp(n,0); 
    temp.user_matrix=matrix_mult(user_matrix,param.user_matrix); 
    return temp;} 

int main(){ 
    int i; 
    /*for(i=1;i<20;i++){ 
    matrix m1(i,2); 
    cout<<i<<" "<<m1.det()<<endl;}*/ 
    matrix m1(16,2),m2(16,2),m3(16,0); 
    for(i=0;i<100000;i++){m3=m1*m2;} 


return 0;} 

代码没有非错误导致的功能。

#include"matrix.h" 
using namespace std; 

matrix::matrix(int n1,int initiate){ 
    srand(3534.34535); 
    n=n1; 
    //float** user_matrix 
    user_matrix=new float* [n]; 
    int i; 
    for(i=0;i<n;i++){ 
    user_matrix[i]=new float [n];} 

    if(initiate==1){ 

    int j; 
    for(i=0;i<n;i++){ 

    for(j=0;j<n;j++){ 

     cout<<"please ["<<i<<"]["<<j<<"]"<<endl; 

     cin>>user_matrix[j][i];} 
    } 
    }else if(initiate==2){ 
    user_matrix=random_matrix(n);} 

} 

float** matrix::random_matrix(int n){ 

    int i,j; 
    float** temp_mat=new float* [n]; 
    for(i=0;i<n;i++){ 
    temp_mat[i]=new float [n];} 

    for(i=0;i<n;i++){ 
    for(j=0;j<n;j++){ 
     temp_mat[i][j]=rand()%10 +1;} 
    } 
    return temp_mat;} 

float** matrix::matrix_mult(float** matrix1,float** matrix2){ 
    int i,j,k; 
    float subresult; 
    float** ret_mat; 
    ret_mat=new float* [n]; 
    for(i=0;i<n;i++){ 
    ret_mat[i]=new float [n];} 
    for(i=0;i<n;i++){ 
    for(j=0;j<n;j++){ 

     for(k=0;k<n;k++){ 
    subresult=subresult + matrix1[k][i]*matrix2[j][k]; 
     } 
     ret_mat[i][j]=subresult;} 
    } 
    return ret_mat;} 
matrix::~matrix(){ 
    int i;for(i=0;i<n;i++){delete[] user_matrix[i];};delete[] user_matrix;} 

matrix matrix::operator* (matrix param) { 
    if(n!=param.n){ 
    matrix blah; 
    return blah;} 

    matrix temp(n,0); 
    temp.user_matrix=matrix_mult(user_matrix,param.user_matrix); 
    return temp;} 

matrix & matrix::operator= (const matrix & param) 
{ 
    int i,j; 
    float** new_array=new float* [param.n]; 
    for(i=0;i<param.n;i++){ 
    new_array[i]=new float [param.n];} 
    for(i=0;i<param.n;i++){ 
    for(j=0;j<param.n;j++){ 
     new_array[i][j]=param.user_matrix[i][j];} 
    } 

    for(i=0;i<param.n;i++){ 
    delete[] user_matrix[i];} 
    delete[] user_matrix; 
    user_matrix=new_array; 
    n=param.n; 
    return *this; 

} 

int main(){ 
    int i; 
    matrix m1(16,2),m2(16,2),m3(16,0); 
    for(i=0;i<100000;i++){m3=m1*m2;} 


return 0;} 

Answer 1

没有看到更多的代码，我只能猜测，你可能

• 没有一个正确的拷贝构造函数，所以user_matrix指针被复制，或者
• 默认构造函数不初始化成员user_matrix或
• matrix_mult返回一个无效指针。

但至少有一件事是错误的。 temp.user_matrix在temp.user_matrix =之前被删除...

Answer 2

当你有一个分配内存的对象时，只需要注意如何复制该对象。你需要定义一个拷贝构造函数和一个赋值运算符，它们根据对象分配的内存来做正确的事情。否则，在删除对象时会出错，通常是因为最终删除了两次相同的内存。

您似乎没有定义复制构造函数或赋值运算符。当你定义一个析构函数时，你几乎总是需要定义一个拷贝构造函数和一个赋值操作符。这被称为“三条规则”。你可以在这里看看如何做到这一点的指导（或者只是阅读一本好的C++书）。

What is The Rule of Three?