snappy 3.1.1__cp311-cp311-win_amd64.whl → 3.2__cp311-cp311-win_amd64.whl

snappy/test.py

@@ -1,26 +1,31 @@
+import os
 import sys
-import getopt
+import glob
+import argparse
 import snappy
 import snappy.snap.test
 import spherogram.test
 import snappy.matrix
 import snappy.verify.test
+import snappy.geometric_structure.test
 import snappy.ptolemy.test
-import snappy.raytracing.cohomology_fractal
-import snappy.raytracing.geodesic_tube_info
-import snappy.raytracing.geodesics
-import snappy.raytracing.ideal_raytracing_data
-import snappy.raytracing.upper_halfspace_utilities
-import snappy.drilling
+import snappy.tiling.test
+import snappy.cusps.test
+import snappy.raytracing.test
+import snappy.len_spec.test
+import snappy.drilling.test
+import snappy.canonical
+import snappy.cusps.cusp_area_matrix
+import snappy.cusps.cusp_areas_from_matrix
+import snappy.isometry_signature
 import snappy.exterior_to_link.test
 import snappy.pari
+import snappy.test_cases
+import snappy.numeric_output_checker
 
-from snappy.sage_helper import (_within_sage, doctest_modules, cyopengl_works,
-                                tk_root, root_is_fake, DocTestParser)
-from snappy import numeric_output_checker
-modules = []
-
-snappy.database.Manifold = snappy.SnapPy.Manifold
+from snappy.sage_helper import _within_sage
+from snappy.testing import (doctest_modules, cyopengl_works,
+                            tk_root, root_is_fake, DocTestParser)
 
 # Augment tests for SnapPy with those that Cython missed
 
