I need to compute **canonical forms** of *many* (~10^6-10^8) vertex-facets incidence graphs of polytope. Two rather big examples I want to consider are

* the [600-cell](https://en.wikipedia.org/wiki/600-cell) with 120 vertices and 600 facets (dimension 4), and
* the smallest known [counter example](https://sites.google.com/site/christopheweibel/research/hirsch-conjecture) to the Hirsch bound with 40 vertices and 36426 facets (dimension 20).

The two options that seem to be best suited to compute these canonical forms are [nauty](http://www3.cs.stonybrook.edu/~algorith/implement/nauty/implement.shtml) and [bliss](http://www.tcs.hut.fi/Software/bliss/).
My questions are

* Is there another option that I have overseen?

* I have not found any benchmark comparisons between the two, so should I prefer one over the other?

* Does the property of being bipartite make a difference in which to choose?

Many Thanks! -- this is not a strictly mathematical question, but I hope it is still suitable here.