netket.sampler.MetropolisLocal

class netket.sampler.MetropolisLocal(machine, n_chains=16, sweep_size=None, batch_size=None)[source]

Sampler acting on one local degree of freedom.

This sampler acts locally only on one local degree of freedom \(s_i\), and proposes a new state: \(s_1 \dots s^\prime_i \dots s_N\), where \(s^\prime_i \neq s_i\).

The transition probability associated to this sampler can be decomposed into two steps:

1. One of the site indices \(i = 1\dots N\) is chosen with uniform probability. 2. Among all the possible (\(m\)) values that \(s_i\) can take, one of them is chosen with uniform probability.

For example, in the case of spin \(1/2\) particles, \(m=2\) and the possible local values are \(s_i = -1,+1\). In this case then MetropolisLocal is equivalent to flipping a random spin.

In the case of bosons, with occupation numbers \(s_i = 0, 1, \dots n_{\mathrm{max}}\), MetropolisLocal would pick a random local occupation number uniformly between \(0\) and \(n_{\mathrm{max}}\).

__init__(machine, n_chains=16, sweep_size=None, batch_size=None)[source]

Constructs a new MetropolisLocal sampler given a machine.

Parameters

  • machine – A machine \(\Psi(s)\) used for the sampling. The probability distribution being sampled from is \(F(\Psi(s))\), where the function \(F(X)\), is arbitrary, by default \(F(X)=|X|^2\).

  • n_chains – The number of Markov Chain to be run in parallel on a single process.

  • sweep_size – The number of exchanges that compose a single sweep. If None, sweep_size is equal to the number of degrees of freedom (n_visible).

  • batch_size – The batch size to be used when calling log_val on the given Machine. If None, batch_size is equal to the number Markov chains (n_chains).

Examples

Sampling from a RBM machine in a 1D lattice of spin 1/2

>>> import netket as nk
>>>
>>> g=nk.graph.Hypercube(length=10,n_dim=2,pbc=True)
>>> hi=nk.hilbert.Spin(s=0.5,graph=g)
>>>
>>> # RBM Spin Machine
>>> ma = nk.machine.RbmSpin(alpha=1, hilbert=hi)
>>>
>>> # Construct a MetropolisLocal Sampler
>>> sa = nk.sampler.MetropolisLocal(machine=ma)
>>> print(sa.machine.hilbert.size)
100

Methods

__init__(machine[, n_chains, sweep_size, …])

Constructs a new MetropolisLocal sampler given a machine.

generate_samples(n_samples[, init_random, …])

reset([init_random])

samples(n_max[, init_random])

Attributes

machine_func