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.


Article metrics loading...

Loading full text...

Full text loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error