The Bareiss algorithm is an efficient method used in numerical linear algebra for the computation of the determinant of a matrix. Developed by Hans Bareiss in the 1960s, this algorithm is particularly notable for its use of rational arithmetic, which helps in reducing numerical errors associated with floating-point computations.
Articles by others on the same topic
The Bareiss algorithm is fraction-free Gaussian elimination for exact determinants. For an integer matrix, put and . At step , choose a nonzero pivot element and update the trailing entries bySet the eliminated entries below the pivot to zero. The determinant identities behind the method make division by exact, so intermediate entries remain integers. Without row exchanges, the last pivot is ; row exchanges change its sign in the usual way. A zero pivot needs an exchange before continuing.
For , the first trailing block becomes . The last update is , giving the determinant without rational intermediate entries. Exact division is different from having no division at all.[ref]
Sylvester's Identity and Multistep Integer-Preserving Gaussian Elimination by Erwin H. Bareiss. Section I, pp. 565–567, equations (3), (6) and (8): bordered minors and exact-divisibility recurrence