你好,我正在使用g++编译器,并且遇到了(我认为)double类型的下溢问题,这是否可能发生,如果是,行为如何定义?
我已经上传了协方差矩阵的csv格式(51x51):http://pastebin.com/r0fx1qsx
这是我使用的计算行列式的代码(需要boost库,现已切换到long double类型,但没有影响):
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
正如您所看到的,其中一个是(稍微)积极的,而另一个是(稍微)消极的。