#include <iostream>
#include <sstream>
#include <fstream>
#include <vector>
#include <cassert>

int main(void)
{
        // path of the input file
        std::string filename = "D:\\Skydrive\\MyCode\\cpp\\TreatString\\TreatString\\src\\test.txt"; // Please edit this line according to you real path of the test.txt file
        std::ifstream infile(filename);

        // check if file exists
        if(!infile) {
                std::cerr << "file "<< filename << " not found." << std::endl;
        }

        std::string line;  // string that stores each line read in
        std::vector<double> vals; // vector that stores values should be read

        // read file line by line
        while(std::getline(infile, line))
        {
                if(line == "Frame"){
                        continue;// if line is "Frame", go to next line
                }
                if(line == ""){
                        continue;// if line is blank, go to next line
                }
                double val1, val2, val3;
                std::istringstream iss(line);
                iss >> val1 >> val2 >> val3;
                vals.push_back(val1);
                vals.push_back(val2);
                vals.push_back(val3);
        }

        // check and make sure 27 values are stored in vector vals.
        assert(vals.size() == 27);

        double a1[3][3], a2[3][3], a3[3][3];

        // pass value from vals to a1, a2, a3
        for(size_t k = 0; k < 3; ++k)
        {
                for(size_t n = 0; n < 3; ++n)
                {
                        a1[k][n] = vals[3*k+n];
                        a2[k][n] = vals[9 + 3*k+n];
                        a3[k][n] = vals[18 + 3*k+n];
                }
        }

        // output a1, a2, a3
        std::cout << "a1:"<<std::endl;
        for(size_t k = 0; k < 3; ++k)
        {
                for(size_t n = 0; n < 3; ++n)
                        std::cout<< a1[k][n] << "\t";
                std::cout << std::endl;
        }

        std::cout << "a2:"<<std::endl;
        for(size_t k = 0; k < 3; ++k)
        {
                for(size_t n = 0; n < 3; ++n)
                        std::cout<< a2[k][n] << "\t";
                std::cout << std::endl;
        }

        std::cout << "a3:"<<std::endl;
        for(size_t k = 0; k < 3; ++k)
        {
                for(size_t n = 0; n < 3; ++n)
                        std::cout<< a3[k][n] << "\t";
                std::cout << std::endl;
        }

        return 0;
}

11楼: Originally posted by zeppe at 2013-06-03 11:58:21
可以试试我这个看看。


#include <iostream>
#include <sstream>
#include <fstream>
#include <vector>
#include <cassert>

