Description Usage Arguments Value References Examples
The algorithm of Swendsen & Wang (1987) forms clusters of neighbouring pixels, then updates all of the labels within a cluster to the same value. When simulating from the prior, such as a Potts model without an external field, this algorithm is very efficient.
1 |
beta |
The inverse temperature parameter of the Potts model. |
k |
The number of unique labels. |
neighbors |
A matrix of all neighbors in the lattice, one row per pixel. |
blocks |
A list of pixel indices, dividing the lattice into independent blocks. |
niter |
The number of iterations of the algorithm to perform. |
random |
Whether to initialize the labels using random or deterministic starting values. |
A list containing the following elements:
alloc
An n by k matrix containing the number of times that pixel i was allocated to label j.
z
An (n+1)
by k matrix containing the final sample from the Potts model after niter iterations of Swendsen-Wang.
sum
An niter
by 1 matrix containing the sum of like neighbors, i.e. the sufficient statistic of the Potts model, at each iteration.
Swendsen, R. H. & Wang, J.-S. (1987) "Nonuniversal critical dynamics in Monte Carlo simulations" Physical Review Letters 58(2), 86–88, DOI: doi: 10.1103/PhysRevLett.58.86
1 2 3 4 5 6 |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.