Dear all, I have just installed ALPS via conda and faced an unexpected problem with the tutorial dwa-01-bosons. While running the program tutorial1a.py on my laptop after thermalization is finished I get a segmentation fault. I changed the code a bit (see below) to run this tutorial on the bigger machine on which there is ALPS also installed via conda. The segmentation fault does not occur, however the program does not produce any output files and the produced graph is empty. I have also tried to run DWA with these parameters without python - using input file, but this have also ended with a segmentation fault.
import pyalps
bin_dir = '/net/people/plgprosniak/miniconda2/bin/'
parms = [] for t in [0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1]: parms.append( { 'LATTICE' : "square lattice", 'MODEL' : "boson Hubbard", 'T' : 0.1, 'L' : 4 , 't' : t , 'mu' : 0.5, 'U' : 1.0 , 'Nmax' : 2 , 'THERMALIZATION' : 100000, 'SWEEPS' : 2000000, 'SKIP' : 500, 'MEASURE[Winding Number]': 1 } )
input_file = pyalps.writeInputFiles('parm1a', parms) res = pyalps.runApplication(bin_dir+'dwa', input_file, Tmin=5, writexml=True)
# Setting graphical output import matplotlib as mpl mpl.use('agg')
# Evaluating the simulation and preparing plots using Python import pyalps import matplotlib.pyplot as plt import pyalps.plot as aplt
results = pyalps.getResultFiles(prefix='parm1a') print(results)
data = pyalps.loadMeasurements(results,'Stiffness') rhos = pyalps.collectXY(data,x='t',y='Stiffness')
fig=plt.figure() aplt.plot(rhos) plt.xlabel('Hopping $t/U$') plt.ylabel('Superfluid density $\rho _s$') fig.savefig('/net/scratch/people/plgprosniak/test/test_dwa.png')
Best regards, Oskar A. Prośniak