summaryrefslogtreecommitdiff
path: root/mixedStokesFromSaddles.cpp
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2021-11-12 16:41:14 +0100
committerJaron Kent-Dobias <jaron@kent-dobias.com>2021-11-12 16:41:14 +0100
commitc9cb8dd499c2b2648a78e63f2a0810451a69f8fd (patch)
tree49f58bf75970ac9b7270f5999a628286bef6f8ab /mixedStokesFromSaddles.cpp
parent646dcb5b35afb2f7dbcdae59fe7280d2de91546c (diff)
downloadcode-c9cb8dd499c2b2648a78e63f2a0810451a69f8fd.tar.gz
code-c9cb8dd499c2b2648a78e63f2a0810451a69f8fd.tar.bz2
code-c9cb8dd499c2b2648a78e63f2a0810451a69f8fd.zip
Some tweaking of the Stokes finder, how the lines are recorded, and fewer executables.
Diffstat (limited to 'mixedStokesFromSaddles.cpp')
-rw-r--r--mixedStokesFromSaddles.cpp54
1 files changed, 0 insertions, 54 deletions
diff --git a/mixedStokesFromSaddles.cpp b/mixedStokesFromSaddles.cpp
deleted file mode 100644
index 2373146..0000000
--- a/mixedStokesFromSaddles.cpp
+++ /dev/null
@@ -1,54 +0,0 @@
-#include <getopt.h>
-#include <chrono>
-#include <fstream>
-
-#include "collectStokesData.hpp"
-
-#include "pcg-cpp/include/pcg_random.hpp"
-#include "randutils/randutils.hpp"
-#include "unsupported/Eigen/CXX11/src/Tensor/TensorFFT.h"
-
-#define PSPIN_P 3
-const int p = PSPIN_P; // polynomial degree of Hamiltonian
-
-using Rng = randutils::random_generator<pcg32>;
-
-int main(int argc, char* argv[]) {
- // model parameters
- unsigned N = 10; // number of spins
- // simulation parameters
- Real ε = 1e-15;
- Real δ = 1;
- unsigned n = 10;
-
- int opt;
-
- while ((opt = getopt(argc, argv, "N:e:d:n:")) != -1) {
- switch (opt) {
- case 'N':
- N = (unsigned)atof(optarg);
- break;
- case 'e':
- ε = atof(optarg);
- break;
- case 'd':
- δ = atof(optarg);
- break;
- case 'n':
- n = atof(optarg);
- break;
- default:
- exit(1);
- }
- }
-
- Rng r;
-
- for (unsigned i = 0; i < n; i++) {
- auto tag = std::chrono::high_resolution_clock::now();
- std::ofstream output("stokes_" + std::to_string(tag.time_since_epoch().count()) + ".dat");
- collectStokesData<2, 4>(output, N, r.engine(), ε, δ, false, 1.0, 0.01);
- }
-
- return 0;
-}