Migrated u1ext-fancy simulation algorithm.

[phys/u1casc.git] / u1casc-flux / u1casc-flux.cpp

diff --git a/u1casc-flux/u1casc-flux.cpp b/u1casc-flux/u1casc-flux.cpp
new file mode 100644 (file)
index 0000000..0061116
--- /dev/null
+++ b/u1casc-flux/u1casc-flux.cpp
@@ -0,0 +1,116 @@
+#include <complex>
+#include <cstring>
+
+#include "latlib/o815/o815.h"
+
+#include "sim.hpp"
+
+o815 *O815;
+sim *Sim;
+
+//const complex<double> _i_ = complex<double>(0.0,1.0);
+
+//#include "obs_phi2.hpp"
+//#include "obs_plaq.hpp"
+//#include "obs_corr.hpp"
+
+o815::comoption specOps[] = {
+  { "kappaone", required_argument, NULL, 'r', "set inverse mass kappa_1", "min:max:inc" },
+  { "kappatwo", required_argument, NULL, 's', "set inverse mass kappa_2", "min:max:inc" },
+  { "beta", required_argument, NULL, 'b', "set inverse gauge coupling beta", "min:max:inc" },
+  { "lambdaone", required_argument, NULL, 'p', "set quartic coupling lambda_1", "min:max:inc" },
+  { "lambdatwo", required_argument, NULL, 'q', "set quartic coupling lambda_2", "min:max:inc" },
+  { "kappa", required_argument, NULL, 'k', "set kappa_1 = kappa_2 = kappa", "min:max:inc" },
+  { "lambda", required_argument, NULL, 'l', "set lambda_1 = lambda_2 = lambda", "min:max:inc" },
+  { "muone", required_argument, NULL, 'n', "set chemical potential mu_1", "min:max:inc" },
+  { "mutwo", required_argument, NULL, 'u', "set chemical potential mu_2", "min:max:inc" },
+  { "mu", required_argument, NULL, 'm', "set mu_1 = mu_2 = mu", "min:max:inc" },
+  { "", 0, NULL, 0, "", "" }
+};
+
+void parseSpecOps() 
+{
+  for (int isopt = 0; isopt < O815->parsedSpecOps.size(); isopt++)
+    switch(O815->parsedSpecOps[isopt].first) {
+    case 'r':
+      O815->paraQ->addRange("kappaone", O815->parsedSpecOps[isopt].second);
+      break;
+    case 's':
+      O815->paraQ->addRange("kappatwo", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'b':
+      O815->paraQ->addRange("beta", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'p':
+      O815->paraQ->addRange("lambdaone", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'q':
+      O815->paraQ->addRange("lambdatwo", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'l':
+      O815->paraQ->addRange("lambdaone", O815->parsedSpecOps[isopt].second);
+      O815->paraQ->linkParas("lambdatwo", "lambdaone");
+      break;
+    case 'k':
+      O815->paraQ->addRange("kappaone", O815->parsedSpecOps[isopt].second);
+      O815->paraQ->linkParas("kappatwo", "kappaone");
+      break;
+    case 'm':
+      O815->paraQ->addRange("muone", O815->parsedSpecOps[isopt].second);
+      O815->paraQ->linkParas("mutwo", "muone");
+      break;
+    case 'n':
+      O815->paraQ->addRange("muone", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'u':
+      O815->paraQ->addRange("mutwo", O815->parsedSpecOps[isopt].second);
+      break;
+    }
+}
+
+void helpHeader() 
+{ 
+  cout << "Usage: ./u1casc-flux [OPTIONS] [obs1] ... [obsN]" << endl << endl;
+}
+
+void parseLonelyArgs()
+{
+  /*
+  for (vector<char*>::iterator lonit = O815->lonelyArgs.begin(); lonit != O815->lonelyArgs.end(); ++lonit) {
+    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, "plaq") == 0 ) {
+      *O815->out->log << "MASTER: registered observable: plaq" << endl << flush;
+      O815->observables.push_back(new obs_plaq(O815));
+    }
+    else if ( strcmp(*lonit, "corr") == 0 ) {
+      *O815->out->log << "MASTER: registered observable: corr" << endl << flush;
+      O815->observables.push_back(new obs_corr(O815));
+    }
+  }
+  */
+}
+
+int main (int argc, char *argv[])
+{
+  O815 = new o815(argc, argv, "u1casc-flux", specOps, &helpHeader);
+
+  O815->addPara("beta", 1);
+  O815->addPara("lambdaone", 1);
+  O815->addPara("lambdatwo", 1);
+  O815->addPara("kappaone", 1);  
+  O815->addPara("kappatwo", 1);
+  O815->addPara("muone", 0);  
+  O815->addPara("mutwo", 0);
+
+  parseSpecOps();
+  O815->postParaInit();
+  O815->Sim = new sim(O815);
+  parseLonelyArgs();
+  O815->mainLoop();
+
+  delete O815;
+  return 0;
+}