ABSTRACT
An iterative refinement approach is taken to rational linear system solving. Such methods produce, for each entry of the solution vector, a rational approximation with denominator a power of 2. From this the correct rational entry can be reconstructed. Our iteration is a numeric-symbolic hybrid in that it uses an approximate numeric solver at each step together with a symbolic (exact arithmetic) residual computation and symbolic rational reconstruction. The rational solution may be checked symbolically (exactly). However, there is some possibility of failure of convergence, usually due to numeric ill-conditioning. Alternatively, the algorithm may be used to obtain an extended precision floating point approximation of any specified precision. In this case we cannot guarantee the result by rational reconstruction and an exact solution check, but the approach gives evidence (not proof) that the probability of error is extremely small. The chief contributions of the method and implementation are (1) confirmed continuation, (2) improved rational reconstruction, and (3) faster and more robust performance.
- Bernhard Beckermann. The condition number of real Vandermonde, Krylov and positive definite Hankel matrices. Numerische Mathematik, 85:553--577, 1997.Google ScholarCross Ref
- C. Bright and A. Storjohan. Vector rational number reconstruction. In ISSAC '11. ACM, 2011. Google ScholarDigital Library
- John Burkardt. TEST MAT Test Matrices. http://people.sc.fsu.edu/~jburkardt/ c_src/test_mat/test_mat.html.Google Scholar
- Stanley Cabay. Exact solution of linear equations. In Proceedings of the second ACM symposium on Symbolic and algebraic manipulation, SYMSAC '71, pages 392--398, New York, NY, USA, 1971. ACM. Google ScholarDigital Library
- Z. Chen and A. Storjohann. A BLAS based C library for exact linear algebra on integer matrices. In Proc. of ISSAC'05, pages 92--99. ACM Press, 2005. Google ScholarDigital Library
- J. D. Dixon. Exact solution of linear equations using p-adic expansion. Numer. Math., pages 137--141, 1982.Google Scholar
- J-G. Dumas, T. Gautier, M. Giesbrecht, P. Giorgi, B. Hovinen, E. Kaltofen, B. D. Saunders, W. Turner, and G. Villard. Linbox: A generic library for exact linear algebra. In ICMS'02, pages 40--50, 2002.Google ScholarCross Ref
- W. Eberly, M. Giesbrecht, P. Giorgi, A. Storjohann, and G. Villard. Solving sparse rational linear systems. In Proc. of ISSAC'06, pages 63--70. ACM Press, 2006. Google ScholarDigital Library
- Joachim Von Zur Gathen and Jurgen Gerhard. Modern Computer Algebra. Cambridge University Press, New York, NY, USA, 2 edition, 2003. Google ScholarDigital Library
- Keith O. Geddes and Wei Wei Zheng. Exploiting fast hardware floating point in high precision computation. In J. Rafael Sendra, editor, ISSAC, pages 111--118. ACM, 2003. Google ScholarDigital Library
- M. J. C. Gover and S. Barnett. Brownian matrices: properties and extensions. International Journal of Systems Science, 17(2):381--386, 1986.Google ScholarCross Ref
- Andrzej Kiełbasinski. Iterative refinement for linear systems in variable-precision arithmetic. BIT, 21(1):97--103, 1981.Google ScholarCross Ref
- E. Kilic and P. Stanica. The Lehmer matrix and its recursive analogue. Journal of Combinatorial Mathematics and Combinatorial Computing, 74(2):195--205, 2010.Google Scholar
- D. H. Lehmer. Solutions to problem E710, proposed by D. H. Lehmer: The inverse of a matrix, November 1946.Google Scholar
- Robert T. Moenck and John H. Carter. Approximate algorithms to derive exact solutions to systems of linear equations. In Proceedings of the International Symposium on on Symbolic and Algebraic Computation, pages 65--73, London, UK, 1979. Springer-Verlag. Google ScholarDigital Library
- Michael B. Monagan. Maximal quotient rational reconstruction: an almost optimal algorithm for rational reconstruction. In Jaime Gutierrez, editor, ISSAC, pages 243--249. ACM, 2004. Google ScholarDigital Library
- Teo Mora, editor. Symbolic and Algebraic Computation, International Symposium ISSAC 2002, Lille, France, July 7--10, 2002, Proceedings. ACM, 2002.Google Scholar
- T. Mulders and A. Storjohann. Certified dense linear system solving. Jounal of symbolic computation, 37(4), 2004.Google Scholar
- Victor Y. Pan and Xinmao Wang. Acceleration of euclidean algorithm and extensions. In Mora {17}, pages 207--213. Google Scholar
- Alicja Smoktunowicz and Jolanta Sokolnicka. Binary cascades iterative refinement in doubled-mantissa arithmetics. BIT, 24(1):123--127, 1984.Google ScholarCross Ref
- Daniel Steffy. Exact solutions to linear systems of equations using output sensitive lifting. ACM Communications in Computer Algebra, 44(4):160--182, 2010. Google ScholarDigital Library
- The LinBox Team. LinBox, a C++ library for exact linear algebra. http://www.linalg.org/.Google Scholar
- John Todd. Basic Numerical Mathematics, Vol. 2: Numerical Algebra. Birkhauser, Basel, and Academic Press, New York, 1977. Google ScholarDigital Library
- Silvio Ursic and Cyro Patarra. Exact solution of systems of linear equations with iterative methods. SIAM Journal on Algebraic and Discrete Methods, 4(1):111--115, 1983.Google ScholarCross Ref
- Mark van Hoeij and Michael B. Monagan. A modular gcd algorithm over number fields presented with multiple extensions. In Mora {17}, pages 109--116. Google Scholar
- Zhengdong Wan. Computing the Smith Forms of Integer Matrices and Solving Related Problems. PhD thesis, University of Delaware, Newark, DE, 2005. Google ScholarDigital Library
- Zhengdong Wan. An algorithm to solve integer linear systems exactly using numerical methods. Journal of Symbolic Computation, 41:621--632, 2006. Google ScholarDigital Library
- James H. Wilkinson. Rounding Errors in Algebraic Processes. Dover Publications, Incorporated, 1994. Google ScholarDigital Library
Index Terms
- Numeric-symbolic exact rational linear system solver
Recommendations
Linear Rational Finite Differences from Derivatives of Barycentric Rational Interpolants
Derivatives of polynomial interpolants lead in a natural way to approximations of derivatives of the interpolated function, e.g., through finite differences. We extend a study of the approximation of derivatives of linear barycentric rational ...
Rational Pythagorean-hodograph space curves
A method for constructing rational Pythagorean-hodograph (PH) curves in R^3 is proposed, based on prescribing a field of rational unit tangent vectors. This tangent field, together with its first derivative, defines the orientation of the curve ...
Comments