code: def test_get_parts_return_index():
      geom = np.array([multi_point, point, multi_polygon])
      expected_parts = []
      expected_index = []
      for i, g in enumerate(geom):
          for j in range(0, shapely.get_num_geometries(g)):
              expected_parts.append(shapely.get_geometry(g, j))
              expected_index.append(i)
  
      parts, index = shapely.get_parts(geom, return_index=True)
      assert len(parts) == len(expected_parts)
      assert_geometries_equal(parts, expected_parts)
      assert np.array_equal(index, expected_index)
  
imports: import warnings
  import numpy
  import pytest
  import shapely
  from shapely import LinearRing, LineString, MultiPolygon, Point, Polygon
  from shapely.testing import assert_geometries_equal
  from shapely.tests.common import all_types
  from shapely.tests.common import empty
  from shapely.tests.common import empty_line_string, empty_line_string_z, empty_point, empty_point_z, empty_polygon, equal_geometries_abnormally_yield_unequal, geometry_collection, geometry_collection_z, ignore_invalid, ignore_warnings, line_string, line_string_nan, line_string_z, linear_ring, multi_line_string, multi_line_string_z, multi_point, multi_point_z, multi_polygon, multi_polygon_z, point, point_z, polygon, polygon_with_hole, polygon_with_hole_z, polygon_z
label_desc: Use image node_id 31 for calling a global function with example usage: test_get_parts_return_index() without return types