On the numerical evaluation of algebro-geometric solutions to integrable equations

Physically meaningful periodic solutions to certain integrable partial differential equations are given in terms of multi-dimensional theta functions associated with real Riemann surfaces. Typical analytical problems in the numerical evaluation of these solutions are studied. In the case of hyperelliptic surfaces efficient algorithms exist even for almost degenerate surfaces. This allows the numerical study of solitonic limits. For general real Riemann surfaces, the choice of a homology basis adapted to the anti-holomorphic involution is important for a convenient formulation of the solutions and smoothness conditions. Since existing algorithms for algebraic curves produce a homology basis not related to automorphisms of the curve, we study symplectic transformations to an adapted basis and give explicit formulae for M-curves. As examples we discuss solutions of the Davey–Stewartson and the multi-component nonlinear Schrödinger equations.