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.