A parallel computing method for the higher order tensor renormalization group

Abstract

In this paper, we propose a parallel computing method for the Higher Order Tensor Renormalization Group (HOTRG) applied to a d-dimensional (d≥2) simple lattice model. Sequential computation of the HOTRG requires O(χ4d−1) computational cost, where χ is bond dimension, in a step to contract indices of tensors. When we simply distribute elements of a local tensor to each process in parallel computing of the HOTRG, frequent communication between processes occurs. The simplest way to avoid such communication is to hold all the tensor elements in each process, however, it requires O(χ2d) memory space. In the presented method, placement of a local tensor element to more than one process is accepted and sufficient local tensor elements are distributed to each process to avoid communication between processes during considering computation step. For the bottleneck part of computational cost, such distribution is achieved by distributing elements of two local tensors to χ2 processes according to one of the indices of each local tensor which are not contracted during considering computation. In the case of d≥3, computational cost in each process is reduced to O(χ4d−3) and memory space requirement in each process is kept to be O(χ2d−1).