This paper describes a general algorithm for finding the commensurator of a nonarithmetic hyperbolic manifold with cusps and for deciding when two such manifolds are commensurable. The method is based on some elementary observations regarding horosphere packings and canonical cell decompositions. For example, we use this to find the commensurators of all nonarithmetic hyperbolic once-punctured torus bundles over the circle.
For hyperbolic 3-manifolds, the algorithm has been implemented using Goodman’s computer program Snap. We use this to determine the commensurability classes of all cusped hyperbolic 3-manifolds triangulated using at most seven ideal tetrahedra, and for the complements of hyperbolic knots and links with up to twelve crossings.
"Commensurators of Cusped Hyperbolic Manifolds." Experiment. Math. 17 (3) 283 - 306, 2008.