0
调用树： 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周围的堆栈已损坏！
我已经看过这个永远不能明白为什么是这种情况。 请帮忙。谢谢。