ffi-geos 2.3.1 → 2.5.0

data/test/geometry_tests.rb

@@ -1,2096 +0,0 @@
-# frozen_string_literal: true
-
-require 'test_helper'
-
-class GeometryTests < Minitest::Test
-  include TestHelper
-
-  def setup
-    super
-    writer.trim = true
-  end
-
-  def test_intersection
-    comparison_tester(
-      :intersection,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((10 10, 10 5, 5 5, 5 10, 10 10))'
-      else
-        'POLYGON ((5 10, 10 10, 10 5, 5 5, 5 10))'
-      end,
-      'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'POLYGON ((5 5, 15 5, 15 15, 5 15, 5 5))'
-    )
-  end
-
-  def test_intersection_with_precision
-    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSIntersectionPrec_r)
-
-    comparison_tester(
-      :intersection,
-      'GEOMETRYCOLLECTION (POLYGON ((1 2, 1 1, 0.5 1, 1 2)), POLYGON ((9.5 1, 2 1, 2 2, 9 2, 9.5 1)), LINESTRING (1 1, 2 1), LINESTRING (2 2, 1 2))',
-      'MULTIPOLYGON(((0 0,5 10,10 0,0 0),(1 1,1 2,2 2,2 1,1 1),(100 100,100 102,102 102,102 100,100 100)))',
-      'POLYGON((0 1,0 2,10 2,10 1,0 1))',
-      precision: 0
-    )
-
-    comparison_tester(
-      :intersection,
-      if Geos::GEOS_NICE_VERSION >= '031000'
-        'GEOMETRYCOLLECTION (LINESTRING (2 0, 4 0), POINT (10 0), POINT (0 0))'
-      else
-        'GEOMETRYCOLLECTION (LINESTRING (2 0, 4 0), POINT (0 0), POINT (10 0))'
-      end,
-      'LINESTRING(0 0, 10 0)',
-      'LINESTRING(9 0, 12 0, 12 20, 4 0, 2 0, 2 10, 0 10, 0 -10)',
-      precision: 2
-    )
-  end
-
-  def test_buffer
-    simple_tester(
-      :buffer,
-      'POLYGON EMPTY',
-      'POINT(0 0)',
-      0
-    )
-
-    snapped_tester(
-      :buffer,
-      'POLYGON ((10 0, 10 -2, 9 -4, 8 -6, 7 -7, 6 -8, 4 -9, 2 -10, 0 -10, -2 -10, -4 -9, -6 -8, -7 -7, -8 -6, -9 -4, -10 -2, -10 0, -10 2, -9 4, -8 6, -7 7, -6 8, -4 9, -2 10, 0 10, 2 10, 4 9, 6 8, 7 7, 8 6, 9 4, 10 2, 10 0))',
-      'POINT(0 0)',
-      10
-    )
-
-    # One segment per quadrant
-    snapped_tester(
-      :buffer,
-      'POLYGON ((10 0, 0 -10, -10 0, 0 10, 10 0))',
-      'POINT(0 0)',
-      10,
-      quad_segs: 1
-    )
-
-    # End cap styles
-    snapped_tester(
-      :buffer,
-      'POLYGON ((100 10, 110 0, 100 -10, 0 -10, -10 0, 0 10, 100 10))',
-      'LINESTRING(0 0, 100 0)',
-      10,
-      quad_segs: 1, endcap: :round
-    )
-
-    snapped_tester(
-      :buffer,
-      'POLYGON ((100 10, 100 -10, 0 -10, 0 10, 100 10))',
-      'LINESTRING(0 0, 100 0)',
-      10,
-      quad_segs: 1, endcap: :flat
-    )
-
-    snapped_tester(
-      :buffer,
-      'POLYGON ((100 10, 110 10, 110 -10, 0 -10, -10 -10, -10 10, 100 10))',
-      'LINESTRING(0 0, 100 0)',
-      10,
-      quad_segs: 1, endcap: :square
-    )
-
-    # Join styles
-    snapped_tester(
-      :buffer,
-      'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 110 0, 107 -7, 100 -10, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))',
-      'LINESTRING(0 0, 100 0, 100 100)',
-      10,
-      quad_segs: 2, join: :round
-    )
-
-    snapped_tester(
-      :buffer,
-      'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 110 0, 100 -10, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))',
-      'LINESTRING(0 0, 100 0, 100 100)',
-      10,
-      quad_segs: 2, join: :bevel
-    )
-
-    snapped_tester(
-      :buffer,
-      'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 110 -10, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))',
-      'LINESTRING(0 0, 100 0, 100 100)',
-      10,
-      quad_segs: 2, join: :mitre
-    )
-
-    snapped_tester(
-      :buffer,
-      if Geos::GEOS_NICE_VERSION >= '031100'
-        'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 110 -4, 104 -10, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))'
-      else
-        'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 109 -5, 105 -9, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))'
-      end,
-      'LINESTRING(0 0, 100 0, 100 100)',
-      10,
-      quad_segs: 2, join: :mitre, mitre_limit: 1.0
-    )
-
-    # Single-sided buffering
-    snapped_tester(
-      :buffer,
-      'POLYGON ((100 0, 0 0, 0 10, 100 10, 100 0))',
-      'LINESTRING(0 0, 100 0)',
-      10,
-      single_sided: true
-    )
-
-    snapped_tester(
-      :buffer,
-      'POLYGON ((0 0, 100 0, 100 -10, 0 -10, 0 0))',
-      'LINESTRING(0 0, 100 0)',
-      -10,
-      single_sided: true
-    )
-  end
-
-  def test_convex_hull
-    geom = read('POINT(0 0)')
-    assert_geom_eql_exact(read('POINT(0 0)'), geom.convex_hull)
-
-    geom = read('LINESTRING(0 0, 10 10)')
-    assert_geom_eql_exact(read('LINESTRING(0 0, 10 10)'), geom.convex_hull)
-
-    geom = read('POLYGON((0 0, 0 10, 5 5, 10 10, 10 0, 0 0))')
-    assert_geom_eql_exact(read('POLYGON((0 0, 0 10, 10 10, 10 0, 0 0))'), geom.convex_hull)
-  end
-
-  def test_difference
-    comparison_tester(
-      :difference,
-      EMPTY_GEOMETRY,
-      'POINT(0 0)',
-      'POINT(0 0)'
-    )
-
-    comparison_tester(
-      :difference,
-      'POINT (0 0)',
-      'POINT(0 0)',
-      'POINT(1 0)'
-    )
-
-    comparison_tester(
-      :difference,
-      'LINESTRING (0 0, 10 0)',
-      'LINESTRING(0 0, 10 0)',
-      'POINT(5 0)'
-    )
-
-    comparison_tester(
-      :difference,
-      EMPTY_GEOMETRY,
-      'POINT(5 0)',
-      'LINESTRING(0 0, 10 0)'
-    )
-
-    comparison_tester(
-      :difference,
-      'POINT (5 0)',
-      'POINT(5 0)',
-      'LINESTRING(0 1, 10 1)'
-    )
-
-    comparison_tester(
-      :difference,
-      'MULTILINESTRING ((0 0, 5 0), (5 0, 10 0))',
-      'LINESTRING(0 0, 10 0)',
-      'LINESTRING(5 -10, 5 10)'
-    )
-
-    comparison_tester(
-      :difference,
-      'LINESTRING (0 0, 5 0)',
-      'LINESTRING(0 0, 10 0)',
-      'LINESTRING(5 0, 20 0)'
-    )
-
-    comparison_tester(
-      :difference,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((0 10, 5 10, 10 10, 10 0, 5 0, 0 0, 0 10))'
-      else
-        'POLYGON ((0 0, 0 10, 5 10, 10 10, 10 0, 5 0, 0 0))'
-      end,
-      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'LINESTRING(5 -10, 5 10)'
-    )
-
-    comparison_tester(
-      :difference,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((0 10, 10 10, 10 0, 0 0, 0 10))'
-      else
-        'POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0))'
-      end,
-      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'LINESTRING(10 0, 20 0)'
-    )
-
-    comparison_tester(
-      :difference,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((0 10, 10 10, 10 5, 5 5, 5 0, 0 0, 0 10))'
-      else
-        'POLYGON ((0 0, 0 10, 10 10, 10 5, 5 5, 5 0, 0 0))'
-      end,
-      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'POLYGON((5 -5, 5 5, 15 5, 15 -5, 5 -5))'
-    )
-  end
-
-  def test_difference_with_precision
-    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSDifferencePrec_r)
-
-    comparison_tester(
-      :difference,
-      'MULTILINESTRING ((2 8, 4 8), (6 8, 10 8))',
-      'LINESTRING (2 8, 10 8)',
-      'LINESTRING (3.9 8.1, 6.1 7.9)',
-      precision: 2
-    )
-  end
-
-  def test_sym_difference
-    %w{ sym_difference symmetric_difference }.each do |method|
-      comparison_tester(
-        method,
-        EMPTY_GEOMETRY,
-        'POINT(0 0)',
-        'POINT(0 0)'
-      )
-
-      comparison_tester(
-        method,
-        'MULTIPOINT (0 0, 1 0)',
-        'POINT(0 0)',
-        'POINT(1 0)'
-      )
-
-      comparison_tester(
-        method,
-        'LINESTRING (0 0, 10 0)',
-        'LINESTRING(0 0, 10 0)',
-        'POINT(5 0)'
-      )
-
-      comparison_tester(
-        method,
-        'LINESTRING (0 0, 10 0)',
-        'POINT(5 0)',
-        'LINESTRING(0 0, 10 0)'
-      )
-
-      comparison_tester(
-        method,
-        'GEOMETRYCOLLECTION (POINT (5 0), LINESTRING (0 1, 10 1))',
-        'POINT(5 0)',
-        'LINESTRING(0 1, 10 1)'
-      )
-
-      comparison_tester(
-        method,
-        'MULTILINESTRING ((0 0, 5 0), (5 0, 10 0), (5 -10, 5 0), (5 0, 5 10))',
-        'LINESTRING(0 0, 10 0)',
-        'LINESTRING(5 -10, 5 10)'
-      )
-
-      comparison_tester(
-        method,
-        'MULTILINESTRING ((0 0, 5 0), (10 0, 20 0))',
-        'LINESTRING(0 0, 10 0)',
-        'LINESTRING(5 0, 20 0)'
-      )
-
-      comparison_tester(
-        method,
-        if Geos::GEOS_NICE_VERSION > '030900'
-          'GEOMETRYCOLLECTION (POLYGON ((0 10, 5 10, 10 10, 10 0, 5 0, 0 0, 0 10)), LINESTRING (5 -10, 5 0))'
-        else
-          'GEOMETRYCOLLECTION (LINESTRING (5 -10, 5 0), POLYGON ((0 0, 0 10, 5 10, 10 10, 10 0, 5 0, 0 0)))'
-        end,
-        'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-        'LINESTRING(5 -10, 5 10)'
-      )
-
-      comparison_tester(
-        method,
-        if Geos::GEOS_NICE_VERSION > '030900'
-          'GEOMETRYCOLLECTION (POLYGON ((0 10, 10 10, 10 0, 0 0, 0 10)), LINESTRING (10 0, 20 0))'
-        else
-          'GEOMETRYCOLLECTION (LINESTRING (10 0, 20 0), POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0)))'
-        end,
-        'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-        'LINESTRING(10 0, 20 0)'
-      )
-
-      comparison_tester(
-        method,
-        if Geos::GEOS_NICE_VERSION > '030900'
-          'MULTIPOLYGON (((0 10, 10 10, 10 5, 5 5, 5 0, 0 0, 0 10)), ((10 0, 10 5, 15 5, 15 -5, 5 -5, 5 0, 10 0)))'
-        else
-          'MULTIPOLYGON (((0 0, 0 10, 10 10, 10 5, 5 5, 5 0, 0 0)), ((5 0, 10 0, 10 5, 15 5, 15 -5, 5 -5, 5 0)))'
-        end,
-        'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-        'POLYGON((5 -5, 5 5, 15 5, 15 -5, 5 -5))'
-      )
-    end
-  end
-
-  def test_sym_difference_with_precision
-    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSSymDifferencePrec_r)
-
-    comparison_tester(
-      :sym_difference,
-      'GEOMETRYCOLLECTION (POLYGON ((0 10, 6 10, 10 10, 10 0, 6 0, 0 0, 0 10)), LINESTRING (6 -10, 6 0))',
-      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'LINESTRING(5 -10, 5 10)',
-      precision: 2
-    )
-  end
-
-  def test_boundary
-    simple_tester(
-      :boundary,
-      'GEOMETRYCOLLECTION EMPTY',
-      'POINT(0 0)'
-    )
-
-    simple_tester(
-      :boundary,
-      'MULTIPOINT (0 0, 10 10)',
-      'LINESTRING(0 0, 10 10)'
-    )
-
-    simple_tester(
-      :boundary,
-      'MULTILINESTRING ((0 0, 10 0, 10 10, 0 10, 0 0), (5 5, 5 6, 6 6, 6 5, 5 5))',
-      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0),( 5 5, 5 6, 6 6, 6 5, 5 5))'
-    )
-  end
-
-  def test_union
-    comparison_tester(
-      :union,
-      'POINT (0 0)',
-      'POINT(0 0)',
-      'POINT(0 0)'
-    )
-
-    comparison_tester(
-      :union,
-      'MULTIPOINT (0 0, 1 0)',
-      'POINT(0 0)',
-      'POINT(1 0)'
-    )
-
-    comparison_tester(
-      :union,
-      'LINESTRING (0 0, 10 0)',
-      'LINESTRING(0 0, 10 0)',
-      'POINT(5 0)'
-    )
-
-    comparison_tester(
-      :union,
-      'LINESTRING (0 0, 10 0)',
-      'POINT(5 0)',
-      'LINESTRING(0 0, 10 0)'
-    )
-
-    comparison_tester(
-      :union,
-      'GEOMETRYCOLLECTION (POINT (5 0), LINESTRING (0 1, 10 1))',
-      'POINT(5 0)',
-      'LINESTRING(0 1, 10 1)'
-    )
-
-    comparison_tester(
-      :union,
-      'MULTILINESTRING ((0 0, 5 0), (5 0, 10 0), (5 -10, 5 0), (5 0, 5 10))',
-      'LINESTRING(0 0, 10 0)',
-      'LINESTRING(5 -10, 5 10)'
-    )
-
-    comparison_tester(
-      :union,
-      'MULTILINESTRING ((0 0, 5 0), (5 0, 10 0), (10 0, 20 0))',
-      'LINESTRING(0 0, 10 0)',
-      'LINESTRING(5 0, 20 0)'
-    )
-
-    comparison_tester(
-      :union,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'GEOMETRYCOLLECTION (POLYGON ((0 10, 5 10, 10 10, 10 0, 5 0, 0 0, 0 10)), LINESTRING (5 -10, 5 0))'
-      else
-        'GEOMETRYCOLLECTION (LINESTRING (5 -10, 5 0), POLYGON ((0 0, 0 10, 5 10, 10 10, 10 0, 5 0, 0 0)))'
-      end,
-      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'LINESTRING(5 -10, 5 10)'
-    )
-
-    comparison_tester(
-      :union,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'GEOMETRYCOLLECTION (POLYGON ((0 10, 10 10, 10 0, 0 0, 0 10)), LINESTRING (10 0, 20 0))'
-      else
-        'GEOMETRYCOLLECTION (LINESTRING (10 0, 20 0), POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0)))'
-      end,
-      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'LINESTRING(10 0, 20 0)'
-    )
-
-    comparison_tester(
-      :union,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((0 10, 10 10, 10 5, 15 5, 15 -5, 5 -5, 5 0, 0 0, 0 10))'
-      else
-        'POLYGON ((0 0, 0 10, 10 10, 10 5, 15 5, 15 -5, 5 -5, 5 0, 0 0))'
-      end,
-      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'POLYGON((5 -5, 5 5, 15 5, 15 -5, 5 -5))'
-    )
-  end
-
-  def test_union_with_precision
-    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSUnionPrec_r)
-
-    geom_a = read('POINT (1.9 8.2)')
-    geom_b = read('POINT (4.1 9.8)')
-
-    result = geom_a.union(geom_b, precision: 2)
-
-    assert_equal('MULTIPOINT (2 8, 4 10)', write(result))
-  end
-
-  def test_union_cascaded
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:union_cascaded)
-
-    simple_tester(
-      :union_cascaded,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0, 0 0), (12 12, 11 12, 11 11, 12 11, 12 12))'
-      else
-        'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))'
-      end,
-      'MULTIPOLYGON(
-        ((0 0, 1 0, 1 1, 0 1, 0 0)),
-        ((10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11)),
-        ((0 0, 11 0, 11 11, 0 11, 0 0))
-      ))'
-    )
-  end
-
-  def test_coverage_union
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:coverage_union)
-
-    simple_tester(
-      :union_cascaded,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((0 1, 1 1, 2 1, 2 0, 1 0, 0 0, 0 1))'
-      else
-        'POLYGON ((0 0, 0 1, 1 1, 2 1, 2 0, 1 0, 0 0))'
-      end,
-      'MULTIPOLYGON(
-        ((0 0, 0 1, 1 1, 1 0, 0 0)),
-        ((1 0, 1 1, 2 1, 2 0, 1 0))
-      ))'
-    )
-  end
-
-  def test_unary_union
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:unary_union)
-
-    simple_tester(
-      :unary_union,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0, 0 0), (12 12, 11 12, 11 11, 12 11, 12 12))'
-      else
-        'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))'
-      end,
-      'MULTIPOLYGON(
-        ((0 0, 1 0, 1 1, 0 1, 0 0)),
-        ((10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11)),
-        ((0 0, 11 0, 11 11, 0 11, 0 0))
-      ))'
-    )
-  end
-
-  def test_unary_union_with_precision
-    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSUnaryUnionPrec_r)
-
-    simple_tester(
-      :unary_union,
-      'POLYGON ((0 0, 0 12, 9 12, 9 15, 15 15, 15 9, 12 9, 12 0, 0 0))',
-      'MULTIPOLYGON(
-        ((0 0, 1 0, 1 1, 0 1, 0 0)),
-        ((10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11)),
-        ((0 0, 11 0, 11 11, 0 11, 0 0))
-      ))',
-      3
-    )
-  end
-
-  def test_node
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:node)
-
-    simple_tester(
-      :node,
-      'MULTILINESTRING ((0 0, 5 0), (5 0, 10 0, 5 -5, 5 0), (5 0, 5 5))',
-      'LINESTRING(0 0, 10 0, 5 -5, 5 5)'
-    )
-  end
-
-  def test_union_without_arguments
-    simple_tester(
-      :union,
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'POLYGON ((0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0, 0 0), (12 12, 11 12, 11 11, 12 11, 12 12))'
-      else
-        'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))'
-      end,
-      'MULTIPOLYGON(
-        ((0 0, 1 0, 1 1, 0 1, 0 0)),
-        ((10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11)),
-        ((0 0, 11 0, 11 11, 0 11, 0 0))
-      ))'
-    )
-  end
-
-  def test_point_on_surface_and_representative_point
-    %w{
-      point_on_surface
-      representative_point
-    }.each do |method|
-      simple_tester(
-        method,
-        'POINT (0 0)',
-        'POINT (0 0)'
-      )
-
-      simple_tester(
-        method,
-        'POINT (5 0)',
-        'LINESTRING(0 0, 5 0, 10 0)'
-      )
-
-      simple_tester(
-        method,
-        'POINT (5 5)',
-        'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))'
-      )
-    end
-  end
-
-  def test_clip_by_rect
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:clip_by_rect)
-
-    %w{
-      clip_by_rect
-      clip_by_rectangle
-    }.each do |method|
-      simple_tester(
-        method,
-        'POINT (0 0)',
-        'POINT (0 0)',
-        -1, -1, 1, 1
-      )
-
-      simple_tester(
-        method,
-        'GEOMETRYCOLLECTION EMPTY',
-        'POINT (0 0)',
-        0, 0, 2, 2
-      )
-
-      simple_tester(
-        method,
-        'LINESTRING (1 0, 2 0)',
-        'LINESTRING (0 0, 10 0)',
-        1, -1, 2, 1
-      )
-
-      simple_tester(
-        method,
-        'POLYGON ((1 1, 1 5, 5 5, 5 1, 1 1))',
-        'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))',
-        1, 1, 5, 5
-      )
-
-      simple_tester(
-        method,
-        'POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))',
-        'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))',
-        -1, -1, 5, 5
-      )
-    end
-  end
-
-  def test_centroid_and_center
-    %w{
-      centroid
-      center
-    }.each do |method|
-      simple_tester(
-        method,
-        'POINT (0 0)',
-        'POINT(0 0)'
-      )
-
-      simple_tester(
-        method,
-        'POINT (5 5)',
-        'LINESTRING(0 0, 10 10)'
-      )
-
-      snapped_tester(
-        method,
-        'POINT (5 4)',
-        'POLYGON((0 0, 0 10, 5 5, 10 10, 10 0, 0 0))'
-      )
-    end
-  end
-
-  def test_minimum_bounding_circle
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:minimum_bounding_circle)
-
-    geom = read('LINESTRING(0 10, 0 20)')
-
-    assert_equal(
-      'POLYGON ((5 15, 5 14, 5 13, 4 12, 4 11, 3 11, 2 10, 1 10, 0 10, -1 10, -2 10, -3 11, -4 11, -4 12, -5 13, -5 14, -5 15, -5 16, -5 17, -4 18, -4 19, -3 19, -2 20, -1 20, 0 20, 1 20, 2 20, 3 19, 4 19, 4 18, 5 17, 5 16, 5 15))',
-      write(geom.minimum_bounding_circle.snap_to_grid(1))
-    )
-  end
-
-  def test_envelope
-    simple_tester(
-      :envelope,
-      'POINT (0 0)',
-      'POINT(0 0)'
-    )
-
-    simple_tester(
-      :envelope,
-      'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))',
-      'LINESTRING(0 0, 10 10)'
-    )
-  end
-
-  def test_relate
-    tester = lambda { |expected, geom_a, geom_b|
-      assert_equal(expected, geom_a.relate(geom_b))
-    }
-
-    geom_a = read('POINT(0 0)')
-    geom_b = read('POINT(0 0)')
-    tester['0FFFFFFF2', geom_a, geom_b]
-
-    geom_a = read('POINT(0 0)')
-    geom_b = read('POINT(1 0)')
-    tester['FF0FFF0F2', geom_a, geom_b]
-
-    geom_a = read('POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))')
-    geom_b = read('POINT(1 0)')
-    tester['FF20F1FF2', geom_a, geom_b]
-  end
-
-  def test_relate_pattern
-    tester = lambda { |pattern, geom_a, geom_b, expected|
-      assert_equal(expected, geom_a.relate_pattern(geom_b, pattern))
-    }
-
-    geom_a = read('POINT(0 0)')
-    geom_b = read('POINT(0 0)')
-    tester['0FFFFFFF2', geom_a, geom_b, true]
-    tester['0*******T', geom_a, geom_b, true]
-    tester['0*******1', geom_a, geom_b, false]
-
-    geom_a = read('POINT(0 0)')
-    geom_b = read('POINT(1 0)')
-    tester['FF0FFF0F2', geom_a, geom_b, true]
-    tester['F*******2', geom_a, geom_b, true]
-    tester['T*******2', geom_a, geom_b, false]
-
-    geom_a = read('POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))')
-    geom_b = read('POINT(1 0)')
-    tester['FF20F1FF2', geom_a, geom_b, true]
-    tester['F****T**T', geom_a, geom_b, true]
-    tester['T*******2', geom_a, geom_b, false]
-  end
-
-  def test_relate_boundary_node_rule
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:relate_boundary_node_rule)
-
-    geom_a = read('LINESTRING(0 0, 2 4, 5 5, 0 0)')
-    geom_b = read('POINT(0 0)')
-
-    ret = geom_a.relate_boundary_node_rule(geom_b, :ogc)
-    assert_equal('0F1FFFFF2', ret)
-
-    ret = geom_a.relate_boundary_node_rule(geom_b, :endpoint)
-    assert_equal('FF10FFFF2', ret)
-
-    assert_raises(TypeError) do
-      geom_a.relate_boundary_node_rule(geom_b, :gibberish)
-    end
-  end
-
-  def test_line_merge
-    simple_tester(
-      :line_merge,
-      'LINESTRING (0 0, 10 10, 10 0, 5 0, 5 -5)',
-      'MULTILINESTRING(
-        (0 0, 10 10),
-        (10 10, 10 0),
-        (5 0, 10 0),
-        (5 -5, 5 0)
-      )'
-    )
-  end
-
-  def test_simplify
-    simple_tester(
-      :simplify,
-      'LINESTRING (0 0, 5 10, 10 0, 10 9, 0 9)',
-      'LINESTRING(0 0, 3 4, 5 10, 10 0, 10 9, 5 11, 0 9)',
-      2
-    )
-  end
-
-  def test_topology_preserve_simplify
-    simple_tester(
-      :topology_preserve_simplify,
-      'LINESTRING (0 0, 5 10, 10 0, 10 9, 5 11, 0 9)',
-      'LINESTRING(0 0, 3 4, 5 10, 10 0, 10 9, 5 11, 0 9)',
-      2
-    )
-  end
-
-  def test_extract_unique_points
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:extract_unique_points)
-
-    geom = read('GEOMETRYCOLLECTION (
-      MULTIPOLYGON (
-        ((0 0, 1 0, 1 1, 0 1, 0 0)),
-        ((10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11))
-      ),
-      POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0)),
-      MULTILINESTRING ((0 0, 2 3), (10 10, 3 4)),
-      LINESTRING (0 0, 2 3),
-      MULTIPOINT (0 0, 2 3),
-      POINT (9 0),
-      POINT(1 0)),
-      LINESTRING EMPTY
-    ')
-
-    simple_tester(
-      :extract_unique_points,
-      'MULTIPOINT (0 0, 1 0, 1 1, 0 1, 10 10, 10 14, 14 14, 14 10, 11 11, 11 12, 12 12, 12 11, 2 3, 3 4, 9 0)',
-      geom.extract_unique_points
-    )
-  end
-
-  def test_relationships
-    tester = lambda { |geom_a, geom_b, tests|
-      tests.each do |test|
-        expected, method, args = test
-        if ENV['FORCE_TESTS'] || geom_a.respond_to?(method)
-          value = geom_a.send(method, *([geom_b] + Array(args)))
-          assert_equal(expected, value)
-        end
-      end
-    }
-
-    tester[read('POINT(0 0)'), read('POINT(0 0)'), [
-      [false, :disjoint?],
-      [false, :touches?],
-      [true, :intersects?],
-      [false, :crosses?],
-      [true, :within?],
-      [true, :contains?],
-      [false, :overlaps?],
-      [true, :eql?],
-      [true, :eql_exact?, TOLERANCE],
-      [true, :covers?],
-      [true, :covered_by?]
-    ]]
-
-    tester[read('POINT(0 0)'), read('LINESTRING(0 0, 10 0)'), [
-      [false, :disjoint?],
-      [true, :touches?],
-      [true, :intersects?],
-      [false, :crosses?],
-      [false, :within?],
-      [false, :contains?],
-      [false, :overlaps?],
-      [false, :eql?],
-      [false, :eql_exact?, TOLERANCE],
-      [false, :covers?],
-      [true, :covered_by?]
-    ]]
-
-    tester[read('POINT(5 0)'), read('LINESTRING(0 0, 10 0)'), [
-      [false, :disjoint?],
-      [false, :touches?],
-      [true, :intersects?],
-      [false, :crosses?],
-      [true, :within?],
-      [false, :contains?],
-      [false, :overlaps?],
-      [false, :eql?],
-      [false, :eql_exact?, TOLERANCE],
-      [false, :covers?],
-      [true, :covered_by?]
-    ]]
-
-    tester[read('LINESTRING(5 -5, 5 5)'), read('LINESTRING(0 0, 10 0)'), [
-      [false, :disjoint?],
-      [false, :touches?],
-      [true, :intersects?],
-      [true, :crosses?],
-      [false, :within?],
-      [false, :contains?],
-      [false, :overlaps?],
-      [false, :eql?],
-      [false, :eql_exact?, TOLERANCE],
-      [false, :covers?],
-      [false, :covered_by?]
-    ]]
-
-    tester[read('LINESTRING(5 0, 15 0)'), read('LINESTRING(0 0, 10 0)'), [
-      [false, :disjoint?],
-      [false, :touches?],
-      [true, :intersects?],
-      [false, :crosses?],
-      [false, :within?],
-      [false, :contains?],
-      [true, :overlaps?],
-      [false, :eql?],
-      [false, :eql_exact?, TOLERANCE],
-      [false, :covers?],
-      [false, :covered_by?]
-    ]]
-
-    tester[read('LINESTRING(0 0, 5 0, 10 0)'), read('LINESTRING(0 0, 10 0)'), [
-      [false, :disjoint?],
-      [false, :touches?],
-      [true, :intersects?],
-      [false, :crosses?],
-      [true, :within?],
-      [true, :contains?],
-      [false, :overlaps?],
-      [true, :eql?],
-      [false, :eql_exact?, TOLERANCE],
-      [true, :covers?],
-      [true, :covered_by?]
-    ]]
-
-    tester[read('POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))'), read('POLYGON((5 -5, 5 5, 15 5, 15 -5, 5 -5))'), [
-      [false, :disjoint?],
-      [false, :touches?],
-      [true, :intersects?],
-      [false, :crosses?],
-      [false, :within?],
-      [false, :contains?],
-      [true, :overlaps?],
-      [false, :eql?],
-      [false, :eql_exact?, TOLERANCE],
-      [false, :covers?],
-      [false, :covered_by?]
-    ]]
-
-    tester[read('POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))'), read('POINT(15 15)'), [
-      [true, :disjoint?],
-      [false, :touches?],
-      [false, :intersects?],
-      [false, :crosses?],
-      [false, :within?],
-      [false, :contains?],
-      [false, :overlaps?],
-      [false, :eql?],
-      [false, :eql_exact?, TOLERANCE],
-      [false, :covers?],
-      [false, :covered_by?]
-    ]]
-  end
-
-  def test_empty
-    refute_geom_empty(read('POINT(0 0)'))
-    assert_geom_empty(read('POINT EMPTY'))
-    refute_geom_empty(read('LINESTRING(0 0, 10 0)'))
-    assert_geom_empty(read('LINESTRING EMPTY'))
-    refute_geom_empty(read('POLYGON((0 0, 10 0, 10 10, 0 0))'))
-    assert_geom_empty(read('POLYGON EMPTY'))
-    refute_geom_empty(read('GEOMETRYCOLLECTION(POINT(0 0))'))
-    assert_geom_empty(read('GEOMETRYCOLLECTION EMPTY'))
-  end
-
-  def test_valid
-    assert_geom_valid(read('POINT(0 0)'))
-    refute_geom_valid(read('POINT(0 NaN)'))
-    refute_geom_valid(read('POINT(0 nan)'))
-  end
-
-  def test_valid_reason
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:valid_reason)
-
-    assert_equal('Valid Geometry', read('POINT(0 0)').valid_reason)
-    assert_equal('Invalid Coordinate[0 nan]', read('POINT(0 NaN)').valid_reason)
-    assert_equal('Invalid Coordinate[0 nan]', read('POINT(0 nan)').valid_reason)
-    assert_equal('Self-intersection[2.5 5]', read('POLYGON((0 0, 0 5, 5 5, 5 10, 0 0))').valid_reason)
-  end
-
-  def test_valid_detail
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:valid_detail)
-
-    tester = lambda { |detail, location, geom, flags|
-      ret = read(geom).valid_detail(flags)
-      assert_equal(detail, ret[:detail])
-      assert_equal(location, write(ret[:location]))
-    }
-
-    assert_nil(read('POINT(0 0)').valid_detail)
-
-    if Geos::GEOS_NICE_VERSION >= '031000'
-      tester['Invalid Coordinate', 'POINT (0 NaN)', 'POINT(0 NaN)', 0]
-    else
-      tester['Invalid Coordinate', 'POINT (0 nan)', 'POINT(0 NaN)', 0]
-    end
-
-    tester['Self-intersection', 'POINT (2.5 5)', 'POLYGON((0 0, 0 5, 5 5, 5 10, 0 0))', 0]
-
-    tester['Ring Self-intersection', 'POINT (0 0)', 'POLYGON((0 0, -10 10, 10 10, 0 0, 4 5, -4 5, 0 0)))', 0]
-
-    assert_nil(
-      read('POLYGON((0 0, -10 10, 10 10, 0 0, 4 5, -4 5, 0 0)))').valid_detail(
-        :allow_selftouching_ring_forming_hole
-      )
-    )
-  end
-
-  def test_simple
-    assert_geom_simple(read('POINT(0 0)'))
-    assert_geom_simple(read('LINESTRING(0 0, 10 0)'))
-    refute_geom_simple(read('LINESTRING(0 0, 10 0, 5 5, 5 -5)'))
-  end
-
-  def test_ring
-    refute_geom_ring(read('POINT(0 0)'))
-    refute_geom_ring(read('LINESTRING(0 0, 10 0, 5 5, 5 -5)'))
-    assert_geom_ring(read('LINESTRING(0 0, 10 0, 5 5, 0 0)'))
-  end
-
-  def test_has_z
-    refute_geom_has_z(read('POINT(0 0)'))
-    assert_geom_has_z(read('POINT(0 0 0)'))
-  end
-
-  def test_num_geometries
-    simple_tester(:num_geometries, 1, 'POINT(0 0)')
-    simple_tester(:num_geometries, 2, 'MULTIPOINT (0 1, 2 3)')
-    simple_tester(:num_geometries, 1, 'LINESTRING (0 0, 2 3)')
-    simple_tester(:num_geometries, 2, 'MULTILINESTRING ((0 1, 2 3), (10 10, 3 4))')
-    simple_tester(:num_geometries, 1, 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
-    simple_tester(:num_geometries, 2, 'MULTIPOLYGON(
-      ((0 0, 1 0, 1 1, 0 1, 0 0)),
-      ((10 10, 10 14, 14 14, 14 10, 10 10),
-      (11 11, 11 12, 12 12, 12 11, 11 11)))')
-    simple_tester(:num_geometries, 6, 'GEOMETRYCOLLECTION (
-      MULTIPOLYGON (
-        ((0 0, 1 0, 1 1, 0 1, 0 0)),
-        ((10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11))
-      ),
-      POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0)),
-      MULTILINESTRING ((0 0, 2 3), (10 10, 3 4)),
-      LINESTRING (0 0, 2 3),
-      MULTIPOINT (0 0, 2 3),
-      POINT (9 0))')
-  end
-
-  # get_geometry_n is segfaulting in the binary GEOS build
-  def test_get_geometry_n
-    skip unless defined?(Geos::FFIGeos)
-
-    simple_tester(:get_geometry_n, 'POINT (0 1)', 'MULTIPOINT (0 1, 2 3)', 0)
-    simple_tester(:get_geometry_n, 'POINT (2 3)', 'MULTIPOINT (0 1, 2 3)', 1)
-    simple_tester(:get_geometry_n, nil, 'MULTIPOINT (0 1, 2 3)', 2)
-  end
-
-  def test_num_interior_rings
-    simple_tester(:num_interior_rings, 0, 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
-    simple_tester(:num_interior_rings, 1, 'POLYGON (
-      (10 10, 10 14, 14 14, 14 10, 10 10),
-      (11 11, 11 12, 12 12, 12 11, 11 11)
-    )')
-    simple_tester(:num_interior_rings, 2, 'POLYGON (
-      (10 10, 10 14, 14 14, 14 10, 10 10),
-      (11 11, 11 12, 12 12, 12 11, 11 11),
-      (13 11, 13 12, 13.5 12, 13.5 11, 13 11))')
-
-    assert_raises(NoMethodError) do
-      read('POINT (0 0)').num_interior_rings
-    end
-  end
-
-  def test_interior_ring_n
-    simple_tester(
-      :interior_ring_n,
-      'LINEARRING (11 11, 11 12, 12 12, 12 11, 11 11)',
-      'POLYGON(
-        (10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11)
-      )',
-      0
-    )
-
-    simple_tester(
-      :interior_ring_n,
-      'LINEARRING (11 11, 11 12, 12 12, 12 11, 11 11)',
-      'POLYGON (
-        (10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11),
-        (13 11, 13 12, 13.5 12, 13.5 11, 13 11)
-      )',
-      0
-    )
-
-    simple_tester(
-      :interior_ring_n,
-      'LINEARRING (13 11, 13 12, 13.5 12, 13.5 11, 13 11)',
-      'POLYGON (
-        (10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11),
-        (13 11, 13 12, 13.5 12, 13.5 11, 13 11)
-      )',
-      1
-    )
-
-    assert_raises(Geos::IndexBoundsError) do
-      simple_tester(
-        :interior_ring_n,
-        nil,
-        'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))',
-        0
-      )
-    end
-
-    assert_raises(NoMethodError) do
-      simple_tester(
-        :interior_ring_n,
-        nil,
-        'POINT (0 0)',
-        0
-      )
-    end
-  end
-
-  def test_exterior_ring
-    simple_tester(
-      :exterior_ring,
-      'LINEARRING (10 10, 10 14, 14 14, 14 10, 10 10)',
-      'POLYGON (
-        (10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11)
-      )'
-    )
-
-    assert_raises(NoMethodError) do
-      read('POINT (0 0)').exterior_ring
-    end
-  end
-
-  def test_interior_rings
-    array_tester(
-      :interior_rings,
-      ['LINEARRING (11 11, 11 12, 12 12, 12 11, 11 11)'],
-      'POLYGON(
-        (10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11)
-      )'
-    )
-
-    array_tester(
-      :interior_rings,
-      [
-        'LINEARRING (11 11, 11 12, 12 12, 12 11, 11 11)',
-        'LINEARRING (13 11, 13 12, 13.5 12, 13.5 11, 13 11)'
-      ],
-      'POLYGON (
-        (10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11),
-        (13 11, 13 12, 13.5 12, 13.5 11, 13 11)
-      )'
-    )
-  end
-
-  def test_num_coordinates
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:num_coordinates)
-
-    simple_tester(:num_coordinates, 1, 'POINT(0 0)')
-    simple_tester(:num_coordinates, 2, 'MULTIPOINT (0 1, 2 3)')
-    simple_tester(:num_coordinates, 2, 'LINESTRING (0 0, 2 3)')
-    simple_tester(:num_coordinates, 4, 'MULTILINESTRING ((0 1, 2 3), (10 10, 3 4))')
-    simple_tester(:num_coordinates, 5, 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
-    simple_tester(:num_coordinates, 15, 'MULTIPOLYGON (
-      ((0 0, 1 0, 1 1, 0 1, 0 0)),
-      ((10 10, 10 14, 14 14, 14 10, 10 10),
-      (11 11, 11 12, 12 12, 12 11, 11 11))
-    )')
-    simple_tester(:num_coordinates, 29, 'GEOMETRYCOLLECTION (
-      MULTIPOLYGON (
-        ((0 0, 1 0, 1 1, 0 1, 0 0)),
-        ((10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11))
-      ),
-      POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0)),
-      MULTILINESTRING ((0 0, 2 3), (10 10, 3 4)),
-      LINESTRING (0 0, 2 3),
-      MULTIPOINT ((0 0), (2 3)),
-      POINT (9 0)
-    )')
-  end
-
-  def test_coord_seq
-    tester = lambda { |expected, g|
-      geom = read(g)
-      cs = geom.coord_seq
-      expected.each_with_index do |c, i|
-        assert_equal(c[0], cs.get_x(i))
-        assert_equal(c[1], cs.get_y(i))
-      end
-    }
-
-    tester[[[0, 0]], 'POINT(0 0)']
-    tester[[[0, 0], [2, 3]], 'LINESTRING (0 0, 2 3)']
-    tester[[[0, 0], [0, 5], [5, 5], [5, 0], [0, 0]], 'LINEARRING(0 0, 0 5, 5 5, 5 0, 0 0)']
-  end
-
-  def test_dimensions
-    types = {
-      dontcare: -3,
-      non_empty: -2,
-      empty: -1,
-      point: 0,
-      curve: 1,
-      surface: 2
-    }
-
-    simple_tester(:dimensions, types[:point], 'POINT(0 0)')
-    simple_tester(:dimensions, types[:point], 'MULTIPOINT (0 1, 2 3)')
-    simple_tester(:dimensions, types[:curve], 'LINESTRING (0 0, 2 3)')
-    simple_tester(:dimensions, types[:curve], 'MULTILINESTRING ((0 1, 2 3), (10 10, 3 4))')
-    simple_tester(:dimensions, types[:surface], 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
-    simple_tester(:dimensions, types[:surface], 'MULTIPOLYGON (
-      ((0 0, 1 0, 1 1, 0 1, 0 0)),
-      ((10 10, 10 14, 14 14, 14 10, 10 10),
-      (11 11, 11 12, 12 12, 12 11, 11 11))
-    )')
-    simple_tester(:dimensions, types[:surface], 'GEOMETRYCOLLECTION (
-      MULTIPOLYGON (
-        ((0 0, 1 0, 1 1, 0 1, 0 0)),
-        ((10 10, 10 14, 14 14, 14 10, 10 10),
-        (11 11, 11 12, 12 12, 12 11, 11 11))
-      ),
-      POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0)),
-      MULTILINESTRING ((0 0, 2 3), (10 10, 3 4)),
-      LINESTRING (0 0, 2 3),
-      MULTIPOINT (0 0, 2 3),
-      POINT (9 0)
-    )')
-  end
-
-  def test_project_and_project_normalized
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:project)
-
-    geom_a = read('POINT(1 2)')
-    geom_b = read('POINT(3 4)')
-
-    # The method only accept lineal geometries
-    assert_raises(Geos::GEOSException) do
-      geom_a.project(geom_b)
-    end
-
-    geom_a = read('LINESTRING(0 0, 10 0)')
-    geom_b = read('POINT(0 0)')
-    assert_equal(0, geom_a.project(geom_b))
-    assert_equal(0, geom_a.project(geom_b, true))
-    assert_equal(0, geom_a.project_normalized(geom_b))
-
-    geom_b = read('POINT(10 0)')
-    assert_equal(10, geom_a.project(geom_b))
-    assert_equal(1, geom_a.project(geom_b, true))
-    assert_equal(1, geom_a.project_normalized(geom_b))
-
-    geom_b = read('POINT(5 0)')
-    assert_equal(5, geom_a.project(geom_b))
-    assert_equal(0.5, geom_a.project(geom_b, true))
-    assert_equal(0.5, geom_a.project_normalized(geom_b))
-
-    geom_a = read('MULTILINESTRING((0 0, 10 0),(20 10, 20 20))')
-    geom_b = read('POINT(20 0)')
-    assert_equal(10, geom_a.project(geom_b))
-    assert_equal(0.5, geom_a.project(geom_b, true))
-    assert_equal(0.5, geom_a.project_normalized(geom_b))
-
-    geom_b = read('POINT(20 5)')
-    assert_equal(10, geom_a.project(geom_b))
-    assert_equal(0.5, geom_a.project(geom_b, true))
-    assert_equal(0.5, geom_a.project_normalized(geom_b))
-  end
-
-  def test_interpolate
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:interpolate)
-
-    simple_tester(:interpolate, 'POINT (0 0)', 'LINESTRING(0 0, 10 0)', 0, false)
-    simple_tester(:interpolate, 'POINT (0 0)', 'LINESTRING(0 0, 10 0)', 0, true)
-
-    simple_tester(:interpolate, 'POINT (5 0)', 'LINESTRING(0 0, 10 0)', 5, false)
-    simple_tester(:interpolate, 'POINT (5 0)', 'LINESTRING(0 0, 10 0)', 0.5, true)
-
-    simple_tester(:interpolate, 'POINT (10 0)', 'LINESTRING(0 0, 10 0)', 20, false)
-    simple_tester(:interpolate, 'POINT (10 0)', 'LINESTRING(0 0, 10 0)', 2, true)
-
-    assert_raises(Geos::GEOSException) do
-      read('POINT(1 2)').interpolate(0)
-    end
-  end
-
-  def test_interpolate_normalized
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:interpolate_normalized)
-
-    tester = lambda { |expected, g, d|
-      geom = read(g)
-      assert_equal(expected, write(geom.interpolate_normalized(d)))
-    }
-
-    writer.trim = true
-
-    tester['POINT (0 0)', 'LINESTRING(0 0, 10 0)', 0]
-    tester['POINT (5 0)', 'LINESTRING(0 0, 10 0)', 0.5]
-    tester['POINT (10 0)', 'LINESTRING(0 0, 10 0)', 2]
-
-    assert_raises(Geos::GEOSException) do
-      read('POINT(1 2)').interpolate_normalized(0)
-    end
-  end
-
-  def test_start_and_end_points
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:start_point)
-
-    geom = read('LINESTRING (10 10, 10 14, 14 14, 14 10)')
-    simple_tester(:start_point, 'POINT (10 10)', geom)
-    simple_tester(:end_point, 'POINT (14 10)', geom)
-
-    geom = read('LINEARRING (11 11, 11 12, 12 11, 11 11)')
-    simple_tester(:start_point, 'POINT (11 11)', geom)
-    simple_tester(:start_point, 'POINT (11 11)', geom)
-  end
-
-  def test_area
-    simple_tester(:area, 1.0, 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
-    simple_tester(:area, 0.0, 'POINT (0 0)')
-    simple_tester(:area, 0.0, 'LINESTRING (0 0 , 10 0)')
-  end
-
-  def test_length
-    simple_tester(:length, 4.0, 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
-    simple_tester(:length, 0.0, 'POINT (0 0)')
-    simple_tester(:length, 10.0, 'LINESTRING (0 0 , 10 0)')
-  end
-
-  def test_distance
-    geom = 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))'
-    simple_tester(:distance, 0.0, geom, read('POINT(0.5 0.5)'))
-    simple_tester(:distance, 1.0, geom, read('POINT (-1 0)'))
-    simple_tester(:distance, 2.0, geom, read('LINESTRING (3 0 , 10 0)'))
-  end
-
-  def test_distance_indexed
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:distance_indexed)
-
-    geom_a = read('POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
-    geom_b = read('POLYGON ((20 30, 10 10, 13 14, 7 8, 20 30))')
-
-    assert_in_delta(9.219544457292887, geom_a.distance_indexed(geom_b), TOLERANCE)
-    assert_in_delta(9.219544457292887, geom_b.distance_indexed(geom_a), TOLERANCE)
-  end
-
-  def test_hausdorff_distance
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hausdorff_distance)
-
-    tester = lambda { |expected, g_1, g_2|
-      geom_1 = read(g_1)
-      geom_2 = read(g_2)
-      assert_in_delta(expected, geom_1.hausdorff_distance(geom_2), TOLERANCE)
-    }
-
-    geom_a = 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))'
-
-    tester[10.0498756211209, geom_a, 'POINT(0 10)']
-    tester[2.23606797749979, geom_a, 'POINT(-1 0)']
-    tester[9.0, geom_a, 'LINESTRING (3 0 , 10 0)']
-  end
-
-  def test_hausdorff_distance_with_densify_fract
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hausdorff_distance)
-
-    tester = lambda { |expected, g_1, g_2|
-      geom_1 = read(g_1)
-      geom_2 = read(g_2)
-      assert_in_delta(expected, geom_1.hausdorff_distance(geom_2, 0.001), TOLERANCE)
-    }
-
-    geom_a = 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))'
-
-    tester[10.0498756211209, geom_a, 'POINT(0 10)']
-    tester[2.23606797749979, geom_a, 'POINT(-1 0)']
-    tester[9.0, geom_a, 'LINESTRING (3 0 , 10 0)']
-  end
-
-  def test_nearest_points
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:nearest_points)
-
-    tester = lambda { |expected, g_1, g_2|
-      geom_1 = read(g_1)
-      geom_2 = read(g_2)
-
-      cs = geom_1.nearest_points(geom_2)
-      result = cs.to_s if cs
-
-      if expected.nil?
-        assert_nil(result)
-      else
-        assert_equal(expected, result)
-      end
-    }
-
-    tester[
-      nil,
-      'POINT EMPTY',
-      'POINT EMPTY'
-    ]
-
-    tester[
-      if Geos::GEOS_NICE_VERSION >= '030800'
-        '5.0 5.0, 8.0 8.0'
-      else
-        '5.0 5.0 NaN, 8.0 8.0 NaN'
-      end,
-      'POLYGON((1 1, 1 5, 5 5, 5 1, 1 1))',
-      'POLYGON((8 8, 9 9, 9 10, 8 8))'
-    ]
-  end
-
-  def test_snap
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:snap)
-
-    geom = read('POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
-    simple_tester(:snap, 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))', geom, read('POINT(0.1 0)'), 0)
-    simple_tester(:snap, 'POLYGON ((0.1 0, 1 0, 1 1, 0 1, 0.1 0))', geom, read('POINT(0.1 0)'), 0.5)
-  end
-
-  def test_polygonize
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:polygonize)
-
-    geom_a = read(
-      'GEOMETRYCOLLECTION(
-        LINESTRING(0 0, 10 10),
-        LINESTRING(185 221, 100 100),
-        LINESTRING(185 221, 88 275, 180 316),
-        LINESTRING(185 221, 292 281, 180 316),
-        LINESTRING(189 98, 83 187, 185 221),
-        LINESTRING(189 98, 325 168, 185 221)
-      )'
-    )
-
-    polygonized = geom_a.polygonize
-    assert_equal(2, polygonized.length)
-    assert_equal(
-      'POLYGON ((185 221, 88 275, 180 316, 292 281, 185 221))',
-      write(polygonized[0].snap_to_grid(0.1))
-    )
-    assert_equal(
-      'POLYGON ((189 98, 83 187, 185 221, 325 168, 189 98))',
-      write(polygonized[1].snap_to_grid(0.1))
-    )
-  end
-
-  def test_polygonize_with_geometry_arguments
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:polygonize)
-
-    geom_a = read('LINESTRING (100 100, 100 300, 300 300, 300 100, 100 100)')
-    geom_b = read('LINESTRING (150 150, 150 250, 250 250, 250 150, 150 150)')
-
-    polygonized = geom_a.polygonize(geom_b)
-    assert_equal(2, polygonized.length)
-    assert_equal(
-      'POLYGON ((100 100, 100 300, 300 300, 300 100, 100 100), (150 150, 250 150, 250 250, 150 250, 150 150))',
-      write(polygonized[0].snap_to_grid(0.1))
-    )
-    assert_equal(
-      'POLYGON ((150 150, 150 250, 250 250, 250 150, 150 150))',
-      write(polygonized[1].snap_to_grid(0.1))
-    )
-  end
-
-  def test_polygonize_valid
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:polygonize_valid)
-
-    geom_a = read(
-      'GEOMETRYCOLLECTION(
-        LINESTRING (100 100, 100 300, 300 300, 300 100, 100 100),
-        LINESTRING (150 150, 150 250, 250 250, 250 150, 150 150)
-      )'
-    )
-
-    polygonized = geom_a.polygonize_valid
-    assert_equal(
-      'POLYGON ((100 100, 100 300, 300 300, 300 100, 100 100), (150 150, 250 150, 250 250, 150 250, 150 150))',
-      write(polygonized.snap_to_grid(0.1))
-    )
-  end
-
-  def test_polygonize_cut_edges
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:polygonize_cut_edges)
-
-    geom_a = read(
-      'GEOMETRYCOLLECTION(
-        LINESTRING(0 0, 10 10),
-        LINESTRING(185 221, 100 100),
-        LINESTRING(185 221, 88 275, 180 316),
-        LINESTRING(185 221, 292 281, 180 316),
-        LINESTRING(189 98, 83 187, 185 221),
-        LINESTRING(189 98, 325 168, 185 221)
-      )'
-    )
-
-    cut_edges = geom_a.polygonize_cut_edges
-    assert_equal(0, cut_edges.length)
-  end
-
-  def test_polygonize_full
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:polygonize_full)
-
-    writer.rounding_precision = if Geos::GEOS_NICE_VERSION >= '031000'
-      0
-    else
-      3
-    end
-
-    geom_a = read(
-      'GEOMETRYCOLLECTION(
-        LINESTRING(0 0, 10 10),
-        LINESTRING(185 221, 100 100),
-        LINESTRING(185 221, 88 275, 180 316),
-        LINESTRING(185 221, 292 281, 180 316),
-        LINESTRING(189 98, 83 187, 185 221),
-        LINESTRING(189 98, 325 168, 185 221)
-      )'
-    )
-
-    polygonized = geom_a.polygonize_full
-
-    assert_kind_of(Array, polygonized[:rings])
-    assert_kind_of(Array, polygonized[:cuts])
-    assert_kind_of(Array, polygonized[:dangles])
-    assert_kind_of(Array, polygonized[:invalid_rings])
-
-    assert_equal(2, polygonized[:rings].length)
-    assert_equal(0, polygonized[:cuts].length)
-    assert_equal(2, polygonized[:dangles].length)
-    assert_equal(0, polygonized[:invalid_rings].length)
-
-    assert_equal(
-      'POLYGON ((185 221, 88 275, 180 316, 292 281, 185 221))',
-      write(polygonized[:rings][0])
-    )
-
-    assert_equal(
-      'POLYGON ((189 98, 83 187, 185 221, 325 168, 189 98))',
-      write(polygonized[:rings][1])
-    )
-
-    assert_equal(
-      'LINESTRING (185 221, 100 100)',
-      write(polygonized[:dangles][0])
-    )
-
-    assert_equal(
-      'LINESTRING (0 0, 10 10)',
-      write(polygonized[:dangles][1])
-    )
-
-    geom_b = geom_a.union(read('POINT(0 0)'))
-    polygonized = geom_b.polygonize_full
-
-    assert_equal(2, polygonized[:dangles].length)
-    assert_equal(0, polygonized[:invalid_rings].length)
-
-    assert_equal(
-      'LINESTRING (132 146, 100 100)',
-      write(polygonized[:dangles][0])
-    )
-
-    assert_equal(
-      'LINESTRING (0 0, 10 10)',
-      write(polygonized[:dangles][1])
-    )
-  end
-
-  def test_polygonize_with_bad_arguments
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:polygonize_full)
-
-    assert_raises(ArgumentError) do
-      geom = read('POINT(0 0)')
-      geom.polygonize(geom, 'gibberish')
-    end
-  end
-
-  def test_build_area
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:build_area)
-
-    geom = read('GEOMETRYCOLLECTION (LINESTRING(0 0, 0 1, 1 1), LINESTRING (1 1, 1 0, 0 0))')
-
-    assert_equal('POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))', write(geom.build_area))
-  end
-
-  def test_make_valid
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:make_valid)
-
-    geom = read('POLYGON((0 0, 1 1, 0 1, 1 0, 0 0))')
-
-    assert_equal(
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'MULTIPOLYGON (((1 0, 0 0, 0.5 0.5, 1 0)), ((1 1, 0.5 0.5, 0 1, 1 1)))'
-      else
-        'MULTIPOLYGON (((0 0, 0.5 0.5, 1 0, 0 0)), ((0.5 0.5, 0 1, 1 1, 0.5 0.5)))'
-      end,
-      write(geom.make_valid)
-    )
-  end
-
-  def test_shared_paths
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:shared_paths)
-
-    geom_a = read('LINESTRING(0 0, 50 0)')
-    geom_b = read('MULTILINESTRING((5 0, 15 0),(40 0, 30 0))')
-
-    paths = geom_a.shared_paths(geom_b)
-    assert_equal(2, paths.length)
-    assert_equal(
-      'MULTILINESTRING ((5 0, 15 0))',
-      write(paths[0])
-    )
-    assert_equal(
-      'MULTILINESTRING ((30 0, 40 0))',
-      write(paths[1])
-    )
-  end
-
-  def test_clone
-    geom_a = read('POINT(0 0)')
-    geom_b = geom_a.clone
-
-    assert_equal(geom_a, geom_b)
-  end
-
-  def test_clone_srid
-    srid = 4326
-    geom_a = read('POINT(0 0)')
-    geom_a.srid = srid
-    geom_b = geom_a.clone
-
-    assert_equal(geom_a, geom_b)
-    assert_equal(srid, geom_b.srid)
-  end
-
-  def test_dup
-    geom_a = read('POINT(0 0)')
-    geom_b = geom_a.dup
-
-    assert_equal(geom_a, geom_b)
-  end
-
-  def test_dup_srid
-    srid = 4326
-    geom_a = read('POINT(0 0)')
-    geom_a.srid = srid
-    geom_b = geom_a.dup
-    assert_equal(geom_a, geom_b)
-    assert_equal(srid, geom_b.srid)
-  end
-
-  def test_line_string_enumerator
-    geom = read('LINESTRING(0 0, 10 10))')
-    assert_kind_of(Enumerable, geom.each)
-    assert_kind_of(Enumerable, geom.to_enum)
-    assert_equal(geom, geom.each(&EMPTY_BLOCK))
-  end
-
-  def test_normalize
-    geom = read('POLYGON((0 0, 5 0, 5 5, 0 5, 0 0))')
-    geom.normalize
-    assert_equal('POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))', write(geom))
-
-    geom = read('POLYGON((0 0, 5 0, 5 5, 0 5, 0 0))').normalize
-    assert_equal('POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))', write(geom))
-  end
-
-  def test_normalize_bang
-    geom = read('POLYGON((0 0, 5 0, 5 5, 0 5, 0 0))')
-    geom.normalize!
-    assert_equal('POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))', write(geom))
-
-    geom = read('POLYGON((0 0, 5 0, 5 5, 0 5, 0 0))').normalize!
-    assert_equal('POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))', write(geom))
-  end
-
-  def test_eql
-    geom_a = read('POINT(1.0 1.0)')
-    geom_b = read('POINT(2.0 2.0)')
-
-    %w{ eql? equals? }.each do |method|
-      assert(geom_a.send(method, geom_a), "Expected geoms to be equal using #{method}")
-      refute(geom_a.send(method, geom_b), "Expected geoms to not be equal using #{method}")
-    end
-  end
-
-  def test_equals_operator
-    geom_a = read('POINT(1.0 1.0)')
-    geom_b = read('POINT(2.0 2.0)')
-
-    assert(geom_a == geom_a, 'Expected geoms to be equal using ==')
-    refute(geom_a == geom_b, 'Expected geoms to not be equal using ==')
-    refute(geom_a == 'test', 'Expected geoms to not be equal using ==')
-  end
-
-  def test_eql_exact
-    geom_a = read('POINT(1.0 1.0)')
-    geom_b = read('POINT(2.0 2.0)')
-
-    %w{ eql_exact? equals_exact? exactly_equals? }.each do |method|
-      refute(geom_a.send(method, geom_b, 0.001), "Expected geoms to not be equal using #{method}")
-    end
-  end
-
-  def test_eql_almost_default
-    geom = read('POINT (1 1)')
-    geom_a = read('POINT (1.0000001 1.0000001)')
-    geom_b = read('POINT (1.000001 1.000001)')
-
-    %w{ eql_almost? equals_almost? almost_equals? }.each do |method|
-      assert(geom.send(method, geom_a), "Expected geoms to be equal using #{method}")
-      refute(geom.send(method, geom_b), "Expected geoms to not be equal using #{method}")
-    end
-  end
-
-  def test_eql_almost
-    geom_a = read('POINT(1.0 1.0)')
-    geom_b = read('POINT(1.1 1.1)')
-
-    refute_equal(geom_a, geom_b)
-
-    %w{ eql_almost? equals_almost? almost_equals? }.each do |method|
-      assert(geom_a.send(method, geom_b, 0), "Expected geoms to be equal using #{method}")
-      refute(geom_a.send(method, geom_b, 1), "Expected geoms to not be equal using #{method}")
-    end
-  end
-
-  def test_srid_copy_policy
-    geom = read('POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))')
-    geom.srid = 4326
-
-    Geos.srid_copy_policy = :zero
-    cloned = geom.clone
-    assert_equal(4326, cloned.srid)
-
-    Geos.srid_copy_policy = :lenient
-    cloned = geom.clone
-    assert_equal(4326, cloned.srid)
-
-    Geos.srid_copy_policy = :strict
-    cloned = geom.clone
-    assert_equal(4326, cloned.srid)
-
-    Geos.srid_copy_policy = :zero
-    geom_b = geom.convex_hull
-    assert_equal(0, geom_b.srid)
-
-    Geos.srid_copy_policy = :lenient
-    geom_b = geom.convex_hull
-    assert_equal(4326, geom_b.srid)
-
-    Geos.srid_copy_policy = :strict
-    geom_b = geom.convex_hull
-    assert_equal(4326, geom_b.srid)
-
-    geom_b = read('POLYGON ((3 3, 3 8, 8 8, 8 3, 3 3))')
-    geom_b.srid = 3875
-
-    Geos.srid_copy_policy = :zero
-    geom_c = geom.intersection(geom_b)
-    assert_equal(0, geom_c.srid)
-
-    Geos.srid_copy_policy = :lenient
-    geom_c = geom.intersection(geom_b)
-    assert_equal(4326, geom_c.srid)
-
-    assert_raises(Geos::MixedSRIDsError) do
-      Geos.srid_copy_policy = :strict
-      geom_c = geom.intersection(geom_b)
-      assert_equal(231_231, geom_c.srid)
-    end
-  ensure
-    Geos.srid_copy_policy = :default
-  end
-
-  def test_bad_srid_copy_policy
-    assert_raises(ArgumentError) do
-      Geos.srid_copy_policy = :blart
-    end
-  end
-
-  def test_srid_copy_policy_default
-    Geos.srid_copy_policy_default = :default
-    assert_equal(:zero, Geos.srid_copy_policy_default)
-
-    Geos.srid_copy_policy_default = :lenient
-    assert_equal(:lenient, Geos.srid_copy_policy_default)
-
-    Geos.srid_copy_policy_default = :strict
-    assert_equal(:strict, Geos.srid_copy_policy_default)
-
-    assert_raises(ArgumentError) do
-      Geos.srid_copy_policy_default = :blart
-    end
-  ensure
-    Geos.srid_copy_policy_default = :default
-  end
-
-  def test_empty_geometry_has_0_area
-    assert_equal(0, read('POLYGON EMPTY').area)
-  end
-
-  def test_empty_geometry_has_0_length
-    assert_equal(0, read('POLYGON EMPTY').length)
-  end
-
-  def test_to_s
-    assert_match(/^\#<Geos::Point: .+>$/, read('POINT(0 0)').to_s)
-  end
-
-  def test_delaunay_triangulation
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:delaunay_triangulation)
-
-    tester = lambda { |expected, geom, *args|
-      geom = read(geom)
-      geom_tri = geom.delaunay_triangulation(*args)
-      geom_tri.normalize!
-
-      assert_equal(expected, write(geom_tri))
-    }
-
-    writer.trim = true
-
-    # empty polygon
-    tester['GEOMETRYCOLLECTION EMPTY', 'POLYGON EMPTY', 0]
-    tester['MULTILINESTRING EMPTY', 'POLYGON EMPTY', 0, only_edges: true]
-
-    # single point
-    tester['GEOMETRYCOLLECTION EMPTY', 'POINT (0 0)', 0]
-    tester['MULTILINESTRING EMPTY', 'POINT (0 0)', 0, only_edges: true]
-
-    # three collinear points
-    tester['GEOMETRYCOLLECTION EMPTY', 'MULTIPOINT(0 0, 5 0, 10 0)', 0]
-    tester['MULTILINESTRING ((5 0, 10 0), (0 0, 5 0))', 'MULTIPOINT(0 0, 5 0, 10 0)', 0, only_edges: true]
-
-    # three points
-    tester['GEOMETRYCOLLECTION (POLYGON ((0 0, 10 10, 5 0, 0 0)))', 'MULTIPOINT(0 0, 5 0, 10 10)', 0]
-    tester['MULTILINESTRING ((5 0, 10 10), (0 0, 10 10), (0 0, 5 0))', 'MULTIPOINT(0 0, 5 0, 10 10)', 0, only_edges: true]
-
-    # polygon with a hole
-    tester[
-      'GEOMETRYCOLLECTION (POLYGON ((8 2, 10 10, 8.5 1, 8 2)), POLYGON ((7 8, 10 10, 8 2, 7 8)), POLYGON ((3 8, 10 10, 7 8, 3 8)), ' \
-      'POLYGON ((2 2, 8 2, 8.5 1, 2 2)), POLYGON ((2 2, 7 8, 8 2, 2 2)), POLYGON ((2 2, 3 8, 7 8, 2 2)), POLYGON ((0.5 9, 10 10, 3 8, 0.5 9)), ' \
-      'POLYGON ((0.5 9, 3 8, 2 2, 0.5 9)), POLYGON ((0 0, 2 2, 8.5 1, 0 0)), POLYGON ((0 0, 0.5 9, 2 2, 0 0)))',
-      'POLYGON((0 0, 8.5 1, 10 10, 0.5 9, 0 0),(2 2, 3 8, 7 8, 8 2, 2 2)))',
-      0
-    ]
-
-    tester[
-      'MULTILINESTRING ((8.5 1, 10 10), (8 2, 10 10), (8 2, 8.5 1), (7 8, 10 10), (7 8, 8 2), (3 8, 10 10), (3 8, 7 8), (2 2, 8.5 1), (2 2, 8 2), (2 2, 7 8), (2 2, 3 8), (0.5 9, 10 10), (0.5 9, 3 8), (0.5 9, 2 2), (0 0, 8.5 1), (0 0, 2 2), (0 0, 0.5 9))',
-      'POLYGON((0 0, 8.5 1, 10 10, 0.5 9, 0 0),(2 2, 3 8, 7 8, 8 2, 2 2)))',
-      0,
-      only_edges: true
-    ]
-
-    # four points with a tolerance making one collapse
-    tester['MULTILINESTRING ((10 0, 10 10), (0 0, 10 10), (0 0, 10 0))', 'MULTIPOINT(0 0, 10 0, 10 10, 11 10)', 2.0, only_edges: true]
-
-    # tolerance as an option
-    tester['MULTILINESTRING ((10 0, 10 10), (0 0, 10 10), (0 0, 10 0))', 'MULTIPOINT(0 0, 10 0, 10 10, 11 10)', tolerance: 2.0, only_edges: true]
-  end
-
-  def test_voronoi_diagram
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:voronoi_diagram)
-
-    tester = lambda { |expected, geom, *args|
-      geom = read(geom)
-      voronoi_diagram = geom.voronoi_diagram(*args)
-
-      assert_equal(expected, write(voronoi_diagram))
-    }
-
-    writer.trim = true
-
-    geom = 'MULTIPOINT(0 0, 100 0, 100 100, 0 100)'
-
-    tester[
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'GEOMETRYCOLLECTION (POLYGON ((200 200, 200 50, 50 50, 50 200, 200 200)), POLYGON ((-100 200, 50 200, 50 50, -100 50, -100 200)), POLYGON ((-100 -100, -100 50, 50 50, 50 -100, -100 -100)), POLYGON ((200 -100, 50 -100, 50 50, 200 50, 200 -100)))'
-      else
-        'GEOMETRYCOLLECTION (POLYGON ((50 200, 200 200, 200 50, 50 50, 50 200)), POLYGON ((-100 50, -100 200, 50 200, 50 50, -100 50)), POLYGON ((50 -100, -100 -100, -100 50, 50 50, 50 -100)), POLYGON ((200 50, 200 -100, 50 -100, 50 50, 200 50)))'
-      end,
-      geom
-    ]
-
-    tester['MULTILINESTRING ((50 50, 50 200), (200 50, 50 50), (50 50, -100 50), (50 50, 50 -100))', geom, tolerance: 0, only_edges: true]
-
-    tester['MULTILINESTRING ((50 50, 50 1100), (1100 50, 50 50), (50 50, -1000 50), (50 50, 50 -1000))', geom,
-      only_edges: true,
-      envelope: read(geom).buffer(1000)
-    ]
-
-    # Allows a tolerance for the first argument
-    writer.rounding_precision = if Geos::GEOS_NICE_VERSION >= '031000'
-      0
-    else
-      3
-    end
-
-    writer.trim = true
-
-    tester[
-      if Geos::GEOS_NICE_VERSION > '030900'
-        'GEOMETRYCOLLECTION (POLYGON ((290 140, 185 140, 185 215, 188 235, 290 252, 290 140)), POLYGON ((80 340, 101 340, 188 235, 185 215, 80 215, 80 340)), POLYGON ((80 140, 80 215, 185 215, 185 140, 80 140)), POLYGON ((290 340, 290 252, 188 235, 101 340, 290 340)))'
-      else
-        'GEOMETRYCOLLECTION (POLYGON ((290 252, 290 140, 185 140, 185 215, 188 235, 290 252)), POLYGON ((80 215, 80 340, 101 340, 188 235, 185 215, 80 215)), POLYGON ((185 140, 80 140, 80 215, 185 215, 185 140)), POLYGON ((101 340, 290 340, 290 252, 188 235, 101 340)))'
-      end,
-      'MULTIPOINT ((150 210), (210 270), (150 220), (220 210), (215 269))',
-      10
-    ]
-  end
-
-  def test_precision
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:precision)
-
-    geom = read('POLYGON EMPTY')
-    scale = geom.precision
-    assert_equal(0.0, scale)
-
-    geom_with_precision = geom.with_precision(2.0)
-
-    assert_equal('POLYGON EMPTY', write(geom_with_precision))
-    scale = geom_with_precision.precision
-    assert_equal(2.0, scale)
-  end
-
-  def test_with_precision
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:with_precision)
-
-    geom = read('LINESTRING(1 0, 2 0)')
-
-    geom_with_precision = geom.with_precision(5.0)
-    assert_equal('LINESTRING EMPTY', write(geom_with_precision))
-
-    geom_with_precision = geom.with_precision(5.0, keep_collapsed: true)
-    assert_equal('LINESTRING (0 0, 0 0)', write(geom_with_precision))
-  end
-
-  def test_minimum_rotated_rectangle
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:minimum_rotated_rectangle)
-
-    geom = read('POLYGON ((1 6, 6 11, 11 6, 6 1, 1 6))')
-    minimum_rotated_rectangle = geom.minimum_rotated_rectangle
-
-    assert_equal('POLYGON ((6 1, 11 6, 6 11, 1 6, 6 1))', write(minimum_rotated_rectangle))
-  end
-
-  def test_minimum_clearance
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:minimum_clearance)
-
-    tester = lambda { |expected_clearance, geom|
-      geom = read(geom)
-      clearance = geom.minimum_clearance
-
-      if expected_clearance.eql?(Float::INFINITY)
-        assert(clearance.infinite?)
-      else
-        assert_in_delta(expected_clearance, clearance, TOLERANCE)
-      end
-    }
-
-    tester[Float::INFINITY, 'LINESTRING EMPTY']
-    tester[20, 'LINESTRING (30 100, 10 100)']
-    tester[100, 'LINESTRING (200 200, 200 100)']
-    tester[3.49284983912134e-05, 'LINESTRING (-112.712119 33.575919, -112.712127 33.575885)']
-  end
-
-  def test_minimum_clearance_line
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:minimum_clearance_line)
-
-    tester = lambda { |expected_geom, geom|
-      geom = read(geom)
-      clearance_geom = geom.minimum_clearance_line
-
-      assert_equal(expected_geom, write(clearance_geom))
-    }
-
-    tester['LINESTRING EMPTY', 'MULTIPOINT ((100 100), (100 100))']
-    tester['LINESTRING (30 100, 10 100)', 'MULTIPOINT ((100 100), (10 100), (30 100))']
-    tester['LINESTRING (200 200, 200 100)', 'POLYGON ((100 100, 300 100, 200 200, 100 100))']
-    tester[
-      'LINESTRING (-112.712119 33.575919, -112.712127 33.575885)',
-      '0106000000010000000103000000010000001a00000035d42824992d5cc01b834e081dca404073b9c150872d5cc03465a71fd4c940400ec00644882d5cc03b8a' \
-      '73d4d1c94040376dc669882d5cc0bf9cd9aed0c940401363997e892d5cc002f4fbfecdc94040ca4e3fa88b2d5cc0a487a1d5c9c940408f1ce90c8c2d5cc06989' \
-      '95d1c8c94040fab836548c2d5cc0bd175fb4c7c940409f1f46088f2d5cc0962023a0c2c940407b15191d902d5cc068041bd7bfc940400397c79a912d5cc0287d' \
-      '21e4bcc940403201bf46922d5cc065e3c116bbc940409d9d0c8e922d5cc0060fd3beb9c940400ef7915b932d5cc09012bbb6b7c940404fe61f7d932d5cc0e4a0' \
-      '8499b6c94040fc71fbe5932d5cc0ea9106b7b5c94040eaec6470942d5cc0c2323674b3c94040601dc70f952d5cc043588d25acc94040aea06989952d5cc03ecf' \
-      '9f36aac94040307f85cc952d5cc0e5eb32fca7c94040dd0a6135962d5cc01b615111a7c9404048a7ae7c962d5cc00a2aaa7ea5c94040f4328ae5962d5cc05eb8' \
-      '7361a4c94040c49448a2972d5cc04d81cccea2c940407c80eecb992d5cc06745d4449fc9404035d42824992d5cc01b834e081dca4040'
-    ]
-    tester['LINESTRING EMPTY', 'POLYGON EMPTY']
-  end
-
-  def test_maximum_inscribed_circle
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:maximum_inscribed_circle)
-
-    geom = read('POLYGON ((100 200, 200 200, 200 100, 100 100, 100 200))')
-    output = geom.maximum_inscribed_circle(0.001)
-    assert_equal('LINESTRING (150 150, 150 200)', write(output))
-  end
-
-  def test_largest_empty_circle
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:largest_empty_circle)
-
-    geom = read('MULTIPOINT ((100 100), (100 200), (200 200), (200 100))')
-    output = geom.largest_empty_circle(0.001)
-    assert_equal('LINESTRING (150 150, 100 100)', write(output))
-
-    geom = read('MULTIPOINT ((100 100), (100 200), (200 200), (200 100))')
-    output = geom.largest_empty_circle(0.001, boundary: read('MULTIPOINT ((100 100), (100 200), (200 200), (200 100))'))
-    assert_equal('LINESTRING (100 100, 100 100)', write(output))
-  end
-
-  def test_minimum_width
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:minimum_width)
-
-    geom = read('POLYGON ((0 0, 0 15, 5 10, 5 0, 0 0))')
-    output = geom.minimum_width
-    assert_equal('LINESTRING (0 0, 5 0)', write(output))
-
-    geom = read('LINESTRING (0 0,0 10, 10 10)')
-    output = geom.minimum_width
-    assert_equal('LINESTRING (5 5, 0 10)', write(output))
-  end
-
-  def test_dump_points
-    geom = read('GEOMETRYCOLLECTION(
-      MULTILINESTRING((0 0, 10 10, 20 20), (100 100, 200 200, 300 300)),
-
-      POINT(10 10),
-
-      POLYGON((0 0, 5 0, 5 5, 0 5, 0 0), (1 1, 4 1, 4 4, 1 4, 1 1))
-    )')
-
-    assert_equal([
-      [
-        [
-          Geos.create_point(0, 0),
-          Geos.create_point(10, 10),
-          Geos.create_point(20, 20)
-        ],
-
-        [
-          Geos.create_point(100, 100),
-          Geos.create_point(200, 200),
-          Geos.create_point(300, 300)
-        ]
-      ],
-
-      [
-        Geos.create_point(10, 10)
-      ],
-
-      [
-        [
-          Geos.create_point(0, 0),
-          Geos.create_point(5, 0),
-          Geos.create_point(5, 5),
-          Geos.create_point(0, 5),
-          Geos.create_point(0, 0)
-        ],
-
-        [
-          Geos.create_point(1, 1),
-          Geos.create_point(4, 1),
-          Geos.create_point(4, 4),
-          Geos.create_point(1, 4),
-          Geos.create_point(1, 1)
-        ]
-      ]
-    ], geom.dump_points)
-  end
-
-  def test_reverse
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:reverse)
-
-    simple_tester(:reverse, 'POINT (3 5)', 'POINT (3 5)')
-    simple_tester(:reverse, 'MULTIPOINT (100 100, 10 100, 30 100)', 'MULTIPOINT (100 100, 10 100, 30 100)')
-    simple_tester(:reverse, 'LINESTRING (200 200, 200 100)', 'LINESTRING (200 100, 200 200)')
-
-    if Geos::GEOS_NICE_VERSION >= '030801'
-      simple_tester(:reverse, 'MULTILINESTRING ((3 3, 4 4), (1 1, 2 2))', 'MULTILINESTRING ((4 4, 3 3), (2 2, 1 1))')
-    else
-      simple_tester(:reverse, 'MULTILINESTRING ((1 1, 2 2), (3 3, 4 4))', 'MULTILINESTRING ((4 4, 3 3), (2 2, 1 1))')
-    end
-
-    simple_tester(:reverse, 'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 2 1, 2 2, 1 2, 1 1))', 'POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0), (1 1, 1 2, 2 2, 2 1, 1 1))')
-    simple_tester(:reverse, 'MULTIPOLYGON (((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 2 1, 2 2, 1 2, 1 1)), ((100 100, 100 200, 200 200, 100 100)))', 'MULTIPOLYGON (((0 0, 0 10, 10 10, 10 0, 0 0), (1 1, 1 2, 2 2, 2 1, 1 1)), ((100 100, 200 200, 100 200, 100 100)))')
-    simple_tester(:reverse, 'GEOMETRYCOLLECTION (LINESTRING (1 1, 2 2), GEOMETRYCOLLECTION (LINESTRING (3 5, 2 9)))', 'GEOMETRYCOLLECTION (LINESTRING (2 2, 1 1), GEOMETRYCOLLECTION(LINESTRING (2 9, 3 5)))')
-    simple_tester(:reverse, 'POINT EMPTY', 'POINT EMPTY')
-    simple_tester(:reverse, 'LINESTRING EMPTY', 'LINESTRING EMPTY')
-    simple_tester(:reverse, 'LINEARRING EMPTY', 'LINEARRING EMPTY')
-    simple_tester(:reverse, 'POLYGON EMPTY', 'POLYGON EMPTY')
-    simple_tester(:reverse, 'MULTIPOINT EMPTY', 'MULTIPOINT EMPTY')
-    simple_tester(:reverse, 'MULTILINESTRING EMPTY', 'MULTILINESTRING EMPTY')
-    simple_tester(:reverse, 'MULTIPOLYGON EMPTY', 'MULTIPOLYGON EMPTY')
-    simple_tester(:reverse, 'GEOMETRYCOLLECTION EMPTY', 'GEOMETRYCOLLECTION EMPTY')
-  end
-
-  def test_frechet_distance
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:frechet_distance)
-
-    assert_in_delta(read('LINESTRING (0 0, 100 0)').frechet_distance(read('LINESTRING (0 0, 50 50, 100 0)')), 70.7106781186548, TOLERANCE)
-  end
-
-  def test_frechet_distance_with_densify
-    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:frechet_distance)
-
-    assert_in_delta(read('LINESTRING (0 0, 100 0)').frechet_distance(read('LINESTRING (0 0, 50 50, 100 0)'), 0.5), 50.0, TOLERANCE)
-  end
-end