code: def test_set_operation_reduce_axis(func, related_func):
      data = [[point] * 2] * 3  # shape = (3, 2)
      actual = func(data, axis=None)  # default
      assert isinstance(actual, Geometry)  # scalar output
      actual = func(data, axis=0)
      assert actual.shape == (2,)
      actual = func(data, axis=1)
      assert actual.shape == (3,)
      actual = func(data, axis=-1)
      assert actual.shape == (3,)
  
imports: import numpy
  import pytest
  import shapely
  from shapely import Geometry, GeometryCollection, Polygon
  from shapely.errors import UnsupportedGEOSVersionError
  from shapely.testing import assert_geometries_equal
  from shapely.tests.common import all_types, empty, ignore_invalid, point, polygon
label_desc: Use image node_id 7 for calling a global function with example usage: test_set_operation_reduce_axis(func, related_func) without return types