简体   繁体   English

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

[英]Reading Fortran datatypes in C++

I have a .txt containing data in Fortran datatypes. 我有一个.txt,其中包含Fortran数据类型的数据。 I am able to read the data using Python in the following manner. 我可以通过以下方式使用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.

I was wondering if there was a way to read Fortran datatypes from C++ (without using the function substring) from the .txt file. 我想知道是否有一种方法可以从.txt文件中从C ++中读取Fortran数据类型(不使用函数子字符串)。

PS. PS。 I have not provided the complete python code because it works properly and I am only using the bit required to explain my question better. 我没有提供完整的python代码,因为它可以正常工作,而我只是在使用所需的位来更好地解释我的问题。

A sample of the data in the .txt file is given below. .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

Using the format provided in the question, I wrote a simple Fortran program to write some data. 使用问题中提供的格式,我编写了一个简单的Fortran程序来写入一些数据。

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

I piped the output of the above program to a file test.dat : 我将上述程序的输出通过管道传输到文件test.dat

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

Then in C++ the data can easily be read in using std::ifstream . 然后在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';
    }
}

Compiling and running outputs 编译并运行输出

1 234 5.67 8.9 0.12346

Answer to the edit: 回答编辑:

If you have to parse such a weird format where spaces are missing, you might want to consider using a proper parser generator such as Boost.Spirit. 如果必须在缺少空格的情况下解析这种奇怪的格式,则可能要考虑使用适当的解析器生成器,例如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';
        }
    }
}

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM