Question

我是C的初学者，我遇到了这个问题我malloc/calloc二维数组，并在使用后尝试free时遇到错误。程序使用随机值填充方形矩阵，并允许用户在不同函数之间进行选择以将它们相乘，以便评估性能。这是代码：

#include <stdio.h>
#include <stdlib.h>
#include <sys/time.h>
#include <cblas.h>
#include <assignment2.h> // declares my non-blas multiplication functions

int main(int argc, const char *argv[]) {

   int n, i, j, blockSize, choice;

   printf("Enter n: ");
   scanf("%d", &n);
   printf("Enter block size: ");
   scanf("%d", &blockSize);

   do {
      printf("Enter method (1 -> naive, 2 -> ijk-blocked, " 
            "3 -> kij-blocked, 4 -> cblas_dgemm): ");
      scanf("%d", &choice);
   } while (choice < 1 || choice > 4);

   double **a, **b, **result;

   /* Init matrices*/
   a = malloc(n * sizeof(double*));
   b = malloc(n * sizeof(double*));
   result = malloc(n * sizeof(double*));
   if (a == NULL || b == NULL || result == NULL) {
      printf("Error.\n");
      return 1;
   }

   for (i = 0; i < n; i++) {
      *(a + i) = malloc(n* sizeof(double));
      *(b + i) = malloc(n* sizeof(double));
      *(result + i) = calloc(n, sizeof(double));
   }

   fillMatrix(n, a);
   fillMatrix(n, b);

   // timing
   struct timeval tv1, tv2;
   struct timezone tz;
   gettimeofday(&tv1, &tz);

   switch(choice) {
      case 1:
         printf("matmul_ijk\n");
         matmul_ijk(n, a, b, result);
         break;
      case 2:
         printf("matmul_ijk_blocked\n");
         matmul_ijk_blocked(n, blockSize, a, b, result);
         break;
      case 3:
         printf("matmul_kij_blocked\n");
         matmul_kij_blocked(n, blockSize, a, b, result);
         break;
      case 4:
         printf("cblas_dgemm\n");
         cblas_dgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans, n, n, n, 1.0,
               a[0], n, b[0], n, 0.0, result[0], n);
         break;
      default:
         puts("Error. Mode not recognized.");
   }
   gettimeofday(&tv2, &tz);

   // Print time
   double elapsed = (double) (tv2.tv_sec-tv1.tv_sec) 
      + (double) (tv2.tv_usec-tv1.tv_usec) * 1.e-6;
   printf("Time elapsed: %lf\n", elapsed);

   // for (i = 0; i < n; i++) {
      // free(a[i]);
      // free(b[i]);
      // free(result[i]);
   // }
   // free(a);
   // free(b);
   // free(result);
   return 0;
}

现在在最后，我评论了我尝试释放我之前分配的内存。因为如果我启用了对free的最后三次调用，我会收到排序错误

2(724,0x7fff7cc76960) malloc: *** error for object 0x7fe62a016000: pointer being freed was not allocated
*** set a breakpoint in malloc_error_break to debug

我不知道为什么，因为我自己分配了内存。但程序运行正常并提供结果，只有在释放错误时才会发生。有趣的是，当n很大时，只会发生。对于例如n=1000没有问题，但n=2000有。{1}}。出于速度原因，我总是选择BLAS例程来测试它，我不知道其他人是否显示相同的行为。

所以我猜错了，因为

，我错过了重要的事情

我记忆malloc似乎不想成为free d
问题似乎与矩阵的大小有关，我不知道怎么做。

任何人都可以指出我的错误吗？

Answer 1

我运行了以下代码而没有发生包含1000和2000的n值的事件：

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

int main(int argc, const char *argv[]) {

    int n, i;

    printf("Enter n: ");
    scanf("%d", &n);

    double **a, **b, **result;

    /* Init matrices*/
    a = malloc(n * sizeof(double*));
    b = malloc(n * sizeof(double*));
    result = malloc(n * sizeof(double*));
    if (a == NULL || b == NULL || result == NULL) {
        printf("Error.\n");
        return 1;
    }

    for (i = 0; i < n; i++) {
        *(a + i) = malloc(n* sizeof(double));
        *(b + i) = malloc(n* sizeof(double));
        *(result + i) = calloc(n, sizeof(double));
    }

    for (i = 0; i < n; i++) {
        free(a[i]);
        free(b[i]);
        free(result[i]);
    }
    free(a);
    free(b);
    free(result);

    return 0;
}

对于操纵a，b和result的函数，可能会发生一些事情。例如，如果要将两个二维数组相乘并将最终结果放入一个名为result的变量中，则乘法函数应该采用指向result的三重指针，例如：

/* allocate space for a, b and result... */
foo(a, b, &result);
/* do something with result... */

您可以仔细检查您的功能签名和代码。

如果您使用gcc，请使用-Wall警告集编译代码：

$ gcc -Wall foo.c -o foo.bin

也许可以考虑使用gdb调试器来确定崩溃发生的确切位置以及当时数据的值。我不是它的专家，我发现RMS's tutorial是一个有用的资源。

无法释放分配的内存

1 个答案: