julia - Julia：BigFloat 正态分布

Question

在 Julia 中，是否有人实现了正态分布 pdf 或 cdf 来支持任意精度的 BigFloats。

例如，此代码返回 0.0，而实际上这些值应该略有不同。

x = parse(BigFloat, "2.1")
x_small = float64(x)
pdf(Normal(), x) - pdf(Normal(), x_small)

Answer 1

不是直接的。我最终的计划是使分布类型参数化，这也将允许Float32参数，但这还需要一段时间。

与此同时，有一个未导出的φ，它给出了你想要的结果：

Distributions.φ(x) - pdf(Normal(), x_small)