  We perform hybrid Monte-Carlo simulation of $N_f=2+1+1 $ lattice QCD with
domain-wall quarks at the physical point. The simulation is carried out on the
$ L^3 \times T = 64^3 \times 64 $ lattice with lattice spacing $a \sim 0.064 $
fm ($ L > 4 $ fm, and $ M_\pi L > 3 $), using the Nvidia DGX-1 (8 V100 GPUs
interconnected by the NVLink). To attain the maximal chiral symmetry for a
finite extent ($N_s=16$) in the fifth dimension, we use the optimal domain-wall
fermion for the quark action, together with the exact one-flavor action for
domain-wall fermion. We outline the salient features of our simulation (without
topology freezing, and small residual mass), together with the preliminary
result of the masses of $ \pi^\pm $, $ K^\pm $, and $ D^\pm $.
