Thanks so much for your reply! I will try it.
By the way, there are may be some missing INTEGER(HID_T) kind definitions of int data type at the Fortran source code about the tebd application, For example,
at tebd/core/Hdf5Interface.f90
-------------------------------------------
INTEGER, INTENT(INOUT) :: groupId --> INTEGER(HID_T), INTENT(INOUT) :: groupId
-------------------------------------------
When I compile the source code, I met some int kind mismatch errors. These errors disappeared after I change the source code. It may be a bug about ALPS.
Best,
Rongyang
-----------------------------------------------------------------
Rongyang Sun
Ph.D. Candidate
Theoretical Quantum Matter Physics
Institute for Advanced Study, Tsinghua University, Beijing 100084, China
E-mail: sunry15@mails.tsinghua.edu.cn; sun-rongyang@outlook.com
On 27 Aug 2017, at 16:52, Sun Rongyang <sun-rongyang@outlook.com> wrote:
Dear ALPS users,
When I follow the ALPS TEBD tutorial about hard-core boson model(http://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:TEBD-01_bhquench) by Python, I measured the average particle number for each site, i.e. 'Avg Number' except the 'Loschmidt Echo' and 'V' by the following python code:
data = pyalps.load.loadTimeEvolution(pyalps.getResultFiles(prefix='1D_HardcoreBosonModel'), measurements=['Loschmidt Echo', 'V', 'Avg Number'])
Then I found the 'Avg Number' equals 0.3! This means that the total particle number is 0.3*10=3, not 5. In the tutorial, the conserved quantum number is 'N_total': 5. When I change the 'N_total' to 1 and 10, I got the same 'Avg Number'=0.3
Can any one give me some help? Thanks a lot!
Best regards,Rongyang Sun
-----------------------------------------------------------------Rongyang SunPh.D. CandidateTheoretical Quantum Matter PhysicsInstitute for Advanced Study, Tsinghua University, Beijing 100084, China
----
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.