src/io/ZConsReader.cpp

#include "io/ZConsReader.hpp"
#include "utils/simError.h"
#include "utils/StringUtils.hpp"
namespace oopse {

ZConsReader::ZConsReader(SimInfo* info) : info_(info){

    std::string zconsFileName_ = getPrefix(info_->getFinalConfigFileName()) + ".fz";
    istream_.open(zconsFileName_.c_str());

    if (!istream_){
        std::cerr << "open " << zconsFileName_ << "error" << std::endl;
        exit(1);
    }

    Globals* simParam = info_->getSimParams();
    int nZconstraints = simParam->getNzConstraints();
    ZconStamp** stamp = simParam->getZconStamp();
    for (int i = 0; i < nZconstraints; i++){
        allZmols_.push_back(stamp[i]->getMolIndex());
    } 
}

ZConsReader::~ZConsReader(){
  istream_.close();
}

void ZConsReader::readNextFrame(){
    char line[MAXBUFFERSIZE];  
    int nFixedZmol;
    int sscanfCount;

    fixedZmolData_.clear();

    istream_.getline(line, MAXBUFFERSIZE);
    sscanfCount = sscanf(line, "%lf", &curTime_);
    if (sscanfCount != 1){
        std::cerr << "ZConsReader Error : reading file error" << std::endl;
        exit(1);
    }

    istream_.getline(line, MAXBUFFERSIZE);
    sscanfCount = sscanf(line, "%d", &nFixedZmol);
    if (sscanfCount != 1){
        std::cerr << "ZConsReader Error : reading file error" << std::endl;
        exit(1);
    }

    ZconsData data;
    for(int i = 0; i < nFixedZmol; i++){
        istream_.getline(line, MAXBUFFERSIZE);
        sscanfCount = sscanf(line, "%d\t%lf\t%lf\t%lf", &data.zmolIndex, &data.zforce, &data.zpos,&data.zconsPos);
        if (sscanfCount != 4){
          std::cerr << "ZConsReader Error : reading file error" << std::endl;
          exit(1);
        }

        fixedZmolData_.push_back(data);
    }

}

bool ZConsReader::hasNextFrame(){
  return istream_.peek() != EOF ? true : false;
}

}
Revision:	1911
Committed:	Mon Jan 10 18:05:45 2005 UTC (20 years, 10 months ago) by tim
File size:	1801 byte(s)
Log Message:	more work in zconstraint
#	User	Rev	Content
1	tim	1492	#include "io/ZConsReader.hpp"
2			#include "utils/simError.h"
3	tim	1911	#include "utils/StringUtils.hpp"
4	tim	1818	namespace oopse {
5	gezelter	1490
6	tim	1911	ZConsReader::ZConsReader(SimInfo* info) : info_(info){
7	gezelter	1490
8	tim	1911	std::string zconsFileName_ = getPrefix(info_->getFinalConfigFileName()) + ".fz";
9			istream_.open(zconsFileName_.c_str());
10	gezelter	1490
11	tim	1911	if (!istream_){
12			std::cerr << "open " << zconsFileName_ << "error" << std::endl;
13			exit(1);
14	gezelter	1490	}
15
16	tim	1911	Globals* simParam = info_->getSimParams();
17			int nZconstraints = simParam->getNzConstraints();
18			ZconStamp** stamp = simParam->getZconStamp();
19			for (int i = 0; i < nZconstraints; i++){
20			allZmols_.push_back(stamp[i]->getMolIndex());
21			}
22	gezelter	1490	}
23
24			ZConsReader::~ZConsReader(){
25	tim	1911	istream_.close();
26	gezelter	1490	}
27
28	tim	1911	void ZConsReader::readNextFrame(){
29			char line[MAXBUFFERSIZE];
30			int nFixedZmol;
31			int sscanfCount;
32	gezelter	1490
33	tim	1911	fixedZmolData_.clear();
34	gezelter	1490
35	tim	1911	istream_.getline(line, MAXBUFFERSIZE);
36			sscanfCount = sscanf(line, "%lf", &curTime_);
37			if (sscanfCount != 1){
38			std::cerr << "ZConsReader Error : reading file error" << std::endl;
39			exit(1);
40			}
41	gezelter	1490
42	tim	1911	istream_.getline(line, MAXBUFFERSIZE);
43			sscanfCount = sscanf(line, "%d", &nFixedZmol);
44			if (sscanfCount != 1){
45			std::cerr << "ZConsReader Error : reading file error" << std::endl;
46			exit(1);
47	gezelter	1490	}
48
49	tim	1911	ZconsData data;
50			for(int i = 0; i < nFixedZmol; i++){
51			istream_.getline(line, MAXBUFFERSIZE);
52			sscanfCount = sscanf(line, "%d\t%lf\t%lf\t%lf", &data.zmolIndex, &data.zforce, &data.zpos,&data.zconsPos);
53			if (sscanfCount != 4){
54			std::cerr << "ZConsReader Error : reading file error" << std::endl;
55			exit(1);
56			}
57	gezelter	1490
58	tim	1911	fixedZmolData_.push_back(data);
59	gezelter	1490	}
60
61			}
62
63			bool ZConsReader::hasNextFrame(){
64	tim	1911	return istream_.peek() != EOF ? true : false;
65	gezelter	1490	}
66
67	tim	1818	}