Dear Wei,
Thanks for the report. Unfortunately I will only be able to look into that in a few weeks. In the meantime I hope you could still obtain your results with the twosite algorithm.
Best regards, Michele
-- ETH Zurich Michele Dolfi Institute for Theoretical Physics HIT G 32.4 Wolfgang-Pauli-Str. 27 8093 Zurich Switzerland
dolfim@phys.ethz.ch www.itp.phys.ethz.ch
+41 44 633 78 56 phone +41 44 633 11 15 fax
On 17 Jul 2015, at 12:22, Wei.Su xichuang_s@163.com wrote:
I occured an error "Intel MKL ERROR: Parameter 10 was incorrect on entry to DGEMM " when using 'singlesite' optimization and 'storagedir' together.
I'm not sure if this is an error of MKL or one of the MPS code in ALPS
However it runs ok if i turn off the "storagedir" option or use the "twosite" optimization schedule.
Here comes the python script
#--------------------------------------------------------------------------------------------#
import pyalps import numpy as np import matplotlib.pyplot as plt import pyalps.plot
#prepare the input parameters parms = [ { 'optimization' : 'singlesite', 'LATTICE' : 'open chain lattice', 'L' : 20, 'MODEL' : 'spin', 'local_S0' : '0.5', 'local_S1' : '1', 'CONSERVED_QUANTUMNUMBERS' : 'N,Sz', 'Sz_total' : 0, 'J' : 1, 'SWEEPS' : 4, 'NUMBER_EIGENVALUES' : 1, 'MAXSTATES' : 50, 'storagedir' : 'storage' } ]
#write the input file and run the simulation input_file = pyalps.writeInputFiles('parm_spin_one',parms) res = pyalps.runApplication('mps_optim',input_file,writexml=True)
#----------------------------------------------------------------------------------------#
-------Su Wei-------
<spin_one.py>