Use scipy.sparse.linalg.bicgstab() instead of own version. Replace trace(dot(a, b)) scalar products with sum(multiply(a, b)) for speed.
Use scipy.sparse.linalg.bicgstab() instead of own version. Replace trace(dot(a, b)) scalar products with sum(multiply(a, b)) for speed.

ftranschel/evoMPS