From d8913e51559169a4954e48217d8a92bc6adf6539 Mon Sep 17 00:00:00 2001 From: Keith Roberts Date: Wed, 14 Apr 2021 15:45:09 -0300 Subject: [PATCH] better threshold --- spyro/solvers/leapfrog_adjoint_level_set.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/spyro/solvers/leapfrog_adjoint_level_set.py b/spyro/solvers/leapfrog_adjoint_level_set.py index bab3353c..47afdf42 100644 --- a/spyro/solvers/leapfrog_adjoint_level_set.py +++ b/spyro/solvers/leapfrog_adjoint_level_set.py @@ -186,7 +186,7 @@ def Leapfrog_adjoint_level_set( # analytical background gradient # Z, _ = SpatialCoordinate(mesh) # c_background = Function(V).interpolate(Min(1.5 + 4.0 * abs(Z), 4.1)) - indicator = Function(V).interpolate(c > 3.9) # 1 inside shape, 0 outside + indicator = Function(V).interpolate(c > 4.3) # 1 inside shape, 0 outside gradc = Function(VF, name="grad_c").interpolate( indicator * grad(c_salt) + (1 - indicator) * grad(c_background)