code: def find_geom(geom, geoms):
      
      for i, g in enumerate(geoms):
          if g is geom:
              return i
docs: Returns the index of a geometry in a list of geometries avoiding
      expensive equality checks of `in` operator.