Crate spin_sim 

Pure-Rust Ising model Monte Carlo on periodic hypercubic lattices.

Algorithms

Move	Function
Metropolis / Gibbs sweep	run_sweep_loop (sweep_mode)
Wolff single-cluster	run_sweep_loop (cluster_mode = "wolff")
Swendsen-Wang	run_sweep_loop (cluster_mode = "sw")
Parallel tempering	run_sweep_loop (pt_interval)
Houdayer ICM	run_sweep_loop (houdayer_interval)

Replicas are parallelized over threads with rayon.

Quick start

use spin_sim::{Lattice, Realization, run_sweep_loop};

let lattice = Lattice::new(vec![16, 16]);
let n_spins = lattice.n_spins;
let n_dims = lattice.n_dims;
let temps = vec![2.0, 2.27, 2.5];

// Uniform ferromagnetic couplings
let couplings = vec![1.0f32; n_spins * n_dims];
let mut real = Realization::new(&lattice, couplings, &temps, 2, 42);

let result = run_sweep_loop(
    &lattice, &mut real,
    2, temps.len(),
    5000, 1000,
    "metropolis",
    Some(1), "wolff",
    Some(1), None,
    &|| {},
);

For a Python interface, see the peapods package.

Re-exports

pub use lattice::Lattice;

pub use simulation::aggregate_results;

pub use simulation::run_sweep_loop;

pub use simulation::Realization;

pub use simulation::SweepResult;

Modules

lattice

simulation