SeismicWaves.jl is a Julia package for acoustic and elastic wave propagation simulations designed to be used in a full-waveform inversion (FWI) framework. It solves different flavours of the wave equation using the finite difference method.
The main features of SeismicWaves.jl are:
- acoustic and elastic simulations
- gradients using the adjoint method
- different misfit/objective functionals provided or the possibility to specifing your own
- device-agnostic backends for parallelization, i.e., CPUs or GPUs (CUDA, AMDGPU, Metal - Apple M-chips) thanks to
ParallelStencil.jl - C-PML boundary conditions
- checkpointing of simulations for the adjoint method
This package integrates well with the inversion framework provided by InverseAlgos, which includes both deterministic (L-BFGS, Quasi-Newton, etc.) and probabilistic (e.g., Hamiltonian Monte Carlo) inverse algorithms, as part of the G⁻¹Lab framework.
More information and an extensive list of features can be found in the documentation, which can be found either online or built locally by running the docs/make.jl file.
Warning
Documentation is currently minimal and work in progress!