Yes, unfortunately `dmplex.closure_ordering()` is still the really expensive operation here. In the next round of optimisations we probably want to push building the `cell_node_list` to C as well, but there is only limited gain in that.