Seismic imaging based on spectral differentiation matrix and GPU implementation

► A new seismic imaging method is presented based on spectral matrix. ► Exponential accuracy and unconditional stability are guaranteed at the same time. ► GPU implementation speeds up the computation.