c - 割り当てられたメモリを解放できません

Question

私は C の初心者で、malloc/calloc2 次元配列を使用した後にそれらを使用しようとするとエラーが発生するという問題があり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;
}

最後に、以前に割り当てたメモリを解放しようとしたことをコメントアウトしました。の最後の 3 つの呼び出しを有効にすると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あります。速度の理由から、私は常にBLASこれをテストするルーチンを選択します。他のルーチンが同じ動作を示すかどうかはわかりません。

だから私はここで何か重要なものを見逃していると思います

1. メモリーmalloced by me したくないfreeらしい
2. 問題はマトリックスのサイズに関連しているように見えますが、その方法がわかりません。

誰かが私の間違いを指摘できますか?