Tensor renormalization group approach to ( 1+1 )-dimensional Hubbard model