8

不知从哪里我得到了这个函数的相当大的结果......它应该很简单，但我现在看不到它。

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

谢谢您的帮助！

4

2 回答 2

14

inhabited_columns是无符号-的，我在它之前看到一个一元： -inhabited_columns.

（请注意，一元-具有非常高的运算符优先级；高于*等）。

那就是你的问题所在！引用 Mike Seymour 的回答：

当你否定它时，结果仍然是无符号的；该值以 2 32为模减少以给出较大的正值。

一种解决方法是写

-(inhabited_columns * log(log_of_sum[j]))

因为那么否定将在浮点中进行

于 2013-10-01T12:51:54.857 回答
6

inhabited_columns是无符号类型。当你否定它时，结果仍然是无符号的；该值以 2 32为模减少以给出较大的正值。

您应该将其更改为足够大的有符号类型（也许int32_t，如果您不会有超过几十亿列），或者可能double因为您将要在双精度算术中使用它。

于 2013-10-01T12:55:18.823 回答