- Comp-phys-alps-users - lists.phys.ethz.ch

s-wave pair field correlation function in fullED
by Vladimir Iglovikov 09 Nov '13

09 Nov '13

[1] I am trying to compute s-wave pair field correlation function, using full diagonalization routine for the fermionic Hubbard model. This function is defined as: P_s = \langle \Delta^{\dagger} \Delta + \Delta \Delta^{\dagger} \rangle, where \Delta^{\dagger} = \frac {1} {\sqrt{N}} \sum_{k} c^{\dagger}_{n \uparrow} c^{\dagger}_{n \downarrow} And I can not figure out how to measure this. [2] How do I measure Kinetic energy? When I am trying use: 'MEASURE_AVERAGE[kinetic-energy]' : 'fermion_kinetic_energy' I get error: Cannot evaluate expression 0.333333 * fermion_kinetic_energy(i) Here is my model definition: ------------------------------ <MODELS> <SITEBASIS name="fermion"> <QUANTUMNUMBER name="Nup" min="0" max="1" type="fermionic"/> <QUANTUMNUMBER name="Ndown" min="0" max="1" type="fermionic"/> <OPERATOR name="Splus" matrixelement="1"> <CHANGE quantumnumber="Ndown" change="-1"/> <CHANGE quantumnumber="Nup" change="1"/> </OPERATOR> <OPERATOR name="Sminus" matrixelement="1"> <CHANGE quantumnumber="Nup" change="-1"/> <CHANGE quantumnumber="Ndown" change="+1"/> </OPERATOR> <OPERATOR name="Sz" matrixelement="(Nup-Ndown)/2"/> <OPERATOR name="Nup" matrixelement="Nup"/> <OPERATOR name="Ndown" matrixelement="Ndown"/> <OPERATOR name="N" matrixelement="Nup+Ndown"/> <OPERATOR name="c_down" matrixelement="1"> <CHANGE quantumnumber="Ndown" change="-1"/> </OPERATOR> <OPERATOR name="cdag_down" matrixelement="1"> <CHANGE quantumnumber="Ndown" change="1"/> </OPERATOR> <OPERATOR name="c_up" matrixelement="1"> <CHANGE quantumnumber="Nup" change="-1"/> </OPERATOR> <OPERATOR name="cdag_up" matrixelement="1"> <CHANGE quantumnumber="Nup" change="1"/> </OPERATOR> <OPERATOR name="n" matrixelement="Nup+Ndown"/> <OPERATOR name="n_up" matrixelement="Nup"/> <OPERATOR name="n_down" matrixelement="Ndown"/> <OPERATOR name="double_occupancy" matrixelement="Nup*Ndown"/> </SITEBASIS> <BONDOPERATOR name="fermion_hop" source="x" target="y"> cdag_up(x)*c_up(y)+cdag_up(y)*c_up(x)+cdag_down(x)*c_down(y)+cdag_down(y)*c_down(x) </BONDOPERATOR> <GLOBALOPERATOR name="fermion_kinetic_energy"> <BONDTERM type="0" source="i" target="j"> -t*fermion_hop(i,j) </BONDTERM> </GLOBALOPERATOR> <BASIS name="fermion"> <SITEBASIS type="1" ref="fermion"/> <SITEBASIS type="2" ref="fermion"/> <SITEBASIS type="3" ref="fermion"/> <CONSTRAINT quantumnumber="Nup" value="Nup_total"/> <CONSTRAINT quantumnumber="Ndown" value="Ndown_total"/> <CONSTRAINT quantumnumber="N" value="N_total"/> <CONSTRAINT quantumnumber="Sz" value="Sz_total"/> </BASIS> <HAMILTONIAN name="fermionic Hubbard"> <PARAMETER name="mu" default="0"/> <PARAMETER name="t" default="1"/> <PARAMETER name="U" default="0"/> <BASIS ref="fermion"/> <SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> +U# * (n_up(i) - 0.5) * (n_down(i) - 0.5) - mu# * n(i) </SITETERM> <BONDTERM source="i" target="j"> -t*fermion_hop(i, j) </BONDTERM> </HAMILTONIAN> </MODELS> ----------------------- Lattice definition: ----------------------- <LATTICES> <LATTICE name="square" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR> a 0 </VECTOR> <VECTOR> 0 a </VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0</VECTOR> <VECTOR>0 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE> <UNITCELL name="Lieb" dimension="2"> <VERTEX id="1" type="1"> <COORDINATE>0 0</COORDINATE> </VERTEX> <VERTEX id="2" type="2"> <COORDINATE>0.5 0</COORDINATE> </VERTEX> <VERTEX id="3" type="3"> <COORDINATE>0 0.5</COORDINATE> </VERTEX> <EDGE> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 0"/> </EDGE> <EDGE> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="3" offset="0 0"/> </EDGE> <EDGE> <SOURCE vertex="1" offset="1 0"/> <TARGET vertex="2" offset="0 0"/> </EDGE> <EDGE> <SOURCE vertex="1" offset="0 1"/> <TARGET vertex="3" offset="0 0"/> </EDGE> </UNITCELL> <LATTICEGRAPH name="Lieb" vt_skip="true"> <FINITELATTICE> <LATTICE ref="square"/> <EXTENT dimension="1" size="nx"/> <EXTENT dimension="2" size="ny"/> <BOUNDARY type="periodic" /> </FINITELATTICE> <UNITCELL ref="Lieb"/> </LATTICEGRAPH> </LATTICES>

2 1

