Symmetric and non-symmetric Abaffy, Broyden and Spedicato (ABS) methods for solving Diophantine systems of equations

In the present paper we describe a new class of algorithms for solving Diophantine systems of equations in integer arithmetic. This algorithm, designated as the integer ABS (iABS) algorithm, is based on the ABS methods in the real space, with extensive modifications to ensure that all calculations remain in the integer space. Importantly, the iABS solves Diophantine systems of equations without determining the Hermite normal form. The algorithm is suitable for solving determined, over- or underdetermined, full rank or rank deficient linear integer equations. We also present a scaled integer ABS system and two special cases for solving general Diophantine systems of equations. In the scaled symmetric iABS (ssiABS), the Abaffian matrix Hi is symmetric, allowing that only half of its elements need to be calculated and stored. The scaled non-symmetric iABS system (snsiABS) provides more freedom in selecting the arbitrary parameters and thus the maximal values of Hi can be maintained at a certainly lower level. In addition to the above theoretical results, we also provide numerical experiments to test the performance of the ssiABS and the snsiABS algorithms. These experiments have confirmed the suitability of the iABS system for practical applications.