Implement Iso Delaunay Tessellation

flatsurf/__init__.py

@@ -25,6 +25,8 @@
 
 from flatsurf.geometry.hyperbolic import HyperbolicPlane
 
+from .geometry.iso_delaunay_tessellation import IsoDelaunayTessellation
+
 from flatsurf.geometry.surface_legacy import (
     Surface_list,
     Surface_dict,

flatsurf/geometry/categories/translation_surfaces.py

@@ -43,6 +43,8 @@
 #  along with sage-flatsurf. If not, see <https://www.gnu.org/licenses/>.
 # ####################################################################
 
+from sage.misc.cachefunc import cached_method
+
 from flatsurf.geometry.categories.surface_category import SurfaceCategoryWithAxiom
 from flatsurf.geometry.categories.half_translation_surfaces import (
     HalfTranslationSurfaces,
@@ -649,6 +651,15 @@ def stratum(self):
 
                     return AbelianStratum([ZZ(a - 1) for a in self.angles()])
 
+                @cached_method
+                def veech_group(self):
+                    from flatsurf.geometry.iso_delaunay_tessellation import IsoDelaunayTessellation
+                    from flatsurf.geometry.veech_group import VeechGroup
+
+                    idt = IsoDelaunayTessellation(self)
+                    idt.explore()
+                    return VeechGroup(idt)
+
                 def canonicalize(self, in_place=None):
                     r"""
                     Return a canonical version of this translation surface.

flatsurf/geometry/hyperbolic.py

@@ -211,6 +211,8 @@
 
 import collections.abc
 
+from flatsurf.geometry.surface import Surface_base
+
 from sage.structure.sage_object import SageObject
 from sage.structure.parent import Parent
 from sage.structure.element import Element