#define _CRT_SECURE_NO_WARNINGS
#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>
void inMatrix(int[][1024], char *fname);
void outMatrix(int[][1024], char *fname);
void Matrixadd(int[][1024], int[][1024], int[][1024]);
void Matrixsub(int[][1024], int[][1024], int[][1024]);
void Matrixmul(int[][1024], int[][1024], int[][1024]); //단순한 알고리즘으로 구하는 행렬의 곱셈 함수
int threshold(int n); //임계값을 구하는 함수
void Strassen(int n, int[][1024], int[][1024], int[][1024]); //쉬트라센 알고리즘으로 구하는 행렬의 곱셈 함수
int result1[1024][1024];
int result2[1024][1024];
int matrix1[1024][1024];
int matrix2[1024][1024];
int a11[1024][1024], a12[1024][1024], a21[1024][1024], a22[1024][1024];
int b11[1024][1024], b12[1024][1024], b21[1024][1024], b22[1024][1024];
int c11[1024][1024], c12[1024][1024], c21[1024][1024], c22[1024][1024];
int m1[1024][1024], m2[1024][1024], m3[1024][1024], m4[1024][1024], m5[1024][1024], m6[1024][1024], m7[1024][1024];
int tmp1[1024][1024], tmp2[1024][1024];
int main() {
int start, end;
char fname1[50], fname2[50];
printf("matrix1 행렬에 입력할 파일명을 입력하세요. : ");
scanf("%s", &fname1);
printf("matrix2 행렬에 입력할 파일명을 입력하세요. : ");
scanf("%s", &fname2);
inMatrix(matrix1, fname1);
inMatrix(matrix2, fname2);
start = clock();
Matrixmul(matrix1, matrix2, result1);
end = clock();
printf("Navie 방식의 결과 시간 : %f ms\n", (double)(end - start));
start = clock();
Strassen(1024, matrix1, matrix2, result2);
end = clock();
printf("Strassen 방식의 결과 시간 : %f ms\n", (double)(end - start));
outMatrix(result1,"Naive result.txt");
outMatrix(result2, "Strassen result.txt");
printf("결과가 Naive result.txt와 Strassen result.txt 파일에 저장되었습니다.\n");
}
void inMatrix(int matrix[][1024], char *fname) {
FILE *fin;
fin = fopen(fname, "rb"); //바이너리 파일로 파일은 연다
char *arr = NULL;
int idx = 0;
if (fin == NULL) {
printf("파일명이 없습니다.\n");
return;
}
for (int i = 0; i < 1024; i++)
for (int j = 0; j < 1024; j++) {
matrix[i][j] = arr[idx];
idx++;
}
fclose(fin);
free(arr);
}
void outMatrix(int matrix[][1024], char *fname) {
FILE *fout;
fout = fopen(fname, "wt");
for (int i = 0; i < 1024; i++) {
for (int j = 0; j < 1024; j++)
fprintf(fout, "%d ", matrix[i][j]);
fprintf(fout, "\n");
}
}
void Matrixadd(int matrix1[][1024], int matrix2[][1024], int result[][1024]) {
for (int i = 0; i < 1024; i++)
for (int j = 0; j < 1024; j++)
result[i][j] = matrix1[i][j] + matrix2[i][j]; //행렬을 더함
}
void Matrixsub(int matrix1[][1024], int matrix2[][1024], int result[][1024]) {
for (int i = 0; i < 1024; i++)
for (int j = 0; j < 1024; j++)
result[i][j] = matrix1[i][j] - matrix2[i][j]; //행렬을 뺌
}
void Matrixmul(int matrix1[][1024], int matrix2[][1024], int result[][1024]) {
for(int i=0;i<1024;i++)
for (int j = 0; j < 1024; j++) {
result[i][j] = 0;
for (int k = 0; k < 1024; k++)
result[i][j] += matrix1[i][k] * matrix2[k][j]; //행렬을 곱함
}
}
int threshold(int n) {
int threshold;
double t = floor(log(n) / log(2) - 6);
threshold = (int)floor(n / pow(2.0, t)) + 1;
return threshold;
}
void Strassen(int n, int matrix1[][1024], int matrix2[][1024], int result[][1024]) {
/*int a11[1024][1024], a12[1024][1024], a21[1024][1024], a22[1024][1024];
int b11[1024][1024], b12[1024][1024], b21[1024][1024], b22[1024][1024];
int c11[1024][1024], c12[1024][1024], c21[1024][1024], c22[1024][1024];
int m1[1024][1024], m2[1024][1024], m3[1024][1024], m4[1024][1024], m5[1024][1024], m6[1024][1024], m7[1024][1024];
int tmp1[1024][1024], tmp2[1024][1024];*/
if (n <= threshold(n)) //임계값을 넘지 않으면
Matrixmul(matrix1, matrix2, result); //일반적인 행렬곱셈으로 계산
else {
for (int i = 0; i < n / 2; i++) { //행렬의 부분을 나눠줌
for (int j = 0; j < n / 2; j++) {
a11[i][j] = matrix1[i][j];
a12[i][j] = matrix1[i][j + n / 2];
a21[i][j] = matrix1[i + n / 2][j];
a22[i][j] = matrix1[i + n / 2][j + n / 2];
b11[i][j] = matrix1[i][j];
b12[i][j] = matrix1[i][j + n / 2];
b21[i][j] = matrix1[i + n / 2][j];
b22[i][j] = matrix1[i + n / 2][j + n / 2];
}
}
//m1=(a11+a22)*(b11+b22)
Matrixadd(a11, a22, tmp1);
Matrixadd(b11, b22, tmp2);
Strassen(n / 2, tmp1, tmp2, m1);
//m2=(a21+a22)*b11
Matrixadd( a21, a22, tmp1);
Strassen(n / 2, tmp1, b11, m2);
//m3=a11*(b12+b22)
Matrixadd(b12, b22, tmp2);
Strassen(n / 2, a11, tmp2, m3);
//m4=a22*(b21+b11)
Matrixadd(b21, b11, tmp2);
Strassen(n / 2, a22, tmp2, m4);
//m5=(a11+a12)*b22
Matrixadd(a11, a12, tmp1);
Strassen(n / 2, tmp1, b22, m5);
//m6=(a21-a11)*(b11+b12)
Matrixsub(a21, a11, tmp1);
Matrixadd(b11, b12, tmp2);
Strassen(n / 2, tmp1, tmp2, m6);
//m7=(a12-a22)*(b21+b22)
Matrixsub(a12, a22, tmp1);
Matrixadd(b21, b22, tmp2);
Strassen(n / 2, tmp1, tmp2, m7);
//c11=m1+m4-m5+m7
Matrixadd(m1, m4, tmp1);
Matrixadd(m5, m7, tmp2);
Matrixsub(tmp1, tmp2, c11);
//c12=m3+m5
Matrixadd(m3, m5, c12);
//c21=m2+m4
Matrixadd(m2, m4, c21);
//c22=m1+m3-m2+m6
Matrixadd(m1, m3, tmp1);
Matrixadd(m2, m6, tmp2);
Matrixsub(tmp1, tmp2, c22);
for (int i = 0; i < n / 2; i++) {
for (int j = 0; j < n / 2; j++) {
result[i][j] = c11[i][j];
result[i][j + n / 2] = c12[i][j];
result[i + n / 2][j] = c21[i][j];
result[i + n / 2][j + n / 2] = c22[i][j];
}
}
}
}
실행을 시키면 파일 이름을 입력받는 부분까지만 실행이 됩니다. 디버깅을 시키면
'예외가 throw됨: 읽기 액세스 위반입니다. arr이(가) 0x1110112였습니다.'
이렇게 뜹니다. 동적할당에 문제가 있는 것 같은데 어떻게 수정을 해야될지 모르겠어요. 도와주세요.