0

我在 MVC++ 2010 中用 C++ 编写了代码。在其中，程序迭代一维指针数组（双 *）的元素。但是，当我将输入（指针数组的大小）设置为非常大（例如 15000）并运行程序时，它会停止工作并显示一个窗口来关闭程序，因为它没有响应！有什么问题？

这是构建我正在谈论的数组的代码的一部分：

map<int, double *> CF;
CoefficientMap(CF);

double *T = new double[I * J];
for (int r = 1; r <= I * J; ++r)
    T[r] = 100;
SOR(T, CF, 1.8);

这是迭代器函数：

void SOR(double *T, map<int, double *> &CF, double w)
{
int iter = 0;
cout << "Stage 2: Solving the linear system of equations using SOR method... ";

const double tol = 0.00001;
double error = tol + 1;

double *TOld = new double[I * J];
for (int i = 1; i <= I * J; ++i)
    TOld[i] = 100;

while (abs(error) > tol)
{
    ++iter;
    for (int i = 1; i <= I * J; ++i)
        T[i] = (CF[i][0] + CF[i][1] * T[i + 1] + CF[i][2] * T[i + J] + CF[i][3] * T[i - J] + CF[i][4] * T[i - 1]) * w + (1 - w) * T[i];

    error = errorCalc(TOld, T, I * J);

    for (int i = 1; i <= I * J; ++i)
        TOld[i] = T[i];

    if (iter % 100 == 0)
    {
        cout << endl << endl;
        cout << "100 iterations done, please wait..." << endl << "Total accumulative error till this point: " << error << endl;
    }

    if (iter > 10000)
        return;
}

cout << "Done!" << endl << endl;
cout << "Converged after " << iter << " iterations!" << endl;
cout << "Final accumulative error: " << error << endl << endl;

}

现在，当 (I * J) 变得足够大（例如 15000）时，程序停止工作！

4

3 回答 3

3

最可能的解释是，您的堆栈空间不足。简单的解决方法是使数组静态或全局。您也可以new从堆中分配它。两者都将数组移出堆栈。

最好的可能是使用智能指针并将其放入堆中：

std::unique_ptr<double[]> arrayOfDoubles(new double[size]);

当智能指针变量超出范围时，这将负责释放内存，无需手动删除。

为了获得更好的答案，请编辑问题以包含代码...

您添加的代码至少在数组索引方面存在问题。索引从 0 开始，数组大小减一。正确循环：

double *T = new double[I * J];
for (int r = 0; r < I * J; ++r)
    T[r] = 100;

你在其他循环中也有同样的错误，同样的修复。

替代解决方案：如果您想从 1 开始索引（例如，因为您以这种方式编写了伪代码算法并且不想更改索引），最简单的方法是分配一个更大的数组而不使用索引 0：

double *T = new double[I * J + 1];

有了它，您可以使用当前的循环。

这种缓冲区溢出一个数组元素是令人讨厌的，因为通常在分配的内存块的末尾可能有未使用的空间，因此在您更改数组大小并且未使用的空间消失之前，错误可能会完全被忽视。即使溢出导致堆损坏，它也可能被忽视，直到您更改代码和损坏更改的效果。因此，例如，如果您不走运，添加调试代码可能会隐藏问题。

于 2013-04-20T18:43:50.353 回答
2

听起来您正在堆栈上分配一个普通数组，如下所示：

void f()
{
    double a[123456];
    ...
}

堆栈的大小是有限的 - 您应该分配new或（更好）使用 a std::vector。

于 2013-04-20T18:44:15.163 回答
2

您在堆栈上分配了太多空间，因此没有足够的内存来满足您的请求。作为替代方案，您可以为对象提供静态存储持续时间，或者您可以将其放在免费存储中new：

std::unique_ptr<int[]> ptr(new int[size]);
于 2013-04-20T18:53:52.353 回答