diff --git a/meshmode/mesh/generation.py b/meshmode/mesh/generation.py index c19c3dd1e6df21f3d197bc4750e616621121c21a..64968fe4ee827d74a47060b166fc2f060badffae 100644 --- a/meshmode/mesh/generation.py +++ b/meshmode/mesh/generation.py @@ -395,12 +395,12 @@ def generate_icosahedron(r, order): top_point = 5 tris = [] - l = len(top_ring) - for i in range(l): - tris.append([top_ring[i], top_ring[(i+1) % l], top_point]) - tris.append([bottom_ring[i], bottom_point, bottom_ring[(i+1) % l], ]) - tris.append([bottom_ring[i], bottom_ring[(i+1) % l], top_ring[i]]) - tris.append([top_ring[i], bottom_ring[(i+1) % l], top_ring[(i+1) % l]]) + m = len(top_ring) + for i in range(m): + tris.append([top_ring[i], top_ring[(i+1) % m], top_point]) + tris.append([bottom_ring[i], bottom_point, bottom_ring[(i+1) % m], ]) + tris.append([bottom_ring[i], bottom_ring[(i+1) % m], top_ring[i]]) + tris.append([top_ring[i], bottom_ring[(i+1) % m], top_ring[(i+1) % m]]) vertices *= r/la.norm(vertices[:, 0])