@@ -43,91 +48,55 @@ browser_tests = [x for x in snappy.SnapPyHP.__test__
 for key in identify_tests + triangulation_tests + browser_tests:
     snappy.SnapPyHP.__test__.pop(key)
 
+def additional_doctests(verbose=False, print_info=True):
+    """
+    I noticed that some of my changes to move code from here into
+    its own files caused some loss of test coverage.
 
-def snap_doctester(verbose):
-    return snappy.snap.test.run_doctests(verbose, print_info=False)
-
-
-snap_doctester.__name__ = 'snappy.snap'
-
-
-def snappy_database_doctester(verbose):
-    # snappy_manifolds's tests is still relying on
-    # SnapPy Number's _accuracy_for_testing.
-    #
-    # Switch to snappy conversion until snappy_manifolds is
-    # is updated.
-    snappy.number.use_field_conversion('snappy')
-    snappy.number.Number._accuracy_for_testing = 8
-    ans = doctest_modules([snappy.database], verbose)
-    snappy.number.Number._accuracy_for_testing = None
-    if _within_sage:
-        snappy.number.use_field_conversion('sage')
-
-    return ans
-
-
-snappy_database_doctester.__name__ = 'snappy.database'
-
-
-def spherogram_doctester(verbose):
-    ans = spherogram.test.run_doctests(verbose, print_info=False)
-
-    # Spherogram's testing is switching to SnapPy numbers and
-    # setting their accuracy.
-    # Switch back to Sage types until Spherogram has been updated.
-    snappy.number.Number._accuracy_for_testing = None
-    if _within_sage:
-        snappy.number.use_field_conversion('sage')
-
-    return ans
-
-
-spherogram_doctester.__name__ = 'spherogram'
-
-
-def ptolemy_doctester(verbose):
-    return snappy.ptolemy.test.run_doctests(verbose, print_info=False)
-
-
-ptolemy_doctester.__name__ = 'snappy.ptolemy'
-
-modules += [numeric_output_checker.run_doctests]
-
-if not _within_sage:
-    modules.append(snappy.number)
-
-modules += [snappy.SnapPy,
-            snappy.SnapPyHP,
-            snappy_database_doctester,
-            snappy,
-            snap_doctester,
-            snappy.matrix,
-            snappy.raytracing.cohomology_fractal,
-            snappy.raytracing.geodesic_tube_info,
-            snappy.raytracing.geodesics,
-            snappy.raytracing.ideal_raytracing_data,
-            snappy.raytracing.upper_halfspace_utilities,
-            snappy.drilling,
-            ptolemy_doctester,
-            spherogram_doctester]
+    In a panic, adding them explicitly here.
 
+    There ought to be a better way to do this...
+    """
 
-def snappy_verify_doctester(verbose):
-    return snappy.verify.test.run_doctests(verbose, print_info=False)
-
-
-snappy_verify_doctester.__name__ = 'snappy.verify'
-modules.append(snappy_verify_doctester)
-
-
-def snappy_exterior_to_link_doctester(verbose):
-    return snappy.exterior_to_link.test.run_doctests(verbose, print_info=False)
-
-
-snappy_exterior_to_link_doctester.__name__ = 'snappy.exterior_to_link'
-modules.insert(0, snappy_exterior_to_link_doctester)
-
+    globs = {'Manifold' : snappy.Manifold,
+             'ManifoldHP' : snappy.ManifoldHP}
+    return doctest_modules(
+        [ snappy.isometry_signature,
+          snappy.canonical,
+          snappy.cusps.cusp_area_matrix,
+          snappy.cusps.cusp_areas_from_matrix
+         ],
+        verbose=verbose,
+        print_info=print_info,
+        extraglobs = globs)
+additional_doctests.__name__ = 'snappy.<HARD TO REACH>'
+
+modules = [
+    snappy.exterior_to_link.test.run_doctests,
+    snappy.numeric_output_checker.run_doctests,
+    snappy.number,
+    snappy.SnapPy,
+    snappy.SnapPyHP,
+    snappy.database,
+    additional_doctests,
+    snappy,
+    snappy.snap.test.run_doctests,
+    snappy.matrix,
+    snappy.geometric_structure.test.run_doctests,
+    snappy.tiling.test.run_doctests,
+    snappy.cusps.test.run_doctests,
+    snappy.raytracing.test.run_doctests,
+    snappy.len_spec.test.run_doctests,
+    snappy.drilling.test.run_doctests,
+    snappy.ptolemy.test.run_doctests,
+    spherogram.test.run_doctests,
+    snappy.verify.test.run_doctests,
+    snappy.test_cases
+]
+
+slow_modules = [
+    snappy.ptolemy.test.run_ptolemy_tests
+]
 
 def graphics_failures(verbose, windows, use_modernopengl):
     if cyopengl_works():
@@ -138,7 +107,7 @@ def graphics_failures(verbose, windows, use_modernopengl):
         snappy.Manifold('m125').cusp_neighborhood().view().test()
         if use_modernopengl:
             snappy.Manifold('m004').inside_view().test()
-        snappy.Manifold('4_1').browse().test()
+        snappy.Manifold('4_1').browse().test(use_modernopengl=use_modernopengl)
         snappy.ManifoldHP('m004').dirichlet_domain().view().test()
         snappy.ManifoldHP('m125').cusp_neighborhood().view().test()
         if use_modernopengl:
@@ -150,77 +119,91 @@ def graphics_failures(verbose, windows, use_modernopengl):
                     print('Close the root window to finish.')
                 else:
                     print('The windows will close in a few seconds.\n'
-                        'Specify -w or --windows to avoid this.')
+                          'Specify -w or --windows to avoid this.')
                     root.after(7000, root.destroy)
                 root.mainloop()
     else:
-        print("***Warning***: CyOpenGL not installed, so not tested")
+        print("***Warning***: Could not test CyOpenGL.")
+        try:
+            import snappy.CyOpenGL
+            print("Reason: Unsuitable Tk configuration for CyOpenGL")
+        except ImportError as e:
+            print("Reason: CyOpenGL could not be imported, %r" % e)
         result = 0
     return result
 
 
+def check_for_docs():
+    path = os.path.dirname(snappy.__file__)
+    match = os.path.join(path, 'doc', '*.html')
+    html_files = [os.path.basename(file) for file in glob.glob(match)]
+    if len(html_files) > 20 and 'index.html' in html_files:
+        print('HTML docs are included.')
+    else:
+        print('WARNING: HTML docs are missing!')
+
+
 def runtests(verbose=False,
              quick=False,
              windows=False,
-             use_modernopengl=True):
+             use_modernopengl=True,
+             graphics=True):
 
     # The default PARI stacksize can (slightly) overflow, causing
     # doctests to fail.
     snappy.pari.allocatemem(2**24, 2**25, silent=True)
 
