A MATLAB Package for Frequency Domain Modeling of Elastic Waves

We developed a MATLAB package for finite difference frequency domain modeling of elastic waveforms. By using array-processing abilities of MATLAB, we efficiently computed the large impedance matrix for realistic model sizes. SuiteSparse library was used to calculate LU factors of complex impedance matrix with reasonable computational costs and METIS library was used to benefit from nested dissection ordering scheme to accelerate LU factorization. Once the factors have been calculated, several seismic sources could be modeled by solving for multiple right hand sides. Perfectly Matched Layers (PML) suppressed the reflections from the boundaries of computational domain. Attenuation effects have been considered by including complex valued velocities. Synthetic examples of Marmousi2 dataset show the performance of the developed package.