Frequency domain wave equation forward modeling using gaussian elimination with static pivoting

Frequency domain wave equation forward modeling is a problem of solving large scale linear sparse systems which is often subject to the limits of computational efficiency and memory storage. Conventional Gaussian elimination cannot resolve the parallel computation of huge data. Therefore, we use the Gaussian elimination with static pivoting (GESP) method for sparse matrix decomposition and multi-source finite-difference modeling. The GESP method does not only improve the computational efficiency but also benefit the distributed parallel computation of matrix decomposition within a single frequency point. We test the proposed method using the classic Marmousi model. Both the single-frequency wave field and time domain seismic section show that the proposed method improves the simulation accuracy and computational efficiency and saves and makes full use of memory. This method can lay the basis for waveform inversion.