code: def test_theano_function_numpy():
      """Test theano_function() vs Numpy implementation."""
      f = theano_function_(
          [x, y], [x + y], dim=1, dtypes={x: "float64", y: "float64"}
      )
      assert (
          np.linalg.norm(f([1, 2], [3, 4]) - np.asarray([4, 6])) < 1e-9
      )
  
      f = theano_function_(
          [x, y], [x + y], dtypes={x: "float64", y: "float64"}, dim=1
      )
      xx = np.arange(3).astype("float64")
      yy = 2 * np.arange(3).astype("float64")
      assert np.linalg.norm(f(xx, yy) - 3 * np.arange(3)) < 1e-9
  
imports: import logging
  from sympy.external import import_module
  from sympy.testing.pytest import raises, SKIP, warns_deprecated_sympy
  import sympy
  from sympy.core.singleton import S
  from sympy.abc import x, y, z, t
  from sympy.printing.theanocode import theano_code, dim_handling, theano_function
label_desc: Use image node_id 26 for calling a global function with example usage: test_theano_function_numpy() without return types