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 #else
  15 #error NO UPDATE-ALGORITHM DEFINED!
  16 #endif
  17
  18 #include "obs_phi2.hpp"
  19 #include "obs_phi2_hist.hpp"
  20
  21 o815 *O815;
  22 sim *Sim;
  23
  24 const complex<double> _i_ = complex<double>(0.0,1.0);
  25
  26 o815::comoption specOps[] = {
  27   { "mass", required_argument, NULL, 'm', "set mass", "min:max:inc" },
  28   { "", 0, NULL, 0, "", "" }
  29 };
  30
  31 void parseSpecOps()
  32 {
  33   for (int isopt = 0; isopt < O815->parsedSpecOps.size(); isopt++)
  34     switch(O815->parsedSpecOps[isopt].first) {
  35     case 'm':
  36       O815->paraQ->addRange("mass", O815->parsedSpecOps[isopt].second);
  37       break;
  38     }
  39 }
  40
  41 void helpHeader()
  42 {
  43   cout << "Usage: ./heatbath-"ALGORITHM" [OPTIONS] [obs1] ... [obsN]" << endl << endl;
  44 }
  45
  46 void parseLonelyArgs()
  47 {
  48   for (vector<char*>::iterator lonit = O815->lonelyArgs.begin(); lonit != O815->lonelyArgs.end(); ++lonit) {
  49     if ( strcmp(*lonit, "phi2") == 0 ) {
  50       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
  51       O815->observables.push_back(new obs_phi2(O815));
  52     }
  53     else if ( strcmp(*lonit, "phi2_hist") == 0 ) {
  54       *O815->out->log << "MASTER: registered observable: phi2_hist" << endl << flush;
  55       O815->observables.push_back(new obs_phi2_hist(O815));
  56     }
  57   }
  58 }
  59
  60 int main (int argc, char *argv[])
  61 {
  62   O815 = new o815(argc, argv, "heatbath-"ALGORITHM, specOps, &helpHeader);
  63
  64   O815->addPara("mass", 1);
  65
  66   parseSpecOps();
  67
  68   O815->postParaInit();
  69
  70   O815->Sim = new sim(O815);
  71
  72   parseLonelyArgs();
  73
  74   O815->mainLoop();
  75
  76   delete O815;
  77   return 0;
  78 }