A polynomial-time algorithm to compute Theta functions of Schottky groups

Contributed by Marc Masdeu ↗

Version 1 / Nov 12, 2026 / CC BY 4.0

Abstract

Abelian integrals, holomorphic differentials, period matrices and the Schottky--Klein prime function of a Riemann surface uniformized by a classical Schottky group are given by Poincar\'e theta series, which are sums over the elements of a free group. The algorithms that evaluate these series sum a truncated, possibly adaptively pruned, portion of the Cayley tree of the group, and the number of terms needed grows exponentially with the number of correct digits. We observe that the sum over all group elements starting with a fixed generator is a holomorphic function on the complement of a disc, and that these functions satisfy an exact linear functional equation coming from the self-similarity of the Cayley tree, given by a transfer operator built from the generators. Representing them by truncated Taylor expansions turns the evaluation of the series into a finite-dimensional linear problem, which we solve iteratively. The resulting algorithm computes the integrals Ωn\Omega_n, the period matrix and the prime function to PP bits of precision; for a fixed Schottky datum the number of arithmetic operations is cubic in PP and in the genus. We prove that it converges geometrically under a standard Schottky-type hypothesis on the group, which is closely related to the absolute convergence of the series. Along the way we give a complete proof of the cross-ratio formula for the period matrix, including a term that is missing from one version of it in the literature. We compare with the existing algorithms, and illustrate the method by reproducing a published value of the prime function and computing it to 3030 digits, and by computing period matrices of genus up to six to about 2525 digits.

Provenance statement

The initial idea and first prototype of the code is original from the author, who planned to develop an complex-analytic analogue to the p-adic construction of Masdeu--Xarles. This paper has been produced after some conversations with Claude's Opus 5.5 during the course of two days. All computations have been carried out on local servers.

Tools used

SageMath
SageMathVersion 10.6
Anthropic
Claude OpusVersion 5.5
Anthropic
Claude SonnetVersion 5.5
Anthropic
Claude HaikuVersion 5.5

References

  1. Morrison, Ralph; Ren, Qingchun. Algorithms for Mumford curves. J. Symbolic Comput., vol. 68, pp. 259–284. 2015DOI
  2. Gerritzen, Lothar; van der Put, Marius. Schottky groups and Mumford curves. Springer, Berlin, vol. 817, pp. viii+317. 1980
  3. The Sage Developers. SageMath, the Sage Mathematics Software System (Version 10.6). 2025link
  4. Fresnel, Jean; van der Put, Marius. Rigid analytic geometry and its applications. Birkhäuser Boston, Inc., Boston, MA, vol. 218, pp. xii+296. 2004DOI
  5. Amorós, Laia; Milione, Piermarco. Mumford curves covering pp-adic Shimura curves and their fundamental domains. Trans. Amer. Math. Soc., vol. 371, no. 2, pp. 1119–1149. 2019DOI
  6. Negrini, Isabella. On the computation of p-adic theta functions arising from the Hurwitz quaternions. 2017
  7. Guitart, Xavier; Masdeu, Marc; Xarles, Xavier. A quaternionic construction of pp-adic singular moduli. Res. Math. Sci., vol. 8, no. 3, pp. Paper No. 45, 20. 2021DOI
  8. Darmon, Henri; Vonk, Jan. Singular moduli for real quadratic fields: a rigid analytic approach. Duke Math. J., vol. 170, no. 1, pp. 23–93. 2021DOI
  9. Schmies, Markus. Computing Poincaré Theta Series for Schottky Groups. Computational Approach to Riemann Surfaces, pp. 165–182. 2011DOI
  10. Bogatyrev, A. B. Effective computation of Chebyshev polynomials for several intervals. Sb. Math., vol. 190, no. 11, pp. 1571–1605. 1999
