Wang and Landau Algorithm - Sample Code

Sample Code

The following is a sample code of the Wang–Landau algorithm in Python, considering that the proposal is such that

The code considers a "system" which is the underlying system being studied.


currentEnergy = system.randomConfiguration # a random initial configuration while (f > epsilon): system.proposeConfiguration # a proposed configuration is proposed proposedEnergy = system.proposedEnergy # the energy of the proposed configuration computed if (random < exp(entropy-entropy)): # if accepted, update the energy and the system: currentEnergy = proposedEnergy system.acceptProposedConfiguration else: # if rejected system.rejectProposedConfiguration H += 1 entropy += f if (isFlat(H)): # isFlat tests whether the histogram is flat (e.g. 95% flatness) H = 0 f *= 0.5 # refine the f parameter

Read more about this topic:  Wang And Landau Algorithm

Famous quotes containing the words sample and/or code:

    “The present war having so long cut off all communication with Great-Britain, we are not able to make a fair estimate of the state of science in that country. The spirit in which she wages war is the only sample before our eyes, and that does not seem the legitimate offspring either of science or of civilization.”
    —Thomas Jefferson (1743–1826)

    “Motion or change, and identity or rest, are the first and second secrets of nature: Motion and Rest. The whole code of her laws may be written on the thumbnail, or the signet of a ring.”
    —Ralph Waldo Emerson (1803–1882)