u1casc-ordinary/u1casc-ordinary.cpp

   1 #include <complex>
   2 #include <cstring>
   3
   4 #include "latlib/o815/o815.h"
   5
   6 #include "sim.hpp"
   7
   8 o815 *O815;
   9 sim *Sim;
  10
  11 //const complex<double> _i_ = complex<double>(0.0,1.0);
  12
  13 //#include "obs_phi2.hpp"
  14
  15 o815::comoption specOps[] = {
  16   { "kappaone", required_argument, NULL, 'k', "set inverse mass kappa_1", "min:max:inc" },
  17   { "kappatwo", required_argument, NULL, 'm', "set inverse mass kappa_2", "min:max:inc" },
  18   { "beta", required_argument, NULL, 'b', "set inverse gauge coupling beta", "min:max:inc" },
  19   { "lambdaone", required_argument, NULL, 'p', "set quartic coupling lambda_1", "min:max:inc" },
  20   { "lambdatwo", required_argument, NULL, 'q', "set quartic coupling lambda_2", "min:max:inc" },
  21   { "", 0, NULL, 0, "", "" }
  22 };
  23
  24 void parseSpecOps()
  25 {
  26   for (int isopt = 0; isopt < O815->parsedSpecOps.size(); isopt++)
  27     switch(O815->parsedSpecOps[isopt].first) {
  28     case 'm':
  29       O815->paraQ->addRange("mass", O815->parsedSpecOps[isopt].second);
  30       break;
  31     }
  32 }
  33
  34 void helpHeader()
  35 {
  36   cout << "Usage: ./u1casc-ordinary [OPTIONS] [obs1] ... [obsN]" << endl << endl;
  37 }
  38
  39 void parseLonelyArgs()
  40 {
  41   for (vector<char*>::iterator lonit = O815->lonelyArgs.begin(); lonit != O815->lonelyArgs.end(); ++lonit) {
  42     /*
  43     if ( strcmp(*lonit, "phi2") == 0 ) {
  44       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
  45       O815->observables.push_back(new obs_phi2(O815));
  46     }
  47     else if ( strcmp(*lonit, "phi2_hist") == 0 ) {
  48       *O815->out->log << "MASTER: registered observable: phi2_hist" << endl << flush;
  49       O815->observables.push_back(new obs_phi2_hist(O815));
  50     }
  51     else if ( strcmp(*lonit, "phip2_hist") == 0 ) {
  52       *O815->out->log << "MASTER: registered observable: phip2_hist" << endl << flush;
  53       O815->observables.push_back(new obs_phip2_hist(O815));
  54     }
  55     */
  56   }
  57 }
  58
  59 int main (int argc, char *argv[])
  60 {
  61   O815 = new o815(argc, argv, "u1casc-ordinary", specOps, &helpHeader);
  62
  63   O815->addPara("beta", 1);
  64   O815->addPara("lambdaone", 1);
  65   O815->addPara("lambdatwo", 1);
  66   O815->addPara("kappaone", 1);
  67   O815->addPara("kappatwo", 1);
  68
  69   parseSpecOps();
  70   O815->postParaInit();
  71   O815->Sim = new sim(O815);
  72   parseLonelyArgs();
  73   O815->mainLoop();
  74
  75   delete O815;
  76   return 0;
  77 }