#include <iostream>
#include <iomanip>

#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

#include "JOscProb/JOscParametersGrid.hh"


/**
 * \file
 * Auxiliary program to test JOscParametersGrid.
 *
 * \author bjung
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;

  JOscParametersGrid      parameters;
  int                     debug;

  try {
    
    JParser<> zap;
    
    zap['@'] = make_field(parameters) = JPARSER::initialised();
    zap['d'] = make_field(debug)      = 3;

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

  NOTICE("User-defined parameters:\n" << parameters);

  ASSERT(parameters.size() > 0);

  JOscParametersGrid parameters1("sinsqTh12", make_grid(5, 0.3,    0.4),
				 "deltaCP",   make_grid(6, 1.5,    1.6),
				 "dM21sq",    make_grid(7, 7.4e-5, 7.5e-5));

  NOTICE("parameters1:\n" << parameters1);

  ASSERT(parameters1.is_valid() && parameters1.size() == 3);

  parameters1.set("sinsqTh23", make_grid(3, -0.6, -0.4));

  NOTICE("parameters1.sinsqTh23 = -0.5 -> parameters1:\n" << parameters1);

  ASSERT(!parameters1.is_valid() && parameters1.size() == 4);

  JOscParametersGrid parameters2("dM31sq",    make_grid(5, -2.6e-3, -2.0e-3),
				 "sinsqTh13", make_grid(2,  0.01,    0.03  ),
				 "sinsqTh23", make_grid(3,  0.57,    0.59 ));

  NOTICE("parameters2:\n" << parameters2);

  ASSERT(parameters2.is_valid() && parameters2.size() == 3);
  ASSERT(!parameters1.contains(parameters2));

  NOTICE("Joining parameters2 with parameters1" << endl);

  parameters1.join(parameters2);

  ASSERT(parameters1.is_valid() && parameters1.size() == 6);
  ASSERT(parameters1.contains(parameters2));
}