Band width reduction as internal step

[brandimarte]

The band width reduction of the adjacency should be made internally when calling, for example, allKekules (if this is really need at all). When doing this the function should return the Kekule structures with respect to the indices of the original adjacency matrix.

[tfrederiksen]

Is there any real gain by reordering when the adjacency matrix is represented with a normal (dense) numpy array? That is, not some sparse format?