HistoryCompare.cc

Go to the documentation of this file.

#include<vector>
#include<string>
#include<fstream>
#include<utils/CmdInterface.hh>
#include<automata/History.hh>

using namespace std;
using namespace Diades::Automata;

typedef enum { HELP =0 } Option;
typedef enum { ARCHV = 0 } FileExtension;
unsigned numberOfOptions = 1;
unsigned numberOfFileExtensions = 1;

vector<string> fileExtensions(numberOfFileExtensions);
vector<string> options(numberOfOptions);
vector<bool> isSet(numberOfOptions,false);

string description="Usage:\t history_compare --help\n\t history_compare hist1.archv hist2.archv\n\n\t Produce a set of data files with the evolution of performance indicators through the histories.";


void initialiseOptions()
{
  options[HELP]               = "--help";
  fileExtensions[ARCHV]   = "archv";
}

void  makeComparison(const History & history1,
                                         const History & history2, const string & prefix);

int main(int argc, char * argv[])
{
  setVerboseLevel(VbOutput);
  initialiseOptions();
  string histFile1;
  string histFile2;
  if(argc < 2)
        {
          printError("Wrong command line. Try option --help.");
        }
  int index = 1;
  while(index < argc)
        {
          Option currentOption;
          if(getOption<Option>(argv[index],options,currentOption))
                {
                  if(isSet[currentOption])
                        {
                          printError("The option '" + options[currentOption] + "' occurs at least twice.");
                        }
                  else
                        {
                          isSet[currentOption] = true;
                        }
                   switch(currentOption)
                        {
                        case HELP:
                          {
                                if((index!=1) || (argc !=2))
                                  {
                                        printError("Incorrect use of option --help.");
                                  }
                                else
                                  {
                                        printUsage(description);
                                  }
                                break;
                          }
                        default:
                          {
                                printError("Unrecognized option: " + string(argv[index]));
                                break;
                          }
                        }
                }
          else
                {
                  // then it should be the name of a file
                  vector<string>::const_iterator it =
                        getFileExtension(string(argv[index]),fileExtensions.begin(),fileExtensions.end());
                  if(it == fileExtensions.end())
                        {
                          printError("Unrecognized file extension in file name: " + string(argv[index]));
                        }
                  if(histFile1.empty())
                        {
                          histFile1 = string(argv[index]);
                        }
                  else
                        {
                          if(histFile2.empty())
                                {
                                  histFile2 = string(argv[index]);
                                }
                          else
                                {
                                  printError("Expecting two history files but get at least three.");
                                }
                        }
                  ++index;
                }
        }
  if(histFile1.empty() || histFile2.empty())
        {
          printError("Cannot get two history file from the command line.");
        }
  
  { verbose<VbOutput>("Loading the history %1%...") % histFile1; }
  History history1;
  {
        std::ifstream ifs(histFile1);
        boost::archive::text_iarchive ia(ifs);
        ia >> history1;
  }

  // cout << "HISTORY 1:" << endl;
  // cout << history1 << endl;

  
  { verbose<VbOutput>("Done.\n"); }
  { verbose<VbOutput>("Loading the history %1%...") % histFile2; }
  History history2;
  {
        std::ifstream ifs(histFile2);
        boost::archive::text_iarchive ia(ifs);
        ia >> history2;
  }

  Diagnosis diag1 = history1.getDiagnosis(history1.lastTimePoint());
  Diagnosis diag2 = history2.getDiagnosis(history2.lastTimePoint());
  
  
  if( (history1.lastTimePoint()-history1.initialTimePoint()) > (history2.lastTimePoint()-history2.initialTimePoint()) )
        {
          history2.publishDiagnosis(history2.initialTimePoint() + (history1.lastTimePoint()-history1.initialTimePoint()),history2.getDiagnosis(history2.lastTimePoint()));
        }
  else
        {
          if ((history2.lastTimePoint()-history2.initialTimePoint()) > (history1.lastTimePoint()-history1.initialTimePoint()))
                {
                  history1.publishDiagnosis(history1.initialTimePoint() + (history2.lastTimePoint()-history2.initialTimePoint()),history1.getDiagnosis(history1.lastTimePoint()));
                }
        }
  // cout << "*********************************" << endl;

  // cout << "*********************************" << endl;
  // cout << "HISTORY 2:" << endl;
  // cout << history2 << endl;

  { verbose<VbOutput>("Done.\n"); }
  
  { verbose<VbOutput>("Comparison... History 1 (%1%) History 2 (%2%)\n") % histFile1 % histFile2 ; }
  makeComparison(history1,history2,"h1h2");

  History perfect1;
  perfect1.publishDiagnosis(perfect1.initialTimePoint(),diag1);
  perfect1.publishDiagnosis(perfect1.initialTimePoint() + (history1.lastTimePoint()-history1.initialTimePoint()), diag1);
  { verbose<VbOutput>("Comparison... History 1 Perfect\n"); }
  makeComparison(history1,perfect1,"h1p");
  { verbose<VbOutput>("Comparison... History 2 Perfect\n"); }

  History perfect2;
  perfect2.publishDiagnosis(perfect2.initialTimePoint(),diag2);
  perfect2.publishDiagnosis(perfect2.initialTimePoint() + (history2.lastTimePoint()-history2.initialTimePoint()), diag2);
  makeComparison(history2,perfect2,"h2p");

  
}





