We present a new package
Theta.jl for computing the Riemann theta function. It is
implemented in Julia and offers accurate numerical evaluation of theta functions with
characteristics and their derivatives of arbitrary order. Our package is optimized for
multiple evaluations of theta functions for the same Riemann matrix, in small
dimensions. As an application, we report on experimental approaches to the Schottky
problem in genus 5.