rgeo 3.0.0.pre.rc.1 → 3.0.0.pre.rc.3

data/ext/geos_c_impl/factory.c

@@ -2,144 +2,91 @@
   Factory and utility functions for GEOS wrapper
 */
 
-
 #include "preface.h"
 
 #ifdef RGEO_GEOS_SUPPORTED
 
-#include <ruby.h>
 #include <geos_c.h>
-#include <ctype.h>
-#include <stdarg.h>
-#include <stdio.h>
-
-#include "globals.h"
+#include <ruby.h>
 
+#include "errors.h"
 #include "factory.h"
 #include "geometry.h"
-#include "point.h"
+#include "geometry_collection.h"
+#include "globals.h"
 #include "line_string.h"
+#include "point.h"
 #include "polygon.h"
-#include "geometry_collection.h"
-#include "errors.h"
+#include "ruby_more.h"
 
 RGEO_BEGIN_C
 
-
 /**** RUBY AND GEOS CALLBACKS ****/
 
-
-// The notice handler is very rarely used by GEOS, only in
-// GEOSIsValid_r (check for NOTICE_MESSAGE in GEOS codebase).
-// We still set it to make sure we do not miss any implementation
-// change. Use `DEBUG=1 rake` to show notice.
-#ifdef RGEO_GEOS_DEBUG
-static void notice_handler(const char* fmt, ...)
-{
-  va_list args;
-  va_start(args, fmt);
-  fprintf(stderr, "GEOS Notice -- ");
-  vfprintf(stderr, fmt, args);
-  fprintf(stderr, "\n");
-  va_end(args);
-}
-#endif
-
-static void error_handler(const char* fmt, ...)
-{
-  // See https://en.cppreference.com/w/c/io/vfprintf
-  va_list args1;
-  va_start(args1, fmt);
-  va_list args2;
-  va_copy(args2, args1);
-  int size = 1+vsnprintf(NULL, 0, fmt, args1);
-  va_end(args1);
-  char geos_full_error[size];
-  vsnprintf(geos_full_error, sizeof geos_full_error, fmt, args2);
-  va_end(args2);
-
-  // NOTE: strtok is destructive, geos_full_error is not to be used afterwards.
-  char *geos_error = strtok(geos_full_error, ":");
-  char *geos_message = strtok(NULL, ":");
-  while(isspace(*geos_message)) geos_message++;
-
-  if (streq(geos_error, "UnsupportedOperationException")) {
-    rb_raise(rb_eRGeoUnsupportedOperation, "%s", geos_message);
-  } else if (streq(geos_error, "IllegalArgumentException")) {
-    rb_raise(rb_eRGeoInvalidGeometry, "%s", geos_message);
-  } else if (geos_message) {
-    rb_raise(rb_eGeosError, "%s: %s", geos_error, geos_message);
-  } else {
-    rb_raise(rb_eGeosError, "%s", geos_error);
-  }
-}
-
 // Destroy function for factory data. We destroy any serialization
-// objects that have been created for the factory, and then destroy
-// the GEOS context, before freeing the factory data itself.
+// objects that have been created for the factory before freeing
+// the factory data itself.
 
-static void destroy_factory_func(void* data)
+static void
+destroy_factory_func(void* data)
 {
   RGeo_FactoryData* factory_data;
-  GEOSContextHandle_t context;
 
   factory_data = (RGeo_FactoryData*)data;
-  context = factory_data->geos_context;
   if (factory_data->wkt_reader) {
-    GEOSWKTReader_destroy_r(context, factory_data->wkt_reader);
+    GEOSWKTReader_destroy(factory_data->wkt_reader);
   }
   if (factory_data->wkb_reader) {
-    GEOSWKBReader_destroy_r(context, factory_data->wkb_reader);
+    GEOSWKBReader_destroy(factory_data->wkb_reader);
   }
   if (factory_data->wkt_writer) {
-    GEOSWKTWriter_destroy_r(context, factory_data->wkt_writer);
+    GEOSWKTWriter_destroy(factory_data->wkt_writer);
   }
   if (factory_data->wkb_writer) {
-    GEOSWKBWriter_destroy_r(context, factory_data->wkb_writer);
+    GEOSWKBWriter_destroy(factory_data->wkb_writer);
   }
   if (factory_data->psych_wkt_reader) {
-    GEOSWKTReader_destroy_r(context, factory_data->psych_wkt_reader);
+    GEOSWKTReader_destroy(factory_data->psych_wkt_reader);
   }
   if (factory_data->marshal_wkb_reader) {
-    GEOSWKBReader_destroy_r(context, factory_data->marshal_wkb_reader);
+    GEOSWKBReader_destroy(factory_data->marshal_wkb_reader);
   }
   if (factory_data->psych_wkt_writer) {
-    GEOSWKTWriter_destroy_r(context, factory_data->psych_wkt_writer);
+    GEOSWKTWriter_destroy(factory_data->psych_wkt_writer);
   }
   if (factory_data->marshal_wkb_writer) {
-    GEOSWKBWriter_destroy_r(context, factory_data->marshal_wkb_writer);
+    GEOSWKBWriter_destroy(factory_data->marshal_wkb_writer);
   }
-  finishGEOS_r(context);
-  free(factory_data);
+  FREE(factory_data);
 }
 
-
 // Destroy function for geometry data. We destroy the internal
 // GEOS geometry (if present) before freeing the data itself.
 
-static void destroy_geometry_func(void* data)
+static void
+destroy_geometry_func(void* data)
 {
   RGeo_GeometryData* geometry_data;
   const GEOSPreparedGeometry* prep;
 
   geometry_data = (RGeo_GeometryData*)data;
   if (geometry_data->geom) {
-    GEOSGeom_destroy_r(geometry_data->geos_context, geometry_data->geom);
+    GEOSGeom_destroy(geometry_data->geom);
   }
   prep = geometry_data->prep;
-  if (prep && prep != (const GEOSPreparedGeometry*)1 && prep != (const GEOSPreparedGeometry*)2 &&
-    prep != (const GEOSPreparedGeometry*)3)
-  {
-    GEOSPreparedGeom_destroy_r(geometry_data->geos_context, prep);
+  if (prep && prep != (const GEOSPreparedGeometry*)1 &&
+      prep != (const GEOSPreparedGeometry*)2 &&
+      prep != (const GEOSPreparedGeometry*)3) {
+    GEOSPreparedGeom_destroy(prep);
   }
-  free(geometry_data);
+  FREE(geometry_data);
 }
 
-
 // Mark function for factory data. This marks the wkt and wkb generator
 // handles so they don't get collected.
 
-static void mark_factory_func(void* data)
+static void
+mark_factory_func(void* data)
 {
   RGeo_FactoryData* factory_data;
 
@@ -156,19 +103,16 @@ static void mark_factory_func(void* data)
   if (!NIL_P(factory_data->wkrep_wkb_parser)) {
     mark(factory_data->wkrep_wkb_parser);
   }
-  if (!NIL_P(factory_data->proj4_obj)) {
-    mark(factory_data->proj4_obj);
-  }
   if (!NIL_P(factory_data->coord_sys_obj)) {
     mark(factory_data->coord_sys_obj);
   }
 }
 
-
 // Mark function for geometry data. This marks the factory and klasses
 // held by the geometry so those don't get collected.
 
-static void mark_geometry_func(void* data)
+static void
+mark_geometry_func(void* data)
 {
   RGeo_GeometryData* geometry_data;
 
@@ -181,35 +125,36 @@ static void mark_geometry_func(void* data)
   }
 }
 
-
 #ifdef HAVE_RB_GC_MARK_MOVABLE
-static void compact_factory_func(void* data)
+static void
+compact_factory_func(void* data)
 {
   RGeo_FactoryData* factory_data;
 
   factory_data = (RGeo_FactoryData*)data;
   if (!NIL_P(factory_data->wkrep_wkt_generator)) {
-    factory_data->wkrep_wkt_generator = rb_gc_location(factory_data->wkrep_wkt_generator);
+    factory_data->wkrep_wkt_generator =
+      rb_gc_location(factory_data->wkrep_wkt_generator);
   }
   if (!NIL_P(factory_data->wkrep_wkb_generator)) {
-    factory_data->wkrep_wkb_generator = rb_gc_location(factory_data->wkrep_wkb_generator);
+    factory_data->wkrep_wkb_generator =
+      rb_gc_location(factory_data->wkrep_wkb_generator);
   }
   if (!NIL_P(factory_data->wkrep_wkt_parser)) {
-    factory_data->wkrep_wkt_parser = rb_gc_location(factory_data->wkrep_wkt_parser);
+    factory_data->wkrep_wkt_parser =
+      rb_gc_location(factory_data->wkrep_wkt_parser);
   }
   if (!NIL_P(factory_data->wkrep_wkb_parser)) {
-    factory_data->wkrep_wkb_parser = rb_gc_location(factory_data->wkrep_wkb_parser);
-  }
-  if (!NIL_P(factory_data->proj4_obj)) {
-    factory_data->proj4_obj = rb_gc_location(factory_data->proj4_obj);
+    factory_data->wkrep_wkb_parser =
+      rb_gc_location(factory_data->wkrep_wkb_parser);
   }
   if (!NIL_P(factory_data->coord_sys_obj)) {
     factory_data->coord_sys_obj = rb_gc_location(factory_data->coord_sys_obj);
   }
 }
 
-
-static void compact_geometry_func(void* data)
+static void
+compact_geometry_func(void* data)
 {
   RGeo_GeometryData* geometry_data;
 
@@ -223,89 +168,95 @@ static void compact_geometry_func(void* data)
 }
 #endif
 
-
-const rb_data_type_t rgeo_factory_type = {
-  .wrap_struct_name = "RGeo/Factory",
-  .function = {
-    .dmark = mark_factory_func,
-    .dfree = destroy_factory_func,
+const rb_data_type_t rgeo_factory_type = { .wrap_struct_name = "RGeo/Factory",
+                                           .function = {
+                                             .dmark = mark_factory_func,
+                                             .dfree = destroy_factory_func,
 #ifdef HAVE_RB_GC_MARK_MOVABLE
-    .dcompact = compact_factory_func,
+                                             .dcompact = compact_factory_func,
 #endif
-  }
-};
+                                           } };
 
-const rb_data_type_t rgeo_geometry_type = {
-  .wrap_struct_name = "RGeo/Geometry",
-  .function = {
-    .dmark = mark_geometry_func,
-    .dfree = destroy_geometry_func,
+const rb_data_type_t rgeo_geometry_type = { .wrap_struct_name = "RGeo/Geometry",
+                                            .function = {
+                                              .dmark = mark_geometry_func,
+                                              .dfree = destroy_geometry_func,
 #ifdef HAVE_RB_GC_MARK_MOVABLE
-    .dcompact = compact_geometry_func,
+                                              .dcompact = compact_geometry_func,
 #endif
-  }
-};
-
+                                            } };
 
 /**** RUBY METHOD DEFINITIONS ****/
 
-
-static VALUE method_factory_srid(VALUE self)
+static VALUE
+method_factory_srid(VALUE self)
 {
   return INT2NUM(RGEO_FACTORY_DATA_PTR(self)->srid);
 }
 
-
-static VALUE method_factory_buffer_resolution(VALUE self)
+static VALUE
+method_factory_buffer_resolution(VALUE self)
 {
   return INT2NUM(RGEO_FACTORY_DATA_PTR(self)->buffer_resolution);
 }
 
-
-static VALUE method_factory_flags(VALUE self)
+static VALUE
+method_factory_flags(VALUE self)
 {
   return INT2NUM(RGEO_FACTORY_DATA_PTR(self)->flags);
 }
 
-VALUE method_factory_supports_z_p(VALUE self)
+VALUE
+method_factory_supports_z_p(VALUE self)
 {
-  return RGEO_FACTORY_DATA_PTR(self)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z ? Qtrue : Qfalse;
+  return RGEO_FACTORY_DATA_PTR(self)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z
+           ? Qtrue
+           : Qfalse;
 }
 
-VALUE method_factory_supports_m_p(VALUE self)
+VALUE
+method_factory_supports_m_p(VALUE self)
 {
-  return RGEO_FACTORY_DATA_PTR(self)->flags & RGEO_FACTORYFLAGS_SUPPORTS_M ? Qtrue : Qfalse;
+  return RGEO_FACTORY_DATA_PTR(self)->flags & RGEO_FACTORYFLAGS_SUPPORTS_M
+           ? Qtrue
+           : Qfalse;
 }
 
-VALUE method_factory_supports_z_or_m_p(VALUE self)
+VALUE
+method_factory_supports_z_or_m_p(VALUE self)
 {
-  return RGEO_FACTORY_DATA_PTR(self)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M ? Qtrue : Qfalse;
+  return RGEO_FACTORY_DATA_PTR(self)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M
+           ? Qtrue
+           : Qfalse;
 }
 
-VALUE method_factory_prepare_heuristic_p(VALUE self)
+VALUE
+method_factory_prepare_heuristic_p(VALUE self)
 {
-  return RGEO_FACTORY_DATA_PTR(self)->flags & RGEO_FACTORYFLAGS_PREPARE_HEURISTIC ? Qtrue : Qfalse;
+  return RGEO_FACTORY_DATA_PTR(self)->flags &
+             RGEO_FACTORYFLAGS_PREPARE_HEURISTIC
+           ? Qtrue
+           : Qfalse;
 }
 
