Dear Prof. Todo, Thank you for your reply! I have added these two lines into my code, and print the measurements out:
set(['Magnetization^2', 'Binder Cumulant', 'Specific Heat', '|Magnetization|', 'Magnetization^4', '|Staggered Magnetization|', 'Staggered Magnetization^2', 'E.Magnetization^4’, 'E.Magnetization^2', 'Energy Density', 'Susceptibility', 'beta * Energy / sqrt(N)', 'Magnetization along Field', 'Magnetization^4 slope', 'Connected Susceptibility', 'Magnetization^2 slope', 'Binder Cumulant U2', 'Cluster size', 'Bond-type Energy', 'Energy', 'Bond-type Energy Density', 'Energy^2', '(beta * Energy)^2 / N', 'Binder Cumulant slope'])
I don’t find an observable like the <\cos\theta>, or which of the measurements above can give me the directions of the magnetization. Please forgive me if this is a silly question, I am a 1st-year graduate student learning how to use ALPS. Thank you! Regards, Wei
On Sep 30, 2020, at 09:45, Synge Todo <wistaria@phys.s.u-tokyo.ac.jp> wrote:
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 https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Falps.comp-phys.org%2F&data=02%7C01%7Cwei.wang%40jhu.edu%7C1fd88b20e2e4436ba74f08d864e2a4ad%7C9fa4f438b1e6473b803f86f8aedf0dec%7C0%7C0%7C637370271890872618&sdata=x2wfhZw4cHs1yHHDs3oc0i61Ht4lqdUvr5igLf%2BnhKI%3D&reserved=0
List info: https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.phys.ethz.ch%2F%2Flistinfo%2Fcomp-phys-alps-users&data=02%7C01%7Cwei.wang%40jhu.edu%7C1fd88b20e2e4436ba74f08d864e2a4ad%7C9fa4f438b1e6473b803f86f8aedf0dec%7C0%7C0%7C637370271890872618&sdata=To%2FPEErYCLzeGTzQ%2F1H0l6oy03PVd4MLhzVak3kctCw%3D&reserved=0 Archive: https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.phys.ethz.ch%2F%2Fpipermail%2Fcomp-phys-alps-users&data=02%7C01%7Cwei.wang%40jhu.edu%7C1fd88b20e2e4436ba74f08d864e2a4ad%7C9fa4f438b1e6473b803f86f8aedf0dec%7C0%7C0%7C637370271890872618&sdata=Nh1D42waDED8xwLau%2BX8OfYK2Fj9sGbDGIFdVOKrmHo%3D&reserved=0
Unsubscribe by writing a mail to comp-phys-alps-users-leave@lists.phys.ethz.ch.
---- Comp-phys-alps-users Mailing List for the ALPS Project https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Falps.comp-phys.org%2F&data=02%7C01%7Cwei.wang%40jhu.edu%7C1fd88b20e2e4436ba74f08d864e2a4ad%7C9fa4f438b1e6473b803f86f8aedf0dec%7C0%7C0%7C637370271890872618&sdata=x2wfhZw4cHs1yHHDs3oc0i61Ht4lqdUvr5igLf%2BnhKI%3D&reserved=0
List info: https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.phys.ethz.ch%2F%2Flistinfo%2Fcomp-phys-alps-users&data=02%7C01%7Cwei.wang%40jhu.edu%7C1fd88b20e2e4436ba74f08d864e2a4ad%7C9fa4f438b1e6473b803f86f8aedf0dec%7C0%7C0%7C637370271890872618&sdata=To%2FPEErYCLzeGTzQ%2F1H0l6oy03PVd4MLhzVak3kctCw%3D&reserved=0 Archive: https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.phys.ethz.ch%2F%2Fpipermail%2Fcomp-phys-alps-users&data=02%7C01%7Cwei.wang%40jhu.edu%7C1fd88b20e2e4436ba74f08d864e2a4ad%7C9fa4f438b1e6473b803f86f8aedf0dec%7C0%7C0%7C637370271890872618&sdata=Nh1D42waDED8xwLau%2BX8OfYK2Fj9sGbDGIFdVOKrmHo%3D&reserved=0
Unsubscribe by writing a mail to comp-phys-alps-users-leave@lists.phys.ethz.ch.