3

嗨，我正在使用 g++ 编译器并且正在经历（我认为）是双打下溢，这可能吗？如果是的话，行为是如何定义的

我在这里上传了协方差矩阵 (51x51) 的 csv 格式：http: //pastebin.com/r0fx1qsx

这是我用来计算行列式的代码（在 c++ 中需要提升）（我已经切换到长双精度并且没有效果）：

int determinant_sign(const boost::numeric::ublas::permutation_matrix<std ::size_t>& pm)
{
    int pm_sign=1;
    std::size_t size = pm.size();
    for (std::size_t i = 0; i < size; ++i)
        if (i != pm(i))
            pm_sign *= -1.0; // swap_rows would swap a pair of rows here, so we change sign
    return pm_sign;
}

long double determinant( boost::numeric::ublas::matrix<long double>& m ) {
    boost::numeric::ublas::permutation_matrix<std ::size_t> pm(m.size1());
    long double det = 1.0;

    if( boost::numeric::ublas::lu_factorize(m,pm) ) {
        det = 0.0;
    } else {
        for(int i = 0; i < (int)m.size1(); i++)
            det *= m(i,i); // multiply by elements on diagonal
        det = det * determinant_sign( pm );
    }
    return det;
}

我为数据给出的结果是-3.59916e-183。

当我运行以下matlab代码时：

M = csvread('path/to/csv');
det(M)
the result I get is:

4.2014e-173

如您所见，一个是（略微）积极的，而一个是（略微）消极的

4

1 回答 1

4

假设浮点单元正确运行，它不会溢出到负值 - 如果值超出有效范围，结果将是“+INF”（正无穷大）。这只会发生在有符号整数上。

当然，完全有可能在计算中出现各种错误，在预期得到肯定答案时给出否定答案。

于 2013-02-19T15:51:42.437 回答