New subject: Time propagation in TEBD stops prematurely/Obtaining wavefunctions in TEBD

19 Sep 2013


      Hi Michael,
Thanks for checking my code, there must be something wrong at my end somewhere. I'll see what I can do in terms of the cluster.
I did have another question, which is whether it is possible to extract the wave function from the HDF5 files produced by TEBD? Can this be done using the pyalps.hdf5 module?
Thanks,
Joseph
On 19 Sep 2013, at 11:00, comp-phys-alps-users-request@lists.phys.ethz.ch
 comp-phys-alps-users-request@lists.phys.ethz.ch wrote:
...
Send Comp-phys-alps-users mailing list submissions to
   comp-phys-alps-users@lists.phys.ethz.ch
To subscribe or unsubscribe via the World Wide Web, visit
   https://lists.phys.ethz.ch/listinfo/comp-phys-alps-users
or, via email, send a message with subject or body 'help' to
   comp-phys-alps-users-request@lists.phys.ethz.ch
You can reach the person managing the list at
   comp-phys-alps-users-owner@lists.phys.ethz.ch
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Comp-phys-alps-users digest..."
Today's Topics:

Re: Time propagation in TEBD stops prematurely (Michael Wall)
Re: Making the scheduler switch tasks,	obtaining MPI rank
(Peter Br?cker)
Re: Making the scheduler switch tasks,	obtaining MPI rank
(Matthias Troyer)


Message: 1
Date: Wed, 18 Sep 2013 08:48:45 -0600
From: Michael Wall mwall.physics@gmail.com
To: comp-phys-alps-users@lists.phys.ethz.ch
Subject: Re: [ALPS-users] Time propagation in TEBD stops prematurely
Message-ID:
   CA+DwVMrz2Y_1Kz7yDSyEqALjskd0iavzFM_9_KTMjPiNzkSzJA@mail.gmail.com
Content-Type: text/plain; charset="iso-8859-1"
Hi Joseph,
Your code runs without stopping for me well past 5 time units.  Is there a
time limit for jobs on your cluster?  You may also check any memory
restrictions and make sure you are not going over them.  Otherwise, I am
not able to reproduce this problem.
-Michael
On Mon, Sep 16, 2013 at 4:10 AM, Joseph Prentice <
Joseph.Prentice@physics.ox.ac.uk> wrote:
...
Hi Michael,
Unfortunately, the fix you suggested, deleting one of my quenches, only
worked for a little bit (although I had to change it a little, by making
POW 1 instead of 0 - otherwise the quench did not occur at all). As I've
increased Chi up to around 800, I've encountered the same issue as before -
namely that the program just stops at around t=5 in the real time
propagation when it should continue to t=15, with no error message at all.
The parameters of my new Python script, incorporating the suggestion from
before is below. Is there anything else I can do? I'm running on a Linux
cluster. Any help would be brilliant.
Thank you very much,
Joseph
Python script parameters:
parms = [ {
       'L'
             : 128,
       'MODEL'
     : 'spin',
       'local_S'
             : 0.5,
       'CONSERVED_QUANTUMNUMBERS'              : 'Sz',
       'Sz'
            : 0,
       'Jxy'
             : 1,
       'Jz'
            : 3.0,
       'ITP_CHIS'
    : [80, 80],
       'ITP_DTS'
     : [0.05, 0.025],
       'ITP_CONVS'                                                     :
[1E-10, 1E-11],
       'INITIAL_STATE'                                                 :
'ground',
       'CHI_LIMIT'
     : 800,
       'TRUNC_LIMIT'                                                   :
1E-12,
       'NUM_THREADS'                                           : 1,
       'TAUS'
    : [15.0],
       'POWS'
    : [1.0],
       'GS'
            : ['Jz'],
       'GIS'
     : [6.0],
       'GFS'
     : [6.0],
       'NUMSTEPS'
    : [625],
       'STEPSFORSTORE'                                         : [10],
       'VERBOSE'
     : 'true'
       } ]
On 22 Aug 2013, at 11:00, <comp-phys-alps-users-request@lists.phys.ethz.ch
...
comp-phys-alps-users-request@lists.phys.ethz.ch wrote:
...
Send Comp-phys-alps-users mailing list submissions to
     comp-phys-alps-users@lists.phys.ethz.ch
