heatbath.cpp

   1 #include <complex>
   2
   3 #include "latlib/o815/o815.h"
   4
   5 #ifdef ALGORITHM_NOR
   6 #define ALGORITHM "nor"
   7 #include "sim-nor.hpp"
   8 #elif ALGORITHM_R
   9 #define ALGORITHM "r"
  10 #include "sim-r.hpp"
  11 #elif ALGORITHM_1MR
  12 #define ALGORITHM "1mr"
  13 #include "sim-1mr.hpp"
  14 #elif ALGORITHM_1MRPLUS
  15 #define ALGORITHM "1mr+"
  16 #include "sim-1mr+.hpp"
  17 #elif ALGORITHM_1MRMINUS
  18 #define ALGORITHM "1mr-"
  19 #include "sim-1mr-.hpp"
  20 #elif ALGORITHM_RMINUS
  21 #define ALGORITHM "r-"
  22 #include "sim-r-.hpp"
  23 #elif ALGORITHM_RPLUS
  24 #define ALGORITHM "r+"
  25 #include "sim-r+.hpp"
  26 #else
  27 #error NO UPDATE-ALGORITHM DEFINED!
  28 #endif
  29
  30 o815 *O815;
  31 sim *Sim;
  32
  33 const complex<double> _i_ = complex<double>(0.0,1.0);
  34
  35 #include "obs_phi2.hpp"
  36 #include "obs_phi2_hist.hpp"
  37 #include "obs_phip2_hist.hpp"
  38
  39 o815::comoption specOps[] = {
  40   { "mass", required_argument, NULL, 'm', "set mass", "min:max:inc" },
  41   { "", 0, NULL, 0, "", "" }
  42 };
  43
  44 void parseSpecOps()
  45 {
  46   for (int isopt = 0; isopt < O815->parsedSpecOps.size(); isopt++)
  47     switch(O815->parsedSpecOps[isopt].first) {
  48     case 'm':
  49       O815->paraQ->addRange("mass", O815->parsedSpecOps[isopt].second);
  50       break;
  51     }
  52 }
  53
  54 void helpHeader()
  55 {
  56   cout << "Usage: ./heatbath-"ALGORITHM" [OPTIONS] [obs1] ... [obsN]" << endl << endl;
  57 }
  58
  59 void parseLonelyArgs()
  60 {
  61   for (vector<char*>::iterator lonit = O815->lonelyArgs.begin(); lonit != O815->lonelyArgs.end(); ++lonit) {
  62     if ( strcmp(*lonit, "phi2") == 0 ) {
  63       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
  64       O815->observables.push_back(new obs_phi2(O815));
  65     }
  66     else if ( strcmp(*lonit, "phi2_hist") == 0 ) {
  67       *O815->out->log << "MASTER: registered observable: phi2_hist" << endl << flush;
  68       O815->observables.push_back(new obs_phi2_hist(O815));
  69     }
  70     else if ( strcmp(*lonit, "phip2_hist") == 0 ) {
  71       *O815->out->log << "MASTER: registered observable: phip2_hist" << endl << flush;
  72       O815->observables.push_back(new obs_phip2_hist(O815));
  73     }
  74   }
  75 }
  76
  77 int main (int argc, char *argv[])
  78 {
  79   O815 = new o815(argc, argv, "heatbath-"ALGORITHM, specOps, &helpHeader);
  80
  81   O815->addPara("mass", 1);
  82
  83   parseSpecOps();
  84
  85   O815->postParaInit();
  86
  87   O815->Sim = new sim(O815);
  88
  89   parseLonelyArgs();
  90
  91   O815->mainLoop();
  92
  93   delete O815;
  94   return 0;
  95 }