Computing Riemann theta functions

Abstract

The Riemann theta function is a complex-valued function of $g$ complex variables. It appears in the construction of many (quasi-)periodic solutions of various equations of mathematical physics. In this paper, algorithms for its computation are given. First, a formula is derived allowing the pointwise approximation of Riemann theta functions, with arbitrary, user-specified precision. This formula is used to construct a uniform approximation formula, again with arbitrary precision.