A divide-and-conquer algorithm for seismic data approximation by the Laguerre series
Abstract An algorithm of the Laguerre transform for approximating functions on large intervals is proposed. The idea of the considered approach is that the calculation of improper integrals of rapidly oscillating functions is replaced by a solution of an initial boundary value problem for the one-dimensional transport equation. It allows one to successfully avoid the problems associated with the stable implementation of the Laguerre transform. A divide-and-conquer algorithm based on shift operations made it possible to significantly reduce the computational cost of the proposed method. Numerical experiments have shown that the methods are economical in the number of operations, stable, and have satisfactory accuracy for seismic data approximation.