Bug#4028

BUG: fixes #4028 scipy.stats.geom.logpmf(1,1) returns 0.0