void  makeComparison(const History & history1,
                                         const History & history2, const string & prefix)
{
  list< pair<ptime,int> > results;
  ofstream file;
  historyCompare(history1,history2,Distance(),results);
  file.open(prefix + "Distance.data");
  for(const pair<ptime,int>  & point : results)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();

  results.clear();
  historyCompare(history1,history2,FcDistance(),results);
  file.open(prefix + "FcDistance.data");
  for(const pair<ptime,int>  & point : results)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();

  results.clear();
  historyCompare(history1,history2,CommonFaults(),results);
  file.open(prefix + "CommonFaults.data");
  for(const pair<ptime,int>  & point : results)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();

  results.clear();
  historyCompare(history1,history2,NonCommonFaults(),results);
  file.open(prefix + "NonCommonFaults.data");
  for(const pair<ptime,int>  & point : results)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();
  
  results.clear();
  historyCompare(history1,history2,BsDistance(),results);
  file.open(prefix + "BsDistance.data");
  for(const  pair<ptime,int>  & point : results)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();

  list< pair<ptime,double> > results2;
  historyCompare(history1,history2,Accuracy(),results2);
  file.open(prefix + "Accuracy.data");
  for(const  pair<ptime,double>  & point : results2)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();
  results2.clear();

  historyCompare(history1,history2,BsAccuracy(),results2);
  file.open(prefix + "BsAccuracy.data");
  for(const  pair<ptime,double>  & point : results2)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();
  results2.clear();
  
  historyCompare(history1,history2,FcAccuracy(),results2);
  file.open(prefix + "FcAccuracy.data");
  for(const  pair<ptime,double>  & point : results2)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();
  results2.clear();

  historyCompare(history1,history2,Precision(),results2);
  file.open(prefix + "Precision.data");
  for(const  pair<ptime,double>  & point : results2)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();
  results2.clear();

  historyCompare(history1,history2,BsPrecision(),results2);
  file.open(prefix + "BsPrecision.data");
  for(const  pair<ptime,double>  & point : results2)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();
  results2.clear();
  
  historyCompare(history1,history2,FcPrecision(),results2);
  file.open(prefix + "FcPrecision.data");
  for(const  pair<ptime,double>  & point : results2)
        {
          file << (point.first - history1.initialTimePoint()).total_seconds() << '\t' << point.second << endl;
        }
  file.close();
  results2.clear();
  
}