Fix #4 avoid OOB in born due to double offset

Project.toml

@@ -1,7 +1,7 @@
 name = "ImageGather"
 uuid = "355d8124-6b2e-49b5-aab5-cdbc0a5fccbe"
 authors = ["mloubout <[email protected]>"]
-version = "0.2.5"
+version = "0.2.6"
 
 [deps]
 JUDI = "f3b833dc-6b2e-5b9c-b940-873ed6319979"

examples/layers_sscig.jl

@@ -49,7 +49,7 @@ wavelet = ricker_wavelet(timeD, dtD, f0)
 q = diff(judiVector(srcGeometry, wavelet))
 
 ###################################################################################################
-opt = Options(space_order=12, isic=false, sum_padding=true)
+opt = Options(space_order=4, isic=false, sum_padding=true)
 # Setup operators
 F = judiModeling(model, srcGeometry, recGeometry; options=opt)
 J0 = judiJacobian(F(model0), q)

src/implementation.py

@@ -96,8 +96,10 @@ def cig_lin(model, src_coords, wavelet, rec_coords, dm_ext, offsets, ic="as", sp
     op.cfunction
 
     # Remove edge for offsets
-    dim_kw = make_kw(oh, DimensionTuple(*u.shape[1:], getters=model.grid.dimensions))
+    dim_kw = make_kw(oh, DimensionTuple(*u.shape[1:], getters=model.grid.dimensions),
+                     born=True)
     op(dt=model.critical_dt, rcvul=rcvl, **dim_kw)
+
     # Output
     return rcvl.data
 
@@ -149,9 +151,10 @@ def _shift(u, oh):
     return uh
 
 
-def make_kw(ohs, shape):
+def make_kw(ohs, shape, born=False):
     kw = dict()
+    scale = 2 if born else 1
     for d, v in ohs.items():
-        kw['%s_m' % d.name] = -np.min(v.data.view(np.ndarray))
-        kw['%s_M' % d.name] = shape[d] - np.max(v.data.view(np.ndarray))
+        kw['%s_m' % d.name] = -scale*np.min(v.data.view(np.ndarray))
+        kw['%s_M' % d.name] = shape[d] - scale*np.max(v.data.view(np.ndarray))
     return kw