Dear developers,
I'm using the dirloop_sse code to simulate the hardcore boson Hubbard model, I met some problems there:
(1) An error occors in many simulations, that is, Internal error: size of prob table is zero. (2) In a ladder model, the stiffness I get from the sse code is much smaller (almost zero) than using other codes (not ALPS code). While in a square lattice, the sse code gives the same stiffness with other codes; BTW, does the definition of stiffness in the ALPS codes miss a factor of 2 in the denominator? (3) When I input some values of chemical potential (mu), I can get the corresponding particle density (rho). When mu is close to zero, the data of rho are not good, while for other mu, the data are nice.
PS: I attach my parameter file below, do I make something wrong in the parameters? LATTICE="ladder"; L=64; MODEL="hardcore boson"; V = 0.0; V' = 0.0; t = 1.0; t' = 1.0; T = 0.03125; SWEEPS=10000000; THERMALIZATION=1000000; { mu=-0.1; }
Best, TaoYing
comp-phys-alps-users@lists.phys.ethz.ch