X-Git-Url: http://git.treefish.org/~alex/phys/heatbath.git/blobdiff_plain/39a5c34637fb0b29f1e9129d6d2095cfb7a7931e..d017fdeec8363b3e871692219bb7eacb614a41f1:/heatbath.cpp

diff --git a/heatbath.cpp b/heatbath.cpp
index de99650..e6a7803 100644
--- a/heatbath.cpp
+++ b/heatbath.cpp
@@ -1,12 +1,27 @@
+#include <complex>
+
 #include "latlib/o815/o815.h"
 
-#include "sim.hpp"
+#ifdef HEATLESS_ALGORITHM
+#define ALGORITHM "heatless"
+#include "sim-heatless.hpp"
+#elif USUALHEAT_ALGORITHM
+#define ALGORITHM "usualheat"
+#include "sim-usualheat.hpp"
+#elif RELAXEDHEAT_ALGORITHM
+#define ALGORITHM "relaxedheat"
+#include "sim-relaxedheat.hpp"
+#else
+#error NO UPDATE-ALGORITHM DEFINED!
+#endif
 
 #include "obs_phi2.hpp"
 
 o815 *O815;
 sim *Sim;
 
+const complex<double> _i_ = complex<double>(0.0,1.0);
+
 o815::comoption specOps[] = {
   { "mass", required_argument, NULL, 'm', "set mass", "min:max:inc" },
   { "", 0, NULL, 0, "", "" }
@@ -23,8 +38,8 @@ void parseSpecOps()
 }
 
 void helpHeader() 
-{
-  cout << "Usage: ./heatbath [OPTIONS] [obs1] ... [obsN]" << endl << endl;
+{ 
+  cout << "Usage: ./heatbath-"ALGORITHM" [OPTIONS] [obs1] ... [obsN]" << endl << endl;
 }
 
 void parseLonelyArgs()
@@ -39,7 +54,7 @@ void parseLonelyArgs()
 
 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);
 
@@ -50,7 +65,7 @@ int main (int argc, char *argv[])
   O815->Sim = new sim(O815);
 
   parseLonelyArgs();
-
+  
   O815->mainLoop();
 
   delete O815;