diff options
| -rw-r--r-- | .gitmodules | 6 | ||||
| m--------- | randutils | 0 | ||||
| -rw-r--r-- | sample.cpp | 117 | ||||
| m--------- | wolff | 0 |
4 files changed, 123 insertions, 0 deletions
diff --git a/.gitmodules b/.gitmodules new file mode 100644 index 0000000..09d63cc --- /dev/null +++ b/.gitmodules @@ -0,0 +1,6 @@ +[submodule "randutils"] + path = randutils + url = https://gist.github.com/imneme/540829265469e673d045 +[submodule "wolff"] + path = wolff + url = https://git.kent-dobias.com/wolff diff --git a/randutils b/randutils new file mode 160000 +Subproject 8486a610a954a8248c12485fb4cfc390a5f5f85 diff --git a/sample.cpp b/sample.cpp new file mode 100644 index 0000000..71ee533 --- /dev/null +++ b/sample.cpp @@ -0,0 +1,117 @@ + +#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; +} + diff --git a/wolff b/wolff new file mode 160000 +Subproject 49ac78a6c04e215950bc9c0f97368605e63da15 |