summaryrefslogtreecommitdiff
path: root/sample.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'sample.cpp')
-rw-r--r--sample.cpp117
1 files changed, 0 insertions, 117 deletions
diff --git a/sample.cpp b/sample.cpp
deleted file mode 100644
index 71ee533..0000000
--- a/sample.cpp
+++ /dev/null
@@ -1,117 +0,0 @@
-
-#include <getopt.h>
-#include <fstream>
-
-#include "randutils/randutils.hpp"
-
-#define WOLFF_NO_FIELD
-#define WOLFF_USE_FINITE_STATES
-#include "wolff/lib/wolff_models/ising.hpp"
-#include "wolff/lib/wolff.hpp"
-
-using namespace wolff;
-
-class sample : public measurement<ising_t, ising_t, graph<>> {
- private:
- typedef struct _dat {
- int e;
- int m;
- } dat;
- std::ofstream e_file;
- dat e;
- unsigned w;
- unsigned N;
- unsigned τ;
- public:
- sample(const wolff::system<ising_t, ising_t, graph<>> &s, unsigned D, unsigned L, double T, unsigned wait) {
- e_file.open("sample_" + std::to_string(D) + "_" + std::to_string(L) + "_" + std::to_string(T) + ".bin" , std::ios::out | std::ios::binary | std::ios::app);
- e.e = - s.ne;
- e.m = s.nv;
- N = 0;
- τ = 2 * (unsigned)ceil(pow(L, 0.3));
- w = wait;
- }
-
- ~sample() {
- e_file.close();
- }
-
- void plain_site_transformed(const wolff::system<ising_t, ising_t, graph<>>&, const typename graph<>::vertex&, const ising_t& si_new) override {
- e.m += (unsigned)2 * si_new;
- }
-
- void plain_bond_visited(const wolff::system<ising_t, ising_t, graph<>>&s, const typename graph<>::halfedge& ed, const ising_t& si_new, double dE) override {
- e.e -= 2 * (si_new * s.s[ed.neighbor.ind]);
- }
-
- void post_cluster(unsigned n, unsigned, const wolff::system<ising_t, ising_t, graph<>>&) override {
- if (N >= w && (N - w) % τ == 0) {
- e_file.write((char *)&e, 2 * sizeof(int));
- }
-
- N++;
- }
-};
-
-int main(int argc, char *argv[]) {
-
- // set defaults
- unsigned N = (unsigned)1e4;
- unsigned w = (unsigned)1e3;
- unsigned D = 2;
- unsigned L = 128;
- double T = 2 / log(1 + sqrt(2));
-
- int opt;
-
- // take command line arguments
- while ((opt = getopt(argc, argv, "N:D:L:T:w:")) != -1) {
- switch (opt) {
- case 'N': // number of steps
- N = (unsigned)atof(optarg);
- break;
- case 'w': // number of steps
- w = (unsigned)atof(optarg);
- break;
- case 'D': // dimension
- D = atoi(optarg);
- break;
- case 'L': // linear size
- L = atoi(optarg);
- break;
- case 'T': // temperature
- T = atof(optarg);
- break;
- default:
- exit(EXIT_FAILURE);
- }
- }
-
- // define the spin-spin coupling
- std::function <double(const ising_t&, const ising_t&)> Z = [] (const ising_t& s1, const ising_t& s2) -> double {
- return (double)(s1 * s2);
- };
-
- // initialize the lattice
- graph<> G(D, L);
-
- // initialize the system
- wolff::system<ising_t, ising_t, graph<>> S(G, T, Z);
-
- // initialize the random number generator
- randutils::auto_seed_128 seeds;
- std::mt19937 rng{seeds};
-
- // define function that generates self-inverse rotations
- std::function <ising_t(std::mt19937&, const wolff::system<ising_t, ising_t, graph<>>&, const graph<>::vertex&)> gen_r = gen_ising<graph<>>;
-
- // initailze the measurement object
- sample A(S, D, L, T, w);
-
- // run wolff N times
- S.run_wolff(N, gen_r, A, rng);
-
- // exit
- return 0;
-}
-