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-------