rgeo 0.5.0 → 0.5.1

data/ext/geos_c_impl/geometry_collection.h

@@ -1,46 +0,0 @@
-/*
-  Geometry collection methods for GEOS wrapper
-*/
-
-
-#ifndef RGEO_GEOS_GEOMETRY_COLLECTION_INCLUDED
-#define RGEO_GEOS_GEOMETRY_COLLECTION_INCLUDED
-
-#include <ruby.h>
-#include <geos_c.h>
-
-#include "factory.h"
-
-RGEO_BEGIN_C
-
-
-/*
-  Initializes the geometry collection module. This should be called after
-  the geometry module is initialized.
-*/
-void rgeo_init_geos_geometry_collection(RGeo_Globals* globals);
-
-/*
-  Comopares the contents of two geometry collections. Does not test the
-  types of the collections themselves, but tests the types, values, and
-  contents of all the contents. The two given geometries MUST be
-  collection types-- i.e. GeometryCollection, MultiPoint, MultiLineString,
-  or MultiPolygon.
-  Returns Qtrue if the contents of the two geometry collections are equal,
-  Qfalse if they are inequal, or Qnil if an error occurs.
-*/
-VALUE rgeo_geos_geometry_collections_eql(GEOSContextHandle_t context, const GEOSGeometry* geom1, const GEOSGeometry* geom2, char check_z);
-
-/*
-  A tool for building up hash values.
-  You must pass in the context, a geos geometry, and a seed hash.
-  Returns an updated hash.
-  This call is useful in sequence, and should be bracketed by calls to
-  rb_hash_start and rb_hash_end.
-*/
-st_index_t rgeo_geos_geometry_collection_hash(GEOSContextHandle_t context, const GEOSGeometry* geom, st_index_t hash);
-
-
-RGEO_END_C
-
-#endif

data/ext/geos_c_impl/line_string.c

@@ -1,675 +0,0 @@
-/*
-  Line string methods for GEOS wrapper
-*/
-
-
-#include "preface.h"
-
-#ifdef RGEO_GEOS_SUPPORTED
-
-#include <string.h>
-#include <ruby.h>
-#include <geos_c.h>
-
-#include "factory.h"
-#include "geometry.h"
-#include "point.h"
-#include "line_string.h"
-
-#include "coordinates.h"
-
-RGEO_BEGIN_C
-
-
-static VALUE method_line_string_geometry_type(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  if (self_data->geom) {
-    result = RGEO_FACTORY_DATA_PTR(self_data->factory)->globals->feature_line_string;
-  }
-  return result;
-}
-
-
-static VALUE method_linear_ring_geometry_type(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  if (self_data->geom) {
-    result = RGEO_FACTORY_DATA_PTR(self_data->factory)->globals->feature_linear_ring;
-  }
-  return result;
-}
-
-
-static VALUE method_line_geometry_type(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  if (self_data->geom) {
-    result = RGEO_FACTORY_DATA_PTR(self_data->factory)->globals->feature_line;
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_length(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-  const GEOSGeometry* self_geom;
-  double len;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_geom = self_data->geom;
-  if (self_geom) {
-    if (GEOSLength_r(self_data->geos_context, self_geom, &len)) {
-      result = rb_float_new(len);
-    }
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_num_points(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-  const GEOSGeometry* self_geom;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_geom = self_data->geom;
-  if (self_geom) {
-    result = INT2NUM(GEOSGetNumCoordinates_r(self_data->geos_context, self_geom));
-  }
-  return result;
-}
-
-
-static void d(VALUE v) {
-    ID sym_puts = rb_intern("puts");
-    ID sym_inspect = rb_intern("inspect");
-    rb_funcall(rb_mKernel, sym_puts, 1,
-        rb_funcall(v, sym_inspect, 0));
-}
-
-static VALUE method_line_string_coordinates(VALUE self)
-{
-  VALUE result;
-  VALUE point_tuple;
-  RGeo_GeometryData* self_data;
-  const GEOSGeometry* self_geom;
-  const GEOSCoordSequence* coord_sequence;
-  int zCoordinate;
-
-  GEOSContextHandle_t context;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_geom = self_data->geom;
-
-  if (self_geom) {
-    zCoordinate = RGEO_FACTORY_DATA_PTR(self_data->factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M;
-    context = self_data->geos_context;
-    coord_sequence = GEOSGeom_getCoordSeq_r(context, self_geom);
-    if(coord_sequence) {
-      result = extract_points_from_coordinate_sequence(context, coord_sequence, zCoordinate);
-    }
-  }
-  return result;
-}
-
-
-
-static VALUE get_point_from_coordseq(VALUE self, const GEOSCoordSequence* coord_seq, unsigned int i, char has_z)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-  GEOSContextHandle_t self_context;
-  double x, y, z;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_context = self_data->geos_context;
-  if (GEOSCoordSeq_getX_r(self_context, coord_seq, i, &x)) {
-    if (GEOSCoordSeq_getY_r(self_context, coord_seq, i, &y)) {
-      if (has_z) {
-        if (!GEOSCoordSeq_getZ_r(self_context, coord_seq, i, &z)) {
-          z = 0.0;
-        }
-      }
-      else {
-        z = 0.0;
-      }
-      result = rgeo_create_geos_point(self_data->factory, x, y, z);
-    }
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_point_n(VALUE self, VALUE n)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-  const GEOSGeometry* self_geom;
-  GEOSContextHandle_t self_context;
-  const GEOSCoordSequence* coord_seq;
-  char has_z;
-  int si;
-  unsigned int i;
-  unsigned int size;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_geom = self_data->geom;
-  if (self_geom) {
-    self_context = self_data->geos_context;
-    coord_seq = GEOSGeom_getCoordSeq_r(self_context, self_geom);
-    if (coord_seq) {
-      has_z = (char)(RGEO_FACTORY_DATA_PTR(self_data->factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
-      si = NUM2INT(n);
-      if (si >= 0) {
-        i = si;
-        if (GEOSCoordSeq_getSize_r(self_context, coord_seq, &size)) {
-          if (i < size) {
-            result = get_point_from_coordseq(self, coord_seq, i, has_z);
-          }
-        }
-      }
-    }
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_points(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-  const GEOSGeometry* self_geom;
-  GEOSContextHandle_t self_context;
-  const GEOSCoordSequence* coord_seq;
-  char has_z;
-  unsigned int size;
-  double x;
-  double y;
-  double z;
-  unsigned int i;
-  VALUE point;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_geom = self_data->geom;
-  if (self_geom) {
-    self_context = self_data->geos_context;
-    coord_seq = GEOSGeom_getCoordSeq_r(self_context, self_geom);
-    if (coord_seq) {
-      has_z = (char)(RGEO_FACTORY_DATA_PTR(self_data->factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
-      if (GEOSCoordSeq_getSize_r(self_context, coord_seq, &size)) {
-        result = rb_ary_new2(size);
-        for (i=0; i<size; ++i) {
-          point = get_point_from_coordseq(self, coord_seq, i, has_z);
-          if (!NIL_P(point)) {
-            rb_ary_store(result, i, point);
-          }
-        }
-      }
-    }
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_start_point(VALUE self)
-{
-  return method_line_string_point_n(self, INT2NUM(0));
-}
-
-
-static VALUE method_line_string_end_point(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-  const GEOSGeometry* self_geom;
-  unsigned int n;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_geom = self_data->geom;
-  if (self_geom) {
-    n = GEOSGetNumCoordinates_r(self_data->geos_context, self_geom);
-    if (n > 0) {
-      result = method_line_string_point_n(self, INT2NUM(n-1));
-    }
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_is_closed(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-  const GEOSGeometry* self_geom;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_geom = self_data->geom;
-  if (self_geom) {
-    result = rgeo_is_geos_line_string_closed(self_data->geos_context, self_geom);
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_is_ring(VALUE self)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-  const GEOSGeometry* self_geom;
-  char val;
-
-  result = Qnil;
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  self_geom = self_data->geom;
-  if (self_geom) {
-    val = GEOSisRing_r(self_data->geos_context, self_geom);
-    if (val == 0) {
-      result = Qfalse;
-    }
-    else if (val == 1) {
-      result = Qtrue;
-    }
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_eql(VALUE self, VALUE rhs)
-{
-  VALUE result;
-  RGeo_GeometryData* self_data;
-
-  result = rgeo_geos_klasses_and_factories_eql(self, rhs);
-  if (RTEST(result)) {
-    self_data = RGEO_GEOMETRY_DATA_PTR(self);
-    result = rgeo_geos_coordseqs_eql(self_data->geos_context, self_data->geom, RGEO_GEOMETRY_DATA_PTR(rhs)->geom, RGEO_FACTORY_DATA_PTR(self_data->factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
-  }
-  return result;
-}
-
-
-static VALUE method_line_string_hash(VALUE self)
-{
-  st_index_t hash;
-  RGeo_GeometryData* self_data;
-  VALUE factory;
-
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  factory = self_data->factory;
-  hash = rb_hash_start(0);
-  hash = rgeo_geos_objbase_hash(factory,
-    RGEO_FACTORY_DATA_PTR(factory)->globals->feature_line_string, hash);
-  hash = rgeo_geos_coordseq_hash(self_data->geos_context, self_data->geom, hash);
-  return LONG2FIX(rb_hash_end(hash));
-}
-
-
-static VALUE method_linear_ring_hash(VALUE self)
-{
-  st_index_t hash;
-  RGeo_GeometryData* self_data;
-  VALUE factory;
-
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  factory = self_data->factory;
-  hash = rb_hash_start(0);
-  hash = rgeo_geos_objbase_hash(factory,
-    RGEO_FACTORY_DATA_PTR(factory)->globals->feature_linear_ring, hash);
-  hash = rgeo_geos_coordseq_hash(self_data->geos_context, self_data->geom, hash);
-  return LONG2FIX(rb_hash_end(hash));
-}
-
-
-static VALUE method_line_hash(VALUE self)
-{
-  st_index_t hash;
-  RGeo_GeometryData* self_data;
-  VALUE factory;
-
-  self_data = RGEO_GEOMETRY_DATA_PTR(self);
-  factory = self_data->factory;
-  hash = rb_hash_start(0);
-  hash = rgeo_geos_objbase_hash(factory,
-    RGEO_FACTORY_DATA_PTR(factory)->globals->feature_line, hash);
-  hash = rgeo_geos_coordseq_hash(self_data->geos_context, self_data->geom, hash);
-  return LONG2FIX(rb_hash_end(hash));
-}
-
-
-static GEOSCoordSequence* coord_seq_from_array(VALUE factory, VALUE array, char close)
-{
-  RGeo_FactoryData* factory_data;
-  VALUE point_type;
-  unsigned int len;
-  char has_z;
-  unsigned int dims;
-  double* coords;
-  GEOSContextHandle_t context;
-  unsigned int i;
-  char good;
-  const GEOSGeometry* entry_geom;
-  const GEOSCoordSequence* entry_cs;
-  double x;
-  GEOSCoordSequence* coord_seq;
-
-  Check_Type(array, T_ARRAY);
-  factory_data = RGEO_FACTORY_DATA_PTR(factory);
-  point_type = factory_data->globals->feature_point;
-  len = (unsigned int)RARRAY_LEN(array);
-  has_z = (char)(RGEO_FACTORY_DATA_PTR(factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
-  dims = has_z ? 3 : 2;
-  coords = ALLOC_N(double, len == 0 ? 1 : len * dims);
-  if (!coords) {
-    return NULL;
-  }
-  context = factory_data->geos_context;
-  for (i=0; i<len; ++i) {
-    good = 0;
-    entry_geom = rgeo_convert_to_geos_geometry(factory, rb_ary_entry(array, i), point_type);
-    if (entry_geom) {
-      entry_cs = GEOSGeom_getCoordSeq_r(context, entry_geom);
-      if (entry_cs) {
-        if (GEOSCoordSeq_getX_r(context, entry_cs, 0, &x)) {
-          coords[i*dims] = x;
-          if (GEOSCoordSeq_getY_r(context, entry_cs, 0, &x)) {
-            coords[i*dims+1] = x;
-            good = 1;
-            if (has_z) {
-              if (GEOSCoordSeq_getZ_r(context, entry_cs, 0, &x)) {
-                coords[i*dims+2] = x;
-              }
-              else {
-                good = 0;
-              }
-            }
-          }
-        }
-      }
-    }
-    if (!good) {
-      free(coords);
-      return NULL;
-    }
-  }
-  if (len > 0 && close) {
-    if (coords[0] == coords[(len-1)*dims] && coords[1] == coords[(len-1)*dims+1]) {
-      close = 0;
-    }
-  }
-  else {
-    close = 0;
-  }
-  coord_seq = GEOSCoordSeq_create_r(context, len + close, 3);
-  if (coord_seq) {
-    for (i=0; i<len; ++i) {
-      GEOSCoordSeq_setX_r(context, coord_seq, i, coords[i*dims]);
-      GEOSCoordSeq_setY_r(context, coord_seq, i, coords[i*dims+1]);
-      GEOSCoordSeq_setZ_r(context, coord_seq, i, has_z ? coords[i*dims+2] : 0);
-    }
-    if (close) {
-      GEOSCoordSeq_setX_r(context, coord_seq, len, coords[0]);
-      GEOSCoordSeq_setY_r(context, coord_seq, len, coords[1]);
-      GEOSCoordSeq_setZ_r(context, coord_seq, len, has_z ? coords[2] : 0);
-    }
-  }
-  free(coords);
-  return coord_seq;
-}
-
-
-static VALUE cmethod_create_line_string(VALUE module, VALUE factory, VALUE array)
-{
-  VALUE result;
-  GEOSCoordSequence* coord_seq;
-  RGeo_FactoryData* factory_data;
-  GEOSGeometry* geom;
-
-  result = Qnil;
-  coord_seq = coord_seq_from_array(factory, array, 0);
-  if (coord_seq) {
-    factory_data = RGEO_FACTORY_DATA_PTR(factory);
-    geom = GEOSGeom_createLineString_r(factory_data->geos_context, coord_seq);
-    if (geom) {
-      result = rgeo_wrap_geos_geometry(factory, geom, factory_data->globals->geos_line_string);
-    }
-  }
-  return result;
-}
-
-
-static VALUE cmethod_create_linear_ring(VALUE module, VALUE factory, VALUE array)
-{
-  VALUE result;
-  GEOSCoordSequence* coord_seq;
-  RGeo_FactoryData* factory_data;
-  GEOSGeometry* geom;
-
-  result = Qnil;
-  coord_seq = coord_seq_from_array(factory, array, 1);
-  if (coord_seq) {
-    factory_data = RGEO_FACTORY_DATA_PTR(factory);
-    geom = GEOSGeom_createLinearRing_r(factory_data->geos_context, coord_seq);
-    if (geom) {
-      result = rgeo_wrap_geos_geometry(factory, geom, factory_data->globals->geos_linear_ring);
-    }
-  }
-  return result;
-}
-
-
-static void populate_geom_into_coord_seq(GEOSContextHandle_t context, const GEOSGeometry* geom, GEOSCoordSequence* coord_seq, unsigned int i, char has_z)
-{
-  const GEOSCoordSequence* cs;
-  double x;
-
-  cs = GEOSGeom_getCoordSeq_r(context, geom);
-  x = 0;
-  if (cs) {
-    GEOSCoordSeq_getX_r(context, cs, 0, &x);
-  }
-  GEOSCoordSeq_setX_r(context, coord_seq, i, x);
-  x = 0;
-  if (cs) {
-    GEOSCoordSeq_getY_r(context, cs, 0, &x);
-  }
-  GEOSCoordSeq_setY_r(context, coord_seq, i, x);
-  x = 0;
-  if (has_z && cs) {
-    GEOSCoordSeq_getZ_r(context, cs, 0, &x);
-  }
-  GEOSCoordSeq_setZ_r(context, coord_seq, i, x);
-}
-
-
-static VALUE cmethod_create_line(VALUE module, VALUE factory, VALUE start, VALUE end)
-{
-  VALUE result;
-  RGeo_FactoryData* factory_data;
-  char has_z;
-  VALUE point_type;
-  GEOSContextHandle_t context;
-  const GEOSGeometry* start_geom;
-  const GEOSGeometry* end_geom;
-  GEOSCoordSequence* coord_seq;
-  GEOSGeometry* geom;
-
-  result = Qnil;
-  factory_data = RGEO_FACTORY_DATA_PTR(factory);
-  has_z = (char)(factory_data->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
-  point_type = factory_data->globals->feature_point;
-  context = factory_data->geos_context;
-
-  start_geom = rgeo_convert_to_geos_geometry(factory, start, point_type);
-  if (start_geom) {
-    end_geom = rgeo_convert_to_geos_geometry(factory, end, point_type);
-    if (end_geom) {
-      coord_seq = GEOSCoordSeq_create_r(context, 2, 3);
-      if (coord_seq) {
-        populate_geom_into_coord_seq(context, start_geom, coord_seq, 0, has_z);
-        populate_geom_into_coord_seq(context, end_geom, coord_seq, 1, has_z);
-        geom = GEOSGeom_createLineString_r(context, coord_seq);
-        if (geom) {
-          result = rgeo_wrap_geos_geometry(factory, geom, factory_data->globals->geos_line);
-        }
-      }
-    }
-  }
-
-  return result;
-}
-
-
-static VALUE impl_copy_from(VALUE klass, VALUE factory, VALUE original, char subtype)
-{
-  VALUE result;
-  const GEOSGeometry* original_geom;
-  GEOSContextHandle_t context;
-  const GEOSCoordSequence* original_coord_seq;
-  GEOSCoordSequence* coord_seq;
-  GEOSGeometry* geom;
-
-  result = Qnil;
-  original_geom = RGEO_GEOMETRY_DATA_PTR(original)->geom;
-  if (original_geom) {
-    context = RGEO_FACTORY_DATA_PTR(factory)->geos_context;
-    if (subtype == 1 && GEOSGetNumCoordinates_r(context, original_geom) != 2) {
-      original_geom = NULL;
-    }
-    if (original_geom) {
-      original_coord_seq = GEOSGeom_getCoordSeq_r(context, original_geom);
-      if (original_coord_seq) {
-        coord_seq = GEOSCoordSeq_clone_r(context, original_coord_seq);
-        if (coord_seq) {
-          geom = subtype == 2 ? GEOSGeom_createLinearRing_r(context, coord_seq) : GEOSGeom_createLineString_r(context, coord_seq);
-          if (geom) {
-            result = rgeo_wrap_geos_geometry(factory, geom, klass);
-          }
-        }
-      }
-    }
-  }
-  return result;
-}
-
-
-static VALUE cmethod_line_string_copy_from(VALUE klass, VALUE factory, VALUE original)
-{
-  return impl_copy_from(klass, factory, original, 0);
-}
-
-
-static VALUE cmethod_line_copy_from(VALUE klass, VALUE factory, VALUE original)
-{
-  return impl_copy_from(klass, factory, original, 1);
-}
-
-
-static VALUE cmethod_linear_ring_copy_from(VALUE klass, VALUE factory, VALUE original)
-{
-  return impl_copy_from(klass, factory, original, 2);
-}
-
-
-void rgeo_init_geos_line_string(RGeo_Globals* globals)
-{
-  VALUE geos_line_string_methods;
-  VALUE geos_linear_ring_methods;
-  VALUE geos_line_methods;
-
-  // Class methods for CAPILineStringImpl
-  rb_define_module_function(globals->geos_line_string, "create", cmethod_create_line_string, 2);
-  rb_define_module_function(globals->geos_line_string, "_copy_from", cmethod_line_string_copy_from, 2);
-
-  // Class methods for CAPILinearRingImpl
-  rb_define_module_function(globals->geos_linear_ring, "create", cmethod_create_linear_ring, 2);
-  rb_define_module_function(globals->geos_linear_ring, "_copy_from", cmethod_linear_ring_copy_from, 2);
-
-  // Class methods for CAPILineImpl
-  rb_define_module_function(globals->geos_line, "create", cmethod_create_line, 3);
-  rb_define_module_function(globals->geos_line, "_copy_from", cmethod_line_copy_from, 2);
-
-  // CAPILineStringMethods module
-  geos_line_string_methods = rb_define_module_under(globals->geos_module, "CAPILineStringMethods");
-  rb_define_method(geos_line_string_methods, "rep_equals?", method_line_string_eql, 1);
-  rb_define_method(geos_line_string_methods, "eql?", method_line_string_eql, 1);
-  rb_define_method(geos_line_string_methods, "hash", method_line_string_hash, 0);
-  rb_define_method(geos_line_string_methods, "geometry_type", method_line_string_geometry_type, 0);
-  rb_define_method(geos_line_string_methods, "length", method_line_string_length, 0);
-  rb_define_method(geos_line_string_methods, "num_points", method_line_string_num_points, 0);
-  rb_define_method(geos_line_string_methods, "point_n", method_line_string_point_n, 1);
-  rb_define_method(geos_line_string_methods, "points", method_line_string_points, 0);
-  rb_define_method(geos_line_string_methods, "start_point", method_line_string_start_point, 0);
-  rb_define_method(geos_line_string_methods, "end_point", method_line_string_end_point, 0);
-  rb_define_method(geos_line_string_methods, "is_closed?", method_line_string_is_closed, 0);
-  rb_define_method(geos_line_string_methods, "is_ring?", method_line_string_is_ring, 0);
-  rb_define_method(geos_line_string_methods, "coordinates", method_line_string_coordinates, 0);
-
-  // CAPILinearRingMethods module
-  geos_linear_ring_methods = rb_define_module_under(globals->geos_module, "CAPILinearRingMethods");
-  rb_define_method(geos_linear_ring_methods, "geometry_type", method_linear_ring_geometry_type, 0);
-  rb_define_method(geos_linear_ring_methods, "hash", method_linear_ring_hash, 0);
-
-  // CAPILineMethods module
-  geos_line_methods = rb_define_module_under(globals->geos_module, "CAPILineMethods");
-  rb_define_method(geos_line_methods, "geometry_type", method_line_geometry_type, 0);
-  rb_define_method(geos_line_methods, "hash", method_line_hash, 0);
-}
-
-
-VALUE rgeo_is_geos_line_string_closed(GEOSContextHandle_t context, const GEOSGeometry* geom)
-{
-  VALUE result;
-  unsigned int n;
-  double x1, x2, y1, y2, z1, z2;
-  const GEOSCoordSequence* coord_seq;
-
-  result = Qnil;
-  n = GEOSGetNumCoordinates_r(context, geom);
-  if (n > 0) {
-    coord_seq = GEOSGeom_getCoordSeq_r(context, geom);
-    if (GEOSCoordSeq_getX_r(context, coord_seq, 0, &x1)) {
-      if (GEOSCoordSeq_getX_r(context, coord_seq, n-1, &x2)) {
-        if (x1 == x2) {
-          if (GEOSCoordSeq_getY_r(context, coord_seq, 0, &y1)) {
-            if (GEOSCoordSeq_getY_r(context, coord_seq, n-1, &y2)) {
-              result = y1 == y2 ? Qtrue : Qfalse;
-            }
-          }
-        }
-        else {
-          result = Qfalse;
-        }
-      }
-    }
-  }
-  return result;
-}
-
-
-RGEO_END_C
-
-#endif