Question

我有一个.txt，其中包含Fortran数据类型的数据。我可以通过以下方式使用Python读取数据。

import  fortranformat as ff

.....#more code in between.
with open(sys.argv[1]) as fh
.....#more code in between.  
    nodeline = ff.FortranRecordReader("(I3,I10,3E12.5)")
    line = fh.readline()
    data = nodeline.read(line)
#note that all the code is not provided.

我想知道是否有一种方法可以从.txt文件中从C ++中读取Fortran数据类型（不使用函数子字符串）。

PS。我没有提供完整的python代码，因为它可以正常工作，而我只是在使用所需的位来更好地解释我的问题。

.txt文件中的数据示例如下。

 -1         1-3.07500E+01-2.96893E+01-1.65000E+01
 -1         2-3.07500E+01 2.96893E+01-1.65000E+01
 -1         3-8.85000E+01 8.74393E+01-1.65000E+01
 -1         4-8.85000E+01-8.74393E+01-1.65000E+01
 -1         5-8.85000E+01 8.74393E+01-2.15000E+01
 -1         6-8.85000E+01-8.74393E+01-2.15000E+01
 -1         7-3.07500E+01 2.96893E+01-2.15000E+01
 -1         8-3.07500E+01-2.96893E+01-2.15000E+01
 -1         9 2.96893E+01-3.07500E+01-1.65000E+01
 -1        10-2.96893E+01-3.07500E+01-1.65000E+01
 -1        11-8.74393E+01-8.85000E+01-1.65000E+01

Answer 1

使用问题中提供的格式，我编写了一个简单的Fortran程序来写入一些数据。

      PROGRAM TEST
      WRITE(*,FMT='(I3,I10,3E12.5)') 1, 234, 5.67, 8.9, 0.123456789
      END PROGRAM

我将上述程序的输出通过管道传输到文件test.dat：

  1       234 0.56700E+01 0.89000E+01 0.12346E+00

然后在C ++中，可以使用std::ifstream轻松读取数据。

#include <fstream>
#include <iostream>

int main() {
    std::ifstream ifs("test.dat");

    int i,j;
    double d,e,f;

    while (ifs >> i >> j >> d >> e >> f) {
        std::cout << i << ' ' << j << ' ' << d << ' ' << e << ' ' << f << '\n';
    }
}

编译和运行输出

1 234 5.67 8.9 0.12346

答案答复：

如果必须在缺少空格的情况下解析这种奇怪的格式，则可能要考虑使用适当的解析器生成器，例如Boost.Spirit。

#include <fstream>
#include <iostream>
#include <tuple>
#include <vector>

#include <boost/fusion/adapted/std_tuple.hpp>
#include <boost/spirit/home/support/iterators/istream_iterator.hpp>
#include <boost/spirit/home/x3.hpp>

int main() {
    std::ifstream input("test.dat");
    input.unsetf(std::ios::skipws);
    std::vector<std::tuple<int, int, double, double, double>> entries;

    boost::spirit::istream_iterator first(input);
    boost::spirit::istream_iterator last;

    using namespace boost::spirit::x3;

    bool r = phrase_parse(first, last,
                          *(int_ >> int_ >> double_ >> double_ >> double_),
                          space, entries);

    if (!r || first != last) {
        std::cerr << "Parsing failed at " << std::string{first, last} << '\n';
    } else {
        for (auto const &entry : entries) {
            int i, j;
            double d, e, f;
            std::tie(i, j, d, e, f) = entry;
            std::cout << i << ' ' << j << ' ' << d << ' ' << e << ' ' << f
                      << '\n';
        }
    }
}

在C ++中读取Fortran数据类型

1 个答案: