f

Author: rafaqz

src/geo_interface.jl

@@ -54,7 +54,7 @@ function GeoInterface.getgeom(::AbstractGeometryTrait, geom::Union{LineString,Li
     return _get_tuple_point(geom, seq, refs, i)
 end
 function GeoInterface.getgeom(::AbstractGeometryTrait, geom::Union{LineString,LinearRing})
-    context = get_global_context()
+    context = get_context(geom)
     seq = getCoordSeq(geom, context)
     n = getSize(seq, context) # Faster thatn `GI.ngeom(geom)` when we already have `seq`
     # Preallocate refse
@@ -65,7 +65,7 @@ function GeoInterface.getgeom(::AbstractGeometryTrait, geom::Union{LineString,Li
     return (_get_tuple_point(geom, seq, refs, i, context, hasz) for i in UInt32(1):UInt32(n))
 end
 
-function _get_tuple_point(geom, seq, (x, y, z), i, context = get_global_context(), hasz=hasZ(geom, context))
+function _get_tuple_point(geom, seq, (x, y, z), i, context=get_context(geom), hasz=hasZ(geom, context))
     if hasz
         GEOSCoordSeq_getXYZ_r(context, seq, i - UInt32(1), x, y, z)
         return x[], y[], Z[]
@@ -74,7 +74,6 @@ function _get_tuple_point(geom, seq, (x, y, z), i, context = get_global_context(
         return x[], y[]
     end
 end
-
 GeoInterface.getgeom(t::AbstractPointTrait, geom::PreparedGeometry) = nothing
 function GeoInterface.getgeom(::PolygonTrait, geom::Polygon, i::Int)
     if i == 1

src/geos_functions.jl

@@ -97,7 +97,7 @@ function cloneCoordSeq(ptr::GEOSCoordSeq, context::GEOSContext = get_global_cont
 end
 
 function destroyCoordSeq(ptr::GEOSCoordSeq, context::GEOSContext = get_global_context())
-    o
+    result = GEOSCoordSeq_destroy_r(context, ptr)
     if result == C_NULL
         error("LibGEOS: Error in GEOSCoordSeq_destroy")
     end
@@ -1306,17 +1306,16 @@ Numbering is one-based.
 For a simple geometry, returns a copy of the input.
 """
 function getGeometry(obj::Geometry, n::Integer, context::GEOSContext = get_context(obj))
-    n in 1:numGeometries(obj, context) || _no_n_geom(obj, context, n)
+    n in 1:numGeometries(obj, context) || error(
+        "GEOSGetGeometryN: $(numGeometries(obj, context)) sub-geometries in geom, therefore n should be in 1:$(numGeometries(obj, context))",
+    )
     result = GEOSGetGeometryN_r(context, obj, n - 1)
     if result == C_NULL
         error("LibGEOS: Error in GEOSGetGeometryN")
     end
     geomFromGEOS(cloneGeom(result, context), context)
 end
 
-@noinline _no_n_geom_error(obj, context, n) =
-    erroe("GEOSGetGeometryN: $(numGeometries(obj, context)) sub-geometries in geom, therefore n: $n should be in 1:$(numGeometries(obj, context))",)
-
 """
     getGeometries(obj::Geometry, context::GEOSContext = get_context(obj))