We develop a state-of-the-art numerical algorithm for solving the relativistic viscous hydrodynamic
equation, using a Riemann solver based on two-shock approximation.
We show validity of our hydrodynamics code from comparison with analytical solutions of
test problems.
Furthermore we perform numerical hydrodynamic calculation with a realistic fluctuating initial condition of high-energy
heavy-ion collisions at LHC produced by TRENTO.
We show time evolution of temperature distributions on $x-y$ and $x-\eta$ planes with lattice QCD equation of state (EoS) and
discuss the EoS dependence of it .