Parallel implementation of the time-evolving block decimation algorithm for the Bose-Hubbard model

Running time: Running time depends on the lattice size, on the number of particles, and on the maximal allowed tensor dimensions. Simulations of simple models take a few seconds on a single CPU while simulations of large and complex models can take up to a week on a cluster with hundreds of CPUs.