Computation of $\wp$-functions on plane algebraic curves
Julia Bernatska
TL;DR
This work develops a rigorous, analytically grounded framework for computing Kleinian $\wp$-functions on plane algebraic curves of low gonality by constructing explicit Riemann-surface models. It combines radical solutions for quadratic, cubic, and quartic equations with carefully designed Abel maps, theta/sigma functions, and Bolza-type relations to obtain direct Abel images and periods, which are then used to evaluate $\wp$-functions and their higher derivatives. The authors provide detailed hyperelliptic and trigonal case studies, including explicit period matrices, $\varkappa$-matrices, and Jacobi-inversion verifications, and they illustrate how to determine the vector of Riemann constants and build paths to points on the curves. The approach enables direct computation of $\wp$-functions without solely relying on theta-function representations, with robust numerical verification against Jacobi inversion and Bolza-type formulas. This has potential impact for applications in integrable systems (e.g., KdV, NLS) where explicit, numerically stable evaluations of multiply periodic functions are required.
Abstract
Numerical tools for computation of $\wp$-functions, also known as Kleinian, or multiply periodic, are proposed. In this connection, computation of periods of the both first and second kinds is reconsidered. An analytical approach to constructing Riemann surfaces of plane algebraic curves of low gonalities is used. The approach is based on explicit radical solutions to quadratic, cubic, and quartic equations, which serve for hyperelliptic, trigonal, and tetragonal curves, respectively. The proposed analytical models of Riemann surfaces give full control over computation of the Abel image of any point or divisor. Therefore, computation of $\wp$-functions at Abel images of given divisors can be done directly. An alternative computation with the help of the Jacobi inversion problem is used for verification. Hyperelliptic and trigonal curves are considered in detail, and illustrated by examples. A method of finding the unique characteristic corresponding to the vector of Riemann constants is suggested for non-hyperelliptic and hyperelliptic curves.
