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:
1. Can you download the latest version from trunk?
2. 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.
----
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.