Fabio Luporini

Reproducer: ``` from devito import * u = Function(name="u", grid=Grid((5, 5))) x, y = u.grid.dimensions ox, oy = u.grid.origin # Operator(Eq(u, exp(-(x-2)**2/.5 - (y-2)**2/.5)))() Operator(Eq(u, x+y))() print(u.data._gather(rank=0)) ```

Producing garbage numbers. At least, after #2136 , we'll get a suitable exception in Python during Operator compilation

Complete discussion here: https://devitocodes.slack.com/archives/C7JMLMSG0/p1673960926948329 No MFE available right now The issue discussion started here: https://github.com/devitocodes/devito/pull/2042/files#r1064672824 The latest suggestion from @mloubout was: ``` `p1-1 - (p1-1 + p0*factor)%factor` this should always...

Introduce devito.Zero to replace sympy.Zero in the FD layer

1

this way we avoid annoying try-excepts such as ``` try: theta = model.theta except AttributeError: theta = 0 ``` or ``` if theta != 0: ```

Dockerfile for ROCm + Omniperf (and more)

8

Hi again, Sorry, first of all, if this is the wrong place to post this. I genuinely wonder whether AMDResearch would be willing to maintain a Dockerfile that ships the...

compiler: Yet another batch of compilation tweaks

1

These are mostly to improve code generation and performance of elastic-like PDEs.

The logic is currently spread over multiple layers (stree, iet) It's also quite chaotic as it currently stands, so we have to rewrite it because in its current form is...