在朱莉娅,有人实现了正态分布pdf或cdf支持任意精度的BigFloats。
例如,这段代码返回0.0,而实际上值应该略有不同。
x = parse(BigFloat, "2.1")
x_small = float64(x)
pdf(Normal(), x) - pdf(Normal(), x_small)发布于 2015-05-26 10:49:31
不是直接的。我的最终计划是使分发类型参数化,这也允许使用Float32参数,但这还需要一段时间。
同时,还有未导出的φ,它提供了所需的结果:
Distributions.φ(x) - pdf(Normal(), x_small)https://stackoverflow.com/questions/30452847
复制相似问题