Dear Wei,
You can obtain a list of measurements by
data = pyalps.loadMeasurements(pyalps.getResultFiles(prefix='data/parmXY')) measurements = { item.props['observable'] for item in pyalps.flatten(data) }
Best, Synge
On Sep 30, 2020, at 9:49, Wei Wang wei.wang@jhu.edu wrote:
Hi,
I am just new to ALPS, and I have a silly question about running MC simulations of XY model:
Following the tutorials, I use python to run the simulations and change the model from “Ising” to “XY”, like:
import pyalps
parms = [ { 'LATTICE': "square lattice", 'T': 1, 'J': 1, 'THERMALIZATION': 10000, 'SWEEPS': 400000, 'UPDATE': "cluster", 'MODEL': "XY", 'L': 8 } ]
input_file = pyalps.writeInputFiles('data/parmXY', parms) pyalps.runApplication('spinmc', input_file, Tmin=5) result_files = pyalps.getResultFiles(prefix='data/parmXY')
data = pyalps.loadMeasurements(result_files)
how can I get some observables like <\cos\theta>, where \theta is the direction of the spin at each site? If such a quantity is not provided directly by ALPS, is there a possible method to get these quantities I want. I just get lost when I see the source code of ALPS.
Thank you in advance!
Best, Wei
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@lists.phys.ethz.ch.