edges_iter = ((i, j, {'boundary-ids': {edge_map[i, j]}})
for i in range(edge_map.shape[0])
for j in edge_map.indices[edge_map.indptr[i]:
edge_map.indptr[i+1]])
or via COO:
coo = edge_map.tocoo()
edges_iter = ((i, j, {'boundary-ids': {k}}) for i, j, k in
zip(coo.row, coo.col, coo.data))
They are unlikely to be equal in performance, but performance might be data-dependent.
We have two options:
They are unlikely to be equal in performance, but performance might be data-dependent.