-static VALUE method_factory_parse_wkt(VALUE self, VALUE str)
+static VALUE
+method_factory_parse_wkt(VALUE self, VALUE str)
 {
   RGeo_FactoryData* self_data;
-  GEOSContextHandle_t self_context;
   GEOSWKTReader* wkt_reader;
   VALUE result;
   GEOSGeometry* geom;
 
   Check_Type(str, T_STRING);
   self_data = RGEO_FACTORY_DATA_PTR(self);
-  self_context = self_data->geos_context;
   wkt_reader = self_data->wkt_reader;
   if (!wkt_reader) {
-    wkt_reader = GEOSWKTReader_create_r(self_context);
+    wkt_reader = GEOSWKTReader_create();
     self_data->wkt_reader = wkt_reader;
   }
   result = Qnil;
   if (wkt_reader) {
-    geom = GEOSWKTReader_read_r(self_context, wkt_reader, RSTRING_PTR(str));
+    geom = GEOSWKTReader_read(wkt_reader, RSTRING_PTR(str));
     if (geom) {
       result = rgeo_wrap_geos_geometry(self, geom, Qnil);
     }
@@ -313,26 +264,31 @@ static VALUE method_factory_parse_wkt(VALUE self, VALUE str)
   return result;
 }
 
-
-static VALUE method_factory_parse_wkb(VALUE self, VALUE str)
+static VALUE
+method_factory_parse_wkb(VALUE self, VALUE str)
 {
   RGeo_FactoryData* self_data;
-  GEOSContextHandle_t self_context;
   GEOSWKBReader* wkb_reader;
   VALUE result;
   GEOSGeometry* geom;
+  char* c_str;
 
   Check_Type(str, T_STRING);
   self_data = RGEO_FACTORY_DATA_PTR(self);
-  self_context = self_data->geos_context;
   wkb_reader = self_data->wkb_reader;
   if (!wkb_reader) {
-    wkb_reader = GEOSWKBReader_create_r(self_context);
+    wkb_reader = GEOSWKBReader_create();
     self_data->wkb_reader = wkb_reader;
   }
   result = Qnil;
   if (wkb_reader) {
-    geom = GEOSWKBReader_read_r(self_context, wkb_reader, (unsigned char*)RSTRING_PTR(str), (size_t)RSTRING_LEN(str));
+    c_str = RSTRING_PTR(str);
+    if (c_str[0] == '\x00' || c_str[0] == '\x01')
+      geom = GEOSWKBReader_read(
+        wkb_reader, (unsigned char*)c_str, (size_t)RSTRING_LEN(str));
+    else
+      geom = GEOSWKBReader_readHEX(
+        wkb_reader, (unsigned char*)c_str, (size_t)RSTRING_LEN(str));
     if (geom) {
       result = rgeo_wrap_geos_geometry(self, geom, Qnil);
     }
@@ -340,26 +296,25 @@ static VALUE method_factory_parse_wkb(VALUE self, VALUE str)
   return result;
 }
 
-
-static VALUE method_factory_read_for_marshal(VALUE self, VALUE str)
+static VALUE
+method_factory_read_for_marshal(VALUE self, VALUE str)
 {
   RGeo_FactoryData* self_data;
-  GEOSContextHandle_t self_context;
   GEOSWKBReader* wkb_reader;
   VALUE result;
   GEOSGeometry* geom;
 
   Check_Type(str, T_STRING);
   self_data = RGEO_FACTORY_DATA_PTR(self);
-  self_context = self_data->geos_context;
   wkb_reader = self_data->marshal_wkb_reader;
   if (!wkb_reader) {
-    wkb_reader = GEOSWKBReader_create_r(self_context);
+    wkb_reader = GEOSWKBReader_create();
     self_data->marshal_wkb_reader = wkb_reader;
   }
   result = Qnil;
   if (wkb_reader) {
-    geom = GEOSWKBReader_read_r(self_context, wkb_reader, (unsigned char*)RSTRING_PTR(str), (size_t)RSTRING_LEN(str));
+    geom = GEOSWKBReader_read(
+      wkb_reader, (unsigned char*)RSTRING_PTR(str), (size_t)RSTRING_LEN(str));
     if (geom) {
       result = rgeo_wrap_geos_geometry(self, geom, Qnil);
     }
@@ -367,25 +322,24 @@ static VALUE method_factory_read_for_marshal(VALUE self, VALUE str)
   return result;
 }
 
-static VALUE method_factory_read_for_psych(VALUE self, VALUE str)
+static VALUE
+method_factory_read_for_psych(VALUE self, VALUE str)
 {
   RGeo_FactoryData* self_data;
-  GEOSContextHandle_t self_context;
   GEOSWKTReader* wkt_reader;
   VALUE result;
   GEOSGeometry* geom;
 
   Check_Type(str, T_STRING);
   self_data = RGEO_FACTORY_DATA_PTR(self);
-  self_context = self_data->geos_context;
   wkt_reader = self_data->psych_wkt_reader;
   if (!wkt_reader) {
-    wkt_reader = GEOSWKTReader_create_r(self_context);
+    wkt_reader = GEOSWKTReader_create();
     self_data->psych_wkt_reader = wkt_reader;
   }
   result = Qnil;
   if (wkt_reader) {
-    geom = GEOSWKTReader_read_r(self_context, wkt_reader, RSTRING_PTR(str));
+    geom = GEOSWKTReader_read(wkt_reader, RSTRING_PTR(str));
     if (geom) {
       result = rgeo_wrap_geos_geometry(self, geom, Qnil);
     }
@@ -397,10 +351,10 @@ static VALUE method_factory_read_for_psych(VALUE self, VALUE str)
 static VALUE marshal_wkb_generator;
 #endif
 
-static VALUE method_factory_write_for_marshal(VALUE self, VALUE obj)
+static VALUE
+method_factory_write_for_marshal(VALUE self, VALUE obj)
 {
   RGeo_FactoryData* self_data;
-  GEOSContextHandle_t self_context;
   GEOSWKBWriter* wkb_writer;
   const GEOSGeometry* geom;
   VALUE result;
@@ -409,23 +363,23 @@ static VALUE method_factory_write_for_marshal(VALUE self, VALUE obj)
   char has_3d;
 
   self_data = RGEO_FACTORY_DATA_PTR(self);
-  self_context = self_data->geos_context;
   has_3d = self_data->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M;
 #ifndef RGEO_GEOS_SUPPORTS_SETOUTPUTDIMENSION
   if (has_3d) {
     if (NIL_P(marshal_wkb_generator)) {
-      marshal_wkb_generator = rb_funcall(
-        rb_const_get_at(rgeo_geos_module, rb_intern("Utils")),
-        rb_intern("marshal_wkb_generator"), 0);
+      marshal_wkb_generator =
+        rb_funcall(rb_const_get_at(rgeo_geos_module, rb_intern("Utils")),
+                   rb_intern("marshal_wkb_generator"),
+                   0);
     }
     return rb_funcall(marshal_wkb_generator, rb_intern("generate"), 1, obj);
   }
 #endif
   wkb_writer = self_data->marshal_wkb_writer;
   if (!wkb_writer) {
-    wkb_writer = GEOSWKBWriter_create_r(self_context);
+    wkb_writer = GEOSWKBWriter_create();
     if (has_3d) {
-      GEOSWKBWriter_setOutputDimension_r(self_context, wkb_writer, 3);
+      GEOSWKBWriter_setOutputDimension(wkb_writer, 3);
     }
     self_data->marshal_wkb_writer = wkb_writer;
   }
@@ -433,10 +387,10 @@ static VALUE method_factory_write_for_marshal(VALUE self, VALUE obj)
   if (wkb_writer) {
     geom = rgeo_get_geos_geometry_safe(obj);
     if (geom) {
-      str = (char*)GEOSWKBWriter_write_r(self_context, wkb_writer, geom, &size);
+      str = (char*)GEOSWKBWriter_write(wkb_writer, geom, &size);
       if (str) {
         result = rb_str_new(str, size);
-        GEOSFree_r(self_context, str);
+        GEOSFree(str);
       }
     }
   }
@@ -447,10 +401,10 @@ static VALUE method_factory_write_for_marshal(VALUE self, VALUE obj)
 static VALUE psych_wkt_generator;
 #endif
 
-static VALUE method_factory_write_for_psych(VALUE self, VALUE obj)
+static VALUE
+method_factory_write_for_psych(VALUE self, VALUE obj)
 {
   RGeo_FactoryData* self_data;
-  GEOSContextHandle_t self_context;
   GEOSWKTWriter* wkt_writer;
   const GEOSGeometry* geom;
   VALUE result;
@@ -458,23 +412,24 @@ static VALUE method_factory_write_for_psych(VALUE self, VALUE obj)
   char has_3d;
 
   self_data = RGEO_FACTORY_DATA_PTR(self);
-  self_context = self_data->geos_context;
   has_3d = self_data->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M;
 #ifndef RGEO_GEOS_SUPPORTS_SETOUTPUTDIMENSION
   if (has_3d) {
     if (NIL_P(psych_wkt_generator)) {
-      psych_wkt_generator = rb_funcall(
-        rb_const_get_at(rgeo_geos_module, rb_intern("Utils")),
-        rb_intern("psych_wkt_generator"), 0);
+      psych_wkt_generator =
+        rb_funcall(rb_const_get_at(rgeo_geos_module, rb_intern("Utils")),
+                   rb_intern("psych_wkt_generator"),
+                   0);
     }
     return rb_funcall(psych_wkt_generator, rb_intern("generate"), 1, obj);
   }
 #endif
   wkt_writer = self_data->psych_wkt_writer;
   if (!wkt_writer) {
-    wkt_writer = GEOSWKTWriter_create_r(self_context);
+    wkt_writer = GEOSWKTWriter_create();
+    GEOSWKTWriter_setTrim(wkt_writer, 1);
     if (has_3d) {
-      GEOSWKTWriter_setOutputDimension_r(self_context, wkt_writer, 3);
+      GEOSWKTWriter_setOutputDimension(wkt_writer, 3);
     }
     self_data->psych_wkt_writer = wkt_writer;
   }
@@ -482,24 +437,24 @@ static VALUE method_factory_write_for_psych(VALUE self, VALUE obj)
   if (wkt_writer) {
     geom = rgeo_get_geos_geometry_safe(obj);
     if (geom) {
-      str = GEOSWKTWriter_write_r(self_context, wkt_writer, geom);
+      str = GEOSWKTWriter_write(wkt_writer, geom);
       if (str) {
         result = rb_str_new2(str);
-        GEOSFree_r(self_context, str);
+        GEOSFree(str);
       }
     }
   }
   return result;
 }
 
-
-static VALUE cmethod_factory_geos_version(VALUE klass)
+static VALUE
+cmethod_factory_geos_version(VALUE klass)
 {
   return rb_str_new2(GEOS_VERSION);
 }
 
-
-static VALUE cmethod_factory_supports_unary_union(VALUE klass)
+static VALUE
+cmethod_factory_supports_unary_union(VALUE klass)
 {
 #ifdef RGEO_GEOS_SUPPORTS_UNARYUNION
   return Qtrue;
@@ -508,103 +463,93 @@ static VALUE cmethod_factory_supports_unary_union(VALUE klass)
 #endif
 }
 
-static VALUE cmethod_factory_create(VALUE klass, VALUE flags, VALUE srid, VALUE buffer_resolution,
-  VALUE wkt_generator, VALUE wkb_generator, VALUE proj4_obj, VALUE coord_sys_obj)
+static VALUE
+cmethod_factory_create(VALUE klass,
+                       VALUE flags,
+                       VALUE srid,
+                       VALUE buffer_resolution,
+                       VALUE wkt_generator,
+                       VALUE wkb_generator,
+                       VALUE coord_sys_obj)
 {
   VALUE result;
   RGeo_FactoryData* data;
-  GEOSContextHandle_t context;
 
   result = Qnil;
   data = ALLOC(RGeo_FactoryData);
   if (data) {
-    context = GEOS_init_r();
-#ifdef RGEO_GEOS_DEBUG
-    GEOSContext_setNoticeHandler_r(context, notice_handler);
-#endif
-    GEOSContext_setErrorHandler_r(context, error_handler);
-
-    if (context) {
-      data->geos_context = context;
-      data->flags = NUM2INT(flags);
-      data->srid = NUM2INT(srid);
-      data->buffer_resolution = NUM2INT(buffer_resolution);
-      data->wkt_reader = NULL;
-      data->wkb_reader = NULL;
-      data->wkt_writer = NULL;
-      data->wkb_writer = NULL;
-      data->psych_wkt_reader = NULL;
-      data->marshal_wkb_reader = NULL;
-      data->psych_wkt_writer = NULL;
-      data->marshal_wkb_writer = NULL;
-      data->wkrep_wkt_generator = wkt_generator;
-      data->wkrep_wkb_generator = wkb_generator;
-      data->wkrep_wkt_parser = Qnil;
-      data->wkrep_wkb_parser = Qnil;
-      data->proj4_obj = proj4_obj;
-      data->coord_sys_obj = coord_sys_obj;
-      result = TypedData_Wrap_Struct(klass, &rgeo_factory_type, data);
-    }
-    else {
-      free(data);
-    }
+    data->flags = RB_NUM2INT(flags);
+    data->srid = RB_NUM2INT(srid);
+    data->buffer_resolution = RB_NUM2INT(buffer_resolution);
+    data->wkt_reader = NULL;
+    data->wkb_reader = NULL;
+    data->wkt_writer = NULL;
+    data->wkb_writer = NULL;
+    data->psych_wkt_reader = NULL;
+    data->marshal_wkb_reader = NULL;
+    data->psych_wkt_writer = NULL;
+    data->marshal_wkb_writer = NULL;
+    data->wkrep_wkt_generator = wkt_generator;
+    data->wkrep_wkb_generator = wkb_generator;
+    data->wkrep_wkt_parser = Qnil;
+    data->wkrep_wkb_parser = Qnil;
+    data->coord_sys_obj = coord_sys_obj;
+    result = TypedData_Wrap_Struct(klass, &rgeo_factory_type, data);
   }
   return result;
 }
 
-
-static VALUE alloc_factory(VALUE klass)
+static VALUE
+alloc_factory(VALUE klass)
 {
-  return cmethod_factory_create(klass, INT2NUM(0), INT2NUM(0), INT2NUM(0), Qnil, Qnil, Qnil, Qnil);
+  return cmethod_factory_create(
+    klass, INT2NUM(0), INT2NUM(0), INT2NUM(0), Qnil, Qnil, Qnil);
 }
 
-
-static VALUE method_factory_initialize_copy(VALUE self, VALUE orig)
+static VALUE
+method_factory_initialize_copy(VALUE self, VALUE orig)
 {
   RGeo_FactoryData* self_data;
   RGeo_FactoryData* orig_data;
-  GEOSContextHandle_t context;
 
   // Clear out existing data
   self_data = RGEO_FACTORY_DATA_PTR(self);
-  context = self_data->geos_context;
   if (self_data->wkt_reader) {
-    GEOSWKTReader_destroy_r(context, self_data->wkt_reader);
+    GEOSWKTReader_destroy(self_data->wkt_reader);
     self_data->wkt_reader = NULL;
   }
   if (self_data->wkb_reader) {
-    GEOSWKBReader_destroy_r(context, self_data->wkb_reader);
+    GEOSWKBReader_destroy(self_data->wkb_reader);
     self_data->wkb_reader = NULL;
   }
   if (self_data->wkt_writer) {
-    GEOSWKTWriter_destroy_r(context, self_data->wkt_writer);
+    GEOSWKTWriter_destroy(self_data->wkt_writer);
     self_data->wkt_writer = NULL;
   }
   if (self_data->wkb_writer) {
-    GEOSWKBWriter_destroy_r(context, self_data->wkb_writer);
+    GEOSWKBWriter_destroy(self_data->wkb_writer);
     self_data->wkb_writer = NULL;
   }
   if (self_data->psych_wkt_reader) {
-    GEOSWKTReader_destroy_r(context, self_data->psych_wkt_reader);
+    GEOSWKTReader_destroy(self_data->psych_wkt_reader);
     self_data->psych_wkt_reader = NULL;
   }
   if (self_data->marshal_wkb_reader) {
-    GEOSWKBReader_destroy_r(context, self_data->marshal_wkb_reader);
+    GEOSWKBReader_destroy(self_data->marshal_wkb_reader);
     self_data->marshal_wkb_reader = NULL;
   }
   if (self_data->psych_wkt_writer) {
-    GEOSWKTWriter_destroy_r(context, self_data->psych_wkt_writer);
+    GEOSWKTWriter_destroy(self_data->psych_wkt_writer);
     self_data->psych_wkt_writer = NULL;
   }
   if (self_data->marshal_wkb_writer) {
-    GEOSWKBWriter_destroy_r(context, self_data->marshal_wkb_writer);
+    GEOSWKBWriter_destroy(self_data->marshal_wkb_writer);
     self_data->marshal_wkb_writer = NULL;
   }
   self_data->wkrep_wkt_generator = Qnil;
   self_data->wkrep_wkb_generator = Qnil;
   self_data->wkrep_wkt_parser = Qnil;
   self_data->wkrep_wkb_parser = Qnil;
-  self_data->proj4_obj = Qnil;
   self_data->coord_sys_obj = Qnil;
 
   // Copy new data from original object
@@ -617,14 +562,13 @@ static VALUE method_factory_initialize_copy(VALUE self, VALUE orig)
     self_data->wkrep_wkb_generator = orig_data->wkrep_wkb_generator;
     self_data->wkrep_wkt_parser = orig_data->wkrep_wkt_parser;
     self_data->wkrep_wkb_parser = orig_data->wkrep_wkb_parser;
-    self_data->proj4_obj = orig_data->proj4_obj;
     self_data->coord_sys_obj = orig_data->coord_sys_obj;
   }
   return self;
 }
 
-
-static VALUE method_set_wkrep_parsers(VALUE self, VALUE wkt_parser, VALUE wkb_parser)
+static VALUE
+method_set_wkrep_parsers(VALUE self, VALUE wkt_parser, VALUE wkb_parser)
 {
   RGeo_FactoryData* self_data;
 
@@ -635,96 +579,126 @@ static VALUE method_set_wkrep_parsers(VALUE self, VALUE wkt_parser, VALUE wkb_pa
   return self;
 }
 
-
-static VALUE method_get_proj4(VALUE self)
-{
-  return RGEO_FACTORY_DATA_PTR(self)->proj4_obj;
-}
-
-
-static VALUE method_get_coord_sys(VALUE self)
+static VALUE
+method_get_coord_sys(VALUE self)
 {
   return RGEO_FACTORY_DATA_PTR(self)->coord_sys_obj;
 }
 
-
-static VALUE method_get_wkt_generator(VALUE self)
+static VALUE
+method_get_wkt_generator(VALUE self)
 {
   return RGEO_FACTORY_DATA_PTR(self)->wkrep_wkt_generator;
 }
 
-
-static VALUE method_get_wkb_generator(VALUE self)
+static VALUE
+method_get_wkb_generator(VALUE self)
 {
   return RGEO_FACTORY_DATA_PTR(self)->wkrep_wkb_generator;
 }
 
-
-static VALUE method_get_wkt_parser(VALUE self)
+static VALUE
+method_get_wkt_parser(VALUE self)
 {
   return RGEO_FACTORY_DATA_PTR(self)->wkrep_wkt_parser;
 }
 
-
-static VALUE method_get_wkb_parser(VALUE self)
+static VALUE
+method_get_wkb_parser(VALUE self)
 {
   return RGEO_FACTORY_DATA_PTR(self)->wkrep_wkb_parser;
 }
 
-
-static VALUE alloc_geometry(VALUE klass)
+static VALUE
+alloc_geometry(VALUE klass)
 {
   return rgeo_wrap_geos_geometry(Qnil, NULL, klass);
 }
 
-
 /**** INITIALIZATION FUNCTION ****/
 
-
-void rgeo_init_geos_factory()
+void
+rgeo_init_geos_factory()
 {
   VALUE geos_factory_class;
 
 #ifndef RGEO_GEOS_SUPPORTS_SETOUTPUTDIMENSION
-  /* We favor rb_gc_register_address over rb_gc_register_mark_object because the value changes at runtime */
+  /* We favor rb_gc_register_address over rb_gc_register_mark_object because
+   * the value changes at runtime */
   psych_wkt_generator = Qnil;
   rb_gc_register_address(&psych_wkt_generator);
   marshal_wkb_generator = Qnil;
   rb_gc_register_address(&marshal_wkb_generator);
 #endif
 
-  geos_factory_class = rb_define_class_under(rgeo_geos_module, "CAPIFactory", rb_cObject);
+  geos_factory_class =
+    rb_define_class_under(rgeo_geos_module, "CAPIFactory", rb_cObject);
   rb_define_alloc_func(geos_factory_class, alloc_factory);
   // Add C constants to the factory.
-  rb_define_const(geos_factory_class, "FLAG_SUPPORTS_Z", INT2FIX(RGEO_FACTORYFLAGS_SUPPORTS_Z));
-  rb_define_const(geos_factory_class, "FLAG_SUPPORTS_M", INT2FIX(RGEO_FACTORYFLAGS_SUPPORTS_M));
-  rb_define_const(geos_factory_class, "FLAG_SUPPORTS_Z_OR_M", INT2FIX(RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M));
-  rb_define_const(geos_factory_class, "FLAG_PREPARE_HEURISTIC", INT2FIX(RGEO_FACTORYFLAGS_PREPARE_HEURISTIC));
+  rb_define_const(geos_factory_class,
+                  "FLAG_SUPPORTS_Z",
+                  INT2FIX(RGEO_FACTORYFLAGS_SUPPORTS_Z));
+  rb_define_const(geos_factory_class,
+                  "FLAG_SUPPORTS_M",
+                  INT2FIX(RGEO_FACTORYFLAGS_SUPPORTS_M));
+  rb_define_const(geos_factory_class,
+                  "FLAG_SUPPORTS_Z_OR_M",
+                  INT2FIX(RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M));
+  rb_define_const(geos_factory_class,
+                  "FLAG_PREPARE_HEURISTIC",
+                  INT2FIX(RGEO_FACTORYFLAGS_PREPARE_HEURISTIC));
   // Add C methods to the factory.
-  rb_define_method(geos_factory_class, "initialize_copy", method_factory_initialize_copy, 1);
-  rb_define_method(geos_factory_class, "_parse_wkt_impl", method_factory_parse_wkt, 1);
-  rb_define_method(geos_factory_class, "_parse_wkb_impl", method_factory_parse_wkb, 1);
+  rb_define_method(
+    geos_factory_class, "initialize_copy", method_factory_initialize_copy, 1);
+  rb_define_method(
+    geos_factory_class, "_parse_wkt_impl", method_factory_parse_wkt, 1);
+  rb_define_method(
+    geos_factory_class, "_parse_wkb_impl", method_factory_parse_wkb, 1);
   rb_define_method(geos_factory_class, "_srid", method_factory_srid, 0);
-  rb_define_method(geos_factory_class, "_buffer_resolution", method_factory_buffer_resolution, 0);
+  rb_define_method(geos_factory_class,
+                   "_buffer_resolution",
+                   method_factory_buffer_resolution,
+                   0);
   rb_define_method(geos_factory_class, "_flags", method_factory_flags, 0);
-  rb_define_method(geos_factory_class, "supports_z?", method_factory_supports_z_p, 0);
-  rb_define_method(geos_factory_class, "supports_m?", method_factory_supports_m_p, 0);
-  rb_define_method(geos_factory_class, "supports_z_or_m?", method_factory_supports_z_or_m_p, 0);
-  rb_define_method(geos_factory_class, "prepare_heuristic?", method_factory_prepare_heuristic_p, 0);
-  rb_define_method(geos_factory_class, "_set_wkrep_parsers", method_set_wkrep_parsers, 2);
-  rb_define_method(geos_factory_class, "_proj4", method_get_proj4, 0);
+  rb_define_method(
+    geos_factory_class, "supports_z?", method_factory_supports_z_p, 0);
+  rb_define_method(
+    geos_factory_class, "supports_m?", method_factory_supports_m_p, 0);
+  rb_define_method(geos_factory_class,
+                   "supports_z_or_m?",
+                   method_factory_supports_z_or_m_p,
+                   0);
+  rb_define_method(geos_factory_class,
+                   "prepare_heuristic?",
+                   method_factory_prepare_heuristic_p,
+                   0);
+  rb_define_method(
+    geos_factory_class, "_set_wkrep_parsers", method_set_wkrep_parsers, 2);
   rb_define_method(geos_factory_class, "_coord_sys", method_get_coord_sys, 0);
-  rb_define_method(geos_factory_class, "_wkt_generator", method_get_wkt_generator, 0);
-  rb_define_method(geos_factory_class, "_wkb_generator", method_get_wkb_generator, 0);
+  rb_define_method(
+    geos_factory_class, "_wkt_generator", method_get_wkt_generator, 0);
+  rb_define_method(
+    geos_factory_class, "_wkb_generator", method_get_wkb_generator, 0);
   rb_define_method(geos_factory_class, "_wkt_parser", method_get_wkt_parser, 0);
   rb_define_method(geos_factory_class, "_wkb_parser", method_get_wkb_parser, 0);
-  rb_define_method(geos_factory_class, "read_for_marshal", method_factory_read_for_marshal, 1);
-  rb_define_method(geos_factory_class, "write_for_marshal", method_factory_write_for_marshal, 1);
-  rb_define_method(geos_factory_class, "read_for_psych", method_factory_read_for_psych, 1);
-  rb_define_method(geos_factory_class, "write_for_psych", method_factory_write_for_psych, 1);
-  rb_define_module_function(geos_factory_class, "_create", cmethod_factory_create, 7);
-  rb_define_module_function(geos_factory_class, "_geos_version", cmethod_factory_geos_version, 0);
-  rb_define_module_function(geos_factory_class, "_supports_unary_union?", cmethod_factory_supports_unary_union, 0);
+  rb_define_method(
+    geos_factory_class, "read_for_marshal", method_factory_read_for_marshal, 1);
+  rb_define_method(geos_factory_class,
+                   "write_for_marshal",
+                   method_factory_write_for_marshal,
+                   1);
+  rb_define_method(
+    geos_factory_class, "read_for_psych", method_factory_read_for_psych, 1);
+  rb_define_method(
+    geos_factory_class, "write_for_psych", method_factory_write_for_psych, 1);
+  rb_define_module_function(
+    geos_factory_class, "_create", cmethod_factory_create, 6);
+  rb_define_module_function(
+    geos_factory_class, "_geos_version", cmethod_factory_geos_version, 0);
+  rb_define_module_function(geos_factory_class,
+                            "_supports_unary_union?",
+                            cmethod_factory_supports_unary_union,
+                            0);
 
   // Define allocation methods for global class types
   rb_define_alloc_func(rgeo_geos_geometry_class, alloc_geometry);
@@ -739,15 +713,13 @@ void rgeo_init_geos_factory()
   rb_define_alloc_func(rgeo_geos_multi_polygon_class, alloc_geometry);
 }
 
-
 /**** OTHER PUBLIC FUNCTIONS ****/
 
-
-VALUE rgeo_wrap_geos_geometry(VALUE factory, GEOSGeometry* geom, VALUE klass)
+VALUE
+rgeo_wrap_geos_geometry(VALUE factory, GEOSGeometry* geom, VALUE klass)
 {
   VALUE result;
   RGeo_FactoryData* factory_data;
-  GEOSContextHandle_t factory_context;
   VALUE klasses;
   VALUE inferred_klass;
   char is_collection;
@@ -756,14 +728,14 @@ VALUE rgeo_wrap_geos_geometry(VALUE factory, GEOSGeometry* geom, VALUE klass)
   result = Qnil;
   if (geom || !NIL_P(klass)) {
     factory_data = NIL_P(factory) ? NULL : RGEO_FACTORY_DATA_PTR(factory);
-    factory_context = factory_data ? factory_data->geos_context : NULL;
 
     // We don't allow "empty" points, so replace such objects with
     // an empty collection.
     if (geom && factory) {
-      if (GEOSGeomTypeId_r(factory_context, geom) == GEOS_POINT && GEOSGetNumCoordinates_r(factory_context, geom) == 0) {
-        GEOSGeom_destroy_r(factory_context, geom);
-        geom = GEOSGeom_createCollection_r(factory_context, GEOS_GEOMETRYCOLLECTION, NULL, 0);
+      if (GEOSGeomTypeId(geom) == GEOS_POINT &&
+          GEOSGetNumCoordinates(geom) == 0) {
+        GEOSGeom_destroy(geom);
+        geom = GEOSGeom_createCollection(GEOS_GEOMETRYCOLLECTION, NULL, 0);
         klass = rgeo_geos_geometry_collection_class;
       }
     }
@@ -772,38 +744,38 @@ VALUE rgeo_wrap_geos_geometry(VALUE factory, GEOSGeometry* geom, VALUE klass)
     if (TYPE(klass) != T_CLASS) {
       inferred_klass = Qnil;
       is_collection = 0;
-      switch (GEOSGeomTypeId_r(factory_context, geom)) {
-      case GEOS_POINT:
-        inferred_klass = rgeo_geos_point_class;
-        break;
-      case GEOS_LINESTRING:
-        inferred_klass = rgeo_geos_line_string_class;
-        break;
-      case GEOS_LINEARRING:
-        inferred_klass = rgeo_geos_linear_ring_class;
-        break;
-      case GEOS_POLYGON:
-        inferred_klass = rgeo_geos_polygon_class;
-        break;
-      case GEOS_MULTIPOINT:
-        inferred_klass = rgeo_geos_multi_point_class;
-        is_collection = 1;
-        break;
-      case GEOS_MULTILINESTRING:
-        inferred_klass = rgeo_geos_multi_line_string_class;
-        is_collection = 1;
-        break;
-      case GEOS_MULTIPOLYGON:
-        inferred_klass = rgeo_geos_multi_polygon_class;
-        is_collection = 1;
-        break;
-      case GEOS_GEOMETRYCOLLECTION:
-        inferred_klass = rgeo_geos_geometry_collection_class;
-        is_collection = 1;
-        break;
-      default:
-        inferred_klass = rgeo_geos_geometry_class;
-        break;
+      switch (GEOSGeomTypeId(geom)) {
+        case GEOS_POINT:
+          inferred_klass = rgeo_geos_point_class;
+          break;
+        case GEOS_LINESTRING:
+          inferred_klass = rgeo_geos_line_string_class;
+          break;
+        case GEOS_LINEARRING:
+          inferred_klass = rgeo_geos_linear_ring_class;
+          break;
+        case GEOS_POLYGON:
+          inferred_klass = rgeo_geos_polygon_class;
+          break;
+        case GEOS_MULTIPOINT:
+          inferred_klass = rgeo_geos_multi_point_class;
+          is_collection = 1;
+          break;
+        case GEOS_MULTILINESTRING:
+          inferred_klass = rgeo_geos_multi_line_string_class;
+          is_collection = 1;
+          break;
+        case GEOS_MULTIPOLYGON:
+          inferred_klass = rgeo_geos_multi_polygon_class;
+          is_collection = 1;
+          break;
+        case GEOS_GEOMETRYCOLLECTION:
+          inferred_klass = rgeo_geos_geometry_collection_class;
+          is_collection = 1;
+          break;
+        default:
+          inferred_klass = rgeo_geos_geometry_class;
+          break;
       }
       if (TYPE(klass) == T_ARRAY && is_collection) {
         klasses = klass;
@@ -813,12 +785,14 @@ VALUE rgeo_wrap_geos_geometry(VALUE factory, GEOSGeometry* geom, VALUE klass)
     data = ALLOC(RGeo_GeometryData);
     if (data) {
       if (geom) {
-        GEOSSetSRID_r(factory_context, geom, factory_data->srid);
+        GEOSSetSRID(geom, factory_data->srid);
       }
-      data->geos_context = factory_context;
       data->geom = geom;
-      data->prep = factory_data && ((factory_data->flags & RGEO_FACTORYFLAGS_PREPARE_HEURISTIC) != 0) ?
-        (GEOSPreparedGeometry*)1 : NULL;
+      data->prep =
+        factory_data &&
+            ((factory_data->flags & RGEO_FACTORYFLAGS_PREPARE_HEURISTIC) != 0)
+          ? (GEOSPreparedGeometry*)1
+          : NULL;
       data->factory = factory;
       data->klasses = klasses;
       result = TypedData_Wrap_Struct(klass, &rgeo_geometry_type, data);
@@ -827,14 +801,17 @@ VALUE rgeo_wrap_geos_geometry(VALUE factory, GEOSGeometry* geom, VALUE klass)
   return result;
 }
 
-VALUE rgeo_wrap_geos_geometry_clone(VALUE factory, const GEOSGeometry* geom, VALUE klass)
+VALUE
+rgeo_wrap_geos_geometry_clone(VALUE factory,
+                              const GEOSGeometry* geom,
+                              VALUE klass)
 {
   VALUE result;
   GEOSGeometry* clone_geom;
 
   result = Qnil;
   if (geom) {
-    clone_geom = GEOSGeom_clone_r(RGEO_FACTORY_DATA_PTR(factory)->geos_context, geom);
+    clone_geom = GEOSGeom_clone(geom);
     if (clone_geom) {
       result = rgeo_wrap_geos_geometry(factory, clone_geom, klass);
     }
@@ -842,16 +819,17 @@ VALUE rgeo_wrap_geos_geometry_clone(VALUE factory, const GEOSGeometry* geom, VAL
   return result;
 }
 
-
-const GEOSGeometry* rgeo_convert_to_geos_geometry(VALUE factory, VALUE obj, VALUE type)
+const GEOSGeometry*
+rgeo_convert_to_geos_geometry(VALUE factory, VALUE obj, VALUE type)
 {
   VALUE object;
 
-  if (NIL_P(type) && RGEO_GEOMETRY_TYPEDDATA_P(obj) && RGEO_GEOMETRY_DATA_PTR(obj)->factory == factory) {
+  if (NIL_P(type) && RGEO_GEOMETRY_TYPEDDATA_P(obj) &&
+      RGEO_GEOMETRY_DATA_PTR(obj)->factory == factory) {
     object = obj;
-  }
-  else {
-    object = rb_funcall(rgeo_feature_module, rb_intern("cast"), 3, obj, factory, type);
+  } else {
+    object =
+      rb_funcall(rgeo_feature_module, rb_intern("cast"), 3, obj, factory, type);
   }
   if (NIL_P(object))
     return NULL;
@@ -860,62 +838,81 @@ const GEOSGeometry* rgeo_convert_to_geos_geometry(VALUE factory, VALUE obj, VALU
   return RGEO_GEOMETRY_DATA_PTR(object)->geom;
 }
 
-
-GEOSGeometry* rgeo_convert_to_detached_geos_geometry(VALUE obj, VALUE factory, VALUE type, VALUE* klasses)
+GEOSGeometry*
+rgeo_convert_to_detached_geos_geometry(VALUE obj,
+                                       VALUE factory,
+                                       VALUE type,
+                                       VALUE* klasses,
+                                       int* state)
 {
   VALUE object;
   GEOSGeometry* geom;
   RGeo_GeometryData* object_data;
   const GEOSPreparedGeometry* prep;
-
   if (klasses) {
     *klasses = Qnil;
   }
-  object = rb_funcall(rgeo_feature_module, rb_intern("cast"), 5, obj, factory, type, ID2SYM(rb_intern("force_new")), ID2SYM(rb_intern("keep_subtype")));
-  geom = NULL;
-  if (!NIL_P(object)) {
-    object_data = RGEO_GEOMETRY_DATA_PTR(object);
-    geom = object_data->geom;
-    if (klasses) {
-      *klasses = object_data->klasses;
-      if (NIL_P(*klasses)) {
-        *klasses = CLASS_OF(object);
-      }
-    }
-    prep = object_data->prep;
-    if (prep && prep != (GEOSPreparedGeometry*)1 && prep != (GEOSPreparedGeometry*)2) {
-      GEOSPreparedGeom_destroy_r(object_data->geos_context, prep);
+
+  object = rb_protect_funcall(rgeo_feature_module,
+                              rb_intern("cast"),
+                              state,
+                              5,
+                              obj,
+                              factory,
+                              type,
+                              ID2SYM(rb_intern("force_new")),
+                              ID2SYM(rb_intern("keep_subtype")));
+  if (*state || NIL_P(object)) {
+    return NULL;
+  }
+
+  object_data = RGEO_GEOMETRY_DATA_PTR(object);
+  geom = object_data->geom;
+  if (klasses) {
+    *klasses = object_data->klasses;
+    if (NIL_P(*klasses)) {
+      *klasses = CLASS_OF(object);
     }
-    object_data->geos_context = NULL;
-    object_data->geom = NULL;
-    object_data->prep = NULL;
-    object_data->factory = Qnil;
-    object_data->klasses = Qnil;
   }
+  prep = object_data->prep;
+  if (prep && prep != (GEOSPreparedGeometry*)1 &&
+      prep != (GEOSPreparedGeometry*)2) {
+    GEOSPreparedGeom_destroy(prep);
+  }
+  object_data->geom = NULL;
+  object_data->prep = NULL;
+  object_data->factory = Qnil;
+  object_data->klasses = Qnil;
+
   return geom;
 }
 
-
-char rgeo_is_geos_object(VALUE obj)
+char
+rgeo_is_geos_object(VALUE obj)
 {
   return RGEO_GEOMETRY_TYPEDDATA_P(obj) ? 1 : 0;
 }
 
-void rgeo_check_geos_object(VALUE obj)
+void
+rgeo_check_geos_object(VALUE obj)
 {
   if (!rgeo_is_geos_object(obj)) {
     rb_raise(rb_eRGeoError, "Not a GEOS Geometry object.");
   }
 }
 
-
-const GEOSGeometry* rgeo_get_geos_geometry_safe(VALUE obj)
+const GEOSGeometry*
+rgeo_get_geos_geometry_safe(VALUE obj)
 {
-  return RGEO_GEOMETRY_TYPEDDATA_P(obj) ? (const GEOSGeometry*)(RGEO_GEOMETRY_DATA_PTR(obj)->geom) : NULL;
+  return RGEO_GEOMETRY_TYPEDDATA_P(obj)
+           ? (const GEOSGeometry*)(RGEO_GEOMETRY_DATA_PTR(obj)->geom)
+           : NULL;
 }
 
-
-VALUE rgeo_geos_coordseqs_eql(GEOSContextHandle_t context, const GEOSGeometry* geom1, const GEOSGeometry* geom2, char check_z)
+VALUE
+rgeo_geos_coordseqs_eql(const GEOSGeometry* geom1,
+                        const GEOSGeometry* geom2,
+                        char check_z)
 {
   VALUE result;
   const GEOSCoordSequence* cs1;
@@ -927,27 +924,30 @@ VALUE rgeo_geos_coordseqs_eql(GEOSContextHandle_t context, const GEOSGeometry* g
 
   result = Qnil;
   if (geom1 && geom2) {
-    cs1 = GEOSGeom_getCoordSeq_r(context, geom1);
-    cs2 = GEOSGeom_getCoordSeq_r(context, geom2);
+    cs1 = GEOSGeom_getCoordSeq(geom1);
+    cs2 = GEOSGeom_getCoordSeq(geom2);
     if (cs1 && cs2) {
       len1 = 0;
       len2 = 0;
-      if (GEOSCoordSeq_getSize_r(context, cs1, &len1) && GEOSCoordSeq_getSize_r(context, cs2, &len2)) {
+      if (GEOSCoordSeq_getSize(cs1, &len1) &&
+          GEOSCoordSeq_getSize(cs2, &len2)) {
         if (len1 == len2) {
           result = Qtrue;
-          for (i=0; i<len1; ++i) {
-            if (GEOSCoordSeq_getX_r(context, cs1, i, &val1) && GEOSCoordSeq_getX_r(context, cs2, i, &val2)) {
+          for (i = 0; i < len1; ++i) {
+            if (GEOSCoordSeq_getX(cs1, i, &val1) &&
+                GEOSCoordSeq_getX(cs2, i, &val2)) {
               if (val1 == val2) {
-                if (GEOSCoordSeq_getY_r(context, cs1, i, &val1) && GEOSCoordSeq_getY_r(context, cs2, i, &val2)) {
+                if (GEOSCoordSeq_getY(cs1, i, &val1) &&
+                    GEOSCoordSeq_getY(cs2, i, &val2)) {
                   if (val1 == val2) {
                     if (check_z) {
                       val1 = 0;
-                      if (!GEOSCoordSeq_getZ_r(context, cs1, i, &val1)) {
+                      if (!GEOSCoordSeq_getZ(cs1, i, &val1)) {
                         result = Qnil;
                         break;
                       }
                       val2 = 0;
-                      if (!GEOSCoordSeq_getZ_r(context, cs2, i, &val2)) {
+                      if (!GEOSCoordSeq_getZ(cs2, i, &val2)) {
                         result = Qnil;
                         break;
                       }
@@ -956,29 +956,24 @@ VALUE rgeo_geos_coordseqs_eql(GEOSContextHandle_t context, const GEOSGeometry* g
                         break;
                       }
                     }
-                  }
-                  else {  // Y coords are different
+                  } else { // Y coords are different
                     result = Qfalse;
                     break;
                   }
-                }
-                else {  // Failed to get Y coords
+                } else { // Failed to get Y coords
                   result = Qnil;
                   break;
                 }
-              }
-              else {  // X coords are different
+              } else { // X coords are different
                 result = Qfalse;
                 break;
               }
-            }
-            else {  // Failed to get X coords
+            } else { // Failed to get X coords
               result = Qnil;
               break;
             }
-          }  // Iteration over coords
-        }
-        else {  // Lengths are different
+          }      // Iteration over coords
+        } else { // Lengths are different
           result = Qfalse;
         }
       }
@@ -987,8 +982,8 @@ VALUE rgeo_geos_coordseqs_eql(GEOSContextHandle_t context, const GEOSGeometry* g
   return result;
 }
 
-
-VALUE rgeo_geos_klasses_and_factories_eql(VALUE obj1, VALUE obj2)
+VALUE
+rgeo_geos_klasses_and_factories_eql(VALUE obj1, VALUE obj2)
 {
   VALUE result;
   VALUE factory;
@@ -996,24 +991,25 @@ VALUE rgeo_geos_klasses_and_factories_eql(VALUE obj1, VALUE obj2)
   result = Qnil;
   if (rb_obj_class(obj1) != rb_obj_class(obj2)) {
     result = Qfalse;
-  }
-  else {
+  } else {
     factory = RGEO_GEOMETRY_DATA_PTR(obj1)->factory;
     /* No need to cache the internal here (https://ips.fastruby.io/4x) */
-    result = rb_funcall(factory, rb_intern("eql?"), 1, RGEO_GEOMETRY_DATA_PTR(obj2)->factory);
+    result = rb_funcall(
+      factory, rb_intern("eql?"), 1, RGEO_GEOMETRY_DATA_PTR(obj2)->factory);
   }
   return result;
 }
 
-
-typedef struct {
+typedef struct
+{
   st_index_t seed_hash;
   double x;
   double y;
   double z;
 } RGeo_Coordseq_Hash_Struct;
 
-st_index_t rgeo_geos_coordseq_hash(GEOSContextHandle_t context, const GEOSGeometry* geom, st_index_t hash)
+st_index_t
+rgeo_geos_coordseq_hash(const GEOSGeometry* geom, st_index_t hash)
 {
   const GEOSCoordSequence* cs;
   unsigned int len;
@@ -1021,17 +1017,18 @@ st_index_t rgeo_geos_coordseq_hash(GEOSContextHandle_t context, const GEOSGeomet
   RGeo_Coordseq_Hash_Struct hash_struct;
 
   if (geom) {
-    cs = GEOSGeom_getCoordSeq_r(context, geom);
+    cs = GEOSGeom_getCoordSeq(geom);
     if (cs) {
-      if (GEOSCoordSeq_getSize_r(context, cs, &len)) {
-        for (i=0; i<len; ++i) {
-          if (GEOSCoordSeq_getX_r(context, cs, i, &hash_struct.x)) {
-            if (GEOSCoordSeq_getY_r(context, cs, i, &hash_struct.y)) {
-              if (!GEOSCoordSeq_getY_r(context, cs, i, &hash_struct.z)) {
+      if (GEOSCoordSeq_getSize(cs, &len)) {
+        for (i = 0; i < len; ++i) {
+          if (GEOSCoordSeq_getX(cs, i, &hash_struct.x)) {
+            if (GEOSCoordSeq_getY(cs, i, &hash_struct.y)) {
+              if (!GEOSCoordSeq_getY(cs, i, &hash_struct.z)) {
                 hash_struct.z = 0;
               }
               hash_struct.seed_hash = hash;
-              hash = rb_memhash(&hash_struct, sizeof(RGeo_Coordseq_Hash_Struct));
+              hash =
+                rb_memhash(&hash_struct, sizeof(RGeo_Coordseq_Hash_Struct));
             }
           }
         }
@@ -1041,14 +1038,15 @@ st_index_t rgeo_geos_coordseq_hash(GEOSContextHandle_t context, const GEOSGeomet
   return hash;
 }
 
-
-typedef struct {
+typedef struct
+{
   st_index_t seed_hash;
   st_index_t h1;
   st_index_t h2;
 } RGeo_Objbase_Hash_Struct;
 
-st_index_t rgeo_geos_objbase_hash(VALUE factory, VALUE type_module, st_index_t hash)
+st_index_t
+rgeo_geos_objbase_hash(VALUE factory, VALUE type_module, st_index_t hash)
 {
   ID hash_method;
   RGeo_Objbase_Hash_Struct hash_struct;
@@ -1060,23 +1058,6 @@ st_index_t rgeo_geos_objbase_hash(VALUE factory, VALUE type_module, st_index_t h
   return rb_memhash(&hash_struct, sizeof(RGeo_Objbase_Hash_Struct));
 }
 
-
-st_index_t rgeo_internal_memhash(const void* ptr, long len)
-{
-  const char* bytes;
-  st_index_t hval;
-  long i;
-
-  bytes = (const char*)ptr;
-  hval = 0x811c9dc5;
-  for (i=0; i<len; ++i) {
-    hval ^= (unsigned int)(*bytes++);
-    hval *= 0x01000193;
-  }
-  return hval;
-}
-
-
 RGEO_END_C
 
 #endif