...

[phys/heatbath.git] / heatbath.cpp

diff --git a/heatbath.cpp b/heatbath.cpp
index de996508e0a5fc6adb056364a0030ea6716fc151..8df557bb3623b84c2e71d9b72227fd32789c5be7 100644 (file)
--- a/heatbath.cpp
+++ b/heatbath.cpp
@@ -1,12 +1,52 @@
+#include <complex>
+#include <cstring>
+

 #include "latlib/o815/o815.h"
 
 #include "latlib/o815/o815.h"
 

-#include "sim.hpp"
+#ifdef ALGORITHM_METRO
+#define ALGORITHM "metro"
+#include "sim-metro.hpp"
+#elif ALGORITHM_AHBMHB
+#include "sim-nonmetro.hpp"
+#define ALGORITHM "AhbMhb"
+#elif ALGORITHM_AHBMOR
+#include "sim-nonmetro.hpp"
+#define ALGORITHM "AhbMor"
+#elif ALGORITHM_AORMHB
+#include "sim-nonmetro.hpp"
+#define ALGORITHM "AorMhb"
+#elif ALGORITHM_AORMOR
+#include "sim-nonmetro.hpp"
+#define ALGORITHM "AorMor"
+#elif ALGORITHM_AORM0
+#include "sim-nonmetro.hpp"
+#define ALGORITHM "AorM0"
+#elif ALGORITHM_AHBM0
+#include "sim-nonmetro.hpp"
+#define ALGORITHM "AhbM0"
+#elif ALGORITHM_A0MOR
+#include "sim-nonmetro.hpp"
+#define ALGORITHM "A0Mor"
+#elif ALGORITHM_A0MHB
+#include "sim-nonmetro.hpp"
+#define ALGORITHM "A0Mhb"
+#else
+#error NO UPDATE-ALGORITHM DEFINED!
+#endif

 
 

-#include "obs_phi2.hpp"
+#ifdef NONMETRO_A_HB
+#error SADSADS
+#endif

 
 o815 *O815;
 sim *Sim;
 
 
 o815 *O815;
 sim *Sim;
 

+const complex<double> _i_ = complex<double>(0.0,1.0);
+
+#include "obs_phi2.hpp"
+#include "obs_phi2_hist.hpp"
+#include "obs_phip2_hist.hpp"
+

 o815::comoption specOps[] = {
   { "mass", required_argument, NULL, 'm', "set mass", "min:max:inc" },
   { "", 0, NULL, 0, "", "" }
 o815::comoption specOps[] = {
   { "mass", required_argument, NULL, 'm', "set mass", "min:max:inc" },
   { "", 0, NULL, 0, "", "" }


@@ -23,8 +63,8 @@ void parseSpecOps()
 }
 
 void helpHeader() 
 }
 
 void helpHeader() 

-{
-  cout << "Usage: ./heatbath [OPTIONS] [obs1] ... [obsN]" << endl << endl;
+{ 
+  cout << "Usage: ./heatbath-"ALGORITHM" [OPTIONS] [obs1] ... [obsN]" << endl << endl;

 }
 
 void parseLonelyArgs()
 }
 
 void parseLonelyArgs()


@@ -33,13 +73,21 @@ void parseLonelyArgs()
     if ( strcmp(*lonit, "phi2") == 0 ) {
       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
       O815->observables.push_back(new obs_phi2(O815));
     if ( strcmp(*lonit, "phi2") == 0 ) {
       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
       O815->observables.push_back(new obs_phi2(O815));

-    } 
+    }
+    else if ( strcmp(*lonit, "phi2_hist") == 0 ) {
+      *O815->out->log << "MASTER: registered observable: phi2_hist" << endl << flush;
+      O815->observables.push_back(new obs_phi2_hist(O815));
+    }
+    else if ( strcmp(*lonit, "phip2_hist") == 0 ) {
+      *O815->out->log << "MASTER: registered observable: phip2_hist" << endl << flush;
+      O815->observables.push_back(new obs_phip2_hist(O815));
+    }

   }
 }
 
 int main (int argc, char *argv[])
 {
   }
 }
 
 int main (int argc, char *argv[])
 {

-  O815 = new o815(argc, argv, "heatbath", specOps, &helpHeader);
+  O815 = new o815(argc, argv, "heatbath-"ALGORITHM, specOps, &helpHeader);

 
   O815->addPara("mass", 1);
 
 
   O815->addPara("mass", 1);
 


@@ -50,7 +98,7 @@ int main (int argc, char *argv[])
   O815->Sim = new sim(O815);
 
   parseLonelyArgs();
   O815->Sim = new sim(O815);
 
   parseLonelyArgs();

-
+  

   O815->mainLoop();
 
   delete O815;
   O815->mainLoop();
 
   delete O815;