Navjot Kukreja

Unexpected aliasing with multiple Grids

9

MFE: ``` from devito import Grid g1 = Grid(shape=(10, 10)) g2 = Grid(shape=(20, 20)) x1, y1 = g1.dimensions x2, y2 = g2.dimensions x1 is x2 # True ``` What happens:...

Absolute value of an integer produces C warning

4

Here is an MFE: ``` import numpy as np from sympy import sign from devito import Function, Grid, Eq, Operator gridR = Grid(shape=(2,)) index = Function(name="index", grid=gridR, dtype=np.int32) R =...

Memory consumption gradually grows when Operator called in loop

11

In current master, there are a few sources of increased memory consumption when operators are created in a loop through the seismic wrapper classes: - memoized methods in Solver -...

A utility for interpolation and resampling

10

Over time the work on resampling input data in time and the spatial interpolation of off-grid points seem to have a lot of overlap and this could be cleanly set...

[WIP] Add gradient subsampling test

1

This PR is the next milestone leading up to the compression followup work. It adds a test that computes the gradient normally first, then does it again for `dt/2` and...

Gradient test in master fails for space order 16

2

Modifying [this](https://github.com/devitocodes/devito/blob/master/tests/test_gradient.py#L58) test to run at space order 16 causes it to fail. Is this expected?

- If `LinearInterpolator._interpolate`/`_inject` could accept `sfunction` as parameters, no state would be required in the `LinearInterpolator` class and it could be singleton-ized (see https://github.com/devitocodes/devito/pull/1025#discussion_r358091918 ) - The initialization of `Interpolator`...