Wrong order of dofs after creating mesh from lists of points and cells

You can use:

o_cell_idx =  mesh.topology.original_cell_index
for ii in range(len(conn1)):
    assert np.allclose(conn1[ii], ee[o_cell_idx[ii]])
2 Likes