Added test

Author: benvanbasten-ns

threedigrid_builder/tests/test_grid.py

@@ -388,6 +388,28 @@ def test_split_partial_cut(cell_polygon):
     assert fragments[0].equals_exact(cell_polygon, tolerance=0.0)
 
 
+def test_split_order(cell_polygon):
+    # This test is to indicate that the order matters
+    cutlines_2 = [
+        shapely.LineString([(5, 0), (5, 6)]),
+        shapely.LineString([(0, 0), (10, 10)]),
+    ]
+    fragments = Grid.split(cell_polygon, cutlines_2)
+    assert len(fragments) == 2
+
+    cutlines_3 = [
+        shapely.LineString([(0, 0), (10, 10)]),
+        shapely.LineString([(5, 0), (5, 5)]),
+    ]
+    fragments = Grid.split(cell_polygon, cutlines_3)
+    assert len(fragments) == 3
+
+    # Combining the lines in a single linestring solves the problem
+    cutlines_3 = [shapely.LineString([(0, 0), (10, 10), (5, 5), (5, 0)])]
+    fragments = Grid.split(cell_polygon, cutlines_3)
+    assert len(fragments) == 3
+
+
 def test_split_double_cut(cell_polygon):
     cutlines = [
         shapely.LineString([(0, 0), (10, 10)]),