Delaunay triangulations of generalized Bolza surfaces

The Bolza surface can be seen as the quotient of the hyperbolic plane, represented by the Poincaré disk model, under the action of the group generated by the hyperbolic isometries identifying opposite sides of a regular octagon centered at the origin. We consider generalized Bolza surfaces Mg, where the octagon is replaced by a regular 4g-gon, leading to a genus g surface. We propose an extension of Bowyer’s algorithm to these surfaces. In particular, we compute the value of the systole of Mg. We also propose algorithms computing small sets of points on Mg that are used to initialize Bowyer’s algorithm.