-    DocTestParser.use_modernopengl = use_modernopengl
+    DocTestParser.use_cymodernopengl = use_modernopengl
+
+    all_modules = modules
+    if not quick:
+        all_modules += slow_modules
+
+    result = doctest_modules(
+        all_modules, verbose=verbose, print_info=True)
 
-    result = doctest_modules(modules, verbose=verbose)
     if not quick:
-        print()
-        # No idea why we mess and set snappy.database.Manifold
-        # to SnapPy.Manifold above... But to make ptolemy work,
-        # temporarily setting it to what it should be.
-        original_db_manifold = snappy.database.Manifold
-        snappy.database.Manifold = snappy.Manifold
-        snappy.ptolemy.test.main(verbose=verbose, doctest=False)
-        snappy.database.Manifold = original_db_manifold
         print()
         spherogram.links.test.run()
     print('\nAll doctests:\n   %s failures out of %s tests.' % result)
 
-    num_graphics_failures = graphics_failures(
-        verbose=verbose,
-        windows=windows,
-        use_modernopengl=use_modernopengl)
+    if graphics:
+        num_graphics_failures = graphics_failures(
+            verbose=verbose,
+            windows=windows,
+            use_modernopengl=use_modernopengl)
+    else:
+        num_graphics_failures = 0
 
     print('Pari stacksize', snappy.pari.stacksize(),
           'max stack size', snappy.pari.stacksizemax())
+
+    check_for_docs()
     return result.failed + num_graphics_failures
 
 
 if __name__ == '__main__':
-
     verbose = False
     quick = False
     windows = False
     use_modernopengl = True
-
-    try:
-        useful_args = [arg for arg in sys.argv[1:] if not arg.startswith('-psn_')]
-        optlist, args = getopt.getopt(
-            useful_args,
-            'ivqws',
-            ['ignore', 'verbose', 'quick', 'windows', 'skip-modern-opengl'])
-        opts = [o[0] for o in optlist]
-        if '-v' in opts or '--verbose' in opts:
-            verbose = True
-        if '-q' in opts or '--quick' in opts:
-            quick = True
-        if '-w' in opts or '--windows' in opts:
-            windows = True
-        if '-s' in opts or '--skip-modern-opengl' in opts:
-            use_modernopengl = False
-
-    except getopt.GetoptError:
-        print("Could not parse arguments")
-
-    sys.exit(runtests(verbose=verbose,
-                      quick=quick,
-                      windows=windows,
-                      use_modernopengl=use_modernopengl))
+    graphics = True
+
+    useful_args = [arg for arg in sys.argv[1:] if not arg.startswith('-psn_')]
+
+    parser = argparse.ArgumentParser()
+    parser.add_argument('-v', '--verbose', action='store_true',
+                        help='show additional information')
+    parser.add_argument('-q', '--quick', action='store_true',
+                        help='skip ptolemy and spherogram.links tests.')
+    parser.add_argument('-w', '--windows', action='store_true',
+                        help='keep windows open until user closes root window.')
+    parser.add_argument('-s', '--skip-modern-opengl', action='store_false',
+                        dest='use_modernopengl',
+                        help='skip tests requiring OpenGL 3.2 or later.')
+    parser.add_argument('-g', '--skip-gui', action='store_false',
+                        dest='graphics',
+                        help='skip tests bringing up GUI windows.')
+    args = parser.parse_args(useful_args)
+
+    sys.exit(runtests(**vars(args)))

snappy/test_cases.py

