使用fscanf()在C语言中读取科学格式数字的问题

时间:2012-06-28 20:34:07

标签: c linux ubuntu scanf

我有多个文本文件,数据如下:

    0.000000E+0 0.000000E+0
    1.800000E-1 0.000000E+0
    4.200000E-1 0.000000E+0
    6.950000E-1 0.000000E+0
    9.450000E-1 0.000000E+0
    1.225000E+0 0.000000E+0

由测试仪器生成。两列数据由制表符分隔。我正在尝试将此数据文件读入程序以进行进一步处理。我写的代码如下:

#include <stdio.h>
#include <stdlib.h>

float **make_array(int size); //function to make 2D array

int main(int argc, char **argv)
{
        FILE *infile;
        infile=fopen(argv[1], "r");

//count number of lines in file
        char c;
        int i=0, numlines=0;
        while((c=fgetc(infile))!=EOF)
                if(c == '\n') numlines++;

        printf("numlines: %d\n", numlines);

        float **entry = make_array(numlines);

//read inputfile
        for(i=0; i< numlines; i++)
                fscanf(infile, " %G\t%G", &entry[0][i], &entry[1][i]);
//verify read values
        printf("\n");
        for(i=0; i<numlines; i++)
                printf("%E\t%E\n", entry[0][i], entry[1][i]);

//clean up before exiting
        free(entry[0]);
        free(entry[1]);
        free(entry);
        fclose(infile);
return 0;
}


//create a 2D array of 2 by {size}
float **make_array(int size){
        int i;
        float **entry = malloc(2*sizeof(float *));
        if(entry==NULL) exit(1);
        entry[0] = malloc(size*sizeof(float *));
        for(i=0;i<size;i++) entry[0][i] == 0.0;
        entry[1] = malloc(size*sizeof(float *));
        for(i=0;i<size;i++) entry[1][i] == 0.0;

return entry;
}

当我编译并运行程序时,它会正确打印出文件中的行数,但是printf语句用于验证读入数组entry的数据只打印出所有零,如下所示:

numlines: 252

0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00
0.000000E+00    0.000000E+00

我为fscanf %g,%e,%f,%G etc尝试了不同的格式说明符,但它们都给出了相同的结果。另外,我正在使用UCCntu 12.04和GCC 4.6。任何指向错误的指针都将非常感激。

哈蒂嘉。

3 个答案:

答案 0 :(得分:4)

您的初始循环:

while((c=fgetc(infile))!=EOF)
    if(c == '\n') numlines++;

读取整个输入文件(以确定数组的大小)。完成后,您对fscanf的调用都会失败,因为您已经在文件的末尾。

假设输入文件是可搜索的(如果它是磁盘文件则应该是这样),您可以使用rewind()fseek()返回到开头并再次阅读。

顺便说一下,您使用%G%E格式说明符是有效的,但%g%e更常见(并且应该以完全相同的方式工作)。

在引用argc之前,您确实应该检查argv[1]的值;如果没有命令行参数调用,程序可能会爆炸。您还应检查fopen()是否成功,fscanf()返回2表示它已成功读取两个数据项。

答案 1 :(得分:1)

所以,问题在于:

while((c=fgetc(infile))!=EOF)
            if(c == '\n') numlines++;

执行while循环后,你的FILE *指向EOF,所以当你到达时:

for(i=0; i< numlines; i++)
            fscanf(infile, " %G\t%G", &entry[0][i], &entry[1][i]);

每个fscanf返回EOF,你根本不扫描值。

答案 2 :(得分:0)

Keith指出了 问题,但你应该考虑其他一些事情:

  • 检查fscanf的结果;它将返回成功转换和分配的次数。如果它回来少于2,那么你就不能正确读取数据;

  • 您在mallocentry[0]的{​​{1}}来电中使用了错误的数据类型;您应该使用entry[1],而不是sizeof (float)。您可以通过编写sizeof (float *)来清理内容 - 这样您就不会出现这种类型的不匹配。

  • 两次读取文件效率低下且不必要;由于您无论如何都在动态分配数组,因此在阅读文件时,可以使用entry[0] = malloc(sizeof *entry[0] * size);扩展entry[0]entry[1]的大小。