Searched refs:matrix_to_invert (Results 1 – 1 of 1) sorted by relevance
105 Matrix<double> GaussianElimination(Matrix<double> matrix_to_invert, in GaussianElimination() argument108 const size_t n = matrix_to_invert.size(); in GaussianElimination()110 for (const std::valarray<double>& column : matrix_to_invert) in GaussianElimination()119 matrix_to_invert = Transpose(matrix_to_invert); in GaussianElimination()130 if (std::abs(matrix_to_invert[column][diagonal_index]) > in GaussianElimination()131 std::abs(matrix_to_invert[diagonal_index][diagonal_index])) { in GaussianElimination()132 std::swap(matrix_to_invert[column], matrix_to_invert[diagonal_index]); in GaussianElimination()140 if (matrix_to_invert[diagonal_index][diagonal_index] == 0.0) { in GaussianElimination()145 matrix_to_invert[diagonal_index][diagonal_index]; in GaussianElimination()146 matrix_to_invert[diagonal_index] /= diagonal_element; in GaussianElimination()[all …]