“…When the wave number k is large, the solution of the problem becomes highly oscillating and efficient numerical methods are required in order to get high performance simulation results. In this topic, various numerical methods were developed in the past decades, such as the finite difference method (see, e.g., [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16]), the finite element method (see, e.g., [17][18][19][20][21][22][23][24][25]), the boundary element method (see, e.g., [26][27][28][29][30]), and other techniques (see, e.g., [17,[31][32][33]). For the finite difference method, two common methods are considered in the literature, namely the parameter method and the high-order method.…”