heatbath.cpp

   1 #include "latlib/o815/o815.h"
   2
   3 #include "sim.hpp"
   4
   5 #include "obs_phi2.hpp"
   6
   7 o815 *O815;
   8 sim *Sim;
   9
  10 o815::comoption specOps[] = {
  11   { "mass", required_argument, NULL, 'm', "set mass", "min:max:inc" },
  12   { "", 0, NULL, 0, "", "" }
  13 };
  14
  15 void parseSpecOps()
  16 {
  17   for (int isopt = 0; isopt < O815->parsedSpecOps.size(); isopt++)
  18     switch(O815->parsedSpecOps[isopt].first) {
  19     case 'm':
  20       O815->paraQ->addRange("mass", O815->parsedSpecOps[isopt].second);
  21       break;
  22     }
  23 }
  24
  25 void helpHeader()
  26 {
  27   cout << "Usage: ./heatbath [OPTIONS] [obs1] ... [obsN]" << endl << endl;
  28 }
  29
  30 void parseLonelyArgs()
  31 {
  32   for (vector<char*>::iterator lonit = O815->lonelyArgs.begin(); lonit != O815->lonelyArgs.end(); ++lonit) {
  33     if ( strcmp(*lonit, "phi2") == 0 ) {
  34       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
  35       O815->observables.push_back(new obs_phi2(O815));
  36     }
  37   }
  38 }
  39
  40 int main (int argc, char *argv[])
  41 {
  42   O815 = new o815(argc, argv, "heatbath", specOps, &helpHeader);
  43
  44   O815->addPara("mass", 1);
  45
  46   parseSpecOps();
  47
  48   O815->postParaInit();
  49
  50   O815->Sim = new sim(O815);
  51
  52   parseLonelyArgs();
  53
  54   O815->mainLoop();
  55
  56   delete O815;
  57   return 0;
  58 }