  The Bolza surface can be seen as the quotient of the hyperbolic plane,
represented by the Poincar\'e 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
$\mathbb{M}_g$, 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 $\mathbb{M}_g$.
We also propose algorithms computing small sets of points on $\mathbb{M}_g$
that are used to initialize Bowyer's algorithm.
