GPU Accelerated Newton for Taylor Series Solutions of Polynomial Homotopies in Multiple Double Precision
Abstract
A polynomial homotopy is a family of polynomial systems, typically in one parameter $t$. Our problem is to compute power series expansions of the coordinates of the solutions in the parameter $t$, accurately, using multiple double arithmetic. One application of this problem is the location of the nearest singular solution in a polynomial homotopy, via the theorem of Fabry. Power series serve as input to construct Padé approximations. Exploiting the massive parallelism of Graphics Processing Units capable of performing several trillions floatingpoint operations per second, the objective is to compensate for the cost overhead caused by arithmetic with power series in multiple double precision. The application of Newton's method for this problem requires the evaluation and differentiation of polynomials, followed by solving a blocked lower triangular linear system. Experimental results are obtained on NVIDIA GPUs, in particular the RTX 2080, P100 and V100. Code generated by the CAMPARY software is used to obtain results in double double, quad double, and octo double precision. The programs in this study are self contained, available in a public github repository under the GPLv3.0 License.
 Publication:

arXiv eprints
 Pub Date:
 January 2023
 DOI:
 10.48550/arXiv.2301.12659
 arXiv:
 arXiv:2301.12659
 Bibcode:
 2023arXiv230112659V
 Keywords:

 Mathematics  Numerical Analysis;
 Computer Science  Distributed;
 Parallel;
 and Cluster Computing;
 Computer Science  Mathematical Software;
 Mathematics  Algebraic Geometry