我毫不费力地获得了此功能的丰硕成果……应该非常简单,但现在看不到它。
double prob_calculator_t::pimpl_t::B_full_term() const
{
double result = 0.0;
for (uint32_t j=0, j_end=U; j<j_end; j++)
{
uint32_t inhabited_columns = doc->row_sums[j];
// DEBUG
cout << "inhabited_columns: " << inhabited_columns << endl;
cout << "log_of_sum[j]: " << log_of_sum[j] << endl;
cout << "sum_of_log[j]: " << sum_of_log[j] << endl;
// end DEBUG
result += ( -inhabited_columns * log( log_of_sum[j] ) + sum_of_log[ j ] );
cout << "result: " << result << endl;
}
return result;
}
跟踪在哪里:
inhabited_columns: 1
log_of_sum[j]: 110.56
sum_of_log[j]: -2.81341
result: 2.02102e+10
inhabited_columns: 42
log_of_sum[j]: 110.56
sum_of_log[j]: -143.064
result: 4.04204e+10
谢谢您的帮助!
最佳答案
inhabited_columns
是未签名的,我在它前面看到一元-
:-inhabited_columns
。
(请注意,一元-
具有很高的运算符优先级;高于*
等)。
那就是你的问题所在!引用迈克·西摩的答案:
一种解决方法是写
-(inhabited_columns * log(log_of_sum[j]))
因为否定将在浮点数中进行
关于c++ - 可怕的数字来自无处,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/19116566/