Check Subdimension.local is honored to avoid OOB accesses

[FabioLuporini]

reproducer:

python tti_example.py

look at the generated code for the padfunc builtin, you'll see something along the lines of

phi[abc_x_l + 4][y + 4][z + 4] = phi[44][y + 4][z + 4];

the 44 should be a symbolic entity, not an hardcoded number

[FabioLuporini]

update: this is due to SubDimension(local=...

[FabioLuporini]

however, the +40 due to nbl should become a symbol otherwise it's gonna be OOB if we supply a (much) smaller Function upon op.apply ?

[mloubout]

gonna be OOB if we supply a (much) smaller Function upon op.apply ?

THis line should prevent it no?

devito/devito/builtins.py

Line 306 in 08aa3f8

if mode == 'reflect' and function.grid.distributor.is_parallel:

[mloubout]

No would OOB for the constant case still.