@@ -0,0 +1,263 @@
+"""
+IMPORTANT: Python only recognises this as a doc string if there is
+nothing before it. In particular, add any includes after the doc string.
+
+Doc tests that we do not want to show up in the documentation.
+
+Decorated isomorphism signature
+-------------------------------
+
+An example where the isosig doesn't pick up the identity perm
+because the lexicographic order of the perm interacts with the
+matrices.
+
+>>> T = Triangulation("L6n1(0,0)(0,0)(0,0)")
+>>> T.triangulation_isosig()
+'gLMzQbcdefffaelaaai_acbBaabCbbabbBC'
+
+Test that slopes are computed correctly.
+
+>>> M=Manifold("L14n63023(-5,1)(5,1)(10,1)")
+>>> M.triangulation_isosig(decorated=False, ignore_orientation=False)
+'vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt'
+
+The canonical orientation (used to compute the unoriented isosig)
+is the reverse of the actual orientation:
+
+>>> M.triangulation_isosig(decorated=False)
+'vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth'
+>>> Mop = M.copy()
+>>> Mop.reverse_orientation()
+>>> Mop.triangulation_isosig(decorated=False, ignore_orientation=False)
+'vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth'
+
+So we expect the oriented isometry signature to flip when neither the isomorphism
+signature nor its decoration capture the orientation.
+
+The following is a baked version of snappy.decorated_isosig.test_slope_transformations().
+The 32 calls to ManifoldHP.isometry_signature are just too expensive for the test suite.
+
+>>> for ignore_cusp_ordering in [False, True]:
+...     for ignore_curves in [False, True]:
+...         for ignore_curve_orientations in [False, True]:
+...             for ignore_filling_orientations in [False, True]:
+...                 for ignore_orientation in [False, True]:
+...                     isosig = M.triangulation_isosig(
+...                         ignore_cusp_ordering = ignore_cusp_ordering,
+...                         ignore_curves = ignore_curves,
+...                         ignore_curve_orientations = ignore_curve_orientations,
+...                         ignore_filling_orientations = ignore_filling_orientations,
+...                         ignore_orientation = ignore_orientation)
+...                     print(isosig)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_dcbaabBBBbBaBbCbBbCb(-5,1)(5,1)(10,1)(0,0)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_cabdBacbbBCbaBBBbabB(-5,1)(5,1)(10,1)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_dcbaabBBBbBaBbCbBbCb(-5,1)(5,1)(10,1)(0,0)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_cabdBacbbBCbaBBBbabB(-5,1)(5,1)(10,1)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_dcbaabBBbBbabBcBbBcB(-5,1)(-5,-1)(-10,-1)(0,0)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_cabdbacbbBcBabBBbaBb(5,1)(5,-1)(-10,1)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_dcbaabBBbBbabBcBbBcB(-5,1)(5,1)(10,1)(0,0)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_cabdbacbbBcBabBBbaBb(5,1)(-5,1)(-10,1)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_dcba(-1,-6)(-6,5)(-12,11)(0,0)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_cabd(7,1)(3,-4)(-1,-11)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_dcba(1,6)(-6,5)(-12,11)(0,0)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_cabd(7,1)(-3,4)(1,11)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_dcba(-1,-6)(-6,5)(-12,11)(0,0)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_cabd(7,1)(3,-4)(-1,-11)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_dcba(1,6)(-6,5)(-12,11)(0,0)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_cabd(7,1)(-3,4)(1,11)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_BbCbBbCbBbBaabBB(0,0)(10,1)(5,1)(-5,1)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_aBBBBacbbBCbbabB(10,1)(-5,1)(5,1)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_BbCbBbCbBbBaabBB(0,0)(10,1)(5,1)(-5,1)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_aBBBBacbbBCbbabB(10,1)(-5,1)(5,1)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_bBcBbBcBbBbaabBB(0,0)(-10,-1)(-5,-1)(-5,1)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_abBBbacbbBcBbaBb(-10,1)(5,1)(5,-1)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt_bBcBbBcBbBbaabBB(0,0)(10,1)(5,1)(-5,1)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth_abBBbacbbBcBbaBb(-10,1)(5,1)(-5,1)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt(0,0)(-12,11)(-6,5)(6,-1)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth(-1,-11)(7,1)(3,-4)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt(0,0)(-12,11)(-6,5)(-6,1)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth(1,11)(7,1)(-3,4)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt(0,0)(-12,11)(-6,5)(6,-1)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth(-1,-11)(7,1)(3,-4)(0,0)
+vLLvvLLMALQQzQQceillmnppqrlmrqtruututiivimllaelaqxrvdoxqltt(0,0)(-12,11)(-6,5)(-6,1)
+vLLvLLPwPQLAMPQcefikkmnplkopqrsttutuuiixvimqlippawidlabavth(1,11)(7,1)(-3,4)(0,0)
+
+The special case where we do not encode the matrices and have the identity perm
+and need it:
+
+>>> T = Triangulation('ovLMvvPQQQccddlmnijklmnmnlgvfamtvfblhaumx_BBbacBbacBbaBCaB(312,241)(332,441)(3322,1441)(726,541)')
+>>> T.homology()
+Z/4684873938228
+>>> isosig = T.triangulation_isosig(ignore_curves=True, ignore_cusp_ordering=False)
+>>> isosig
+'ovLMvvPQQQccddlmnijklmnmnlgvfamtvfblhaumx_abcd(71,312)(-1105,332)(-8085,3322)(726,1993)'
+>>> Triangulation(isosig).homology()
+Z/4684873938228
+
+
+isometry_signature
+------------------
+
+>>> M = Manifold('m004(1,2)')
+>>> M.isometry_signature()
+'cPcbbbiht(3,2)'
+
+sage: M.isometry_signature(verified=True)
+'cPcbbbiht(3,2)'
+
+>>> M = ManifoldHP('m004(1,2)')
+>>> M.isometry_signature()
+'cPcbbbiht(3,2)'
+
+sage: M.isometry_signature(verified=True) # Test bug reported by Nathan
+'cPcbbbiht(3,2)'
+
+# Cases where the drilled manifold had a canonical cell decomposition
+# with non-tetrahedral cells.
+
+>>> M = Manifold('m137(3,2)')
+>>> M.isometry_signature()
+'sLLvwzvQPAQPQccghmiljkpmqnoorqrrqfafaoaqoofaoooqqaf(3,2)'
+
+sage: M.isometry_signature(verified=True)
+'sLLvwzvQPAQPQccghmiljkpmqnoorqrrqfafaoaqoofaoooqqaf(3,2)'
+
+>>> M = ManifoldHP('m137(3,2)')
+>>> M.isometry_signature()
+'sLLvwzvQPAQPQccghmiljkpmqnoorqrrqfafaoaqoofaoooqqaf(3,2)'
+
+sage: M.isometry_signature(verified=True)
+'sLLvwzvQPAQPQccghmiljkpmqnoorqrrqfafaoaqoofaoooqqaf(3,2)'
+
+sage: M.isometry_signature(verified=True, exact_bits_prec_and_degrees=[]) # doctest: +ELLIPSIS +IGNORE_EXCEPTION_DETAIL
+Traceback (most recent call last):
+...
+RuntimeError: Could not compute or verify canonical retriangulation of drilled manifold. Geodesic was: abCDaDAd.
+
+Test isometry_signature's ignore_orientation
+
+>>> M = Manifold("m006")
+>>> M.isometry_signature(ignore_orientation=False)
+'eLMkaccddjgbaj'
+>>> M.isometry_signature(ignore_orientation=True)
+'eLAkaccddngbak'
+>>> M.isometry_signature() # default value
+'eLAkaccddngbak'
+>>> M.dehn_fill((3,4))
+>>> M.isometry_signature(ignore_orientation=False)
+'eLMkaccddjgbaj(-1,4)'
+>>> M.isometry_signature(ignore_orientation=True)
+'eLAkaccddngbak(-3,4)'
+
+Test isometry_signature's of_link and filling.
+
+>>> Manifold("o9_44206(2,3)").isometry_signature(of_link=True)
+'jLLvMQQacggfiihhijkkjkehhtb_abBabBbabaab'
+
+Class hierarchy
+---------------
+
+>>> isinstance(Manifold("m004"), Triangulation)
+True
+
+>>> isinstance(ManifoldHP("m004"), TriangulationHP)
+True
+
+Low precision and high precision comparisons
+--------------------------------------------
+>>> M   = Manifold("m004")
+>>> Mhp = Manifold("m004")
+>>> N   = Manifold("m003")
+>>> Nhp = Manifold("m003")
+>>> M.is_isometric_to(M)
+True
+>>> M.is_isometric_to(Mhp)
+True
+>>> Mhp.is_isometric_to(M)
+True
+>>> Mhp.is_isometric_to(Mhp)
+True
+>>> M.is_isometric_to(N)
+False
+>>> M.is_isometric_to(Nhp)
+False
+>>> Mhp.is_isometric_to(N)
+False
+>>> Mhp.is_isometric_to(Nhp)
+False
+
+>>> O = Triangulation("mvvLALQQQhfghjjlilkjklaaaaaffffffff",
+... remove_finite_vertices = False)
+>>> O.has_finite_vertices()
+True
+>>> Ohp = TriangulationHP("mvvLALQQQhfghjjlilkjklaaaaaffffffff",
+... remove_finite_vertices = False)
+>>> Ohp.has_finite_vertices()
+True
+
+>>> len(O.isomorphisms_to(O))
+8
+>>> len(O.isomorphisms_to(Ohp))
+8
+>>> len(Ohp.isomorphisms_to(O))
+8
+>>> len(Ohp.isomorphisms_to(Ohp))
+8
+>>> len(M.isomorphisms_to(O))
+0
+>>> len(M.isomorphisms_to(Ohp))
+0
+>>> len(Mhp.isomorphisms_to(O))
+0
+>>> len(Mhp.isomorphisms_to(Ohp))
+0
+
+Canonical retriangulation
+-------------------------
+
+Some cases that should be rejected
+
+>>> M = Manifold("m004(3,4)")
+>>> M.canonical_retriangulation() # doctest: +ELLIPSIS
+Traceback (most recent call last):
+...
+ValueError: Canonical retriangulation needs all cusps to be complete.
+
+sage: M.canonical_retriangulation(verified=True) # doctest: +ELLIPSIS
+Traceback (most recent call last):
+...
+ValueError: Canonical retriangulation needs all cusps to be complete.
+
+Cusp areas
+----------
+
+>>> M = Manifold('o9_44210')
+>>> M.cusp_areas(policy='greedy') # doctest: +NUMERIC9
+[7.053940530873898, 3.2712450270, 2.7091590087]
+>>> M.cusp_areas(policy='greedy', first_cusps=[]) # doctest: +NUMERIC9
+[7.053940530873898, 3.2712450270, 2.7091590087]
+>>> M.cusp_areas(policy='greedy', first_cusps=[0,]) # doctest: +NUMERIC9
+[7.053940530873898, 3.2712450270, 2.7091590087]
+>>> M.cusp_areas(policy='greedy', first_cusps=[0,1]) # doctest: +NUMERIC9
+[7.053940530873898, 3.2712450270, 2.7091590087]
+>>> M.cusp_areas(policy='greedy', first_cusps=[0,1,2]) # doctest: +NUMERIC9
+[7.053940530873898, 3.2712450270, 2.7091590087]
+>>> M.cusp_areas(policy='greedy', first_cusps=[0,2,1]) # doctest: +NUMERIC9
+[7.053940530873898, 2.3513135103, 3.7690945490]
+>>> M.cusp_areas(policy='greedy', first_cusps=[1,]) # doctest: +NUMERIC9
+[2.30025338030798, 10.0315765558665, 0.883442685721903]
+
+Cusp translations
+-----------------
+
+>>> M = Manifold("s776")
+>>> M.cusp_translations(policy = 'greedy', first_cusps = [], bits_prec = 100) # doctest: +NUMERIC21
+[(0.70710678118654752440084436210 + 1.8708286933869706927918743662*I, 2.8284271247461900976033774484), (0.35355339059327376220042218105 + 0.93541434669348534639593718308*I, 1.4142135623730950488016887242), (0.35355339059327376220042218105 + 0.93541434669348534639593718308*I, 1.4142135623730950488016887242)]
+
+
+"""
+
+if not __doc__:
+    raise Exception("doc string with tests was not recognized.")
+
+from . import Manifold, ManifoldHP, Triangulation, TriangulationHP