Bill Allombert on Wed, 15 Feb 2012 16:09:59 +0100 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Complex AGM |
Hello PARI developers, I noticed the complex AGM function in PARI do not always return the so-called 'optimal' AGM. However this is documented: ??agm ... In the case of complex or negative numbers, the principal square root is always chosen. ... So I am not sure whether this should be fixed. However, do we have a proof of that PARI agm never diverges to zero ? Example of bad case are: agm(1,-.00001) which gives -0.0403257332 - 0.00327376459*I. but the optimal AGM is 0.114955789 + 0.0279973724*I Cheers, Bill.