Our paper was accepted by Commun. Comput. Phys. This is a joint work with Prof. Zhennan Zhou (Westlake University) and Shengbo Ma (Zhejiang University).