LMCAS::anonymous_namespace{symbolic_algebra_adapter.cpp} 命名空间参考

类型定义
using	ResultantMatrix = std::vector<std::vector<Rational>>

函数
void	eliminate_resultant_column (ResultantMatrix &matrix, int pivot, const Rational &previous)
Rational	resultant_determinant (ResultantMatrix &matrix)
ResultantMatrix	sylvester_matrix (const Polynomial< Rational > &left, const Polynomial< Rational > &right)

类型定义说明

ResultantMatrix

using LMCAS::anonymous_namespace{symbolic_algebra_adapter.cpp}::ResultantMatrix = std::vector<std::vector<Rational>>

在文件 symbolic_algebra_adapter.cpp 第 7 行定义.

函数说明

eliminate_resultant_column()

void LMCAS::anonymous_namespace{symbolic_algebra_adapter.cpp}::eliminate_resultant_column	(	ResultantMatrix &	matrix,
		int	pivot,
		const Rational &	previous )

在文件 symbolic_algebra_adapter.cpp 第 29 行定义.

                                                          {
    const int size = static_cast<int>(matrix.size());
    for (int row = pivot + 1; row < size; ++row) {
        for (int column = pivot + 1; column < size; ++column) {
            matrix[row][column] =
                (matrix[row][column] * matrix[pivot][pivot] -
                 matrix[row][pivot] * matrix[pivot][column]) / previous;
        }
        matrix[row][pivot] = Rational(0);
    }
}

引用了 eliminate_resultant_column().

被这些函数引用 eliminate_resultant_column() , 以及 resultant_determinant().

函数调用图:

这是这个函数的调用关系图:

resultant_determinant()

Rational LMCAS::anonymous_namespace{symbolic_algebra_adapter.cpp}::resultant_determinant

(

ResultantMatrix &

matrix

)

在文件 symbolic_algebra_adapter.cpp 第 42 行定义.

                                                        {
    Rational previous(1);
    int sign = 1;
    const int size = static_cast<int>(matrix.size());
    for (int pivot = 0; pivot < size; ++pivot) {
        if (matrix[pivot][pivot] == Rational(0)) {
            int replacement = pivot + 1;
            while (replacement < size &&
                   matrix[replacement][pivot] == Rational(0)) {
                ++replacement;
            }
            if (replacement == size) { return Rational(0); }
            std::swap(matrix[pivot], matrix[replacement]);
            sign = -sign;
        }
        eliminate_resultant_column(matrix, pivot, previous);
        previous = matrix[pivot][pivot];
    }
    auto determinant = matrix.back().back();
    if (sign < 0) { determinant = Rational(0) - determinant; }
    return determinant;
}

引用了 eliminate_resultant_column() , 以及 resultant_determinant().

被这些函数引用 resultant_determinant().

函数调用图:

这是这个函数的调用关系图:

sylvester_matrix()

ResultantMatrix LMCAS::anonymous_namespace{symbolic_algebra_adapter.cpp}::sylvester_matrix	(	const Polynomial< Rational > &	left,
		const Polynomial< Rational > &	right )

在文件 symbolic_algebra_adapter.cpp 第 9 行定义.

                                                                    {
    const int left_degree = left.degree();
    const int right_degree = right.degree();
    const int size = left_degree + right_degree;
    ResultantMatrix matrix(size, std::vector<Rational>(size, Rational(0)));
    for (int row = 0; row < right_degree; ++row) {
        for (int column = 0; column <= left_degree; ++column) {
            matrix[row][row + column] = left.coeffs[left_degree - column];
        }
    }
    for (int row = 0; row < left_degree; ++row) {
        for (int column = 0; column <= right_degree; ++column) {
            matrix[right_degree + row][row + column] =
                right.coeffs[right_degree - column];
        }
    }
    return matrix;
}

引用了 LMCAS::Polynomial< CoeffType >::coeffs, LMCAS::Polynomial< CoeffType >::degree(), LMCAS::right , 以及 sylvester_matrix().

被这些函数引用 sylvester_matrix().

函数调用图:

这是这个函数的调用关系图: