Computing Real Roots of Real Polynomials
|
|
- Oliver Bridges
- 5 years ago
- Views:
Transcription
1 Computing Real Roots of Real Polynomials and now For Real! Alexander Kobel Max-Planck-Institute for Informatics, Saarbrücken, Germany Fabrice Rouillier INRIA & Université Pierre et Marie Curie, Paris, France Michael Sagraloff Max-Planck-Institute for Informatics, Saarbrücken, Germany
2 Computing Real Roots A classical approach
3 Computing Real Roots of Real Polynomials Given a square-free polynomial P(x) = p n x n + p n 1 x n p 1 x + p 0 R[x] and an open interval I = (a 0, b 0 ) R, find isolating intervals I 1,..., I r I such that each I j contains exactly one real root of P, the I j are pairwise disjoint, and j I j covers all real roots of P in I. 1
4 Subdivision approach Queue {I} Isol while Queue is not empty: pop I = (a, b) from Queue if I contains no root: (exclusion predicate) discard I else if I contains exactly one root: (inclusion-and-isolation predicate) add I to Isol otherwise (or we-don t-know): split I at m = a+b (or some other point) if P(m) = 0: add {m} to Isol add (a, m) and (m, b) to Queue
5 Subdivision approach Queue {I} Isol while Queue is not empty: pop I = (a, b) from Queue if I contains no root: (exclusion predicate) discard I else if I contains exactly one root: (inclusion-and-isolation predicate) add I to Isol otherwise (or we-don t-know): split I at m = a+b (or some other point) if P(m) = 0: add {m} to Isol add (a, m) and (m, b) to Queue
6 Descartes method Descartes Rule of Signs r = # positive real roots of P v = # sign changes in coeffs of P r v and r v (mod ) 3
7 Descartes method Descartes Rule of Signs r = # positive real roots of P v = # sign changes in coeffs of P r v and r v (mod ) localized version r I = # real roots of P in I = (a, b) v I = # sign changes of (x + 1) n P ( ) ax+b x+1 r I v I and r I v I (mod ) 3
8 Descartes method Descartes Rule of Signs r = # positive real roots of P v = # sign changes in coeffs of P r v and r v (mod ) localized version r I = # real roots of P in I = (a, b) v I = # sign changes of (x + 1) n P ( ) ax+b x+1 r I v I and r I v I (mod ) exclusion predicate v I = 0 I contains no root of P inclusion-and-isolation predicate v I = 1 I contains exactly one root of P 3
9 Descartes method Descartes Rule of Signs r = # positive real roots of P v = # sign changes in coeffs of P r v and r v (mod ) localized version r I = # real roots of P in I = (a, b) v I = # sign changes of (x + 1) n P ( ) ax+b x+1 r I v I and r I v I (mod ) exclusion predicate v I = 0 I contains no root of P inclusion-and-isolation predicate v I = 1 I contains exactly one root of P r I = v I {0, 1} if other (complex) roots well-separated from I 3
10 Descartes method v = up to O(nτ) steps
11 Descartes method v = up to O(nτ) steps
12 Descartes method v = up to O(nτ) steps
13 Descartes method v = up to O(nτ) steps
14 Descartes method v = up to O(nτ) steps
15 Descartes method v = up to O(nτ) steps
16 Descartes method v = up to O(nτ) steps
17 Descartes method v = up to O(nτ) steps
18 Descartes method v = up to O(nτ) steps
19 Descartes method v = up to O(nτ) steps
20 Descartes method v = up to O(nτ) steps
21 Complexity ( benchmark problem ) isolate all roots of an integer polynomial of degree n and coefficient bitsize τ # roots / tree width: n min. root separation: Õ(nτ) 5
22 Complexity ( benchmark problem ) isolate all roots of an integer polynomial of degree n and coefficient bitsize τ # roots / tree width: n min. root separation: Õ(nτ) subdivision rule and subdivision precision overall bit model of computation tree size demand complexity bisection, exact over Z / Q Õ(nτ) 1 Õ(n τ) Õ(n 4 τ ) 1 1 [Eigenwillig, Sharma, Yap: ISSAC 006] [Collins: JSC 016] 5
23 of Real Polynomials Deficiencies and remedies
24 Beyond the RealRAM-model P(x) = πx (π + 4 )x + = (πx )(x 1) x x R[x]? v =?? P( 1 ) = 0? 1 π 1 6
25 Beyond the RealRAM-model P(x) = πx (π + 4 )x + = (πx )(x 1) 6. x 9. x + 3. R[x]? v =?? P( 1 ) = 0? 1 π 1 6
26 Beyond the RealRAM-model P(x) = πx (π + 4 )x + = (πx )(x 1) 6.3 x 8.8 x +.8 R[x]? v =?? P( 1 ) = 0? 1 π 1 6
27 Beyond the RealRAM-model P(x) = πx (π + 4 )x + = (πx )(x 1) 6.8 x 8.79 x +.83 R[x]? v =?? P( 1 ) = 0? 1 π 1 6
28 Beyond the RealRAM-model P(x) = πx (π + 4 )x + = (πx )(x 1) 6.83 x x +.88 R[x]? v =?? P( 1 ) = 0? 1 π 1 6
29 Beyond the RealRAM-model P(x) = πx (π + 4 )x + = (πx )(x 1) 6.83 x x R[x]? v =?? P( 1 ) = 0? 1 π 1 6
30 Beyond the RealRAM-model P(x) = πx (π + 4 )x + = (πx )(x 1) 6.83 x x R[x]? v =?? P( 1 ) = 0? 1 π 1 problems with arbitrarily approximable ( bitstream ) coefficients no termination on subdivision on an exact root precision demand for processing I = (a, b) depends on P(a) and P(b) 6
31 Admissible points ADsc (Approximate Descartes) Instead of subdivision at m I: sample suitable points near m, choose an admissible point m with large value P(m ) among the samples and subdivide at m. 7
32 Admissible points ADsc (Approximate Descartes) Instead of subdivision at m I: sample suitable points near m, choose an admissible point m with large value P(m ) among the samples and subdivide at m. In particular: P(m ) = 0! guaranteed termination for arbitrary (square-free) inputs keeps precision demand near theoretical optimum 7
33 Complexity ( benchmark problem ) isolate all roots of an integer polynomial of degree n and coefficient bitsize τ # roots / tree width: n min. root separation: Õ(nτ) subdivision rule and subdivision precision overall bit model of computation tree size demand complexity bisection, exact over Z / Q Õ(nτ) Õ(n τ) Õ(n 4 τ ) bisection, approximate Õ(nτ) Õ(nτ) 1 Õ(n 3 τ ) 1 1 [Sagraloff: JSC 014] 8
34 Clustered roots v = up to O(nτ) steps
35 Clustered roots v = up to O(nτ) steps
36 Clustered roots v = up to O(nτ) steps
37 Combining Newton s and Descartes method inspired by (Approximate) Quadratic Interval Refinement and the Brent-Dekker method 3 [Abbott: ISSAC 006; ACM CCA 014], [Kerber, Sagraloff: JCAM 015] 3 [Brent: 1973] 10
38 Combining Newton s and Descartes method inspired by (Approximate) Quadratic Interval Refinement and the Brent-Dekker method 3 certified trial-and-error approach: multiplicity-agnostic variant of Newton iteration determine candidate subinterval of I containing a cluster verify using Descartes tests; on failure, resort to bisection [Abbott: ISSAC 006; ACM CCA 014], [Kerber, Sagraloff: JCAM 015] 3 [Brent: 1973] 10
39 Combining Newton s and Descartes method inspired by (Approximate) Quadratic Interval Refinement and the Brent-Dekker method 3 certified trial-and-error approach: multiplicity-agnostic variant of Newton iteration determine candidate subinterval of I containing a cluster verify using Descartes tests; on failure, resort to bisection successful in almost all steps; reduces chains of subdivisions around a cluster to poly-logarithmic length [Abbott: ISSAC 006; ACM CCA 014], [Kerber, Sagraloff: JCAM 015] 3 [Brent: 1973] 10
40 Combining Newton s and Descartes method inspired by (Approximate) Quadratic Interval Refinement and the Brent-Dekker method 3 certified trial-and-error approach: multiplicity-agnostic variant of Newton iteration determine candidate subinterval of I containing a cluster verify using Descartes tests; on failure, resort to bisection successful in almost all steps; reduces chains of subdivisions around a cluster to poly-logarithmic length similar techniques in complex domain: upcoming talk by Juan Xu [Abbott: ISSAC 006; ACM CCA 014], [Kerber, Sagraloff: JCAM 015] 3 [Brent: 1973] 10
41 Complexity ( benchmark problem ) isolate all roots of an integer polynomial of degree n and coefficient bitsize τ # roots / tree width: n min. root separation: Õ(nτ) subdivision rule and subdivision precision overall bit model of computation tree size demand complexity bisection, exact over Z / Q Õ(nτ) Õ(n τ) Õ(n 4 τ ) bisection, approximate Õ(nτ) Õ(nτ) Õ(n 3 τ ) Newton, exact over Z / Q Õ(n) 1 Õ(n τ) Õ(n 3 τ) 1 amortized over entire tree: Õ(nτ) 1 [Sagraloff: ISSAC 01] 11
42 Complexity ( benchmark problem ) isolate all roots of an integer polynomial of degree n and coefficient bitsize τ # roots / tree width: n min. root separation: Õ(nτ) subdivision rule and subdivision precision overall bit model of computation tree size demand complexity bisection, exact over Z / Q Õ(nτ) Õ(n τ) Õ(n 4 τ ) bisection, approximate Õ(nτ) Õ(nτ) Õ(n 3 τ ) Newton, exact over Z / Q Õ(n) Õ(n τ) Õ(n 3 τ) Newton, approximate Õ(n) Õ(nτ) Õ(n 3 + n τ) 1 amortized over entire tree: Õ(nτ) amortized over entire tree: Õ(n + τ) 1 [Sagraloff, Mehlhorn: JSC 016] best known: Õ(n τ) [Pan: JSC 00] [Mehlhorn, Sagraloff, Wang: JSC 015] 11
43 and now For Real! Implementation results
44 Implementation RS C library for real root solving, refinement & more sophisticated implementation of classical Descartes tailored for approximate arithmetic high-performance general purpose solver default real root solver in Maple since version 11 [Rouillier, Zimmermann: JCAM 003] 1
45 Implementation ANewDsc (Approximate Arithmetic Newton-Descartes) implemented on top of RS merged admissible point selection & Newton-Descartes heuristics to reduce / eliminate overhead on easy instances additional optimizations certified output (based on interval arithmetic using MPFI) matches theoretical worst-case bit complexity (Las Vegas) (assuming asymptotically fast polynomial arithmetic) to be integrated in Maple 01x? [Sagraloff, Mehlhorn: JSC 016] 1
46 Benchmark excerpt: tame (well-separated) instances instance degree bitsize RS ANewDsc time [s] nodes time [s] nodes speedup random Hermite Legendre Wilkinson
47 Benchmark excerpt: hard (clustered) instances instance degree bitsize RS ANewDsc time [s] nodes time [s] nodes speedup monic random random
48 Benchmark excerpt: hard (clustered) instances instance degree bitsize Mignotte polynomials: x n (ax 1) RS ANewDsc time [s] nodes time [s] nodes speedup a = a = a = two clustered roots near 1/a with separation of appx. a n/ 14
49 Comparison to other solvers instance n τ MPSolve CF Sage RS ANewDsc random Hermite Legendre Wilkinson random > random monic > 700 > x n ( a x 1) x n (ax 1) x n (ax 1) > nested 4-fold
50 Bit complexity estimation: tame instances 6k time [s] 0 time [s] 4k k 16K 3K 48K n 10 τ 16K 3K 48K 64K estimated exponent for degree: 1.97 estimated exponent for bitsize: 0.03 theoretical vs. estimated bit complexity: Õ(n 3 + n τ) vs. O(n 1.97 τ 0.03 ) degree n, integer coefficients in range ( τ, τ ) chosen uniformly at random left: τ = 104; right: n =
51 Bit complexity estimation: hard instances 1k time [s] 300 time [s] K K 3K n 1.5K 3K 4.5K τ estimated exponent for degree:.37 estimated exponent for bitsize: 1.45 theoretical vs. estimated bit complexity: Õ(n 3 + n τ) vs. O(n.37 τ 1.45 ) Mignotte-like: x n (ax 1) 3 (cluster of multiplicity 3 around irrational center) left: a = 56 1; right: n =
52 ANewDsc: Approximate Arithmetic Newton-Descartes highly efficient general-purpose solver 18
53 ANewDsc: Approximate Arithmetic Newton-Descartes highly efficient general-purpose solver tremendous speedups in degenerate situations without tailored tweaks (e.g., does not (yet) exploit sparsity) 18
54 ANewDsc: Approximate Arithmetic Newton-Descartes highly efficient general-purpose solver tremendous speedups in degenerate situations without tailored tweaks (e.g., does not (yet) exploit sparsity) no 18
55 ANewDsc: Approximate Arithmetic Newton-Descartes highly efficient general-purpose solver tremendous speedups in degenerate situations without tailored tweaks (e.g., does not (yet) exploit sparsity) easily solves previously infeasible instances no 18
56 ANewDsc: Approximate Arithmetic Newton-Descartes highly efficient general-purpose solver tremendous speedups in degenerate situations without tailored tweaks (e.g., does not (yet) exploit sparsity) easily solves previously infeasible instances first available solver with expected performance near theoretical optimum native support for inputs with arbitrary real coefficients no 18
57 ANewDsc: Approximate Arithmetic Newton-Descartes highly efficient general-purpose solver tremendous speedups in degenerate situations without tailored tweaks (e.g., does not (yet) exploit sparsity) easily solves previously infeasible instances first available solver with expected performance near theoretical optimum native support for inputs with arbitrary real coefficients implementation & benchmarks available at ANewDsc.mpi-inf.mpg.de no 18
58 ANewDsc: Approximate Arithmetic Newton-Descartes highly efficient general-purpose solver tremendous speedups in degenerate situations without tailored tweaks (e.g., does not (yet) exploit sparsity) easily solves previously infeasible instances first available solver with expected performance near theoretical optimum native support for inputs with arbitrary real coefficients implementation & benchmarks available at ANewDsc.mpi-inf.mpg.de no Thank you for your attention! 18
59
60 Newton-Descartes 0 0 choose samples x i
61 Newton-Descartes 0 0 choose samples x i compute Newton correction terms N(x i ) = P(x i )/P (x i )
62 Newton-Descartes 0 0 choose samples x i compute Newton correction terms N(x i ) = P(x i )/P (x i ) guess multiplicity k s.t. x i k N(x i ) is approximately the same
63 Newton-Descartes 0 0 choose samples x i compute Newton correction terms N(x i ) = P(x i )/P (x i ) guess multiplicity k s.t. x i k N(x i ) is approximately the same verify that no roots are lost; otherwise, resort to bisection
64
65 Comparison to other solvers instance n τ MPSolve CF Sage SLV RS ANewDsc random Hermite Legendre Wilkinson (80.9) random > > random monic > > 700 > 700 > x n ( a x 1) x n (ax 1) x n (ax 1) > > nested 4-fold
66
Computing Real Roots of Real Polynomials
Computing Real Roots of Real Polynomials and now For Real! Alexander Kobel, 1,2 Fabrice Rouillier 3 and Michael Sagraloff 1 alexander.kobel@mpi-inf.mpg.de Fabrice.Rouillier@inria.fr michael.sagraloff@mpi-inf.mpg.de
More informationRoot Clustering for Analytic and their Complexity
Root Clustering for Analytic and their Complexity Chee Yap Courant Institute, NYU August 5, 2015 Mini-Symposium: Algorithms and Complexity in Polynomial System Solving SIAM AG15, Daejeon, Korea Joint work
More informationCertified Complex Numerical Root Finding
Intro Subdivision Numerical Comparison Goals Max Planck Institute for Informatics Dep. 1: Algorithms and Complexity Seminar on Computational Geometry and Geometric Computing May 6, 2010 Intro Subdivision
More informationComputing Real Roots of Real Polynomials
Computing Real Roots of Real Polynomials Michael Sagraloff Kurt Mehlhorn October 18, 2013 Abstract Computing the roots of a univariate polynomial is a fundamental and long-studied problem of computational
More informationComputing real roots of real polynomials... and now for real!
Computing real roots of real polynomials...... and now for real! Alexander Kobel Max-Planck-Institut für Informatik Campus E1 4 66123 Saarbrücken, Germany alexander.kobel@mpiinf.mpg.de Fabrice Rouillier
More informationA General Approach to Isolating Roots of a Bitstream
A General Approach to Isolating Roots of a Bitstream Polynomial Michael Sagraloff Abstract. We describe a new approach to isolate the roots (either real or complex) of a square-free polynomial F with real
More informationImplementation of a Near-Optimal Complex Root Clustering Algorithm
Implementation of a Near-Optimal Complex Root Clustering Algorithm Rémi Imbach 1, Victor Y. Pan 2, and Chee Yap 3 1 TU Kaiserslautern Email: imbach@mathematik.uni-kl.de www.mathematik.uni-kl.de/en/agag/members/
More informationNew Bounds in the Analysis of Subdivision Algorithms for Real Root Isolation. Chee Yap Courant Institute of Mathematical Sciences New York University
New Bounds in the Analysis of Subdivision Algorithms for Real Root Isolation 1 Chee Yap Courant Institute of Mathematical Sciences New York University Joint work with Vikram Sharma, Zilin Du, Chris Wu,
More informationQuadratic Interval Refinement
(QIR) Analysis of Seminar on Computational Geometry and Geometric Computing Outline Introduction (QIR) Analysis of 1 Introduction in Short Main Goal of this Work 2 3 (QIR) Bisection Method Algorithm 4
More informationReliable and Efficient Geometric Computation p.1/27
Reliable and Efficient Geometric Computation Kurt Mehlhorn Max-Planck-Institut für Informatik slides and papers are available at my home page partially supported by EU-projects ECG and ACS Effective Computational
More informationA New Method for Real Root Isolation of Univariate
A New Method for Real Root Isolation of Univariate Polynomials Ting Zhang and Bican Xia Abstract. A new algorithm for real root isolation of univariate polynomials is proposed, which is mainly based on
More informationarxiv: v1 [cs.na] 30 Apr 2013
Fast Approximate Polynomial Multipoint Evaluation and Applications Alexander Kobel 1 3 Michael Sagraloff 1 1 Max-Planck-Institut für Informatik 2 International Max Planck Research School for Computer Science
More informationOn multiple roots in Descartes Rule and their distance to roots of higher derivatives
Journal of Computational and Applied Mathematics 200 (2007) 226 230 www.elsevier.com/locate/cam Short Communication On multiple roots in Descartes Rule and their distance to roots of higher derivatives
More informationThe Bernstein Basis and Real Root Isolation
Combinatorial and Computational Geometry MSRI Publications Volume 52, 2005 The Bernstein Basis and Real Root Isolation BERNARD MOURRAIN, FABRICE ROUILLIER, AND MARIE-FRANÇOISE ROY Abstract. In this mostly
More informationComplexity Analysis of Root Clustering for a Complex Polynomial
Complexity Analysis of Root Clustering for a Complex Polynomial Ruben Becker MPI for Informatics, Saarbrücken Graduate School of Computer Science Saarbrücken, Germany ruben@mpi-inf.mpg.de Juan Xu Courant
More informationImproved Budan-Fourier Count for Root Finding
Improved Budan-Fourier Count for Root Finding André Galligo To cite this version: André Galligo. Improved Budan-Fourier Count for Root Finding. 2011. HAL Id: hal-00653762 https://hal.inria.fr/hal-00653762
More informationACS. Algorithms for Complex Shapes with Certified Numerics and Topology
ACS Algorithms for Complex Shapes with Certified Numerics and Topology Experimental implementation of more operations on algebraic numbers, possibly with the addition of numeric filters, and of robust
More informationExperimental evaluation and cross-benchmarking of univariate real solvers
Experimental evaluation and cross-benchmarking of univariate real solvers Michael Hemmer Max-Planck-Institut für Informatik Saarbrücken, Germany Ioannis Z. Emiris University of Athens Athens, Greece Elias
More informationReal Solving on Algebraic Systems of Small Dimension
Real Solving on Algebraic Systems of Small Dimension Master s Thesis Presentation Dimitrios I. Diochnos University of Athens March 8, 2007 D. I. Diochnos (Univ. of Athens, µ Q λ ) Real Solving on Bivariate
More informationComputing Rational Points in Convex Semi-algebraic Sets and Sums-of-Squares Decompositions
Computing Rational Points in Convex Semi-algebraic Sets and Sums-of-Squares Decompositions Mohab Safey El Din 1 Lihong Zhi 2 1 University Pierre et Marie Curie, Paris 6, France INRIA Paris-Rocquencourt,
More informationReal Root Isolation of Regular Chains.
Real Root Isolation of Regular Chains. François Boulier 1, Changbo Chen 2, François Lemaire 1, Marc Moreno Maza 2 1 University of Lille I (France) 2 University of London, Ontario (Canada) ASCM 2009 (Boulier,
More informationMath Models of OR: Branch-and-Bound
Math Models of OR: Branch-and-Bound John E. Mitchell Department of Mathematical Sciences RPI, Troy, NY 12180 USA November 2018 Mitchell Branch-and-Bound 1 / 15 Branch-and-Bound Outline 1 Branch-and-Bound
More informationExperimental Evaluation and Cross-Benchmarking of Univariate Real Solvers
Experimental Evaluation and Cross-Benchmarking of Univariate Real Solvers Michael Hemmer MPI for Informatics Saarbrücken, Germany hemmer@mpi-inf.mpg.de Ioannis Z. Emiris University of Athens Athens, Greece
More informationx 2 + 6x 18 x + 2 Name: Class: Date: 1. Find the coordinates of the local extreme of the function y = x 2 4 x.
1. Find the coordinates of the local extreme of the function y = x 2 4 x. 2. How many local maxima and minima does the polynomial y = 8 x 2 + 7 x + 7 have? 3. How many local maxima and minima does the
More informationReal Solving on Bivariate Systems with Sturm Sequences and SLV Maple TM library
Real Solving on Bivariate Systems with Sturm Sequences and SLV Maple TM library Dimitris Diochnos University of Illinois at Chicago Dept. of Mathematics, Statistics, and Computer Science September 27,
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationCylindrical Algebraic Decomposition in Coq
Cylindrical Algebraic Decomposition in Coq MAP 2010 - Logroño 13-16 November 2010 Assia Mahboubi INRIA Microsoft Research Joint Centre (France) INRIA Saclay Île-de-France École Polytechnique, Palaiseau
More informationSPECTRA - a Maple library for solving linear matrix inequalities in exact arithmetic
SPECTRA - a Maple library for solving linear matrix inequalities in exact arithmetic Didier Henrion Simone Naldi Mohab Safey El Din Version 1.0 of November 5, 2016 Abstract This document briefly describes
More informationComparing real algebraic numbers of small degree
Comparing real algebraic numbers of small degree Ioannis Z. Emiris and Elias P. Tsigaridas Department of Informatics and Telecommunications National Kapodistrian University of Athens, Greece {emiris,et}@di.uoa.gr
More informationRational Exponents. Polynomial function of degree n: with leading coefficient,, with maximum number of turning points is given by (n-1)
Useful Fact Sheet Final Exam Interval, Set Builder Notation (a,b) = {x a
More informationImproved algorithms for solving bivariate systems via Rational Univariate Representations
Improved algorithms for solving bivariate systems via Rational Univariate Representations Yacine Bouzidi, Sylvain Lazard, Guillaume Moroz, Marc Pouget, Fabrice Rouillier, Michael Sagraloff To cite this
More informationOn Computation of Positive Roots of Polynomials and Applications to Orthogonal Polynomials. University of Bucharest CADE 2007.
On Computation of Positive Roots of Polynomials and Applications to Orthogonal Polynomials Doru Ştefănescu University of Bucharest CADE 2007 21 February 2007 Contents Approximation of the Real Roots Bounds
More informationSolving bivariate systems using Rational Univariate Representations
Solving bivariate systems using Rational Univariate Representations Yacine Bouzidi, Sylvain Lazard, Guillaume Moroz, Marc Pouget, Fabrice Rouillier, Michael Sagraloff To cite this version: Yacine Bouzidi,
More informationComputation of the error functions erf and erfc in arbitrary precision with correct rounding
Computation of the error functions erf and erfc in arbitrary precision with correct rounding Sylvain Chevillard Arenaire, LIP, ENS-Lyon, France Sylvain.Chevillard@ens-lyon.fr Nathalie Revol INRIA, Arenaire,
More informationZeroes of Transcendental and Polynomial Equations. Bisection method, Regula-falsi method and Newton-Raphson method
Zeroes of Transcendental and Polynomial Equations Bisection method, Regula-falsi method and Newton-Raphson method PRELIMINARIES Solution of equation f (x) = 0 A number (real or complex) is a root of the
More informationON THE VARIOUS BISECTION METHODS DERIVED FROM VINCENT S THEOREM. Alkiviadis G. Akritas, Adam W. Strzeboński, Panagiotis S. Vigklas
Serdica J. Computing (008), 89 104 ON THE VARIOUS BISECTION METHODS DERIVED FROM VINCENT S THEOREM Alkiviadis G. Akritas, Adam W. Strzeboński, Panagiotis S. Vigklas Dedicated to Professors Alberto Alesina
More informationReal Algebraic Numbers: Complexity Analysis and Experimentation
Real Algebraic Numbers: Complexity Analysis and Experimentation Ioannis Z. Emiris 1 Bernard Mourrain and Elias P. Tsigaridas 1 1 Department of Informatics and Telecommunications National Kapodistrian University
More informationAlgebraic algorithms and applications to geometry
Algebraic algorithms and applications to geometry Elias P. Tsigaridas Department of Informatics and Telecommunications National Kapodistrian University of Athens, HELLAS et@di.uoa.gr Abstract. Real algebraic
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationWe consider the problem of finding a polynomial that interpolates a given set of values:
Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.
More information6.1 Polynomial Functions
6.1 Polynomial Functions Definition. A polynomial function is any function p(x) of the form p(x) = p n x n + p n 1 x n 1 + + p 2 x 2 + p 1 x + p 0 where all of the exponents are non-negative integers and
More informationNumerical Methods in Informatics
Numerical Methods in Informatics Lecture 2, 30.09.2016: Nonlinear Equations in One Variable http://www.math.uzh.ch/binf4232 Tulin Kaman Institute of Mathematics, University of Zurich E-mail: tulin.kaman@math.uzh.ch
More informationReal Root Isolation of Polynomial Equations Ba. Equations Based on Hybrid Computation
Real Root Isolation of Polynomial Equations Based on Hybrid Computation Fei Shen 1 Wenyuan Wu 2 Bican Xia 1 LMAM & School of Mathematical Sciences, Peking University Chongqing Institute of Green and Intelligent
More informationComplete Numerical Isolation of Real Zeros in General Triangular Systems
Complete Numerical Isolation of Real Zeros in General Triangular Systems Jin-San Cheng 1, Xiao-Shan Gao 1 and Chee-Keng Yap 2,3 1 Key Lab of Mathematics Mechanization Institute of Systems Science, AMSS,
More informationSolving Sparse Rational Linear Systems. Pascal Giorgi. University of Waterloo (Canada) / University of Perpignan (France) joint work with
Solving Sparse Rational Linear Systems Pascal Giorgi University of Waterloo (Canada) / University of Perpignan (France) joint work with A. Storjohann, M. Giesbrecht (University of Waterloo), W. Eberly
More informationUniversity of Bristol - Explore Bristol Research. Peer reviewed version. Link to publication record in Explore Bristol Research PDF-document
Huelse, D., & Hemmer, M. (29). Generic implementation of a modular gcd over algebraic extension fields. Paper presented at 25th European Workshop on Computational Geometry, Brussels, Belgium. Peer reviewed
More informationHierarchical Matrices. Jon Cockayne April 18, 2017
Hierarchical Matrices Jon Cockayne April 18, 2017 1 Sources Introduction to Hierarchical Matrices with Applications [Börm et al., 2003] 2 Sources Introduction to Hierarchical Matrices with Applications
More informationThere are self-avoiding walks of steps on Z 3
There are 7 10 26 018 276 self-avoiding walks of 38 797 311 steps on Z 3 Nathan Clisby MASCOS, The University of Melbourne Institut für Theoretische Physik Universität Leipzig November 9, 2012 1 / 37 Outline
More informationNear Optimal Tree Size Bounds on a Simple Real Root Isolation Algorithm
Near Optimal Tree Size Bounds on a Simple Real Root Isolation Algorithm ABSTRACT Vikram Sharma Institute of Mathematical Sciences Chennai, India 6003 vikram@imsc.res.in The problem of isolating all real
More informationAsymptotic redundancy and prolixity
Asymptotic redundancy and prolixity Yuval Dagan, Yuval Filmus, and Shay Moran April 6, 2017 Abstract Gallager (1978) considered the worst-case redundancy of Huffman codes as the maximum probability tends
More informationStatements, Implication, Equivalence
Part 1: Formal Logic Statements, Implication, Equivalence Martin Licht, Ph.D. January 10, 2018 UC San Diego Department of Mathematics Math 109 A statement is either true or false. We also call true or
More informationCSE 206A: Lattice Algorithms and Applications Spring Basic Algorithms. Instructor: Daniele Micciancio
CSE 206A: Lattice Algorithms and Applications Spring 2014 Basic Algorithms Instructor: Daniele Micciancio UCSD CSE We have already seen an algorithm to compute the Gram-Schmidt orthogonalization of a lattice
More informationAdvanced Optimization
Advanced Optimization Lecture 3: 1: Randomized Algorithms for for Continuous Discrete Problems Problems November 22, 2016 Master AIC Université Paris-Saclay, Orsay, France Anne Auger INRIA Saclay Ile-de-France
More informationarxiv: v1 [math.ag] 27 Jul 2018
Noname manuscript No. (will be inserted by the editor) Bounds for polynomials on algebraic numbers and application to curve topology Daouda Niang Diatta Sény Diatta Fabrice Rouillier Marie-Françoise Roy
More informationComplete Numerical Isolation of Real Zeros in General Triangular Systems 1)
MM Research Preprints, 1 31 KLMM, AMSS, Academia Sinica Vol. 25, December 2006 1 Complete Numerical Isolation of Real Zeros in General Triangular Systems 1) Jin-San Cheng 1, Xiao-Shan Gao 1 and Chee-Keng
More information) = nlog b ( m) ( m) log b ( ) ( ) = log a b ( ) Algebra 2 (1) Semester 2. Exponents and Logarithmic Functions
Exponents and Logarithmic Functions Algebra 2 (1) Semester 2! a. Graph exponential growth functions!!!!!! [7.1]!! - y = ab x for b > 0!! - y = ab x h + k for b > 0!! - exponential growth models:! y = a(
More informationNumerical integration in arbitrary-precision ball arithmetic
/ 24 Numerical integration in arbitrary-precision ball arithmetic Fredrik Johansson (LFANT, Bordeaux) Journées FastRelax, Inria Sophia Antipolis 7 June 208 Numerical integration in Arb 2 / 24 New code
More informationSieving for Shortest Vectors in Ideal Lattices:
Sieving for Shortest Vectors in Ideal Lattices: a Practical Perspective Joppe W. Bos Microsoft Research LACAL@RISC Seminar on Cryptologic Algorithms CWI, Amsterdam, Netherlands Joint work with Michael
More informationComputing the RSA Secret Key is Deterministic Polynomial Time Equivalent to Factoring
Computing the RSA Secret Key is Deterministic Polynomial Time Equivalent to Factoring Alexander May Faculty of Computer Science, Electrical Engineering and Mathematics University of Paderborn 33102 Paderborn,
More informationAccurate Multiple-Precision Gauss-Legendre Quadrature
Accurate Multiple-Precision Gauss-Legendre Quadrature Laurent Fousse Université Henri-Poincaré Nancy 1 laurent@komite.net Abstract Numerical integration is an operation that is frequently available in
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More information23. Cutting planes and branch & bound
CS/ECE/ISyE 524 Introduction to Optimization Spring 207 8 23. Cutting planes and branch & bound ˆ Algorithms for solving MIPs ˆ Cutting plane methods ˆ Branch and bound methods Laurent Lessard (www.laurentlessard.com)
More informationLecture 5: Random numbers and Monte Carlo (Numerical Recipes, Chapter 7) Motivations for generating random numbers
Lecture 5: Random numbers and Monte Carlo (Numerical Recipes, Chapter 7) Motivations for generating random numbers To sample a function in a statistically controlled manner (i.e. for Monte Carlo integration)
More informationContinued fractions and number systems: applications to correctly-rounded implementations of elementary functions and modular arithmetic.
Continued fractions and number systems: applications to correctly-rounded implementations of elementary functions and modular arithmetic. Mourad Gouicem PEQUAN Team, LIP6/UPMC Nancy, France May 28 th 2013
More informationAn Improved Approximation Algorithm for Virtual Private Network Design
An Improved Approximation Algorithm for Virtual Private Network Design Friedrich Eisenbrand Fabrizio Grandoni Abstract Virtual private network design deals with the reservation of capacities in a network,
More informationPartitions and Covers
University of California, Los Angeles CS 289A Communication Complexity Instructor: Alexander Sherstov Scribe: Dong Wang Date: January 2, 2012 LECTURE 4 Partitions and Covers In previous lectures, we saw
More informationA strongly polynomial algorithm for linear systems having a binary solution
A strongly polynomial algorithm for linear systems having a binary solution Sergei Chubanov Institute of Information Systems at the University of Siegen, Germany e-mail: sergei.chubanov@uni-siegen.de 7th
More informationHonors Advanced Mathematics November 4, /2.6 summary and extra problems page 1 Recap: complex numbers
November 4, 013.5/.6 summary and extra problems page 1 Recap: complex numbers Number system The complex number system consists of a + bi where a and b are real numbers, with various arithmetic operations.
More informationLecture: Local Spectral Methods (1 of 4)
Stat260/CS294: Spectral Graph Methods Lecture 18-03/31/2015 Lecture: Local Spectral Methods (1 of 4) Lecturer: Michael Mahoney Scribe: Michael Mahoney Warning: these notes are still very rough. They provide
More information1 The Fundamental Theorem of Arithmetic. A positive integer N has a unique prime power decomposition. Primality Testing. and. Integer Factorisation
1 The Fundamental Theorem of Arithmetic A positive integer N has a unique prime power decomposition 2 Primality Testing Integer Factorisation (Gauss 1801, but probably known to Euclid) The Computational
More information1. Definition of a Polynomial
1. Definition of a Polynomial What is a polynomial? A polynomial P(x) is an algebraic expression of the form Degree P(x) = a n x n + a n 1 x n 1 + a n 2 x n 2 + + a 3 x 3 + a 2 x 2 + a 1 x + a 0 Leading
More informationOn solving zero-dimensional and parametric systems
On solving zero-dimensional and parametric systems Fabrice Rouillier Fabrice.Rouillier@inria.fr - http://fgbrs.lip6.fr/~rouillie SALSA (INRIA) project and CALFOR (LIP6) team Paris, France General Objectives
More informationCOMP 382: Reasoning about algorithms
Fall 2014 Unit 4: Basics of complexity analysis Correctness and efficiency So far, we have talked about correctness and termination of algorithms What about efficiency? Running time of an algorithm For
More informationShort Division of Long Integers. (joint work with David Harvey)
Short Division of Long Integers (joint work with David Harvey) Paul Zimmermann October 6, 2011 The problem to be solved Divide efficiently a p-bit floating-point number by another p-bit f-p number in the
More informationCOMP 355 Advanced Algorithms
COMP 355 Advanced Algorithms Algorithm Design Review: Mathematical Background 1 Polynomial Running Time Brute force. For many non-trivial problems, there is a natural brute force search algorithm that
More informationAccelerating the Numerical Computation of Positive Roots of Polynomials using Improved Bounds
Accelerating the Numerical Computation of Positive Roots of Polynomials using Improved Bounds Kinji Kimura 1, Takuto Akiyama, Hiroyuki Ishigami 3, Masami Takata 4, and Yoshimasa Nakamura 5 1,,3,5 Graduate
More informationMath /Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined
Math 400-001/Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined using limits. As a particular case, the derivative of f(x)
More informationSCALED REMAINDER TREES
Draft. Aimed at Math. Comp. SCALED REMAINDER TREES DANIEL J. BERNSTEIN Abstract. It is well known that one can compute U mod p 1, U mod p 2,... in time n(lg n) 2+o(1) where n is the number of bits in U,
More informationEquations and Inequalities
Chapter 3 Equations and Inequalities Josef Leydold Bridging Course Mathematics WS 2018/19 3 Equations and Inequalities 1 / 61 Equation We get an equation by equating two terms. l.h.s. = r.h.s. Domain:
More informationFinding a Heaviest Triangle is not Harder than Matrix Multiplication
Finding a Heaviest Triangle is not Harder than Matrix Multiplication Artur Czumaj Department of Computer Science New Jersey Institute of Technology aczumaj@acm.org Andrzej Lingas Department of Computer
More informationCSE548, AMS542: Analysis of Algorithms, Fall 2016 Date: Nov 30. Final In-Class Exam. ( 7:05 PM 8:20 PM : 75 Minutes )
CSE548, AMS542: Analysis of Algorithms, Fall 2016 Date: Nov 30 Final In-Class Exam ( 7:05 PM 8:20 PM : 75 Minutes ) This exam will account for either 15% or 30% of your overall grade depending on your
More informationCSE373: Data Structures and Algorithms Lecture 2: Math Review; Algorithm Analysis. Hunter Zahn Summer 2016
CSE373: Data Structures and Algorithms Lecture 2: Math Review; Algorithm Analysis Hunter Zahn Summer 2016 Today Finish discussing stacks and queues Review math essential to algorithm analysis Proof by
More informationACS. Benchmarks and evaluation of experimental algebraic kernels. Ioannis Z. Emiris. ACS Technical Report No.: ACS-TR
ACS Algorithms for Complex Shapes with Certified Numerics and Topology Benchmarks and evaluation of experimental algebraic kernels Dimitrios I. Diochnos Ioannis Z. Emiris Elias P. Tsigaridas ACS Technical
More informationA Worst-case Bound for Topology Computation of Algebraic Curves
arxiv:1104.1510v1 [cs.sc] 8 Apr 2011 A Worst-case Bound for Topology Computation of Algebraic Curves Michael Kerber Institute of Science and Technology (IST) Austria Klosterneuburg, Austria mkerber@ist.ac.at
More information1 Computational Problems
Stanford University CS254: Computational Complexity Handout 2 Luca Trevisan March 31, 2010 Last revised 4/29/2010 In this lecture we define NP, we state the P versus NP problem, we prove that its formulation
More informationOn Newton-Raphson iteration for multiplicative inverses modulo prime powers
On Newton-Raphson iteration for multiplicative inverses modulo prime powers Jean-Guillaume Dumas To cite this version: Jean-Guillaume Dumas. On Newton-Raphson iteration for multiplicative inverses modulo
More informationSo far we have implemented the search for a key by carefully choosing split-elements.
7.7 Hashing Dictionary: S. insert(x): Insert an element x. S. delete(x): Delete the element pointed to by x. S. search(k): Return a pointer to an element e with key[e] = k in S if it exists; otherwise
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationarxiv: v5 [cs.sc] 15 May 2018
On Newton-Raphson iteration for multiplicative inverses modulo prime powers arxiv:109.666v5 [cs.sc] 15 May 018 Jean-Guillaume Dumas May 16, 018 Abstract We study algorithms for the fast computation of
More informationDeterministic Global Optimization for Dynamic Systems Using Interval Analysis
Deterministic Global Optimization for Dynamic Systems Using Interval Analysis Youdong Lin and Mark A. Stadtherr Department of Chemical and Biomolecular Engineering University of Notre Dame, Notre Dame,
More informationarxiv: v1 [cs.sc] 24 Mar 2014
Model-based construction of Open Non-uniform Cylindrical Algebraic Decompositions arxiv:1403.6487v1 [cs.sc] 24 Mar 2014 Christopher W. Brown United States Naval Academy Annapolis, Maryland 21402 wcbrown@usna.edu
More informationSolving fixed-point equations over semirings
Solving fixed-point equations over semirings Javier Esparza Technische Universität München Joint work with Michael Luttenberger and Maximilian Schlund Fixed-point equations We study systems of equations
More informationA note on the complexity of univariate root isolation
A note on the complexity of univariate root isolation Ioannis Emiris, Elias P. Tsigaridas To cite this version: Ioannis Emiris, Elias P. Tsigaridas. A note on the complexity of univariate root isolation.
More informationCOMP 355 Advanced Algorithms Algorithm Design Review: Mathematical Background
COMP 355 Advanced Algorithms Algorithm Design Review: Mathematical Background 1 Polynomial Time Brute force. For many non-trivial problems, there is a natural brute force search algorithm that checks every
More informationRoot Isolation. Lecture 9
Lecture 9 Root Isolation The determination of the roots of a univariate real polynomial is ubiquitous in geometric computing. Let f be such a polynomial and let n be its degree, i.e., f = f i x i R[x].
More informationGoals for This Lecture:
Goals for This Lecture: Learn the Newton-Raphson method for finding real roots of real functions Learn the Bisection method for finding real roots of a real function Look at efficient implementations of
More informationDense Arithmetic over Finite Fields with CUMODP
Dense Arithmetic over Finite Fields with CUMODP Sardar Anisul Haque 1 Xin Li 2 Farnam Mansouri 1 Marc Moreno Maza 1 Wei Pan 3 Ning Xie 1 1 University of Western Ontario, Canada 2 Universidad Carlos III,
More informationWhat s the best data structure for multivariate polynomials in a world of 64 bit multicore computers?
What s the best data structure for multivariate polynomials in a world of 64 bit multicore computers? Michael Monagan Center for Experimental and Constructive Mathematics Simon Fraser University British
More informationA Branch-and-Bound Algorithm for Unconstrained Global Optimization
SCAN 2010, Lyon, September 27 30, 2010 1 / 18 A Branch-and-Bound Algorithm for Unconstrained Global Optimization Laurent Granvilliers and Alexandre Goldsztejn Université de Nantes LINA CNRS Interval-based
More information