To subscribe or unsubscribe via the World Wide Web, visit
     https://lists.phys.ethz.ch/listinfo/comp-phys-alps-users
or, via email, send a message with subject or body 'help' to
     comp-phys-alps-users-request@lists.phys.ethz.ch
You can reach the person managing the list at
     comp-phys-alps-users-owner@lists.phys.ethz.ch
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Comp-phys-alps-users digest..."
Today's Topics:

Alps 2.1.2 r6963 hybridization expansion: No parameter
available (backes@th.physik.uni-frankfurt.de)
Time propagation in TEBD stops prematurely (Joseph Prentice)
Re: Alps 2.1.2 r6963 hybridization expansion: No parameter
available (Emanuel Gull)
Re: Alps 2.1.2 r6963 hybridization expansion: No parameter
available (Hartmut Hafermann)
Re: Time propagation in TEBD stops prematurely (Michael Wall)
Re: Alps 2.1.2 r6963 hybridization expansion: No parameter
available (backes@th.physik.uni-frankfurt.de)


Message: 1
Date: Wed, 21 Aug 2013 13:51:25 +0200
From: backes@th.physik.uni-frankfurt.de
To: comp-phys-alps-users@lists.phys.ethz.ch
Subject: [ALPS-users] Alps 2.1.2 r6963 hybridization expansion: No
     parameter       available
Message-ID:
     <
70b90575915bdffa61ec46f7481c0ecb.squirrel@th.physik.uni-frankfurt.de>
...
Content-Type: text/plain;charset=utf-8
Dear all,
I want to use the newest version of the ALPS Hybridization Expansion
Impurity Solver from the nightly build 2.1.2-r6963.
Compilation is successful, but when running the solver with a simple test
case I get the following Error (I shortened the paths for readability):
...
hybridization parm.h5
terminate called after throwing an instance of 'std::runtime_error'
what():  No parameter available
In [..]/alps/src/alps/ngs/detail/paramproxy.hpp on 78 in cast
[..]applications/dmft/qmc/hybridization(_ZNK4alps6detail10paramproxy4castIiEET_v+0x13d)
...
[0x42dbad]
[..]applications/dmft/qmc/hybridization(_ZN13hybridizationC1ERKN4alps6paramsEi+0x10b8)
...
[0x43f968]
[..]applications/dmft/qmc/hybridization(main+0x103) [0x41dfe3]
/lib64/libc.so.6(__libc_start_main+0xfd) [0x7ffff415fcdd]
[..]applications/dmft/qmc/hybridization() [0x41d409]
I tested a previous version 2.1.1-r6670 with identical compilation
process
...
on the same input files and it works flawlessly.
There seem to be some significant changes in this version: There is a new
hybridization2 folder in applications/dmft/qmc/ and also the
Documentation
...
inside this directory vanished. Do I have to change something in my
parm.h5 input file for this new version?
Thanks for your help.
Best wishes,
Steffen

Message: 2
Date: Wed, 21 Aug 2013 12:05:05 +0000
From: Joseph Prentice Joseph.Prentice@physics.ox.ac.uk
To: "comp-phys-alps-users@lists.phys.ethz.ch"
     comp-phys-alps-users@lists.phys.ethz.ch
Subject: [ALPS-users] Time propagation in TEBD stops prematurely
Message-ID:
     747861816F6A5C47BE50FCA07F5B2B4A0E25BB@EXCHNG16.physics.ox.ac.uk
Content-Type: text/plain; charset="us-ascii"
Hi all,
I'm having a bit of trouble with getting a quench using the XXZ model to
run in TEBD. I'm running the job on a Linux cluster, quenching Delta
instantaneously from 3 to 6, and the program stops earlier than it should.
I turned on the Verbose option, and it seems as though the imaginary time
propagation works fine, but the real time propagation stops prematurely,
with no error message. Does anyone have any idea why this has happened? It
seems very odd. I will attach my Python script in case there is something
wrong with it, and the output I obtained.
...
Thanks very much,
Joseph Prentice

Re: [ALPS-users] Time propagation in TEBD stops prematurely/Obtaining wavefunctions in TEBD