我在C中编写了一些代码来读取包含复数的二进制文件。它有效,但我对我需要表演的演员感到不舒服。有没有更好的办法?速度在我的程序中至关重要(通过将代码从C ++ iostreams更改为C stdio函数,我将执行时间增加了一倍)。还可以加快速度吗?
这是我的代码:
#include<complex.h>
#include<errno.h>
#define spaceh 6912
#define Nc 3
#define dirac 4
... ...
typedef double complex dcomplex;
long size;
size_t result;
char filename[84];
char* buffer;
dcomplex* zbuff;
int i, j, k, srccol, srcdir;
srcdir = 1;
srccol = 2;
/* allocate array dcomplex QM[dirac][Nc][space] on the heap */
sprintf(filename, "/<path>/file.%d.%d.bin", srcdir, srccol);
FILE* input;
input = fopen(filename, "rb");
if(readfile)
{
fseek(input, 0, SEEK_END);
size = ftell(input);
rewind(input);
buffer = (char*)malloc(sizeof(char)*size);
if(buffer == NULL)
{
fputs("Buffer allocation failed.", stderr);
exit(1);
}
result = fread(buffer, 1, size, input);
if(result != size)
{
fputs("File reading failed.", stderr);
exit(2);
}
/* The cast I'm referring to */
zbuff = (dcomplex*)buffer;
}
else
{
printf("File was not successfully opened: %s\n", strerror(errno));
}
count = 0;
for(k = 0; k < space; k++)
{
for(j = 0; j < Nc; j++)
{
for(i = 0; i < dirac; i++)
{
QM[i][j][k] = convert_complex(zbuff(count));
count++;
}
}
}
free(buffer);
fclose(input);
convert_complex函数反转单个复数的字节顺序。我对此更不舒服,但我不希望我的问题变得太大。
答案 0 :(得分:2)
直接声明zbuff,无需中间缓冲区。为此,您需要在适当的位置进行以下更改。在fread中,不是读取大小1,而是读取sizeof(dcomplex)。这应该做到。
//buffer = (char*)malloc(sizeof(char)*size);
zbuff = (dcomplex*)malloc(sizeof(dcomplex)*size);
if(zbuff == NULL)
{
fputs("Buffer allocation failed.", stderr);
exit(1);
}
result = fread(zbuff, sizeof(dcomplex), size, input);
if(result != size)
{
fputs("File reading failed.", stderr);
exit(2);
}
/* The cast I'm referring to */
//zbuff = (dcomplex*)buffer;
.......
free(zbuff);
将所有出现的'buffer'替换为'zbuff'。