We can integrate arbitrary rational functions of the trigonometric functions because of the tangent half-angle substitution (https://en.m.wikipedia.org/wiki/Tangent_half-angle_substitution). This led me to the following question:
Is every integral of a rational function of Jacobian elliptic functions (https://dlmf.nist.gov/22) of a fixed modulus expressible as a finite composition of elementary functions, Jacobian elliptic functions, elliptic integrals and the Jacobi epsilon (https://dlmf.nist.gov/22.16#ii)? If so, does there exist an algorithm for it?