JMultiMap.cc File Reference

Example program to test JTOOLS::JMultiMap. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <cmath>
#include "JTools/JMultiMap.hh"
#include "JTools/JMap.hh"
#include "JLang/JManip.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

Example program to test JTOOLS::JMultiMap.

Author

mdejong

Definition in file JMultiMap.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 52 of file JMultiMap.cc.

{
  using namespace std;
  using namespace JPP;

  double      precision;
  int         debug;

  try {

    JParser<> zap("Example program to test multi-dimensional map.");

    zap['e'] = make_field(precision)        = 1.0e-10;
    zap['d'] = make_field(debug)            = 3;

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


  const double xmin = -1.0;
  const double xmax = +1.0;
  const double nx   = 3;
   
  typedef JMAPLIST<JMap,
                   JMap,
                   JMap>::maplist                JMaplist_t;

  typedef JMultiMap<double, double, JMaplist_t>  JMultimap_t;

  JMultimap_t buffer;

  for (double x = xmin; x <= xmax; x += (xmax - xmin)/(nx - 1)) {
    for (double y = xmin; y <= xmax; y += (xmax - xmin)/(nx - 1)) {
      for (double z = xmin; z <= xmax; z += (xmax - xmin)/(nx - 1)) {
        buffer[x][y][z] = f1(x, y, z);
      }
    }
  }

  const JFormat_t format(6, 1, std::ios::fixed);

  setFormat< JMultiKey<3, const double> >(format);
  setFormat< JMultiKey<2, const double> >(format);
  setFormat< JMultiKey<1, const double> >(format);

  DEBUG("i->[second]*->(first|second)" << endl);
    
  for (JMultimap_t::super_const_iterator i = buffer.super_begin(); i != buffer.super_end(); ++i) {
    DEBUG(""
          << format << i->first                  << ' ' 
          << format << i->second->first          << ' ' 
          << format << i->second->second->first  << ' ' 
          << format << i->second->second->second << endl);

    ASSERT(fabs(f1(i->first,
                   i->second->first,
                   i->second->second->first) - i->second->second->second) <= precision, "Test iterator equality");
  }

  DEBUG("i->[second]*->(first|second)" << endl);

  for (JMultimap_t::super_const_reverse_iterator i = buffer.super_rbegin(); i != buffer.super_rend(); ++i) {
    DEBUG(""
          << format << i->first                  << ' ' 
          << format << i->second->first          << ' ' 
          << format << i->second->second->first  << ' ' 
          << format << i->second->second->second << endl);

    ASSERT(fabs(f1(i->first,
                   i->second->first,
                   i->second->second->first) - i->second->second->second) <= precision, "Test iterator equality");
  }

  DEBUG("*i.[second]*.(first|second)" << endl);

  for (JMultimap_t::super_const_reverse_iterator i = buffer.super_rbegin(); i != buffer.super_rend(); ++i) {
    DEBUG(""
          << format << (*i).first                << ' ' 
          << format << (*i).second.first         << ' ' 
          << format << (*i).second.second.first  << ' ' 
          << format << (*i).second.second.second << endl);

    ASSERT(fabs(f1(i->first,
                   i->second->first,
                   i->second->second->first) - i->second->second->second) <= precision, "Test iterator equality");
  }

  DEBUG("i.getKey() i.getValue()" << endl);

  for (JMultimap_t::super_const_iterator i = buffer.super_begin(); i != buffer.super_end(); ++i) {

    DEBUG(i.getKey() << ' ' << format << i.getValue() << endl);

    ASSERT(fabs(f1(i.getKey()) - i.getValue()) <= precision, "Test iterator equality");
  }

  return 0;
}