Question

我正在用C ++编写代码，用于高级消除教学方法。编译时没有错误，但在执行程序时出现堆栈粉碎错误。

这是我写的：

#include <iostream>
#include <stdlib.h>
#include <math.h>

using namespace std;

int main(){
    int i,j,l,n;
    double *c;
    int *indice;
    double a[2][2];

    n=2;
    a[1][1]=1;
    a[1][2]=2;
    a[2][1]=3;
    a[2][2]=4;

    c=new double[n];
    indice=new int[n];
    /*Inizialize indice*/
    for(i=0;i<n;i++){indice[i]=i;}
    /*find scaling factor*/
    for (i=0;i<n;i++){
        double c1=0;
        for (j=0;j<n;j++){
            double c0=abs(a[i][j]);
            if(c0>c1) c1=c0;}
        c[i]=c1;}
    /*find pivot*/
    int k=0;
    for(i=0;i<n-1;i++){
        double pi1=0;
        for(j=i;j<n;j++){
            double pi0=abs(a[indice[j]][i]);
            pi0/=c[indice[j]];
            if(pi0>pi1){
                pi1=pi0;
                k=j;}}
    /*interchange rows according to pivoting order*/
        int itemp=indice[i];
        indice[i]=indice[k];
        indice[k]=itemp;
        for(j=i+1;j<n;j++){
            double pj=a[indice[j]][i]/(a[indice[i]][i]);
    /*recording pivot ratio below diagonal*/
            a[indice[j]][i]=pj;
    /*modify other elements accordingly*/
            for(l=i+1;l<n;l++){
                a[indice[j]][l]-=pj*a[indice[i]][l];}}}

    delete c;
    delete indice;
    return 0;}

这是什么意思，我的错误在哪里？

Answer 1

您在此处有错误：

double a[2][2];

a[1][1]=1;
a[1][2]=2;
a[2][1]=3;
a[2][2]=4;

有效数组索引的范围是0..1，但是您正在访问索引2处的元素，这些元素超出范围。

你可能想要的是：

a[0][0]=1;
a[0][1]=2;
a[1][0]=3;
a[1][1]=4;

或更简洁：

double a[2][2] = { { 1, 2 }, { 3, 4 } };

Answer 2

正如在另一个答案中已经指出的那样，您的索引超出了范围，例如：在您的a矩阵中：

double a[2][2];

a[1][1]=1;
a[1][2]=2;
a[2][1]=3;
a[2][2]=4;

在这种情况下，有效索引为0和1，不 2（因为在C ++索引中从0开始）。

我想补充说，使用某些矩阵C ++ 类代替原始C风格数组可能会派上用场。您的矩阵类可以定义自定义索引访问方法（即使是operator()(int row, int column)的重载也可以），并且在 debug 构建中，您可以检查输入索引是否在有效范围内，并且抛出一个异常，甚至只是断言它们是否超出范围。通过这种方式，您可以更轻松地调试代码，因为您可以立即发现无效索引。

如果您发现索引检查明显减慢了生产代码的速度，您可以在发布版本中禁用索引检查，并将它们保留在调试版本中。

即使使用std::array也比原始C风格数组更好，因为您可以使用其at()方法通过边界检查来访问元素。（而且我也认为std::array重载operator[]也会在调试版本中检查一些边界。）

作为旁注，当您使用new[]动态分配数组内存时，应使用delete[]而不是delete发布它：

// Your code:
// delete c; 
// delete indice

// Fix:
delete[] c;
delete[] indice;

更好的是，如果您需要动态内存分配，请考虑使用std::vector而不是原始new[] / delete[]。

堆栈粉碎检测：高斯消除

2 个答案:

***堆栈粉碎检测***：高斯消除

2 个答案:

堆栈粉碎检测：高斯消除