lu_factorize данные о повреждениях (определитель)

Question

Folks.Когда я пытаюсь вычислить определитель, у меня возникают проблемы с моей матрицей , иногда . На этом снимке вы видите проблему с элементом (2, 2). Это настолько редкий случай. Моя функция:

    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;
    }

Как это исправить?Благодарю.