source: molecuilder/src/Parser/TremoloParser.cpp@ 2704e2

/*
 * TremoloParser.cpp
 *
 *  Created on: Mar 2, 2010
 *      Author: metzler
 */

#include "TremoloParser.hpp"
#include "World.hpp"
#include "atom.hpp"
#include "element.hpp"
#include "periodentafel.hpp"
#include <map>
#include <vector>

using namespace std;

/**
 * Constructor.
 */
TremoloParser::TremoloParser() {
  knownKeys["x"] = x;
  knownKeys["u"] = u;
  knownKeys["F"] = F;
  knownKeys["stress"] = stress;
  knownKeys["Id"] = Id;
  knownKeys["neighbors"] = neighbors;
  knownKeys["imprData"] = imprData;
  knownKeys["GroupMeasureTypeNo"] = GroupMeasureTypeNo;
  knownKeys["Type"] = Type;
  knownKeys["extType"] = extType;
  knownKeys["name"] = name;
  knownKeys["resName"] = resName;
  knownKeys["chainID"] = chainID;
  knownKeys["resSeq"] = resSeq;
  knownKeys["occupancy"] = occupancy;
  knownKeys["tempFactor"] = tempFactor;
  knownKeys["segID"] = segID;
  knownKeys["Charge"] = Charge;
  knownKeys["charge"] = charge;
  knownKeys["GrpTypeNo"] = GrpTypeNo;
}

/**
 * Destructor.
 */
TremoloParser::~TremoloParser() {
}

/**
 * Stores keys from the ATOMDATA line.
 *
 * \param line to parse the keys from
 * \param with which offset the keys begin within the line
 */
void TremoloParser::parseAtomDataKeysLine(string line, int offset) {
  string keyword;
  stringstream lineStream;

  lineStream << line.substr(offset);
  while (lineStream.good()) {
    lineStream >> keyword;
    usedFields.push_back(keyword);
  }
}

/**
 * Reads one data line of a tremolo file and interprets it according to the keys
 * obtained from the ATOMDATA line.
 *
 * \param line to parse as an atom
 */
void TremoloParser::readAtomDataLine(string line) {
  vector<string>::iterator it;
  stringstream lineStream;
  string word;

  lineStream << line;
  for (it=usedFields.begin(); it < usedFields.end(); it++) {
    cout << *it << " -- " << it->substr(0, it->find("=")) << " -- " << knownKeys[it->substr(0, it->find("="))] << endl;
    switch (knownKeys[it->substr(0, it->find("="))]) {
      case x :
        lineStream >> word;
        cout<< "Found an x: word: " << word << endl;
        break;
      default :
        lineStream >> word;
        cout << "Unknown key: " << *it << ", word: " << word << endl;
        break;
    }
  }
}


/**
 * Loads atoms from a tremolo-formatted file.
 *
 * \param tremolo file
 */
void TremoloParser::load(istream* file) {
  string line;
  string::size_type location;

  usedFields.clear();
  while (file->good()) {
    std::getline(*file, line, '\n');
    if (usedFields.empty()) {
      location = line.find("ATOMDATA", 0);
      if (location != string::npos) {
       parseAtomDataKeysLine(line, location + 8);
      }
    }
    if (line.length() > 0 && line.at(0) != '#') {
      readAtomDataLine(line);
    }
  }
}

/*
#ATOMDATA  <record_entry_1> ... <record_entry_n>
# <record_entry>: <dataname>[=<n>]
# <dataname>    : x | u | F | stress | Id | neighbors | imprData
#               | GroupMeasureTypeNo | Type | extType
#               | name | resName | chainID | resSeq
#               | occupancy | tempFactor | segID | Charge
#               | charge
ATOMDATA name Id x=3 mass charge epsilon sigma eps14 sig14 name type protein protno neighbors=4
*/

//MatrixContainer* data = readData(fileName, getHeaderSize('#'), 0);


/**
 * Saves the World's current state into as a tremolo file.
 *
 * \param file where to save the state
 */
void TremoloParser::save(ostream* file) {
/*
write header
for (each atom in world) {
  write coordinates and additional data
}
*/
}