whoops

test/truncated/lognormal.jl

@@ -27,7 +27,7 @@ naive_moment(d, n, μ, σ²) = (σ = sqrt(σ²); expectation(x -> ((x - μ) / σ
         # Without a `max(_, 0)`, this would be within machine precision of 0 (as above) but
         # numerically negative, which could cause downstream issues that assume a nonnegative
         # variance
-        @test v1 > 0
+        @test v1 >= 0
         # Compare results with not specifying a lower bound at all
         d2 = truncated(LogNormal(1, 5); upper=1e5)
         @test mean(d1) == mean(d2)