gpt4 book ai didi

c++ - 此方法中堆栈溢出的原因( float 学)

转载 作者:可可西里 更新时间:2023-11-01 16:52:12 27 4
gpt4 key购买 nike

我偶尔会在这个方法中遇到 stackoverflow 异常。

double norm_cdf(const double x) {
double k = 1.0/(1.0 + 0.2316419*x);
double k_sum = k*(0.319381530 + k*(-0.356563782 + k*(1.781477937 + k*(-1.821255978 + 1.330274429*k))));

if (x >= 0.0) {
return (1.0 - (1.0/(pow(2*M_PI,0.5)))*exp(-0.5*x*x) * k_sum);
} else {
return 1.0 - norm_cdf(-x);
}
}

关于我为什么会得到它有什么建议吗?我可以采取任何措施来纠正错误吗?

最佳答案

您的问题是 x 不是数字。 NAN >= 0.0 为假,-NAN >= 0.0 也为假。

您可以像其他人建议的那样专门检查 NAN,但我建议简化事情:

static double norm_cdf_positive(const double x) {
double k = 1.0/(1.0 + 0.2316419*x);
double k_sum = k*(0.319381530 + k*(-0.356563782 + k*(1.781477937 + k*(-1.821255978 + 1.330274429*k))));

return (1.0 - (1.0/(pow(2*M_PI,0.5)))*exp(-0.5*x*x) * k_sum);
}

double norm_cdf(const double x) {
if (x >= 0.0) {
return norm_cdf_positive(x);
} else {
return 1.0 - norm_cdf_positive(-x);
}
}

这样做的好处是编译器可以对其行为做出更明智的假设。请注意,我已将“内部”函数标记为静态(这会将其范围限制在当前编译单元)。您还可以使用未命名的 namespace 。 (编辑:实际上 Timothy Shields 有一种更简单的方法来删除递归,它将所有内容保存在一个函数中)

关于c++ - 此方法中堆栈溢出的原因( float 学),我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/18154909/

27 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com