import mir.test: should;

0.geometricInvCDF(0.5).should == -1;
0.5.geometricInvCDF(0.5).should == 0;
0.75.geometricInvCDF(0.5).should == 1;
0.875.geometricInvCDF(0.5).should == 2;
0.95.geometricInvCDF(0.5).should == 4;

0.geometricInvCDF(0.25).should == -1;
0.25.geometricInvCDF(0.25).should == 0;
0.4375.geometricInvCDF(0.25).should == 1;
0.578125.geometricInvCDF(0.25).should == 2;
0.95.geometricInvCDF(0.25).should == 10;

0.5.geometricInvCDF(1).should == 0;
1.geometricInvCDF(0.5).should == double.infinity;