Compute pi to base 12 using Python?

[mensanator at aol.com]

Nick Craig-Wood wrote:
> mensanator at aol.com <mensanator at aol.com> wrote:
> >  I'm using GMPY (see code).
> [snip]
>
> If you are using gmpy you might as well do it like this.....
>
> gmpy.pi() uses the Brent-Salamin Arithmetic-Geometric Mean formula
for
> pi IIRC.  This converges quadratically, and it will calculate you a
> million places without breaking a sweat.

It would be nice if that were documented. What do I have to do, go get
the documentation for the original GMP to find out what else is in GMPY
that they didn't include in the doc file?


>
> >>> import gmpy
> >>> from math import log
> >>> bits = int(3003*log(12)/log(2))
> >>> pi=gmpy.pi(bits+100)
> >>> gmpy.fdigits(pi, 12, 3003)
>
'3.184809493b918664573a6211bb151551a05729290a7809a492742140a60a55256a0661a03753a3aa54805646880181a3683083272bbba0a370b12265529a828903b4b256b8403759a71626b8a54687621849b849a8225616b442796a31737b229b2391489853943b8763725616447236b027a421aa17a38b52a18a838b01514a51144a23315a3009a8906b61b8b48a62253a88a50a43ba0944572315933664476b3aabb77583975120683526b75b462060bb03b432551913772729a2147553531793848a0402b999b5058535374465a68806716644039539a8431935198527b9399b112990abb0383b107645424577a51601b3624a88b7a676a3992912121a213887b92873946a61332242217aa7354115357744939112602ba4b888818a3269222b528487747839994ab223b65b8762695422822669ba00a586097842a51750362073b5a768363b21bb1a97a4a194447749399804922175a068a46739461990a2065bb0a30bbab7024a585b1a84428195489784a07a331a7b0a1574565b373b05b03a5a80a13ab87857734679985558a5373178a7b28271992a3894a5776085083b9b238b2220542462888641a2bab8b3083ab49659172a312b78518654494a068662586a181835a64440b2970a122813975898815367208905801032881449223841428763329617531239b9a657405584014534390b587625606bb80923795944b43757a431b039556282978a6a49590553490ba1844947175637a908247b50127722464441380a852b0847b5813019bb70a67663b426565434069884476132193344ba55a2128a03838974606b851b2979321a408067225a5aa4b3464a1a17473595333909ab9127079655b3164b68b9b28a9b818a220a025ab0934203995b7a62a7aa739355340539ba3182905b193905603a43b660b9426a92294697144a896a5b2339358bb2b7294bb89635b071a6351211360b820b1882ab8433b54757b87a373284b1ba182a10326476b369a4a6365b58b8018994bb152556765475a704bb94b6b2a39458971a8b90512786b5029404818644323552916170b3abb7363496427b088b68725a68570040617949289077b278069a09b559324b8a66828b40549b0296065b2300330592569a7b76b92ba1293585b6a9b604567a0901362856373b4b56897946256b4172b1b50474351364749a33996a81ba8847347a8411b850b79a03018291672aa0945656a159aa6aa0a845531a592005b8a34366b882257107b190969a846474836a9800750778920ba797297a2791101b0685a86bb704b9baa17b055293679843b35215b0a8b1182b611953b080aa5431b219907a8448a81b1a9493245676b88013b470335240859594158621014216619553246570601967448b470174b9244892444817453865a4003b5aa7176451aab90681a949786154aa040477382ba69371041710b8728458a23979252b254236753a44a1900aa283536a227648812525743868b410a567794663359a6726a5286783328135114789b7645505b047848020a730a9557b206776aa56a19682744107901306b29008808619866b4911a05264b872a46b5376383932699531b449195640b62a63622830886247a47b3957169861239358041aa281333622aa15912b0a636047a489bb0726282a78b96671b27305a9652496b9b999011a7ba36898891665b1a6009058978850a21b01a158a1473b84a192b8672542a2a7056581995207a436a5b3ba2824637a3112abb57176468206a071200a327b3216425148100786502aa21236abb35499277670a1269730583403b1922a483856007301983989159bb688a58b602339806b63002a339a50b0ba533b84827793913081070a32595a101803a9a20234691b1a0b623274b69b0b44688195169461059543a252bb05208720ba13118266a872b26b9b584959b45179519534b221a335a2bb6971b3276b3a63a5b791723109b176529bb90651584279b7825712521b8269800738b07a62b14547884414451224092937165625696557a78799a82126613535a36b410309b759976119777b895801074b9779b9b5137538b2799951012273399bb818b722 at 0'
>
> --
> Nick Craig-Wood <nick at craig-wood.com> --
http://www.craig-wood.com/nick