ci: add stagg test for inner param

Author: mloubout

devito/finite_differences/differentiable.py

@@ -1134,7 +1134,7 @@ def _(expr, x0, **kwargs):
     from devito.finite_differences.derivative import Derivative
     x0_expr = {d: v for d, v in x0.items() if v is not expr.indices_ref[d]}
     if expr.is_parameter:
-        return expr._evaluate(**kwargs)
+        return expr
     elif x0_expr:
         dims = tuple((d, 0) for d in x0_expr)
         fd_o = tuple([expr.interp_order]*len(dims))

tests/test_derivatives.py

@@ -762,6 +762,24 @@ def test_nested_call(self):
         # Should be commutative
         assert simplify(deriv.evaluate - derivc.evaluate) == 0
 
+    def test_param_stagg_inner(self):
+        space_order = 2
+        nx, ny = 5, 5
+
+        grid = Grid((nx, ny))
+
+        x, y = grid.dimensions
+        yp = y + y.spacing / 2
+        xp = x + x.spacing / 2
+
+        f = TimeFunction(name="f", grid=grid, space_order=space_order, staggered=y)
+        p = Function(name="p", grid=grid, space_order=space_order, parameter=True)
+        g = TimeFunction(name="g", grid=grid, space_order=space_order, staggered=(x, y))
+
+        eqn = Eq(g, (p * f).dx)
+        eqne = eqn.evaluate.rhs
+        assert simplify(eqne - (p._subs(y, yp).evaluate * f).dx(x0=xp).evaluate) == 0
+
 
 class TestTwoStageEvaluation: