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 #include "obs_plaq.hpp"
  15
  16 o815::comoption specOps[] = {
  17   { "kappaone", required_argument, NULL, 'r', "set inverse mass kappa_1", "min:max:inc" },
  18   { "kappatwo", required_argument, NULL, 's', "set inverse mass kappa_2", "min:max:inc" },
  19   { "beta", required_argument, NULL, 'b', "set inverse gauge coupling beta", "min:max:inc" },
  20   { "lambdaone", required_argument, NULL, 'p', "set quartic coupling lambda_1", "min:max:inc" },
  21   { "lambdatwo", required_argument, NULL, 'q', "set quartic coupling lambda_2", "min:max:inc" },
  22   { "kappa", required_argument, NULL, 'k', "set kappa_1 = kappa_2 = kappa", "min:max:inc" },
  23   { "lambda", required_argument, NULL, 'l', "set lambda_1 = lambda_2 = lambda", "min:max:inc" },
  24   { "", 0, NULL, 0, "", "" }
  25 };
  26
  27 void parseSpecOps()
  28 {
  29   for (int isopt = 0; isopt < O815->parsedSpecOps.size(); isopt++)
  30     switch(O815->parsedSpecOps[isopt].first) {
  31     case 'r':
  32       O815->paraQ->addRange("kappaone", O815->parsedSpecOps[isopt].second);
  33       break;
  34     case 's':
  35       O815->paraQ->addRange("kappatwo", O815->parsedSpecOps[isopt].second);
  36       break;
  37     case 'b':
  38       O815->paraQ->addRange("beta", O815->parsedSpecOps[isopt].second);
  39       break;
  40     case 'p':
  41       O815->paraQ->addRange("lambdaone", O815->parsedSpecOps[isopt].second);
  42       break;
  43     case 'q':
  44       O815->paraQ->addRange("lambdatwo", O815->parsedSpecOps[isopt].second);
  45       break;
  46     case 'l':
  47       O815->paraQ->addRange("lambdaone", O815->parsedSpecOps[isopt].second);
  48       O815->paraQ->linkParas("lambdatwo", "lambdaone");
  49       break;
  50     case 'k':
  51       O815->paraQ->addRange("kappaone", O815->parsedSpecOps[isopt].second);
  52       O815->paraQ->linkParas("kappatwo", "kappaone");
  53       break;
  54     }
  55 }
  56
  57 void helpHeader()
  58 {
  59   cout << "Usage: ./u1casc-ordinary [OPTIONS] [obs1] ... [obsN]" << endl << endl;
  60 }
  61
  62 void parseLonelyArgs()
  63 {
  64   for (vector<char*>::iterator lonit = O815->lonelyArgs.begin(); lonit != O815->lonelyArgs.end(); ++lonit) {
  65     if ( strcmp(*lonit, "phi2") == 0 ) {
  66       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
  67       O815->observables.push_back(new obs_phi2(O815));
  68     }
  69     else if ( strcmp(*lonit, "plaq") == 0 ) {
  70       *O815->out->log << "MASTER: registered observable: plaq" << endl << flush;
  71       O815->observables.push_back(new obs_plaq(O815));
  72     }
  73   }
  74 }
  75
  76 int main (int argc, char *argv[])
  77 {
  78   O815 = new o815(argc, argv, "u1casc-ordinary", specOps, &helpHeader);
  79
  80   O815->addPara("beta", 1);
  81   O815->addPara("lambdaone", 1);
  82   O815->addPara("lambdatwo", 1);
  83   O815->addPara("kappaone", 1);
  84   O815->addPara("kappatwo", 1);
  85
  86   parseSpecOps();
  87   O815->postParaInit();
  88   O815->Sim = new sim(O815);
  89   parseLonelyArgs();
  90   O815->mainLoop();
  91
  92   delete O815;
  93   return 0;
  94 }