Show all 60 referencesShow fewer references
  1. Bogatyrev, Andrei. Extremal Polynomials and Riemann Surfaces. Springer, Berlin. 2012
  2. Burnside, W. On a class of automorphic functions. Proc. London Math. Soc., vol. 23, pp. 49–88. 1891
  3. Schmies, Markus. Computational Methods for Riemann Surfaces and Helicoids with Handles. 2005
  4. Crowdy, D. G.; Marshall, J. S. Computing the Schottky–Klein prime function on the Schottky double of planar domains. Comput. Methods Funct. Theory, vol. 7, no. 1, pp. 293–308. 2007
  5. Lyamaev, S. Yu. Summation of Poincaré theta series in the Schottky model. Comput. Math. Math. Phys., vol. 62, no. 7, pp. 1059–1073. 2022DOI
  6. Lyamaev, S. Yu. On approximate summation of Poincaré series in the Schottky model. Dokl. Math., vol. 106, no. 1, pp. 247–250. 2022DOI
  7. Fairchild, Samantha; Ríos Ortiz, Ángel David. Crossing the transcendental divide: from Schottky groups to algebraic curves. 2024arXiv
  8. Deconinck, Bernard; Heil, Matthias; Bobenko, Alexander; van Hoeij, Mark; Schmies, Marcus. Computing Riemann theta functions. Math. Comp., vol. 73, no. 247, pp. 1417–1442. 2004
  9. Seppälä, M. Myrberg's numerical uniformization of hyperelliptic curves. Ann. Acad. Sci. Fenn. Math., vol. 29, pp. 3–20. 2004
  10. Masdeu, Marc; Xarles, Xavier. Efficient computation of non-archimedean theta functions. Math. Comp., vol. 95, no. 357. 2026arXivDOI
  11. Bobenko, Alexander I. Introduction to compact Riemann surfaces. Computational Approach to Riemann Surfaces, vol. 2013, pp. 3–64. 2011
  12. Crowdy, Darren; Marshall, Jonathan S. Uniformizing real hyperelliptic MM-curves using the Schottky–Klein prime function. Computational Approach to Riemann Surfaces, vol. 2013, pp. 183–193. 2011
  13. Hidalgo, Rubén A.; Seppälä, Mika. Numerical Schottky uniformizations: Myrberg's opening process. Computational Approach to Riemann Surfaces, vol. 2013, pp. 195–209. 2011
  14. Hidalgo, Rubén A. Numerical uniformization of hyperelliptic MM-symmetric Riemann surfaces. Proyecciones, vol. 20, no. 3, pp. 351–365. 2001
  15. Hidalgo, Rubén A.; Figueroa, J. Numerical Schottky uniformizations. Geom. Dedicata, vol. 111, pp. 125–157. 2005
  16. Bowen, Rufus. Hausdorff dimension of quasi-circles. Inst. Hautes Études Sci. Publ. Math., vol. 50, pp. 11–25. 1979
  17. Sullivan, Dennis. The density at infinity of a discrete group of hyperbolic motions. Inst. Hautes Études Sci. Publ. Math., vol. 50, pp. 171–202. 1979DOI
  18. McMullen, Curtis T. Hausdorff dimension and conformal dynamics. III. Computation of dimension. Amer. J. Math., vol. 120, no. 4, pp. 691–721. 1998DOI
  19. Pollicott, Mark; Vytnova, Polina. Hausdorff dimension estimates applied to Lagrange and Markov spectra, Zaremba theory, and limit sets of Fuchsian groups. Trans. Amer. Math. Soc. Ser. B, vol. 9, pp. 1102–1159. 2022
  20. Bandtlow, Oscar F.; Jenkinson, Oliver. Explicit eigenvalue estimates for transfer operators acting on spaces of holomorphic functions. Adv. Math., vol. 218, pp. 902–925. 2008
  21. Crowdy, Darren G.; Kropf, Everett H.; Green, Christopher C.; Nasser, Mohamed M. S. The Schottky–Klein prime function: a theoretical and computational tool for applications. IMA J. Appl. Math., vol. 81, no. 3, pp. 589–628. 2016DOI
  22. Crowdy, Darren. Solving Problems in Multiply Connected Domains. SIAM, Philadelphia. 2020
  23. Yamamoto, Shuntaro; Miyoshi, Hiroyuki. Fast computation and convergence analysis of the infinite-product representation of the Schottky–Klein prime function. 2026arXiv
  24. Baker, H. F. Abelian Functions: Abel's Theorem and the Allied Theory of Theta Functions. Cambridge University Press. 1897
  25. Hejhal, Dennis A. Theta Functions, Kernel Functions, and Abelian Integrals. Amer. Math. Soc., Providence, RI, vol. 129. 1972
  26. Bobenko, A. I. Schottky uniformization and finite-gap integration. Soviet Math. Dokl., vol. 36, no. 1, pp. 38–42. 1988
  27. Belokolos, E. D.; Bobenko, A. I.; Enolskii, V. Z.; Its, A. R.; Matveev, V. B. Algebro-Geometric Approach to Nonlinear Integrable Equations. Springer, Berlin. 1994
  28. Liu, Haoxuan; You, Zuhong; Yuan, Xiaoping. A parameterization of small amplitude KP finite gap solutions via classical Schottky uniformization and persistence of one and two gap solutions. 2026arXiv
  29. Baune, Konstantin; Broedel, Johannes; Im, Egor; Lisitsyn, Artyom; Zerbini, Federico. Schottky–Kronecker forms and hyperelliptic polylogarithms. J. Phys. A, vol. 57, pp. 445202. 2024arXiv
  30. Berger, Manuel; Broedel, Johannes. Translating auxiliary symmetries between Schottky uniformization and Jacobi parametrization. J. Phys. A, vol. 59, pp. 385206. 2026arXiv
  31. Bobenko, Alexander I.; Bobenko, Nikolai; Suris, Yuri B. Dimers and MM-curves. 2024arXiv
  32. Bobenko, Alexander I.; Bobenko, Nikolai. Dimers and MM-curves: limit shapes from Riemann surfaces. 2024arXiv
  33. Molin, Pascal; Neurohr, Christian. Computing period matrices and the Abel–Jacobi map of superelliptic curves. Math. Comp., vol. 88, pp. 847–888. 2019
  34. Kaya, Enis; Masdeu, Marc; Müller, J. Steffen; van der Put, Marius. Algorithms for hyperelliptic Mumford curves: pp-adic uniformization, pp-adic integrals and pp-adic heights. 2026arXiv
  35. Bogatyrev, A. B.; Goreinov, S. A.; Lyamaev, S. Yu. Analytical approach to multiband filter synthesis and comparison to other approaches. Probl. Inf. Transm., vol. 53, no. 3, pp. 260–273. 2017
  36. Bezrodnykh, S.; Bogatyrev, A.; Goreinov, S.; Grigoriev, O.; Hakula, H.; Vuorinen, M. On capacity computation for symmetric polygonal condensers. J. Comput. Appl. Math., vol. 361, pp. 271–282. 2019
  37. Bogatyrev, A. B.; Grigor'ev, O. A. Filtration under a stepped dam and Riemann theta functions. Proc. Steklov Inst. Math., vol. 311, pp. 10–21. 2020
  38. Bogatyrev, A. B. Real meromorphic differentials: a language for describing meron configurations in planar magnetic nanoelements. Theoret. and Math. Phys., vol. 193, no. 1, pp. 1547–1559. 2017
  39. Barré-Sirieix, K.; Diaz, G.; Gramain, F.; Philibert, G. Une preuve de la conjecture de Mahler–Manin. Invent. Math., vol. 124, pp. 1–9. 1996
  40. Hidalgo, Rubén A. Towards a proof of the classical Schottky uniformization conjecture. 2017. Revised 2025arXiv
  41. Shimura, Goro. On the real points of an arithmetic quotient of a bounded symmetric domain. Math. Ann., vol. 215, pp. 135–164. 1975DOI
  42. Čerednik, I. V. Uniformization of algebraic curves by discrete arithmetic subgroups of PGL2(kw){\rm PGL}_2(k_w) with compact quotient spaces. Math. USSR Sb., vol. 29, pp. 55–78. 1976DOI
  43. Drinfeld, V. G. Coverings of pp-adic symmetric domains. Funct. Anal. Appl., vol. 10, pp. 107–115. 1976DOI
  44. Manin, Yu. I.; Drinfeld, V. G. Periods of pp-adic Schottky groups. J. Reine Angew. Math., vol. 262/263, pp. 239–247. 1973
  45. Bruin, Nils; Sijsling, Jeroen; Zotine, Alexandre. Numerical computation of endomorphism rings. Open Book Series, vol. 2, pp. 155–171. 2019arXivDOI
  46. Mityushev, Vladimir. Convergence of the Poincaré series for some classical Schottky groups. Proc. Amer. Math. Soc., vol. 126, no. 8, pp. 2399–2406. 1998DOI
  47. Mityushev, Vladimir V.; Rogosin, Sergei V. Constructive Methods for Linear and Nonlinear Boundary Value Problems for Analytic Functions: Theory and Applications. Chapman & Hall/CRC, Boca Raton, FL, vol. 108. 2000
  48. Mityushev, Vladimir; Rylko, Natalia. A fast algorithm for computing the flux around non-overlapping disks on the plane. Math. Comput. Modelling, vol. 57, pp. 1350–1359. 2013DOI
  49. Deligne, P.; Rapoport, M. Les schémas de modules de courbes elliptiques. Modular Functions of One Variable, II, vol. 349, pp. 143–316. 1973DOI
  50. The PARI Group. PARI/GPlink

Version history

  1. v1 · 24 pagesSubmitted by Marc MasdeuInitial depositCurrentNov 05, 2026