int main(void)
{
        // path of the input ...

13楼: Originally posted by mljphy at 2013-06-03 17:54:42
谢谢你的回复。你的代码我仔细检验了。对于我原来的输入文件，在修改了路径那行：
std::string filename = "D:\\Skydrive\\MyCode\\cpp\\TreatString\\TreatString\\src\\test.txt"; //
改为
const c ...

12楼: Originally posted by zhanglinfeng at 2013-06-03 12:16:29
我机子上测试了是对的，只要filename是你文件所在的绝对路径，那就OK。

std::vector<double> vt;
bool bItem = false;
int nrow = 0;
while(std::getline(infile, line)
{
nrow += 1; //记录行号
// 其他代码
if (line == "ITEM: TIMESTEP")
{
     bItem = true; // turn on bItem
    nItem = nrow; //
}
if (bItem && nrow > nItem)
{
     if(line == "Atoms")
    {
          bItem = false;  // turn off bItem
    }
   else
   {
          std::istringstream iss1(line);
          double vab;
          iss1 >> vab;
         vt.push_back(vab);
   }
}
// 其他代码
}

16楼: Originally posted by zeppe at 2013-06-03 18:22:02
1.我发现如果输入文件(test.txt)中只加入第1，8，15，22行（对应ITEM: TIMESTEP），并在代码中加入下列代码行：      
if(line == "ITEM: TIMESTEP" { continue;  }
程序能够正常运行，并输出正确结果。 ...

#include <iostream>
#include <sstream>
#include <fstream>
#include <vector>
#include <cassert>

int main(void)
{
        // path of the input file ###### please edit this line according to your real path
        std::string filename = "D:\\Skydrive\\MyCode\\cpp\\TreatString\\TreatString\\src\\test_new.txt";
        std::ifstream infile(filename);

        // check if file exists
        if(!infile) {
                std::cerr << "file "<< filename << " not found." << std::endl;
        }

        // variables for the timesteps
        std::vector<double> timesteps;
        bool bTimeStep = false;
        int nTimeStep = 0;

        // variables for the atoms
        std::vector<int> atoms;  
        bool bAtoms = false;
        int nAtoms = 0;

        // variables for the coordinates
        std::vector<double> coordinates;
        bool bCoordinates = false;
        int nCoordinates = 0;

        // variables for the frames
        std::vector<double> frames;   
        bool bFrame = false;
        int nFrame = 0;

        std::string line;  // string that stores each line read in
        int nLine = 0;

        // read file line by line
        while(std::getline(infile, line))
        {
                nLine += 1;
                // read in TimeStep
                if(line == "ITEM: TIMESTEP"){
                        bTimeStep = true;
                        nTimeStep = nLine;
                }
                if(nLine > nTimeStep && bTimeStep){
                        if(line == "Atoms"){
                                bTimeStep = false;
                        }
                        else {
                                double val = std::numeric_limits<double>::quiet_NaN();
                                std::istringstream iss(line);
                                iss >> val;
                                timesteps.push_back(val);
                        }
                }
                // read in Atoms
                if(line == "Atoms"){
                        bAtoms = true;
                        nAtoms = nLine;
                }
                if(nLine > nAtoms && bAtoms){
                        if(line == "Coordinates"){
                                bAtoms = false;
                        }
                        else {
                                int val = 0;
                                std::istringstream iss(line);
                                iss >> val;
                                atoms.push_back(val);
                        }
                }
                //read in coordinates
                if(line == "Coordinates"){
                        bCoordinates = true;
                        nCoordinates = nLine;
                }
                if(nLine > nCoordinates && bCoordinates){
                        if(line == "Frame"){
                                bCoordinates = false;
                        }
                        else {
                                double val1, val2;
                                std::istringstream iss(line);
                                iss >> val1 >> val2;
                                coordinates.push_back(val1);
                                coordinates.push_back(val2);
                        }
                }
                //read in frame
                if(line == "Frame"){
                        bFrame = true;
                        nFrame = nLine;
                }
                if(nLine > nFrame && bFrame){
                        // ############ should be modified when file changes
                        if(line == "......共171行......"){
                                bFrame = false;
                        }
                        else if(line == "ITEM: TIMESTEP") {
                                bFrame = false;
                        }
                        else {
                                double val1, val2, val3, val4, val5;
                                std::istringstream iss(line);
                                iss >> val1 >> val2 >> val3 >> val4 >> val5;
                                frames.push_back(val1);
                                frames.push_back(val2);
                                frames.push_back(val3);
                                frames.push_back(val4);
                                frames.push_back(val5);
                        }
                }
        }

        // check if successfully read
        if(timesteps.empty()){
                std::cerr << "ERROR: failed to read timesteps." << std::endl;
        }
        if(atoms.empty()){
                std::cerr << "ERROR: failed to read atoms." << std::endl;
        }
        if(coordinates.empty()){
                std::cerr << "ERROR: failed to read coordinates." << std::endl;
        }
        if(frames.empty()){
                std::cerr << "ERROR: failed to read frames." << std::endl;
        }
        // ############ should be modified when file changes
        assert(timesteps.size() == atoms.size());
        assert(timesteps.size()*6 == coordinates.size());
        assert(timesteps.size()*5*4 == frames.size());

        //output values
        for(size_t i = 0; i < timesteps.size(); ++i)
        {
                std::cout << "TimeSteps" << std::endl;
                std::cout << timesteps[i] << std::endl;

                std::cout << "Atoms" << std::endl;
                std::cout << atoms[i] << std::endl;

                std::cout << "Coordinates" << std::endl;
                for(size_t j = 0; j < 3; ++j)
                {
                        for(size_t k = 0; k < 2; ++k)
                        {
                                std::cout<<coordinates[6*i + j*2 + k] << "\t";
                        }
                        std::cout << std::endl;
                }

                std::cout << "Frames" << std::endl;
                // ############# should be modified when file changes
                for(size_t m = 0; m < 4; ++m)
                {
                        for(size_t n = 0; n < 5; ++n)
                        {
                                std::cout << frames[20*i + m*5 + n] << "\t";
                        }
                        std::cout<< std::endl;
                }
        }

        return 0;
}

18楼: Originally posted by zeppe at 2013-06-04 11:04:19
楼主你好，请试用下面这些可以来实现读取17楼提供的文件，假设存为"test_new.txt"
在我电脑上是"D:\\Skydrive\\MyCode\\cpp\\TreatString\\TreatString\\src\\test_new.txt"
请根据你自己的情 ...