code: def test_blockwise_stacked_new_axes(concatenate):
      def f(x):
          return x[..., None] * np.ones((1, 7))
  
      x = da.ones(5, chunks=2)
      y = da.blockwise(
          f,
          "aq",
          x,
          "a",
          new_axes={"q": 7},
          concatenate=concatenate,
          dtype=x.dtype,
      )
      z = da.blockwise(
          f,
          "abq",
          y,
          "ab",
          new_axes={"q": 7},
          concatenate=concatenate,
          dtype=x.dtype,
      )
      assert z.chunks == ((2, 2, 1), (7,), (7,))
      assert_eq(z, np.ones((5, 7, 7)))
  
imports: from __future__ import annotations
  import collections
  from operator import add
  import numpy
  import pytest
  import dask
  import dask.array
  from dask.array.utils import assert_eq
  from dask.blockwise import _BLOCKWISE_DEFAULT_PREFIX, Blockwise, _unique_dep, index_subs, optimize_blockwise, rewrite_blockwise
  from dask.highlevelgraph import HighLevelGraph
  from dask.utils_test import dec, hlg_layer_topological, inc
label_desc: Use image node_id 15 for calling a global function with example usage: test_blockwise_stacked_new_axes(concatenate) and returns: unknown