dmrg with complex hopping amplitudes
by benoît grémaud 28 Oct '13

28 Oct '13

Dear All, I'm using DMRG (bose-hubbard) with complex tunneling amplitudes, the hamiltonian being hermitian, of course. The results in simple situation is reasonable, i.e. the energy and local density are invariant by local U(1) transformation. In addition, I am measuring the ODLRO, i.e. I have the following line in the input file: 'MEASURE_CORRELATIONS[ODLRO]' : "bdag:b", However, the results are always real numbers, whereas I would expect complex numbers. Do they correspond to the real part of the correlation function? Is there a way to get the imaginary part? thank you for you help. my version of alps is: alps-2.1.1-r6176 regards, Benoît Grémaud -- -------------------- Researcher at CNRS, Laboratoire Kastler Brossel (Paris, France) Visiting Associate Professor at the Centre for Quantum Technologies (NUS, Singapore)

1 0

Error Installation in Ubuntu 12.04 Lts(64 bit)
by Javad Vahedi 26 Oct '13

26 Oct '13

Dear Michele Many thanks for your consideration. Heres my output ######################################################## maha@maha-pc:~/Downloads/build$ cmake . -- Compiler version: c++ (Ubuntu/Linaro 4.6.3-1ubuntu5) 4.6.3 -- Build type: Release -- Python interpreter /usr/bin/python -- Python interpreter ok : version 2.7.3 -- PYTHON_INCLUDE_DIRS = /usr/include/python2.7 -- PYTHON_NUMPY_INCLUDE_DIR = /usr/lib/python2.7/dist -packages/numpy/core/include -- PYTHON_SITE_PKG = /usr/lib/python2.7/dist-packages -- PYTHON_LIBRARY = /usr/lib/python2.7/config/libpython2.7.so -- PYTHON_EXTRA_LIBS =-lssl -lcrypto -lssl -lcrypto -L/usr/lib -lz -lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -Wl,-O1 -Wl, -Bsymbolic-functions -- ALPS version: 2.2.b1 -- Looking for Boost Source -- Found Boost Source: /home/maha/Downloads/alps-2.2.b1-r7195-src -with-boost/boost -- Boost Version: 1_54_0 -- Adding Boost dir: /home/maha/Downloads/alps-2.2.b1-r7195-src -with-boost/boost -- MPI compiler was /usr/bin/mpicxx -- LAPACK DEBUG::Compiler id GNU -- LAPACK DEBUG::Compiler version -- LAPACK DEBUG::Processor x86_64 -- LAPACK DEBUG::ENV{MKL} -- LAPACK DEBUG::ENV{MKLROOT} /opt/intel/composer_xe_2013.5.192/mkl -- LAPACK DEBUG::ENV{MKL_HOME} -- LAPACK DEBUG::MKL_VERSION (raw) 2013.5.192 -- LAPACK DEBUG::MKL_VERSION 10.3.1 -- Found intel/mkl library -- LAPACK DEBUG::LAPACK_LIBRARY = -L/opt/intel/composer_xe_2013.5.192/mkl/lib/intel64;-lmkl_intel_lp64;-lmkl_sequential;-lmkl_core;-lpthread;-lm -- SQLite Library: not found -- LPSolve Library: not found -- Could NOT find SZIP (missing: SZIP_LIBRARIES SZIP_INCLUDE_DIRS) -- HDF5: /usr/include debug;/usr/lib/libhdf5_hl.so;/usr/lib/libhdf5.so;optimized;/usr/lib/libhdf5_hl.so;/usr/lib/libhdf5.so HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND -- Python interpreter /usr/bin/python -- Python interpreter ok : version 2.7.3 -- PYTHON_INCLUDE_DIRS = /usr/include/python2.7 -- PYTHON_NUMPY_INCLUDE_DIR = /usr/lib/python2.7/dist -packages/numpy/core/include -- PYTHON_SITE_PKG = /usr/lib/python2.7/dist-packages -- PYTHON_LIBRARY = /usr/lib/python2.7/config/libpython2.7.so -- PYTHON_EXTRA_LIBS =-lssl -lcrypto -lssl -lcrypto -L/usr/lib -lz -lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -Wl,-O1 -Wl, -Bsymbolic-functions -- Numpy include in /usr/lib/python2.7/dist-packages/numpy/core/include -- ALPS XML dir is /opt/alps/lib/xml -- HDF5: /usr/include debug;/usr/lib/libhdf5_hl.so;/usr/lib/libhdf5.so;optimized;/usr/lib/libhdf5_hl.so;/usr/lib/libhdf5.so HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND;HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND -- HDF5: /usr/include debug;/usr/lib/libhdf5_hl.so;/usr/lib/libhdf5.so;optimized;/usr/lib/libhdf5_hl.so;/usr/lib/libhdf5.so HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND;HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND -- tebd will not be built -- HDF5: /usr/include debug;/usr/lib/libhdf5_hl.so;/usr/lib/libhdf5.so;optimized;/usr/lib/libhdf5_hl.so;/usr/lib/libhdf5.so HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND;HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND -- Configuring done -- Generating done -- Build files have been written to: /home/maha/Downloads/build #################################################################### Your sincerely Javad

1 0

dmrg with complex tunneling
by benoît grémaud 26 Oct '13

26 Oct '13

