Repository files navigation

GeoPhyInv provides architecture-agnostic elastic and acoustic wave equation solvers using eitherCUDA.jl forBase.Threads for high-performance computations on GPUs and CPUs, respectively. For GPU computations, a performance similar to CUDA C is achieved, thanks toParallelStencil.jl.The finite-difference simulations are performed in both 2-D and 3-Dusing a staggered-grid velocity-stress formulation.Finally,distributed computing shipped with Julia ensures that the modelling of the super-sources can be parallelized.

For complete installation, enter these package manager commands in the REPL:

using PkgPkg.add(PackageSpec(name="GeoPhyInv",url="https://github.com/pawbz/GeoPhyInv.jl.git"))

• Some implementation ideas are borrowed from Jan Thorbecke'sfdelmodc software.
• [paper] P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method.
• Work ofKomatitsch and Martin (2007) on convolutional perfectly matched layers for seismic wave equation.
• Charles Clerget@chclerget tested some methods of this package.
• The Poisson solver (GeoPhyInv.Poisson) was developed by Niels Grobbe, after adapting scripts fromAime Fournier.
• Thanks toEarth Resources Laboratory, MIT. A few developments of this project were supported by Aime Fournier via research funds from Equinor.