/**********************************************

  Class for parsing ECI file (multicomponent)

  Written by k-fleak

 *********************************************/

#include <iostream>
#include <fstream>
#include <sstream>
#include <iterator>
#include <boost/lexical_cast.hpp>
#include <string>
#include <vector>
#include <map>

#ifndef __PARSEECI_CPP
#define __PARSEECI_CPP




class ParseECI{
private:
  std::ifstream input;
  std::map<int, std::map<std::vector<int>, double> > eciMap;
public:
  ParseECI(char* filename);
  ~ParseECI();
  std::map<int, std::map<std::vector<int>, double> > getECI();
};

ParseECI::ParseECI(char* filename){
  input.open(filename);
  std::string line;

  while(!input.eof()){
    std::getline(input,line);
    std::stringstream str;
    std::string stmp;
    str << line;
    std::istream_iterator<std::string> iss_iter (str);
    std::istream_iterator<std::string> iss_end;
    if ((*iss_iter).size() < 1){
      break;
    }
    int clusterIndex = boost::lexical_cast<int>(*iss_iter);
    std::advance(iss_iter, 1);
    std::vector<int> basis;
    double eci;
    while (iss_iter != iss_end){
      stmp = *iss_iter;
      for (int i = 0; i < stmp.size(); ++i){
        basis.push_back(boost::lexical_cast<int>(stmp[i]));
      }
      ++iss_iter;
      eci = boost::lexical_cast<double>(*iss_iter);
      ++iss_iter;
      eciMap[clusterIndex][basis] = eci;
      basis.clear();
    }
  }

  input.close();
}

ParseECI::~ParseECI(){
}

std::map<int, std::map<std::vector<int>, double> > ParseECI::getECI(){
  return eciMap;
}


#endif

/*
int main(){
  ParseECI pe("ECI");
  std::map<int, std::map<std::vector<int>, double> > eci = pe.getECI();
  std::map<int, std::map<std::vector<int>, double> >::iterator it;
  std::map<std::vector<int>, double>::iterator it2;
  it = eci.begin();
  while (it != eci.end()){
    int cluster = (*it).first;
    it2 = (*it).second.begin();
    std::cout << "cluster: " << cluster << std::endl;
    while (it2 != (*it).second.end()){
      std::vector<int> basis;
      double ecis;
      basis = (*it2).first;
      ecis  = (*it2).second;
      for (int i = 0; i < basis.size(); ++i){
        std::cout << basis[i] << " ";
      }
      std::cout << " :  " << ecis << std::endl;
      ++it2;
    }
    ++it;
  }
  
  
};
*/