How do I get the number of cells around each node

[zyhsjtu]

How do I get the number of cells around each node?

[KnutAM]

You can use the topology features.

grid = generate_grid(Quadrilateral, (3,3));
top = ExclusiveTopology(grid);
cells_connected_to_node_2 = top.vertex_to_cell[2]

Note that this only works for nodes at the vertices of the cell (i.e. not for nodes on faces / edges / internal in cells existing for higher-order geometric interpolations)

[zyhsjtu]

Thanks for your quick response!