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