  We investigate the phase transition of the four-dimensional Ising model with
two types of tensor network scheme, one is the higher-order tensor
renormalization group and the other is the anisotropic tensor renormalization
group. The results for the internal energy and magnetization obtained by the
former algorithm with the impure tensor method, enlarging the lattice volume up
to $1024^4$, are consistent with the weak first-order phase transition. For the
later algorithm, our implementation successfully reduces the execution time
thanks to the parallel computation and the results provided by ATRG seems
comparable to those with HOTRG.
