运行时检查失败#2 - 围绕bMatix堆栈被损坏
问题描述:
调用树: populateMatrix(_系数,bMatrix,行的cols); solveMatrix(coeff,bMatrix,cols + 1);运行时检查失败#2 - 围绕bMatix堆栈被损坏
void solveMatrix(float aMatrix[][DEFCOLS+1],float bMatrix[DEFCOLS+1], size_t cols){
std::cout << "\nInside solveMatrix...: " << endl;
size_t N2 = cols;
std::cout << "\N2 is...: " << N2 << endl;
for(size_t p=0; p<N2; p++){
//std::cout << "\nInside 1st for loop...: " << endl;
// find pivot row and swap
int max = p;
for(size_t i=p+1; i<N2; i++){
//std::cout << "\nInside 2nd for loop...: " << endl;
if (abs(aMatrix[i][p]) > abs(aMatrix[max][p])){
max = i;
}
}
//std::cout << "\nJust b4 all the swapping...: " << endl;
float temp[] = { *aMatrix[p] };
*aMatrix[p] = *aMatrix[max];
*aMatrix[max] = *temp;
float t = bMatrix[p];
bMatrix[p] = bMatrix[max];
bMatrix[max] = t;
//std::cout << "\nDone all the swapping...: " << endl;
if (abs(aMatrix[p][p]) <= MINISCULE) {
std::cout << "***** Error matrix value too small. Matrix is singular" << endl;
return;
}
// Pivot /in A and b
for(size_t i=p+1; i<N2; i++){
float alpha = aMatrix[i][p]/aMatrix[p][p];
bMatrix[i] = alpha * bMatrix[p];
for(size_t j=p; j<N2; j++){
aMatrix[i][j] -= alpha * aMatrix[p][j];
}
}
std::cout << "\nAbout to do the back subst..: " << endl;
// back subst.
std::vector<float> outMatrix(N2, 0.0);
int i =0;
for(i=N2-1; i>=0; i--){
float sum = 0.0;
int j=0;
for(j=i+1; j<N2; j++){
sum += aMatrix[i][j] * outMatrix[j];
}
if (aMatrix[i][i] > 0){
outMatrix[i] = (bMatrix[i] - sum)/aMatrix[i][i];
}else {
outMatrix[i] = 0.0;
}
}
int g = 0;
std::cout << "\nSolution: " << endl;
for (;g < N2-1; g++) {
cout << outMatrix[g] << " ";
}
cout << endl;
return;
}
} //结束solveMatrix()
这高斯消去具有部分枢转的实施方式。应用程序 运行并给出结果 - 不正确的重试,但在返回后,我得到下面的错误。
只要它返回到通话顶部,我就会得到: 运行时检查失败#2 - bMatrix周围的堆栈已损坏!
我已经看过这个永远不能明白为什么是这种情况。 请帮忙。谢谢。
答
问题发生是因为您传递给solveMatrix()(即cols + 1)的实际矩阵大小大于实际数组的大小(coeff或bMatrix)。
调试最简单的方法是:你的函数
- 彻底清除体内。确保您不会再收到错误消息。
- 恢复()循环的第一个。看看是否会导致错误出现。
- ...
- 重复,直到你找到它的功能的确切部分导致错误。再仔细查看一下,确保您不超过你的数组的限制条件(即,如果bMatrix为float [16],你不写bMatrix [16],[17]以上。