def test_density(): h = StreamHist() for p in [1., 2., 2., 3.]: h.update(p) assert about(0.0, h.density(0.0), 1e-10) assert about(0.0, h.density(0.5), 1e-10) assert about(0.5, h.density(1.0), 1e-10) assert about(1.5, h.density(1.5), 1e-10) assert about(2.0, h.density(2.0), 1e-10) assert about(1.5, h.density(2.5), 1e-10) assert about(0.5, h.density(3.0), 1e-10) assert about(0.0, h.density(3.5), 1e-10) assert about(0.0, h.density(4.0), 1e-10)
def test_point_density_at_zero(): h = StreamHist().update(-1).update(0).update(1) assert h.density(0) == 1 h = StreamHist().update(0) assert h.density(0) == float("inf")