Ising rule with Metropolis–Hastings single-spin-flip dynamics. More...

#include <tpl_ca_stochastic_rules.H>

Public Types
using	engine_type = Engine

Public Member Functions
	Ising_Metropolis_Rule (const double T, const double J=1.0, const double H=0.0, const std::uint64_t master_seed=0)
	Build a Metropolis Ising rule.

double	temperature () const noexcept

double	coupling () const noexcept

double	field () const noexcept

std::uint64_t	master_seed () const noexcept

void	set_master_seed (std::uint64_t s) noexcept

template<typename State , std::size_t Rank>
State	operator() (const State &current, Neighbor_View< State > neighbours, const Cell_Context< Rank > &ctx) const
	Compute the next Ising state for cell `current` (Metropolis dynamics).

Private Attributes
double	T_

double	J_

double	H_

std::uint64_t	master_seed_

Detailed Description

template<typename Engine = std::mt19937_64>
class Aleph::CA::Ising_Metropolis_Rule< Engine >

Ising rule with Metropolis–Hastings single-spin-flip dynamics.

Always flip when deltaE <= 0; otherwise flip with probability exp(-deltaE / T). Like Glauber dynamics, Metropolis dynamics satisfy detailed balance and target the Boltzmann distribution of the Ising Hamiltonian.

Template Parameters

Engine RNG engine; defaults to std::mt19937_64.

Definition at line 489 of file tpl_ca_stochastic_rules.H.

Member Typedef Documentation

◆ engine_type

template<typename Engine = std::mt19937_64>

using Aleph::CA::Ising_Metropolis_Rule< Engine >::engine_type = Engine

Definition at line 497 of file tpl_ca_stochastic_rules.H.

Constructor & Destructor Documentation

◆ Ising_Metropolis_Rule()

template<typename Engine = std::mt19937_64>

Aleph::CA::Ising_Metropolis_Rule< Engine >::Ising_Metropolis_Rule	(	const double	T,
		const double	J = `1.0`,
		const double	H = `0.0`,
		const std::uint64_t	master_seed = `0`
	)

inline

Build a Metropolis Ising rule.

Parameters

[in]	T	temperature; must be strictly positive.
[in]	J	coupling constant (default `+1`).
[in]	H	external magnetic field (default `0`).
[in]	master_seed	master seed forwarded to `cell_seed()`.

Exceptions

std::domain_error if T <= 0.

Definition at line 507 of file tpl_ca_stochastic_rules.H.

References ah_domain_error_if.

Member Function Documentation

◆ coupling()

template<typename Engine = std::mt19937_64>

double Aleph::CA::Ising_Metropolis_Rule< Engine >::coupling ( ) const

inlinenoexcept

Definition at line 521 of file tpl_ca_stochastic_rules.H.

References Aleph::CA::Ising_Metropolis_Rule< Engine >::J_.

◆ field()

template<typename Engine = std::mt19937_64>

double Aleph::CA::Ising_Metropolis_Rule< Engine >::field ( ) const

inlinenoexcept

Definition at line 526 of file tpl_ca_stochastic_rules.H.

References Aleph::CA::Ising_Metropolis_Rule< Engine >::H_.

◆ master_seed()

template<typename Engine = std::mt19937_64>

std::uint64_t Aleph::CA::Ising_Metropolis_Rule< Engine >::master_seed ( ) const

inlinenoexcept

Definition at line 531 of file tpl_ca_stochastic_rules.H.

References Aleph::CA::Ising_Metropolis_Rule< Engine >::master_seed_.

◆ operator()()

template<typename Engine = std::mt19937_64>

template<typename State , std::size_t Rank>

State Aleph::CA::Ising_Metropolis_Rule< Engine >::operator()	(	const State &	current,
		Neighbor_View< State >	neighbours,
		const Cell_Context< Rank > &	ctx
	)		const