agm(1,-1) = 0

This commit is contained in:
Fredrik Johansson 2016-03-08 17:19:54 +01:00
parent b9aab31961
commit 69564af91c

View file

@ -351,6 +351,10 @@ acb_agm1(acb_t m, const acb_t z, slong prec)
{
acb_agm1_basecase(m, z, prec);
}
else if (acb_equal_si(z, -1))
{
acb_zero(m);
}
else
{
/* use M(z) = (z+1)/2 * M(2 sqrt(z) / (z+1)) */