Indeed it was a matter of old binaries. Anyway I get the segmentation fault even in the case of V=0 which is a BH model, and this used to work before. I will try the older versions.
The input that gives me a segmentation fault is just (version I have just downloaded from svn):
LATTICE="open chain lattice" L=30
MODEL="boson Hubbard" Nmax=5 N_total=5
t=.5 U=10 mu=0.2
THERMALIZATION=15000 SWEEPS=200000 SKIP=50
MEASURE[Local Density]=1 tof_phase=0
{ T=0.08 }
On 11 Oct 2016, at 17:59, Tama Ma tamama@yotcopi.com wrote:
Hi Mateusz,
Extended BHM is supported in the later version of DWA.
Several checks:
- Can you download the latest version from trunk?
- Can you check whether +V is supported also?
Best regards, Tamama
Sent from my iPhone
On 11 Oct 2016, at 11:13 PM, Mateusz Łącki mateusz.lacki@gmail.com wrote:
Dear All, could I ask if the DWA code actually supports the extended BH model (as defined in “boson Hubbard” hamiltonian in a standard issue of models.xml)? the results I obtain seem not to depend at all on the value of the V parameter. With sparsediag I observe no such behaviour.
thebondterm contains “hopping” (proportional to t#) and interaction (proportional to V#). If I wanted to split (as the hopping is treated separately in worm-type algorithms I guess…) it into two bond-terms it seems that only the last bondterm is taken into account:
this results in t=0 result
<SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2 </SITETERM> <BONDTERM source="i" target="j"> <PARAMETER name="t#" default="0"/> -t#*(bdag(i)*b(j)+bdag(j)*b(i)) </BONDTERM> <BONDTERM source="i" target="j"> <PARAMETER name="V#" default="0"/> -V#*(n(j)*n(i)) </BONDTERM>
this in V=0 result:
<SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2 </SITETERM> <BONDTERM source="i" target="j"> <PARAMETER name="V#" default="0"/> -V#*(n(j)*n(i)) </BONDTERM> <BONDTERM source="i" target="j"> <PARAMETER name="t#" default="0"/> -t#*(bdag(i)*b(j)+bdag(j)*b(i)) </BONDTERM>
Best, Mateusz Łącki
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.