Elastic 2D P-SV adjoints and gradients

Project.toml

@@ -5,6 +5,7 @@ version = "0.7.0"
 
 [deps]
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
+Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 Logging = "56ddb016-857b-54e1-b83d-db4d58db5568"
 ParallelStencil = "94395366-693c-11ea-3b26-d9b7aac5d958"

examples/gradient_check_fd/constant_density/gradient_check_fd.jl → examples/gradient_check_fd/acoustic/constant_density/gradient_check_fd.jl

@@ -14,9 +14,9 @@ warn_logger = ConsoleLogger(stderr, Logging.Warn)
 error_logger = ConsoleLogger(stderr, Logging.Error)
 debug_logger = ConsoleLogger(stderr, Logging.Debug)
 
-include("../../models.jl")
-include("../../geometries.jl")
-include("../../plotting_utils.jl")
+include("../../../models.jl")
+include("../../../geometries.jl")
+include("../../../plotting_utils.jl")
 
 function setup(nt, c0, c0max, r, dx, dy, dt, halo, rcoef, nx, ny, parall)
     ##========================================
@@ -81,7 +81,7 @@ function gradient_fd_check(wavesim, shots, matprop_const, matprop_gauss; compute
             copy(shots[i].recs.positions),
             nt;
             observed=copy(shots[i].recs.seismograms),
-            invcov=Diagonal(ones(nt))
+            invcov=1.0 * I(nt)
         )
         # add pair as shot
         push!(shots_obs, ScalarShot(; srcs=shots[i].srcs, recs=recs)) # srcs => recs)

examples/gradient_check_fd/variable_density/gradient_check_fd.jl → examples/gradient_check_fd/acoustic/variable_density/gradient_check_fd.jl

@@ -14,9 +14,9 @@ warn_logger = ConsoleLogger(stderr, Logging.Warn)
 error_logger = ConsoleLogger(stderr, Logging.Error)
 debug_logger = ConsoleLogger(stderr, Logging.Debug)
 
-include("../../models.jl")
-include("../../geometries.jl")
-include("../../plotting_utils.jl")
+include("../../../models.jl")
+include("../../../geometries.jl")
+include("../../../plotting_utils.jl")
 
 function setup(nt, c0, c0max, rho0, rho0max, r, dx, dy, dt, halo, rcoef, nx, ny, parall)
     ##========================================
@@ -96,7 +96,7 @@ function gradient_fd_check(wavesim, shots, matprop_const, matprop_gauss)
             copy(shots[i].recs.positions),
             nt;
             observed=copy(shots[i].recs.seismograms),
-            invcov=Diagonal(ones(nt))
+            invcov=1.0 * I(nt)
         )
         # add pair as shot
         push!(shots_obs, ScalarShot(; srcs=shots[i].srcs, recs=recs)) # srcs => recs)

examples/gradient_check_fd/variable_density/plot_gradient.jl → examples/gradient_check_fd/acoustic/variable_density/plot_gradient.jl

@@ -1,8 +1,8 @@
 using Serialization
 
-include("../../models.jl")
-include("../../geometries.jl")
-include("../../plotting_utils.jl")
+include("../../../models.jl")
+include("../../../geometries.jl")
+include("../../../plotting_utils.jl")
 
 # Numerical parameters
 nt = 1000