We describe algorithms to compute elliptic functions and their relatives (Jacobi theta functions, modular forms, elliptic integrals, and the arithmeticgeometric mean) numerically to arbitrary precision with rigorous error bounds for arbitrary complex variables. Implementations in ball arithmetic are available in the open source Arb library. We discuss the algorithms from a concrete implementation point of view, with focus on performance at tens to thousands of digits of precision. 1 Available at http://arblib.org. The functionality for modular forms and elliptic functions can be found in the acb_modular (http://arblib.org/acb_modular.html) and acb_elliptic (http://arblib.org/acb_elliptic.html) modules. 2 Of course, for applications that do not require rigorous error bounds, all the algorithms can just as well be implemented in ordinary floating-point arithmetic.