2

複素数を含むバイナリ ファイルを読み取るコードを C で記述しました。それはうまくいきますが、私が演じる必要があるキャストに不快感を覚えます. より良い方法はありますか？私のプログラムでは速度が重要です (コードを C++ iostream から C stdio 関数に変更することで、実行時間を 2 倍にしました)。また、高速化することはできますか？

これが私のコードです：

#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 関数は、単一の複素数のバイト順を逆にします。それにはさらに不快ですが、質問が大きくなりすぎないようにします。

4

1 に答える 1

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」に置き換えます。

于 2011-02-06T23:01:07.600 に答える