import mir.math.common: exp;
import mir.test: shouldApprox;

0.geometricLPMF(0.5).exp.shouldApprox == 0.geometricPMF(0.5);
1.geometricLPMF(0.5).exp.shouldApprox == 1.geometricPMF(0.5);
2.geometricLPMF(0.25).exp.shouldApprox == 2.geometricPMF(0.25);