Modified Gaussian Elimination without Division Operations

Abstract

A new modified method based on the Gaussian elimination method for solution of linear system of equations in the projective space is formulated. It is based on application of projective extension of the Euclidean space and use of homogeneous coordinates. It leads to an elimination of division operation and higher precision due to division operation elimination. The approach is based on understanding that a solution of the linear system is equivalent to the extended cross-product, i.e. . As it can be seen there no division is needed. Use of the projective representation enables to avoid division operation and use advantages of the matrix-vector architectures. Division operations have to be used only if the final result of computation has to be in the Euclidean representation. The proposed method was implemented in C# and C++ and experimentally verified. It is especially convenient for computations on GPUs based architectures.