슈트라센 알고리즘을 작성 중인데, 알고리즘은 다 작성했는데, int **c로 넘겨받은 이차원 포인터 배열에 값이 대입이 안돼요.

읽기 엑세스 위반이라고 뜨는데 왜 뜨는지 모르겠어요

#define _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <string.h>
void read(int **A, int **B, int **C);
void strassen(int row, int **a, int **b, int **c);
void freeArr(int **a, int **b, int **c, int row);

int size;

void main() {
    int** a = NULL;
    int** b = NULL;
    int** c = NULL;

    read(a, b, c);
    strassen(size, a, b, c);

    int i, j;
    printf("Strassen 행렬곱의 결과 : \n");
    for (i = 0;i < 2;i++) {
        for (j = 0;j < 2;j++) {
            printf("%d ", c[i][j]);
        }
        printf("\n");
    }

    free(a, b, c, size);
}

void read(int **A, int **B, int **C) {
    FILE *fp;
    fp = fopen("matrix.txt", "r");

    char buffer[50];
    int dim = 0;
    fgets(buffer, sizeof(buffer), fp);      //buffer에 읽은 문자열 입력
    char *ptr = strtok(buffer, " ");
    dim = atoi(ptr); //행렬의 차원 수

    A = (int**)malloc(sizeof(int*)*dim);
    B = (int**)malloc(sizeof(int*)*dim);
    C = (int**)malloc(sizeof(int*)*dim);

    int i, j;
    for (i = 0;i < dim;i++) {
        A[i] = (int*)malloc(sizeof(int)*dim);
        B[i] = (int*)malloc(sizeof(int)*dim);
        C[i] = (int*)malloc(sizeof(int)*dim);

    }
    for (i = 0;i < dim;i++) {
        for (j = 0;j < dim;j++) {
            C[i][j] = 0;
        }
    }

    //행렬A 입력받음  
    printf("행렬 A:\n");
    fgets(buffer, sizeof(buffer), fp);
    ptr = strtok(buffer, " ");          //빈칸을 기준으로 끊는다
    for (i = 0;i < dim;i++) {
        for (j = 0;j < dim;j++) {
            A[i][j] = atoi(ptr);
            ptr = strtok(NULL, " ");
            printf("%d ", A[i][j]);
        }
        printf("\n");
    }

    //행렬B 입력받음
    printf("행렬 B:\n");
    fgets(buffer, sizeof(buffer), fp);
    ptr = strtok(buffer, " ");          //빈칸을 기준으로 끊는다
    for (i = 0;i < dim;i++) {
        for (j = 0;j < dim;j++) {
            B[i][j] = atoi(ptr);
            ptr = strtok(NULL, " ");
            printf("%d ", B[i][j]);
        }
        printf("\n");
    }

    fclose(fp);
    size = dim;
}

void strassen(int row, int **a, int **b, int **c) {
    int i, j, q;

    if (row <= 2) {
        int i, j, k;
        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++) {
                //c[i][j] = 0;
                *(*(c + i) + j) = 0;        ///////////////////////////////읽기 엑세스 위반////////////////////////////
                for (k = 0;k < row;k++) {
                    c[i][j] += a[i][k] * b[k][j];
                }
            }
        }
    }
    else {          //배열 동적할당
        int row = size / 2;
        int** p1 = (int**)malloc(sizeof(int*)*row);
        int** p2 = (int**)malloc(sizeof(int*)*row);
        int** p3 = (int**)malloc(sizeof(int*)*row);
        int** p4 = (int**)malloc(sizeof(int*)*row);
        int** p5 = (int**)malloc(sizeof(int*)*row);
        int** p6 = (int**)malloc(sizeof(int*)*row);
        int** p7 = (int**)malloc(sizeof(int*)*row);
        int** temp1 = (int**)malloc(sizeof(int*)*row);
        int** temp2 = (int**)malloc(sizeof(int*)*row);


        for (q = 0;q < row;q++) {       //2차원 배열 동적할당
            int row = (size) / 2;
            p1[q] = (int*)malloc(sizeof(int)*row);
            p2[q] = (int*)malloc(sizeof(int)*row);
            p3[q] = (int*)malloc(sizeof(int)*row);
            p4[q] = (int*)malloc(sizeof(int)*row);
            p5[q] = (int*)malloc(sizeof(int)*row);
            p6[q] = (int*)malloc(sizeof(int)*row);
            p7[q] = (int*)malloc(sizeof(int)*row);
            temp1[q] = (int*)malloc(sizeof(int)*row);
            temp2[q] = (int*)malloc(sizeof(int)*row);

        }

        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++) {
                temp1[i][j] = 0;
                temp1[i][j] = a[i][j] + a[i + row][j + row];
                temp2[i][j] = b[i][j] + b[i + row][j + row];
            }
        }
        strassen(row, temp1, temp2, p1);        //p1

        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++) {
                temp1[i][j] = a[i + row][j] + a[i + row][j + row];
                temp2[i][j] = b[i][j];

            }
        }
        strassen(row, temp1, temp2, p2);        //p2

        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++)
            {
                temp1[i][j] = a[i][j];  // a11
                temp2[i][j] = b[i][j + row] - b[i + row][j + row]; // b12 - b22
            }
        }
        strassen(row, temp1, temp2, p3);  // p3

        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++) {
                temp1[i][j] = a[i + row][j + row];
                temp2[i][j] = b[i + row][j] - b[i][j]; //b21-b11
            }
        }
        strassen(row, temp1, temp2, p4);    //p4

        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++) {
                temp1[i][j] = a[i][j] + a[i][j + row];
                temp2[i][j] = b[i + row][j + row];
            }
        }
        strassen(row, temp1, temp2, p5);    //p5

        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++) {
                temp1[i][j] = a[i + row][j] - a[i][j];
                temp2[i][j] = b[i][j] + b[i][j + row];

            }
        }
        strassen(row, temp1, temp2, p6);

        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++) {
                temp1[i][j] = a[i][j + row] = a[i + row][j + row];
                temp2[i][j] = b[i + row][j] + b[i + row][j + row];
            }
        }
        strassen(row, temp1, temp2, p7);

        for (i = 0;i < row;i++) {
            for (j = 0;j < row;j++) {
                c[i][j] = p1[i][j] + p4[i][j] - p5[i][j] + p7[i][j]; //c11
                c[i][j + row] = p3[i][j] + p5[i][j]; //c12
                c[i + row][j] = p2[i][j] + p4[i][j]; //c21
                c[i + row][j + row] = p1[i][j] + p3[i][j] - p2[i][j] + p6[i][j]; //c22
            }
        }

        for (j = 0;j < row;j++) // 동적 배열 해제
        {
            free(p7[j]);
            free(p6[j]);
            free(p5[j]);
            free(p4[j]);
            free(p3[j]);
            free(p2[j]);
            free(p1[j]);
            free(temp2[j]);
            free(temp1[j]);
        }

        free(temp2);
        free(temp1);
        free(p7);
        free(p6);
        free(p5);
        free(p4);
        free(p3);
        free(p2);
        free(p1);

    } //크기가 2보다 큰 경우 - else
}

void freeArr(int **a, int **b, int **c, int row) {
    for (int j = 0;j < row;j++) // 동적 배열 해제
    {
        free(a[j]);
        free(b[j]);
        free(c[j]);
    }
    free(a);
    free(b);
    free(c);
}