Files:

[PDF] 262.1kB  [gzipped postscript] 330.4kB  

Abstract:

We present algorithmic, complexity and implementation results concerning real root isolation of integer univariate polynomials using the continued fraction expansion of real algebraic numbers. One motivation is to explain the method's good performance in practice. %% We derive an expected complexity bound of $\sOB(d^6 + d^4 \tau^2)$, where $d$ is the polynomial degree and $\tau$ bounds the coefficient bit size, using a standard bound on the expected bit size of the integers in the continued fraction expansion, thus matching the current worst-case complexity bound for real root isolation by exact methods (Sturm, Descartes and Bernstein subdivision). %% Moreover, using a homothetic transformation we improve the expected complexity bound to $\sOB( d^3 \tau)$. We compute the multiplicities within the same complexity and extend the algorithm to non square-free polynomials. %% Finally, we present an open-source \textttC++ implementation in the algebraic library \synaps, and illustrate its completeness and efficiency as compared to some other available software. For this we use polynomials with coefficient bit size up to 25000 bits and degree up to 2000.

BibTeX:
@Article{te-tcs-2008,
  author =       {Elias~P. Tsigaridas and Ioannis~Z. Emiris},
  title =        {{On the complexity of real root isolation using
                  Continued Fractions}},
  journal =      TCS,
  volume =       392,
  pages =        {158--173},
  year =         2008,
  editor =       {L. Bus{\'e} and M. Elkadi and B. Mourrain},
  abstract =     " We present algorithmic, complexity and
                  implementation results concerning real root
                  isolation of integer univariate polynomials using
                  the continued fraction expansion of real algebraic
                  numbers.  One motivation is to explain the method's
                  good performance in practice.  %% We derive an
                  expected complexity bound of $\sOB(d^6 + d^4
                  \tau^2)$, where $d$ is the polynomial degree and
                  $\tau$ bounds the coefficient bit size, using a
                  standard bound on the expected bit size of the
                  integers in the continued fraction expansion, thus
                  matching the current worst-case complexity bound for
                  real root isolation by exact methods (Sturm,
                  Descartes and Bernstein subdivision).  %% Moreover,
                  using a homothetic transformation we improve the
                  expected complexity bound to $\sOB( d^3 \tau)$.  We
                  compute the multiplicities within the same
                  complexity and extend the algorithm to non
                  square-free polynomials.  %% Finally, we present an
                  open-source \texttt{C++} implementation in the
                  algebraic library \synaps, and illustrate its
                  completeness and efficiency as compared to some
                  other available software.  For this we use
                  polynomials with coefficient bit size up to 25\,000
                  bits and degree up to 2\,000.",
}

Generated by bib2html.pl (written by Patrick Riley , modified by Elias ) on Wed Oct 23, 2019 21:41:02