Make Reshape work for Numba scalars

aesara/link/numba/dispatch.py

@@ -1154,7 +1154,10 @@ def numba_funcify_Reshape(op, **kwargs):
 
         @numba.njit(inline="always")
         def reshape(x, shape):
-            return np.reshape(np.ascontiguousarray(x), to_fixed_tuple(shape, ndim))
+            # TODO: Use this until https://github.com/numba/numba/issues/7353 is closed.
+            return np.reshape(
+                np.ascontiguousarray(np.asarray(x)), to_fixed_tuple(shape, ndim)
+            )
 
     return reshape
 

tests/link/test_numba.py

@@ -817,6 +817,21 @@ def test_Reshape(v, shape, ndim):
     )
 
 
+def test_Reshape_scalar():
+    v = aet.vector()
+    v.tag.test_value = np.array([1.0], dtype=config.floatX)
+    g = Reshape(1)(v[0], (1,))
+    g_fg = FunctionGraph(outputs=[g])
+    compare_numba_and_py(
+        g_fg,
+        [
+            i.tag.test_value
+            for i in g_fg.inputs
+            if not isinstance(i, (SharedVariable, Constant))
+        ],
+    )
+
+
 @pytest.mark.parametrize(
     "v, shape, fails",
     [