It's worth noting that for practical clock pendulum amplitudes, say up to 10 degrees, the difference between the AGM formula and the usual "A^2/16" approximation to the circular deviation is very very small: 3.3 usec for a 1 second period and 10 degrees. I've done pendulum simulations for the order of amplitude used in "Clock B" and compared them to simple formulas based on the same approximation with virtually identical results.
Where the AGM really is useful is for methods where you need to quickly compute elliptic integrals. One example is designing very selective LC filters, where there are formulas for the component values that incorporate elliptic integrals. Once you needed to have extensive tables of these, but now they can be calculated more accurately in a few lines of code.