***检测到堆栈崩溃***：高斯消去法 我在C++中编写了一个用于高斯消去的教学方法。编译时间没有错误，但在执行程序时出现了堆栈崩溃错误
以下是我写的：***检测到堆栈崩溃***：高斯消去法 我在C++中编写了一个用于高斯消去的教学方法。编译时间没有错误，但在执行程序时出现了堆栈崩溃错误,c++,runtime-error,C++,Runtime Error,以下是我写的： #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
#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;}
#包括
#包括
#包括
使用名称空间std；
int main（）{
int i，j，l，n；
双*c；
int*标记；
双a[2][2]；
n=2；
a[1][1]=1；
a[1][2]=2；
a[2][1]=3；
a[2][2]=4；
c=新的双精度[n]；
indice=新整数[n]；
/*印度化*/
对于（i=0；i您在此处有一个错误：
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 } };
正如在另一个答案中已经指出的，您的索引超出了范围，例如在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样式数组可能会很方便。您的矩阵类可以定义自定义索引访问方法（即使是代码<操作程序（）的超载（int行，int列）会这样做）。，并且，在调试版本中，您可以检查输入索引是否在有效范围内，如果超出范围，则引发异常，甚至只是断言。这样，调试代码变得更容易，因为您可以立即发现无效索引
如果您发现索引检查明显降低了生产代码的速度，那么可以在发布版本中禁用索引检查，而只将其保留在调试版本中
甚至使用也比原始的C样式数组好，因为您可以使用它的at（）
方法通过边界检查来访问元素。（我还认为std:：array
的重载操作符[]
也会在调试构建中进行一些边界检查。）
作为旁注，当您使用new[]
动态分配阵列内存时，您应该使用delete[]
释放它，而不是delete
：
// Your code:
// delete c;
// delete indice
// Fix:
delete[] c;
delete[] indice;
更好的是，如果需要动态内存分配，请考虑使用<代码> STD:：vector <代码>而不是原始代码> NeX[] <代码> >代码>删除[]/COD> .谢谢。最愚蠢的错误……我猜你来自MATLAB背景？C++中的数组索引从0开始，而不是从1开始。此外，如果使用std:：vector
而不是原始指针和数组，则可能更容易检测到错误（如果随后使用vector:：at（）
访问元素）。还存在错误的delete
（应该是delete[]
）使用<代码> STD:：数组和/或<代码> STD:：vector ，然后不必担心这些事情。请考虑使用<代码> STD:：数组< /代码>代替固定长度的C样式数组。前者同样快，但更安全。请注意<代码>（）的用法。new[]delete[]