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