source: util/src/averager.cpp@ e34098

/*
 * averager.cpp
 *
 *  Created on: Apr 15, 2010
 *      Author: heber
 */

using namespace std;

#include <iostream>
#include <iomanip>
#include <fstream>
#include <sstream>
#include <string>


#include "average.hpp"
#include "version.h"

int main(int argc, char **argv) {
  MeanErrorMap *Results = NULL;
  IndexSet Columns;
  int col;
        int skiplines=0;
        string zeile;
        stringstream line;
        
        cout << ESPACKVersion << endl;

        if (argc <= 2) {
                cout << "Usage: " << argv[0] << " <column (zero-based)> <data file> [#skip initial lines]\n";
                return(1);
        }
        {
    line.str(argv[1]);
    line >> col;
    cout << "Using column " << col << "." << endl;
    Columns.insert(col);
    line.clear();
        }
        if (argc > 3) {
    line.str(argv[3]);
    line >> skiplines;
    cout << "Will skip first " << skiplines << " lines." << endl;
    line.clear();
        }
  ifstream test(argv[2]);
  if (test == NULL) {
    cout << "Can't open File " << argv[2] << "\n";
    return(255);
  } else {
    // skip initial lines
    for (;skiplines>0;skiplines--)
      getline(test, zeile, '\n');
    // call average function and print result
    if (!test.eof()) {
      Results = AverageColumns(test, Columns);
      cout << "Mean\tStd_deviation" << endl;
      for (MeanErrorMap::const_iterator Runner = Results->begin(); Runner != Results->end(); ++Runner)
        cout << Runner->second.first << "\t" << Runner->second.second << endl;
      delete(Results);
    }
  }

  return(0);
}