如何在C ++中读取FORTRAN格式的数字

时间:2013-08-29 20:34:03

标签: c++

我有一个包含数字的文件,如下所示:

1.000000-5 2.436700+0 2.530000-2 2.436700+0 5.000000-2 2.436700+0
1.000000+1 2.436700+0 1.000000+2 2.433800+0 1.000000+3 2.433800+0

我需要用C ++来读取数字

1.0E-5 2.4367E0 2.53E-2 2.4367E0 5.0E-2 2.4367E0
1.0E1 2.4367E0 1.0E2 2.4338E0 1.0E3 2.4338E0

挑战在于文件的编号没有E; E表示指数表示法。

我如何将这样的东西读成浮点数?它需要非常高效,因为我必须为每个文件读取数十万或数百万次这样的数字。

有关如何实现这一目标的任何建议?

4 个答案:

答案 0 :(得分:3)

我只需在每次 - 或+之前添加E,然后在其上使用atofstrtof

如果这对你来说不够快,请查看你最喜欢的atof实现(我找不到一个简单的搜索,但它不应该太难)并更改它以便它不会寻找E,而只是寻找 - 或+。

答案 1 :(得分:3)

你应该写一个解析器。一个简单的实现是:

#include <iostream>
#include <sstream>
#include <vector>

int main() {

    // Signed data set
    std::istringstream input(
        "+1.000000-5 -2.436700+0 +2.530000-2 -2.436700+0 +5.000000-2 -2.436700+0\n"
        "+1.000000+1 -2.436700+0 +1.000000+2 -2.433800+0 +1.000000+3 -2.433800+0");

    std::vector<double> result;

    std::string in;
    while(input >> in) {
        auto e = in.find_last_of("+-");
        if(0 < e && e != std::string::npos) {
            in.replace(e, 1, std::string("E") + in[e]); // this might get a optimization
        }
        std::istringstream number(in);
        double d;
        number >> d;
        result.push_back(d);
    }
    for(auto d: result) std::cout << std::fixed << d << std::endl;
    return 0;
}

答案 2 :(得分:0)

将字符串扫描到+或 - (或空格)字符并存储。用null替换+或 - 并将字符串传递给atof以获取第一部分。然后,如果你在指数中进行了+或a - 扫描。一旦你有了int或double,使用value * pow(10.0, exponent)。或者你可以使用strtod并避免空值修改。

或如上所述,修改现有的atof实现。

答案 3 :(得分:0)

这是我刚刚快速入侵的内容,您需要测试它是否适用于您的所有情况:

ifstream file;
file.open("f:\\stackoverflow\\fortranfloat\\fortranfloats.txt");

string line;
if (file.is_open())
{
    while (file.good())
    {
        getline(file, line);
        for (int i = 0; i < line.length(); i++)
        {
            int j = 0;
            char buf[10];
            while (i < line.length() && line[i] != ' ')
            {
                if (line[i] == '-' || line[i] == '+')
                    buf[j++] = 'e';

                buf[j++] = line[i++];
            }
            float number = atof(buf);
            cout << "Number: " << number << endl;
        }
    }
    file.close();
}
else
{
    cout << "Failed to open file" << endl;
}