В настоящее время я делаю небольшую тестовую программу для простой проверки файлов. Программа записывает две небольшие матрицы (A и B) в файлы, закрывает и снова открывает их, считывает матрицы из файлов, умножает их и записывает полученную матрицу (C) в новый файл. Затем он закрывает и снова открывает этот файл, содержащий ответ, и распечатывает его для меня, чтобы проверить, правильно ли прошла операция ввода-вывода.
Моя проблема в том, что матрица результатов читается не так, как ожидалось.
Я считаю себя новичком в Си и в операциях ввода / вывода файлов, и этот код доставляет мне неприятности. Я использую WinXP, Codeblocks и Mingw.
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define bufferA(i,k) (bufferA[i*cols+k])
#define bufferB(k,j) (bufferB[k*cols+j])
#define bufferC(i,j) (bufferC[i*cols+j])
void printMatrix(int *nMatrixToPrint, int nNumberOfElements, int nDimension) {
// This function prints out the element of an Array. This array represents a matrix in memory.
int nIndex;
printf("\n");
for (nIndex = 0; nIndex < nNumberOfElements; nIndex++) {
if (nIndex % nDimension == 0)
printf("\n");
printf("%d,",nMatrixToPrint[nIndex]);
}
return;
}
int main(int argc, char *argv[]) {
int nElements = 16, nDim = 4;
int A[4][4] = {{1,2,3,1},{2,2,1,2},{4,2,3,1},{5,1,1,3}};
int B[4][4] = {{3,2,1,4},{2,2,3,3},{4,1,3,2},{2,2,5,1}};
// Create files of A and B, delete old ones if present
FILE *fpA = fopen("A.dat", "w+");
FILE *fpB = fopen("B.dat", "w+");
// Write data to them
fwrite((int*)A, sizeof(*A), nElements, fpA);
fwrite((int*)B, sizeof(*B), nElements, fpB);
// and close them
fclose(fpA);
fclose(fpB);
// Reopen files
fpA = fopen("A.dat", "r");
fpB = fopen("B.dat", "r");
// Allocate memory
int *bufferA = (int*)malloc(nElements * sizeof(*bufferA));
int *bufferB = (int*)malloc(nElements * sizeof(*bufferB));
int *bufferC = (int*)calloc(nElements, sizeof(*bufferC));
// Read files
fread(bufferA, sizeof(int), nElements, fpA);
fread(bufferB, sizeof(int), nElements, fpB);
printf("\nA");
printMatrix(bufferA, nElements, nDim);
printf("\n\nB");
printMatrix(bufferB, nElements, nDim);
// Matrix multiplication
// Calculate and write to C
int i,j,k = 0; // Loop indices
int n = nDim,l = nDim, m = nDim, cols = nDim;
// multiply
for (i = 0; i < n; i++) { // Columns
for (j = 0; j < m; j++) { // Rows
//C(i,j) = 0;
for (k = 0; k < l; k++) {
bufferC(i,j) += bufferA(i,k) * bufferB(k,j);
}
}
}
printf("\n\nC_buffer");
printMatrix(bufferC, nElements, nDim);
// Create C and write to it
FILE* Cfile = fopen("C.dat", "w");
fwrite(bufferC, sizeof(*bufferC), nElements, Cfile);
// Close files
fclose(fpA);
fclose(fpB);
fclose(Cfile);
// reopen C for reading
Cfile = fopen("C.dat", "r");
// Obtain file size
fseek(Cfile , 0 , SEEK_END);
long lSize = ftell(Cfile);
rewind(Cfile);
printf("\nC file length is: %ld", lSize);
// read data into bufferA
fread(bufferA, sizeof(int), lSize, Cfile);
fclose(Cfile);
printf("\n\nC_file");
printMatrix(bufferA, nElements, nDim);
// Free allocated memory and remove dangling pointers
free(bufferA); bufferA = NULL;
free(bufferB); bufferB = NULL;
free(bufferC); bufferC = NULL;
exit(0);
}
Что дает мне следующий вывод:
A
1,2,3,1
2,2,1,2,
4,2,3,1,
5,1,1,3
B
3,2,1,4
2,2,3,3,
4,1,3,2,
2,2,5,1
C_buffer
21,11,21,17
18,13,21,18
30,17,24,29
27,19,26,28
Длина файла C: 64
C_File
21,11,21,17
18,13,21,18
30,17,24,29
27,19,1,3
Как видите, последние два элемента в C_file неверны, вместо этого выходные данные показывают последние два элемента в A, когда я записывал содержимое файла в bufferA. Переключение на bufferB поменяет местами последние два символа с последними элементами в B, что все еще ошибочно. Копия файла в другой проект выдаст два последних целых числа, как и в ram по этому адресу malloc.
Мой вопрос таков: почему fwrite не записывает нужные данные в файл. Почему он управляет первыми 14 элементами, а не последними двумя? И как это отличается от моего предыдущего правильного использования fwrite и fread, когда я писал и получал элементы A и B?