Bug: #4028 scipy.stats.geom.logpmf(1,1) returns 0.0. #4028 fixed.

BUG: #4028 fixed. Modified scipy.stats.geom.logpmf() to:
           def _logpmf(self, k, p):
                    return special.xlog1py(k - 1, -p) + log(p)