JCombinatorics.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>

#include "JTools/JCombinatorics.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"


/**
 * \file
 * Example program to test combinatorics of pairs of indices (class JTOOLS::JCombinatorics).
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;

  int     numberOfElements;
  int     debug;

  try { 

    JParser<> zap("Example program to test combinatorics of pairs of indices.");
    
    zap['n'] = make_field(numberOfElements);
    zap['d'] = make_field(debug)           = 1;

    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }


  using namespace JPP;

  if (numberOfElements <= 0) {
    FATAL("Number of elements " << numberOfElements << " <= 0" << endl);
  }


  JCombinatorics combinatorics(numberOfElements);

  int n = 0;

  for (int i = 0; i != (int) combinatorics.getNumberOfPairs(); ++i) {

    STATUS("event: " << setw(8) << i << '\r'); DEBUG(endl);

    JCombinatorics::pair_type pair = combinatorics.getPair(i);

    if (i != combinatorics.getIndex(pair.first, pair.second)) {

      ++n;

      ERROR(endl                   << ' ' << 
            setw(3) << pair.first  << ' ' <<
            setw(3) << pair.second << ' ' <<
            setw(4) << combinatorics.getIndex(pair.first, pair.second) << endl);
    }
  }
  DEBUG(endl); 

  NOTICE("Number of elements " << setw(8) << right << numberOfElements                 << endl);
  NOTICE("Number of pairs    " << setw(8) << right << combinatorics.getNumberOfPairs() << endl);
  NOTICE("Number of errors   " << setw(8) << right << n                                << endl);

  ASSERT(n == 0);
}