Dear All, I'm using DMRG (bose-hubbard) with complex tunneling amplitudes, the hamiltonian being hermitian, of course. The results in simple situation is reasonable, i.e. the energy and local density are invariant by local U(1) transformation. In addition, I am measuring the ODLRO, i.e. I have the following line in the input file: 'MEASURE_CORRELATIONS[ODLRO]' : "bdag:b", However, the results are always real numbers, whereas I would expect complex numbers. Do they correspond to the real part of the correlation function? Is there a way to get the imaginary part? thank you for you help. my version of alps is: alps-2.1.1-r6176 regards, Benoît Grémaud -- -------------------- Researcher at CNRS, Laboratoire Kastler Brossel (Paris, France) Visiting Associate Professor at the Centre for Quantum Technologies (NUS, Singapore)

1 0

Error Installation in Ubuntu 12.04 Lts(64 bit)
by Javad Vahedi 25 Oct '13

25 Oct '13

Dear ALPS user I have been trying to install 'alps-2.2.b1-r7195-src-with-boost'' version on my pc. I went through the ALPS building from source page prescription, but unfortunatly something goes wrong and I don't know how to fix it. Any help will be kindly appreciated After running "make" I faced the following error: . . . . Linking CXX executable jacobidavidson1 /opt/intel/composer_xe_2013.5.192/mkl/lib/intel64/libmkl_core.so: undefined reference to `pthread_mutex_trylock' collect2: ld returned 1 exit status make[2]: *** [example/ietl/jacobidavidson1] Error 1 make[1]: *** [example/ietl/CMakeFiles/jacobidavidson1.dir/all] Error 2 make: *** [all] Error 2 ####################################################### Then after running " make test " maha@maha-pc:~/Downloads/build$ make test Running tests... Test project /home/maha/Downloads/build Start 1: ising_single 1/141 Test #1: ising_single ......................... Passed 0.23 sec Start 2: ising 2/141 Test #2: ising ................................ Passed 0.42 sec Start 3: heisenberg 3/141 Test #3: heisenberg ........................... Passed 0.57 sec Start 4: loop_single 4/141 Test #4: loop_single .......................... Passed 0.14 sec Start 5: exchange 5/141 Test #5: exchange ............................. Passed 0.22 sec Start 6: wanglandau 6/141 Test #6: wanglandau ........................... Passed 0.77 sec Start 7: enum_as_class 7/141 Test #7: enum_as_class ........................***Failed 0.02 sec Start 8: enum_vectorizable 8/141 Test #8: enum_vectorizable ....................***Failed 0.02 sec Start 9: pair_int_vectorizable 9/141 Test #9: pair_int_vectorizable ................***Failed 0.02 sec Start 10: complexobservable 10/141 Test #10: complexobservable ....................***Failed 0.01 sec Start 11: detailedbinning 11/141 Test #11: detailedbinning ......................***Failed 0.01 sec Start 12: histogram 12/141 Test #12: histogram ............................***Failed 0.01 sec Start 13: signed 13/141 Test #13: signed ...............................***Failed 0.02 sec Start 14: simpleobseval 14/141 Test #14: simpleobseval ........................***Failed 0.01 sec Start 15: testobservableset 15/141 Test #15: testobservableset ....................***Failed 0.01 sec Start 16: observableset_hdf5 16/141 Test #16: observableset_hdf5 ...................***Failed 0.01 sec Start 17: observableset_xml 17/141 Test #17: observableset_xml ....................***Failed 0.02 sec Start 18: vectorobseval 18/141 Test #18: vectorobseval ........................***Failed 0.02 sec Start 19: mcdata 19/141 Test #19: mcdata ...............................***Failed 0.02 sec Start 20: mcdata2 20/141 Test #20: mcdata2 ..............................***Failed 0.02 sec Start 21: mcanalyze 21/141 Test #21: mcanalyze ............................***Failed 0.02 sec Start 22: fixed_capacity_deque 22/141 Test #22: fixed_capacity_deque .................***Failed 0.02 sec Start 23: fixed_capacity_traits 23/141 Test #23: fixed_capacity_traits ................***Failed 0.02 sec Start 24: fixed_capacity_vector 24/141 Test #24: fixed_capacity_vector ................***Failed 0.02 sec Start 25: test_deque 25/141 Test #25: test_deque ...........................***Failed 0.01 sec Start 26: test_vector 26/141 Test #26: test_vector ..........................***Failed 0.01 sec Start 27: hdf5_complex 27/141 Test #27: hdf5_complex .........................***Failed 0.02 sec Start 28: hdf5_copy 28/141 Test #28: hdf5_copy ............................***Failed 0.02 sec Start 29: hdf5_real_complex_vec 29/141 Test #29: hdf5_real_complex_vec ................***Failed 0.01 sec Start 30: hdf5_real_complex_matrix 30/141 Test #30: hdf5_real_complex_matrix .............***Failed 0.01 sec Start 31: hdf5_bool 31/141 Test #31: hdf5_bool ............................***Failed 0.01 sec Start 32: hdf5_replace 32/141 Test #32: hdf5_replace .........................***Failed 0.02 sec Start 33: hdf5_pair 33/141 Test #33: hdf5_pair ............................***Failed 0.01 sec Start 34: hdf5_parms 34/141 Test #34: hdf5_parms ...........................***Failed 0.02 sec Start 35: hdf5_ising 35/141 Test #35: hdf5_ising ...........................***Failed 0.02 sec Start 36: hdf5_family 36/141 Test #36: hdf5_family ..........................***Failed 0.02 sec Start 37: hdf5_valgrind 37/141 Test #37: hdf5_valgrind ........................***Failed 0.02 sec Start 38: hdf5_observableset 38/141 Test #38: hdf5_observableset ...................***Failed 0.02 sec Start 39: hdf5_multi_array 39/141 Test #39: hdf5_multi_array .....................***Failed 0.02 sec Start 40: hdf5_memory 40/141 Test #40: hdf5_memory ..........................***Failed 0.02 sec Start 41: hdf5_misc 41/141 Test #41: hdf5_misc ............................***Failed 0.02 sec Start 42: hdf5_vecveccplx 42/141 Test #42: hdf5_vecveccplx ......................***Failed 0.01 sec Start 43: hdf5_exceptions 43/141 Test #43: hdf5_exceptions ......................***Failed 0.02 sec Start 44: hdf5_multiarchive 44/141 Test #44: hdf5_multiarchive ....................***Failed 0.02 sec Start 45: hdf5_fortran_string 45/141 Test #45: hdf5_fortran_string ..................***Failed 0.02 sec Start 46: lattice_example1 46/141 Test #46: lattice_example1 .....................***Failed 0.01 sec Start 47: lattice_example2 47/141 Test #47: lattice_example2 .....................***Failed 0.02 sec Start 48: lattice_example3 48/141 Test #48: lattice_example3 .....................***Failed 0.02 sec Start 49: lattice_example4 49/141 Test #49: lattice_example4 .....................***Failed 0.02 sec Start 50: lattice_example5 50/141 Test #50: lattice_example5 .....................***Failed 0.01 sec Start 51: lattice_example6 51/141 Test #51: lattice_example6 .....................***Failed 0.01 sec Start 52: lattice_example7 52/141 Test #52: lattice_example7 .....................***Failed 0.01 sec Start 53: lattice_example8 53/141 Test #53: lattice_example8 .....................***Failed 0.02 sec Start 54: lattice_example9 54/141 Test #54: lattice_example9 .....................***Failed 0.02 sec Start 55: lattice_example10 55/141 Test #55: lattice_example10 ....................***Failed 0.02 sec Start 56: lattice_example11 56/141 Test #56: lattice_example11 ....................***Failed 0.01 sec Start 57: lattice_label 57/141 Test #57: lattice_label ........................***Failed 0.01 sec Start 58: coloring 58/141 Test #58: coloring .............................***Failed 0.01 sec Start 59: parity 59/141 Test #59: parity ...............................***Failed 0.02 sec Start 60: model_example1 60/141 Test #60: model_example1 .......................***Failed 0.02 sec Start 61: model_example2 61/141 Test #61: model_example2 .......................***Failed 0.02 sec Start 62: model_example3 62/141 Test #62: model_example3 .......................***Failed 0.01 sec Start 63: model_example4 63/141 Test #63: model_example4 .......................***Failed 0.01 sec Start 64: model_example5 64/141 Test #64: model_example5 .......................***Failed 0.02 sec Start 65: model_example6 65/141 Test #65: model_example6 .......................***Failed 0.02 sec Start 66: model_example7 66/141 Test #66: model_example7 .......................***Failed 0.01 sec Start 67: model_example8 67/141 Test #67: model_example8 .......................***Failed 0.01 sec Start 68: model_example9 68/141 Test #68: model_example9 .......................***Failed 0.02 sec Start 69: model_example10 69/141 Test #69: model_example10 ......................***Failed 0.01 sec Start 70: model_example11 70/141 Test #70: model_example11 ......................***Failed 0.01 sec Start 71: model_example12 71/141 Test #71: model_example12 ......................***Failed 0.01 sec Start 72: model_example13 72/141 Test #72: model_example13 ......................***Failed 0.01 sec Start 73: model_example14 73/141 Test #73: model_example14 ......................***Failed 0.02 sec Start 74: model_example15 74/141 Test #74: model_example15 ......................***Failed 0.01 sec Start 75: model_example16 75/141 Test #75: model_example16 ......................***Failed 0.01 sec Start 76: model_example17 76/141 Test #76: model_example17 ......................***Failed 0.01 sec Start 77: model_example18 77/141 Test #77: model_example18 ......................***Failed 0.02 sec Start 78: ngs_hdf5 78/141 Test #78: ngs_hdf5 .............................***Failed 0.02 sec Start 79: ngs_hash 79/141 Test #79: ngs_hash .............................***Failed 0.02 sec Start 80: param_default 80/141 Test #80: param_default ........................***Failed 0.02 sec Start 81: param_ordering 81/141 Test #81: param_ordering .......................***Failed 0.01 sec Start 82: param_stream 82/141 Test #82: param_stream .........................***Failed 0.02 sec Start 83: param_assign 83/141 Test #83: param_assign .........................***Failed 0.02 sec Start 84: param_not_found 84/141 Test #84: param_not_found ......................***Failed 0.02 sec Start 85: matrix_unit_tests 85/141 Test #85: matrix_unit_tests ....................***Failed 0.02 sec Start 86: matrix_column_view 86/141 Test #86: matrix_column_view ...................***Failed 0.02 sec Start 87: matrix_transpose_view 87/141 Test #87: matrix_transpose_view ................***Failed 0.02 sec Start 88: matrix_hdf5 88/141 Test #88: matrix_hdf5 ..........................***Failed 0.02 sec Start 89: matrix_deprecated_hdf5_format_test 89/141 Test #89: matrix_deprecated_hdf5_format_test ...***Failed 0.02 sec Start 90: matrix_ublas_sparse_functions 90/141 Test #90: matrix_ublas_sparse_functions ........***Failed 0.01 sec Start 91: accumulate_if 91/141 Test #91: accumulate_if ........................***Failed 0.02 sec Start 92: vector_functions 92/141 Test #92: vector_functions .....................***Failed 0.02 sec Start 93: vector_valarray_conversion 93/141 Test #93: vector_valarray_conversion ...........***Failed 0.02 sec Start 94: matrix_algorithms 94/141 Test #94: matrix_algorithms ....................***Failed 0.02 sec Start 95: real_tests 95/141 Test #95: real_tests ...........................***Failed 0.01 sec Start 96: os 96/141 Test #96: os ...................................***Failed 0.02 sec Start 97: sizeof 97/141 Test #97: sizeof ...............................***Failed 0.02 sec Start 98: xdrdump 98/141 Test #98: xdrdump ..............................***Failed 0.02 sec Start 99: xdrdump2 99/141 Test #99: xdrdump2 .............................***Failed 0.02 sec Start 100: boostdump 100/141 Test #100: boostdump ............................***Failed 0.02 sec Start 101: boostdump2 101/141 Test #101: boostdump2 ...........................***Failed 0.02 sec Start 102: boostdump3 102/141 Test #102: boostdump3 ...........................***Failed 0.02 sec Start 103: boostdump4 103/141 Test #103: boostdump4 ...........................***Failed 0.01 sec Start 104: expression 104/141 Test #104: expression ...........................***Failed 0.02 sec Start 105: expression2 105/141 Test #105: expression2 ..........................***Failed 0.02 sec Start 106: flatten 106/141 Test #106: flatten ..............................***Failed 0.01 sec Start 107: parameter 107/141 Test #107: parameter ............................***Failed 0.02 sec Start 108: parameterlist 108/141 Test #108: parameterlist ........................***Failed 0.02 sec Start 109: parameterlist_xml 109/141 Test #109: parameterlist_xml ....................***Failed 0.02 sec Start 110: parameters 110/141 Test #110: parameters ...........................***Failed 0.02 sec Start 111: parameters_xml 111/141 Test #111: parameters_xml .......................***Failed 0.02 sec Start 112: parameters_hdf5 112/141 Test #112: parameters_hdf5 ......................***Failed 0.02 sec Start 113: clone_info 113/141 Test #113: clone_info ...........................***Failed 0.01 sec Start 114: clone_phase 114/141 Test #114: clone_phase ..........................***Failed 0.01 sec Start 115: clone_timer 115/141 Test #115: clone_timer ..........................***Failed 0.01 sec Start 116: exmc_optimize 116/141 Test #116: exmc_optimize ........................***Failed 0.01 sec Start 117: exp_number 117/141 Test #117: exp_number ...........................***Failed 0.02 sec Start 118: footprint 118/141 Test #118: footprint ............................***Failed 0.01 sec Start 119: info_test 119/141 Test #119: info_test ............................***Failed 0.01 sec Start 120: integer_range 120/141 Test #120: integer_range ........................***Failed 0.02 sec Start 121: linear_regression 121/141 Test #121: linear_regression ....................***Failed 0.02 sec Start 122: merge 122/141 Test #122: merge ................................***Failed 0.02 sec Start 123: percentage 123/141 Test #123: percentage ...........................***Failed 0.02 sec Start 124: temperature_scan 124/141 Test #124: temperature_scan .....................***Failed 0.02 sec Start 125: time 125/141 Test #125: time .................................***Failed 0.03 sec Start 126: version 126/141 Test #126: version ..............................***Failed 0.02 sec Start 127: wl_weight 127/141 Test #127: wl_weight ............................***Failed 0.02 sec Start 128: xmlhandler 128/141 Test #128: xmlhandler ...........................***Failed 0.01 sec Start 129: xmlparser 129/141 Test #129: xmlparser ............................***Failed 0.02 sec Start 130: xmlstream 130/141 Test #130: xmlstream ............................***Failed 0.02 sec Start 131: python_pyioarchive 131/141 Test #131: python_pyioarchive ................... Passed 2.94 sec Start 132: python_pyhdf5io 132/141 Test #132: python_pyhdf5io ...................... Passed 0.67 sec Start 133: python_pyparams 133/141 Test #133: python_pyparams ...................... Passed 0.61 sec Start 134: python_hlist_test 134/141 Test #134: python_hlist_test .................... Passed 0.58 sec Start 135: python_mcdata 135/141 Test #135: python_mcdata ........................ Passed 0.58 sec Start 136: python_pyhdf5 136/141 Test #136: python_pyhdf5 ........................ Passed 0.62 sec Start 137: python_mcanalyze 137/141 Test #137: python_mcanalyze ..................... Passed 1.22 sec Start 138: random_choice 138/141 Test #138: random_choice ........................***Failed 0.02 sec Start 139: uniform_on_sphere_n 139/141 Test #139: uniform_on_sphere_n ..................***Failed 0.02 sec Start 140: vmusage 140/141 Test #140: vmusage ..............................***Failed 0.02 sec Start 141: bitops 141/141 Test #141: bitops ...............................***Failed 0.02 sec 9% tests passed, 128 tests failed out of 141 Label Time Summary: fixed_capacity = 0.08 sec hdf5 = 0.30 sec lattice = 0.21 sec model = 0.27 sec numeric = 0.18 sec osiris = 0.13 sec parameter = 0.14 sec parapack = 0.24 sec parser = 0.05 sec pyalps = 7.21 sec random = 0.03 sec utility = 0.03 sec Total Test time (real) = 11.65 sec The following tests FAILED: 7 - enum_as_class (Failed) 8 - enum_vectorizable (Failed) 9 - pair_int_vectorizable (Failed) 10 - complexobservable (Failed) 11 - detailedbinning (Failed) 12 - histogram (Failed) 13 - signed (Failed) 14 - simpleobseval (Failed) 15 - testobservableset (Failed) 16 - observableset_hdf5 (Failed) 17 - observableset_xml (Failed) 18 - vectorobseval (Failed) 19 - mcdata (Failed) 20 - mcdata2 (Failed) 21 - mcanalyze (Failed) 22 - fixed_capacity_deque (Failed) 23 - fixed_capacity_traits (Failed) 24 - fixed_capacity_vector (Failed) 25 - test_deque (Failed) 26 - test_vector (Failed) 27 - hdf5_complex (Failed) 28 - hdf5_copy (Failed) 29 - hdf5_real_complex_vec (Failed) 30 - hdf5_real_complex_matrix (Failed) 31 - hdf5_bool (Failed) 32 - hdf5_replace (Failed) 33 - hdf5_pair (Failed) 34 - hdf5_parms (Failed) 35 - hdf5_ising (Failed) 36 - hdf5_family (Failed) 37 - hdf5_valgrind (Failed) 38 - hdf5_observableset (Failed) 39 - hdf5_multi_array (Failed) 40 - hdf5_memory (Failed) 41 - hdf5_misc (Failed) 42 - hdf5_vecveccplx (Failed) 43 - hdf5_exceptions (Failed) 44 - hdf5_multiarchive (Failed) 45 - hdf5_fortran_string (Failed) 46 - lattice_example1 (Failed) 47 - lattice_example2 (Failed) 48 - lattice_example3 (Failed) 49 - lattice_example4 (Failed) 50 - lattice_example5 (Failed) 51 - lattice_example6 (Failed) 52 - lattice_example7 (Failed) 53 - lattice_example8 (Failed) 54 - lattice_example9 (Failed) 55 - lattice_example10 (Failed) 56 - lattice_example11 (Failed) 57 - lattice_label (Failed) 58 - coloring (Failed) 59 - parity (Failed) 60 - model_example1 (Failed) 61 - model_example2 (Failed) 62 - model_example3 (Failed) 63 - model_example4 (Failed) 64 - model_example5 (Failed) 65 - model_example6 (Failed) 66 - model_example7 (Failed) 67 - model_example8 (Failed) 68 - model_example9 (Failed) 69 - model_example10 (Failed) 70 - model_example11 (Failed) 71 - model_example12 (Failed) 72 - model_example13 (Failed) 73 - model_example14 (Failed) 74 - model_example15 (Failed) 75 - model_example16 (Failed) 76 - model_example17 (Failed) 77 - model_example18 (Failed) 78 - ngs_hdf5 (Failed) 79 - ngs_hash (Failed) 80 - param_default (Failed) 81 - param_ordering (Failed) 82 - param_stream (Failed) 83 - param_assign (Failed) 84 - param_not_found (Failed) 85 - matrix_unit_tests (Failed) 86 - matrix_column_view (Failed) 87 - matrix_transpose_view (Failed) 88 - matrix_hdf5 (Failed) 89 - matrix_deprecated_hdf5_format_test (Failed) 90 - matrix_ublas_sparse_functions (Failed) 91 - accumulate_if (Failed) 92 - vector_functions (Failed) 93 - vector_valarray_conversion (Failed) 94 - matrix_algorithms (Failed) 95 - real_tests (Failed) 96 - os (Failed) 97 - sizeof (Failed) 98 - xdrdump (Failed) 99 - xdrdump2 (Failed) 100 - boostdump (Failed) 101 - boostdump2 (Failed) 102 - boostdump3 (Failed) 103 - boostdump4 (Failed) 104 - expression (Failed) 105 - expression2 (Failed) 106 - flatten (Failed) 107 - parameter (Failed) 108 - parameterlist (Failed) 109 - parameterlist_xml (Failed) 110 - parameters (Failed) 111 - parameters_xml (Failed) 112 - parameters_hdf5 (Failed) 113 - clone_info (Failed) 114 - clone_phase (Failed) 115 - clone_timer (Failed) 116 - exmc_optimize (Failed) 117 - exp_number (Failed) 118 - footprint (Failed) 119 - info_test (Failed) 120 - integer_range (Failed) 121 - linear_regression (Failed) 122 - merge (Failed) 123 - percentage (Failed) 124 - temperature_scan (Failed) 125 - time (Failed) 126 - version (Failed) 127 - wl_weight (Failed) 128 - xmlhandler (Failed) 129 - xmlparser (Failed) 130 - xmlstream (Failed) 138 - random_choice (Failed) 139 - uniform_on_sphere_n (Failed) 140 - vmusage (Failed) 141 - bitops (Failed) Errors while running CTest make: *** [test] Error 8

2 1

alps installation on ubuntu 13.10
by Rafał 23 Oct '13

23 Oct '13

Hi all, For a long time I used alps on ubuntu 12.04 and I had no problem with the installation. Right know I wanted to try installing alps on ubuntu 13.10 and I encounter an problem that I can't solve. I have extracted alps and boost (from alps-2.1.1-r6176-src-with-boost.tar.gz) into separated folders in ~/src/. I made an folder for building. I installed dependencies with sudo apt-get install cmake-curses-gui libboost-all-dev libhdf5-serial-dev libfftw3-dev gfortran python-matplotlib python-scipy python-h5py liblapack-dev xsltproc python-dev I found some problem with hdf5 when building alps, that does not existed in ubuntu 12.04, so I installed hdf5 from sources like in: http://www.hdfgroup.org/ftp/HDF5/current/src/unpacked/release_docs/INSTALL adding --enable-fortran configure option. then I tried to execute in the build folder following command: sudo cmake -D BLAS_LIBRARY=/usr/lib/libblas.so -D LAPACK_LIBRARY=/usr/lib/liblapack.so -D PYTHON_LIBRARY=/usr/lib/ libpython2.7.so -D Boost_ROOT_DIR:PATH=/home/arymanus/src/alps2/boost HDF5_LIBRARIES=/usr/local/hdf5/lib -D HDF5_INCLUDE_DIR=/usr/local/hdf5/include/ /home/arymanus/src/alps2/alps and result is below. -- Build type: Release -- Python interpreter /usr/bin/python -- Python interpreter ok : version 2.7.5+ -- PYTHON_INCLUDE_DIRS = /usr/include/python2.7 -- PYTHON_NUMPY_INCLUDE_DIR = /usr/lib/python2.7/dist-packages/numpy/core/include -- PYTHON_SITE_PKG = /usr/lib/python2.7/dist-packages -- PYTHON_LIBRARY = /usr/lib/libpython2.7.so -- PYTHON_EXTRA_LIBS =-L/usr/lib -lz -lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -Wl,-O1 -Wl,-Bsymbolic-functions -- ALPS version: 2.1.1 -- Looking for Boost Source -- Found Boost Source: /home/arymanus/src/alps2/boost -- Boost Version: 1_49_0 -- MPI compiler was /usr/bin/mpicxx -- SQLite Library: not found -- LPSolve Library: not found -- Could NOT find SZIP (missing: SZIP_LIBRARIES SZIP_INCLUDE_DIRS) -- Found HDF5: debug;/usr/lib/x86_64-linux-gnu/libhdf5_hl.so;/usr/lib/x86_64-linux-gnu/libhdf5.so;/usr/lib/x86_64-linux-gnu/libpthread.so;/usr/lib/x86_64-linux-gnu/libz.so;/usr/lib/x86_64-linux-gnu/libdl.so;/usr/lib/x86_64-linux-gnu/libm.so;optimized;/usr/lib/x86_64-linux-gnu/libhdf5_hl.so;/usr/lib/x86_64-linux-gnu/libhdf5.so;/usr/lib/x86_64-linux-gnu/libpthread.so;/usr/lib/x86_64-linux-gnu/libz.so;/usr/lib/x86_64-linux-gnu/libdl.so;/usr/lib/x86_64-linux-gnu/libm.so -- HDF5: /usr/local/hdf5/include/;/usr/include debug;/usr/lib/x86_64-linux-gnu/libhdf5_hl.so;/usr/lib/x86_64-linux-gnu/libhdf5.so;/usr/lib/x86_64-linux-gnu/libpthread.so;/usr/lib/x86_64-linux-gnu/libz.so;/usr/lib/x86_64-linux-gnu/libdl.so;/usr/lib/x86_64-linux-gnu/libm.so;optimized;/usr/lib/x86_64-linux-gnu/libhdf5_hl.so;/usr/lib/x86_64-linux-gnu/libhdf5.so;/usr/lib/x86_64-linux-gnu/libpthread.so;/usr/lib/x86_64-linux-gnu/libz.so;/usr/lib/x86_64-linux-gnu/libdl.so;/usr/lib/x86_64-linux-gnu/libm.so HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND -- Python interpreter /usr/bin/python -- Python interpreter ok : version 2.7.5+ -- PYTHON_INCLUDE_DIRS = /usr/include/python2.7 -- PYTHON_NUMPY_INCLUDE_DIR = /usr/lib/python2.7/dist-packages/numpy/core/include -- PYTHON_SITE_PKG = /usr/lib/python2.7/dist-packages -- PYTHON_LIBRARY = /usr/lib/libpython2.7.so -- PYTHON_EXTRA_LIBS =-L/usr/lib -lz -lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -Wl,-O1 -Wl,-Bsymbolic-functions -- Numpy include in /usr/lib/python2.7/dist-packages/numpy/core/include -- ALPS XML dir is /opt/alps/lib/xml -- Adding Boost dir: /home/arymanus/src/alps2/boost -- HDF5: /usr/local/hdf5/include/;/usr/include debug;/usr/lib/x86_64-linux-gnu/libhdf5_hl.so;/usr/lib/x86_64-linux-gnu/libhdf5.so;/usr/lib/x86_64-linux-gnu/libpthread.so;/usr/lib/x86_64-linux-gnu/libz.so;/usr/lib/x86_64-linux-gnu/libdl.so;/usr/lib/x86_64-linux-gnu/libm.so;optimized;/usr/lib/x86_64-linux-gnu/libhdf5_hl.so;/usr/lib/x86_64-linux-gnu/libhdf5.so;/usr/lib/x86_64-linux-gnu/libpthread.so;/usr/lib/x86_64-linux-gnu/libz.so;/usr/lib/x86_64-linux-gnu/libdl.so;/usr/lib/x86_64-linux-gnu/libm.so HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND;HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND -- HDF5: /usr/local/hdf5/include/;/usr/include debug;/usr/lib/x86_64-linux-gnu/libhdf5_hl.so;/usr/lib/x86_64-linux-gnu/libhdf5.so;/usr/lib/x86_64-linux-gnu/libpthread.so;/usr/lib/x86_64-linux-gnu/libz.so;/usr/lib/x86_64-linux-gnu/libdl.so;/usr/lib/x86_64-linux-gnu/libm.so;optimized;/usr/lib/x86_64-linux-gnu/libhdf5_hl.so;/usr/lib/x86_64-linux-gnu/libhdf5.so;/usr/lib/x86_64-linux-gnu/libpthread.so;/usr/lib/x86_64-linux-gnu/libz.so;/usr/lib/x86_64-linux-gnu/libdl.so;/usr/lib/x86_64-linux-gnu/libm.so HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND;HDF5_DLL-NOTFOUND;HDF5_HL_DLL-NOTFOUND -- tebd will not be built -- Configuring done -- Generating done -- Build files have been written to: /home/arymanus/src I have no idea how to solve this problem I will be very glad for your help. Best regards, Rafał Skolasinski

2 4

Why the default confidence interval 63%
by Ryo IGARASHI 01 Oct '13

01 Oct '13

Dear ALPS developers, According to the new Tutorial(MC-01)[1], the default configuration interval is set to 0.63. Why do you choose this value as default? If this value is meant to be inside $\pm \sigma$, it should be around 0.6827 and not 0.63. Is there a specific reason you choose this value? [1] https://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:MC-01_Equil… Best regards, -- Ryo IGARASHI, Ph.D. rigarash(a)issp.u-tokyo.ac.jp OpenPGP fingerprint: BAD9 71E3 28F3 8952 5640 6A53 EC79 A280 6A19 2319

2 1

Complex site- and bond-operators
by pradeep thakur 30 Sep '13

30 Sep '13

Dear All, How can I introduce "i" (square-root of -1) into the model library (models.xml file)? I'm unable to define new site- and bond-operators which involve "i". Any help will be very much appreciated. With regards, Pradeep Thakur Dept. of Physics., University of Pune, Pune, Maharashtra, INDIA.

3 13

error bars in dirloop_sse code
by Tao Ying 28 Sep '13

28 Sep '13

Dear developers, I'm using the ALPS dirloop_sse code to simulate the hard core boson. I have one question about the error bars from the code. How the error bars are obtained in this code? Does it assume the autocorrelation time is short? I mean, is it assuming every measurement is statistically independent? Or is the code doing different binnings? Best, TaoYing -- Tao Ying, Ph.D candidate, Physics department Harbin Institute of Technology

2 1

Re: [ALPS-users] Time propagation in TEBD stops prematurely/Obtaining wavefunctions in TEBD
by Joseph Prentice 19 Sep '13

19 Sep '13

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(a)lists.phys.ethz.ch> <comp-phys-alps-users-request(a)lists.phys.ethz.ch> wrote: > Send Comp-phys-alps-users mailing list submissions to > comp-phys-alps-users(a)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(a)lists.phys.ethz.ch > > You can reach the person managing the list at > comp-phys-alps-users-owner(a)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: > > 1. Re: Time propagation in TEBD stops prematurely (Michael Wall) > 2. Re: Making the scheduler switch tasks, obtaining MPI rank > (Peter Br?cker) > 3. 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(a)gmail.com> > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Subject: Re: [ALPS-users] Time propagation in TEBD stops prematurely > Message-ID: > <CA+DwVMrz2Y_1Kz7yDSyEqALjskd0iavzFM_9_KTMjPiNzkSzJA(a)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(a)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(a)lists.phys.ethz.ch >>> >> <comp-phys-alps-users-request(a)lists.phys.ethz.ch> wrote: >> >>> Send Comp-phys-alps-users mailing list submissions to >>> comp-phys-alps-users(a)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(a)lists.phys.ethz.ch >>> >>> You can reach the person managing the list at >>> comp-phys-alps-users-owner(a)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: >>> >>> 1. Alps 2.1.2 r6963 hybridization expansion: No parameter >>> available (backes(a)th.physik.uni-frankfurt.de) >>> 2. Time propagation in TEBD stops prematurely (Joseph Prentice) >>> 3. Re: Alps 2.1.2 r6963 hybridization expansion: No parameter >>> available (Emanuel Gull) >>> 4. Re: Alps 2.1.2 r6963 hybridization expansion: No parameter >>> available (Hartmut Hafermann) >>> 5. Re: Time propagation in TEBD stops prematurely (Michael Wall) >>> 6. Re: Alps 2.1.2 r6963 hybridization expansion: No parameter >>> available (backes(a)th.physik.uni-frankfurt.de) >>> >>> >>> ---------------------------------------------------------------------- >>> >>> Message: 1 >>> Date: Wed, 21 Aug 2013 13:51:25 +0200 >>> From: backes(a)th.physik.uni-frankfurt.de >>> To: comp-phys-alps-users(a)lists.phys.ethz.ch >>> Subject: [ALPS-users] Alps 2.1.2 r6963 hybridization expansion: No >>> parameter available >>> Message-ID: >>> < >> 70b90575915bdffa61ec46f7481c0ecb.squirrel(a)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(a)physics.ox.ac.uk> >>> To: "comp-phys-alps-users(a)lists.phys.ethz.ch" >>> <comp-phys-alps-users(a)lists.phys.ethz.ch> >>> Subject: [ALPS-users] Time propagation in TEBD stops prematurely >>> Message-ID: >>> <747861816F6A5C47BE50FCA07F5B2B4A0E25BB(a)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 >>>

2 1