Fix tiling in trirefine.py

rth · rth · commit fec1fa075fa4 · 2019-01-22T16:49:49.000+01:00
diff --git a/lib/matplotlib/tri/trirefine.py b/lib/matplotlib/tri/trirefine.py
@@ -238,7 +238,7 @@ def _refine_triangulation_once(triangulation, ancestors=None):
         # (can be -1 if border edge)
         # For slave and master we will identify the apex pointing to the edge
         # start
-        edge_elems = np.tile(np.arange(3, dtype=np.int32), ntri)
+        edge_elems = np.tile(np.arange(ntri, dtype=np.int32), 3)
         edge_apexes = np.repeat(np.arange(3, dtype=np.int32), ntri)
         edge_neighbors = neighbors[edge_elems, edge_apexes]
         mask_masters = (edge_elems > edge_neighbors)
