A Landen-type method for computation of Weierstrass functions

We establish a version of the Landen's transformation for Weierstrass functions and invariants that is applicable to general lattices in complex plane. Using it we present an effective method for computing Weierstrass functions, their periods, and elliptic integral in Weierstrass form given Weierstrass invariants $g_2$ and $g_3$ of an elliptic curve. Similarly to the classical Landen's method our algorithm has quadratic rate of convergence.