View source: R/isingLibWrapper.R
transitionProbability1D | R Documentation |
Given a vector of flip sites, 1s or -1s, representing up and down spins respectively, and an other flip sites, return the transition probability, applying periodic boundary conditions, i.e., cyclic. This function calls the C function 'transitionProbability1D'.
transitionProbability1D(ikBT, x, xflip, J, H, probSel)
ikBT |
1/kB*T (Boltzmann factor) |
x |
1D Spin sites on the lattice. |
xflip |
1D Spin sites on the lattice: after a flip. |
J |
Interaction strength |
H |
External field |
probSel |
Which transition probability to use. 1 for Metropolis 2 for Glauber |
Returns transition probability.
Mehmet Suzen <mehmet.suzen@physics.org>
n <- 10 # 10 spin sites
mySites <- genConfig1D(n) # Generate sites
mySitesNew <- flipConfig1D(mySites)
# only short-range part
transitionProbability1D(1.0, mySites, mySitesNew, 1.0, 0.0, 1) # Metropolis
transitionProbability1D(1.0, mySites, mySitesNew, 1.0, 0.0, 2) # Glauber
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.