lu_factorize损害赔偿数据（行列式）

Question

伙计。当我尝试计算行列式时，有时遇到矩阵有时的麻烦。 At this picture you can see the problem with the (2, 2) member. 那是非常罕见的情况。 我的功能：

    double_t Determinant(const ublas::matrix<double_t> & matrix)
    {
        // boost::numeric::ublas
        ublas::matrix<double_t> mLu(matrix);
        ublas::permutation_matrix<std::size_t> pivots(matrix.size1());

        std::cout << "BEFORE FACTORIZE" << std::endl;
        PrintMatrix(mLu); // just printing the matrix

        if (ublas::lu_factorize(mLu, pivots)) // there are some problems here
        {
            return 0.0;
        }

        std::cout << std::endl << "AFTER FACTORIZE" << std::endl;
        PrintMatrix(mLu); // just printing the matrix

        double_t det = 1.0;
        for (std::size_t i = 0; i < pivots.size(); ++i)
        {
            det *= (pivots(i) != i ? -1. : 1.) * mLu(i, i);
        }

        return det;
    }

该如何解决？谢谢。