Update gmpy2 compatibility to 2.0.0
summations.py: return Piecewise for symbolic geometric sums in finite case (issue 3703)
Merge pull request #1941 from shikil/gmpy2

Update gmpy2 compatibility to 2.0.0
Add Stephen Loo to the AUTHORS/aboutus. Welcome to SymPy!
Merge pull request #1910 from skirpichev/3703-finite-geom-sum

summations.py: return Piecewise for symbolic geometric sums in finite case (issue 3703)
Don't do costly expand() on top of Order.__new__() method (issue 3654)

    Also, change compute_leading_term -> as_leading_term, as
    compute_leading_term(x) for x*(1 + x**2)**10000 calls expand() too.
    This seems to be related to Add._eval_nseries method.  Sometimes,
    it doesn't return order term (e.g. (1 + x**2).nseries(x)) - that's
    why Pow._eval_nseries() does actually expand_multinomial() for
    (1 + x**2)**1000.
Add _eval_as_leading_term() for polygamma

    fix test failure O(1/digamma(1/x)) == O(1/log(x))
    in test_order.py
Added test for issue 3654
Cleanup around expr.as_leading_term calculation