Arrays 为什么阵列在我需要足够的空间后会出现泄漏？

我想用mkl来计算一个方程。但阵列a[]似乎一直在泄漏，如下所示：

  44.62  -0.09 -6277438562204192487878988888393020692503707483087375482269988814848.00 -6277438562204192487878988888393020692503707483087375482269988814848.00 -6277438562204192487878988888393020692503707483087375482269988814848.00
  -0.09  11.29  -0.09 -6277438562204192487878988888393020692503707483087375482269988814848.00 -6277438562204192487878988888393020692503707483087375482269988814848.00
 -6277438562204192487878988888393020692503707483087375482269988814848.00  -0.09   0.18  -0.09 -6277438562204192487878988888393020692503707483087375482269988814848.00
 -6277438562204192487878988888393020692503707483087375482269988814848.00 -6277438562204192487878988888393020692503707483087375482269988814848.00  -0.09  11.29  -0.09
 -6277438562204192487878988888393020692503707483087375482269988814848.00 -6277438562204192487878988888393020692503707483087375482269988814848.00 -6277438562204192487878988888393020692503707483087375482269988814848.00  -0.09  44.62

我的代码是：

#include <stdlib.h>
#include <stdio.h>
#include "mkl_lapacke.h"

/* Auxiliary routines prototypes */
extern void my_print_matrix(char* desc, MKL_INT m, MKL_INT n, double* a, MKL_INT lda, FILE* fpWrite);
extern void print_matrix(char* desc, MKL_INT m, MKL_INT n, double* a, MKL_INT lda);

/* Parameters */
#define N 5//nstep
#define LDA N
#define RMIN -10.0
#define RMAX 10.0
/* Main program */
int main() {
    /* Locals */
    MKL_INT n = N, lda = LDA, info;

    /* Local arrays */
    double h = (RMAX - RMIN) / (double(N) + 1.0);;
    double xi;
    double *w;
    double *a;

    w= (double*)malloc(sizeof(double) * N);
    a = (double*)malloc(sizeof(double) * N*LDA);

    for (int i = 0; i < N; i++) {
        xi = RMIN + double(1.0+i) * h;
        a[i*(N+1)] = 2.0 / h / h+xi * xi;
        if (i==0) {
            a[1] = -1.0 / h / h;
        }
        else if (i == N - 1) {
            a[LDA * N-2] =- 1.0 / h / h;
        }
        else {
            a[i *(N + 1)+1] = -1.0/h/h;
            a[i * (N + 1) - 1] = -1.0/h/h;
        }
    }
    print_matrix("Matrix", n, n, a, lda);
    /* Executable statements */
    printf("LAPACKE_dsyev (row-major, high-level) Example Program Results\n");
    /* Solve eigenproblem */
    info = LAPACKE_dsyev(LAPACK_ROW_MAJOR, 'V', 'U', n, a, lda, w);
    /* Check for convergence */
    if (info > 0) {
        printf("The algorithm failed to compute eigenvalues.\n");
        exit(1);
    }
    exit(0);
} /* End of LAPACKE_dsyev Example */

/* Auxiliary routine: printing a matrix */
void print_matrix(char* desc, MKL_INT m, MKL_INT n, double* a, MKL_INT lda) {
    MKL_INT i, j;
    printf("\n %s\n", desc);
    for (i = 0; i < m; i++) {
        for (j = 0; j < n; j++) printf(" %6.2f", a[i * lda + j]);
        printf("\n");
    }
}

#包括
#包括
#包括“mkl_lapacke.h”
/*辅助例程原型*/
外部无效my_print_矩阵（char*desc、MKL_INT m、MKL_INT n、double*a、MKL_INT lda、FILE*fpWrite）；
外部无效打印矩阵（char*desc、MKL\u INT m、MKL\u INT n、double*a、MKL\u INT lda）；
/*参数*/
#定义N 5//N步骤
#定义LDA N
#定义RMIN-10.0
#定义RMAX 10.0
/*主程序*/
int main（）{
/*当地人*/
MKL_INT n=n，lda=lda，info；
/*本地阵列*/
双h=（RMAX-RMIN）/（双（N）+1.0）；；
双席；
双*w；
双*a；
w=（双*）malloc（sizeof（双）*N）；
a=（双*）malloc（双*N*LDA）；
对于（int i=0；i0）{
printf（“算法无法计算特征值。\n”）；
出口（1）；
}
出口（0）；
}/*LAPACKE_dsyev示例结束*/
/*辅助例程：打印矩阵*/
无效打印矩阵（char*desc、MKL\u INT m、MKL\u INT n、double*a、MKL\u INT lda）{
MKL_INT i，j；
printf（“\n%s\n”，desc）；
对于（i=0；i

前两个元素看起来是对的，但数字是错的。。。最后的数字是一个相同但非常大的数字。我认为是阵列泄漏，但我不知道如何处理。所以我请求帮助

---

a[]显示只是初始化，而不是结果。我的问题是初始化错误在哪里？

让我们仔细看看如何初始化

a

：

a[i*(N+1)] = 2.0 / h / h+xi * xi;
if (i==0) {
    a[1] = -1.0 / h / h;
}
else if (i == N - 1) {
    a[LDA * N-2] =- 1.0 / h / h;
}
else {
    a[i *(N + 1)+1] = -1.0/h/h;
    a[i * (N + 1) - 1] = -1.0/h/h;
}

让我们以

i==0

和

i==1

为例：

i==0

在这里，您首先执行无条件初始化

a[i*(N+1)] = 2.0 / h / h+xi * xi;

如果我们计算指数

i*（N+1）

它是

0*（N+1）

，即

0

。因此，您将初始化

a[0]

然后有

if（i==0）

在这里初始化

a[1]

i==1

首先是无条件初始化索引

i*（N+1）

，然后是

1*（50+1）

，它等于

51

。因此，在这里初始化

a[51]

然后，条件

i==1

和

i==N-1

都是假的，因此我们最终进入最后的

else

子句：

a[i *(N + 1)+1] = -1.0/h/h;
a[i * (N + 1) - 1] = -1.0/h/h;

第一个索引

i*（N+1）+1

将是

1*（50+1）+1

，即

52

。因此，您可以初始化

a[52]

下一个索引

i*（N+1）-1

将是

1*（50+1）-1

，即

50

。因此，您可以初始化

a[50]

这种模式在整个循环中重复，索引越来越高，但从不降低

这意味着您永远不会将索引

2

初始化为

49

。这些元素将具有不确定的值，如果您运气不好，其中一个值可能是陷阱值，这将导致在使用它们时出现未定义的行为

---

关于内存泄漏，您需要重新编写算法来初始化数组的所有元素。

a

请不要忘记通过调用free（a）和free（w）函数来释放所有已分配给*w和*数组的内存。

在一个不相关的注释中，

1.0+i

的结果已经是一个

double

，你不需要演员。请假设我是盲人，我的屏幕阅读器无法解读你的照片。它是否显示文本的图片？

w=malloc（sizeof（*w）*N）对于您的问题，从一个较小的输入集开始，使用调试器在监视变量及其值的同时逐条检查代码语句。我还建议您将所有表达式的结果存储在临时变量中，以便更容易查看它们的值。然后，在执行步骤时，注意实际初始化的
a
元素。我敢打赌，你没有正确地初始化所有元素。我不喜欢<代码> A[i*（n+1）]＝2 /h /h +席*席；代码>