Solving Sparse Rational Linear Systems. Pascal Giorgi. University of Waterloo (Canada) / University of Perpignan (France) joint work with
|
|
- Fay Hamilton
- 5 years ago
- Views:
Transcription
1 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 (University of Calgary), G. Villard (ENS Lyon) ISSAC'26, Genova - July, 26
2 Problem Let A a non-singular matrix and b a vector dened over Z. Problem : Compute x = A b over the rational numbers A = B x = A b = B C A, b = C A B C A. Main diculty : expression swell
3 Problem Let A a non-singular matrix and b a vector dened over Z. Problem : Compute x = A b over the rational numbers A = B x = A b = B C A, b = C A B C A. Main diculty : expression swell and take advantage of sparsity
4 Motivations Large linear systems are involved in many mathematical applications Over a nite eld : integers factorization [Odlyzko 999], discrete logarithm [Odlyzko 999 ; Thome 23]. Over the integers : number theory [Cohen 993], group theory [Newman 972], integer programming [Aardal, Hurkens, Lenstra 999]. Rational linear systems are central in recent linear algebra algorithms I Determinant [Abbott, Bronstein, Mulders 999 ; Storjohann 25] I Smith form [Giesbrecht 995 ; Eberly, Giesbrecht, Villard 2] I Nullspace, Kernel [Chen, Storjohann 25]
5 Algorithms for non-singular system solving Dense matrices : I Gaussian elimination and CRA,! O~(n!+ log jjajj) bit operations I P-adic lifting [Monck, Carter 979 ; Dixon 982],! O~(n 3 log jjajj) bit operations I High order lifting [Storjohann 25],! O~(n! log jjajj) bit operations Sparse matrices : I P-adic lifting or CRA [Kaltofen, Saunders 99],! O~(n 2 log jjajj) bit operations with non-zero elts.
6 P-adic algorithm with matrix inversion Scheme to compute A b [Dixon 982] : (-) B := A mod p (-2) r := b for i := to k (2-) x i := B:r mod p (2-2) r := (=p)(r A:x i ) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x
7 P-adic algorithm with matrix inversion Scheme to compute A b [Dixon 982] : (-) B := A mod p O~(n 3 log jjajj) (-2) r := b for i := to k k = O~(n) (2-) x i := B:r mod p O~(n 2 log jjajj) (2-2) r := (=p)(r A:x i ) O~(n 2 log jjajj) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x
8 P-adic algorithm with matrix inversion Scheme to compute A b [Dixon 982] : (-) B := A mod p O~(n 3 log jjajj) (-2) r := b for i := to k k = O~(n) (2-) x i := B:r mod p O~(n 2 log jjajj) (2-2) r := (=p)(r A:x i ) O~(n 2 log jjajj) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x Main operations : matrix inversion and matrix-vector products
9 Dense linear system solving in practice Ecient implementations are available : LinBox. [ Use tuned BLAS oating-point library for exact arithmetic matrix inversion over prime eld [Dumas, Giorgi, Pernet 24] BLAS matrix-vector product with CRT over the integers Rational number reconstruction half GCD [Schonage 97] heuristic using integer multiplication [NTL library] random dense linear system with 3 bits entries (P4-3.4Ghz) n Time :6s 4:3s 3:s 99:6s 236:8s 449:2s performances improvement of a factor over NTL's tuned implementation
10 What does happen when matrices are sparse? We consider sparse matrices with O(n) non zero elements,! matrix-vector product needs only O(n) operations.
11 Sparse linear system and P-adic lifting Computing the modular inverse is proscribed due to ll-in Solution [Kaltofen, Saunders 99] :,! use modular minimal polynomial instead of inverse
12 Sparse linear system and P-adic lifting Computing the modular inverse is proscribed due to ll-in Solution [Kaltofen, Saunders 99] :,! use modular minimal polynomial instead of inverse Wiedemann approach (986) Let A 2 IF n n of full rank and b 2 IF n. Then x = A b can be expressed as a linear combination of the Krylov subspace fb; Ab; :::; A n bg Let f A () = f + f + ::: + f d d 2 IF[] be the minimal polynomial of A
13 Sparse linear system and P-adic lifting Computing the modular inverse is proscribed due to ll-in Solution [Kaltofen, Saunders 99] :,! use modular minimal polynomial instead of inverse Wiedemann approach (986) Let A 2 IF n n of full rank and b 2 IF n. Then x = A b can be expressed as a linear combination of the Krylov subspace fb; Ab; :::; A n bg Let f A () = f + f + ::: + f d d 2 IF[] be the minimal polynomial of A A b = f (f b + f 2 Ab + ::: + f d A d b) {z }
14 Sparse linear system and P-adic lifting Computing the modular inverse is proscribed due to ll-in Solution [Kaltofen, Saunders 99] :,! use modular minimal polynomial instead of inverse Wiedemann approach (986) Let A 2 IF n n of full rank and b 2 IF n. Then x = A b can be expressed as a linear combination of the Krylov subspace fb; Ab; :::; A n bg Let f A () = f + f + ::: + f d d 2 IF[] be the minimal polynomial of A A b = f (f b + f 2 Ab + ::: + f d A d b) {z } x Applying minpoly in each lifting steps cost O~(nd) eld operations, then giving a worst case complexity of O~(n 3 log jjajj) bit operations.
15 Sparse linear system solving in practice use of LinBox library on Itanium II -.3Ghz, 28Gb RAM random systems with 3 bits entries and elts/row (plus identity) system order Maple 64:7s 849s 98s CRA-Wied 4:8s 68s 7s 3857s 452s P-adic-Wied :2s 3s 693s 2629s 834s Dixon :9s s 42s 78s 429s
16 Sparse linear system solving in practice use of LinBox library on Itanium II -.3Ghz, 28Gb RAM random systems with 3 bits entries and elts/row (plus identity) system order Maple 64:7s 849s 98s CRA-Wied 4:8s 68s 7s 3857s 452s P-adic-Wied :2s 3s 693s 2629s 834s Dixon :9s s 42s 78s 429s main dierence : (2-) x i = B:r mod p (Dixon) P (2-) x i := deg A f f i= f i A i r mod p (P-adic-Wied) Remark : n sparse matrix applications is far from level 2 BLAS in practice.
17 Our objectives In practice : Integrate level 2 and 3 BLAS in integer sparse solver In theory : Improve bit complexity of sparse linear system solving =) O~(n ) bits operations with < 3?
18 Our alternative to Block Wiedemann Express the inverse of the sparse matrix through a structured form,! block Hankel/Toeplitz structures Let u 2 IF s n and v 2 IF n s s.t. following matrices are non-singular U = B u ua. ua m C B C A, V = v Av : : : Am v A 2 IFn n
19 Our alternative to Block Wiedemann Express the inverse of the sparse matrix through a structured form,! block Hankel/Toeplitz structures Let u 2 IF s n and v 2 IF n s s.t. following matrices are non-singular U = B u ua. ua m C B C A, V = v Av : : : Am v A 2 IFn n then we can dene the block Hankel matrix H = UAV = B 2 m 2 3 m+ C. A, m m 2m C i = ua i v 2 IF s s and thus we have A = VH U
20 Alternative to Block Wiedemann Nice property on block Hankel matrix inverse [Gohberg, Krupnik 972, Labahn, Choi, Cabay 99] H = B : : :.... C A B : : : C A {z } {z } B : : :.... C A B : : : C A {z } {z } H T H 2 T 2 where H, H 2 are block Hankel matrices and T, T 2 are block Toeplitz matrices
21 Alternative to Block Wiedemann Nice property on block Hankel matrix inverse [Gohberg, Krupnik 972, Labahn, Choi, Cabay 99] H = B : : :.... C A B : : : C A {z } {z } B : : :.... C A B : : : C A {z } {z } H T H 2 T 2 where H, H 2 are block Hankel matrices and T, T 2 are block Toeplitz matrices Block coecients in H, H 2, T, T 2 come from Hermite Pade approximants of H(z) = + 2 z + : : : + 2m z 2m 2 [Labahn, Choi, Cabay 99]. Complexity of H reduces to polynomial matrix multiplication [Giorgi, Jeannerod, Villard 23].
22 Alternative to Block Wiedemann Scheme to compute A b : (-) H(z) := 2m X i= ua i v :z i mod p (-2) compute H mod p from H(z) (-3) r := b for i := to k (2-) x i := V H U:r mod p (2-2) r := (=p)(r A:x i ) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x
23 Alternative to Block Wiedemann Scheme to compute A b : (-) H(z) := 2m X i= ua i v :z i mod p O~(sn 2 log jjajj) (-2) compute H mod p from H(z) O~(s 2 n log jjajj) (-3) r := b for i := to k k = O~(n) (2-) x i := V H U:r mod p O~((n 2 + sn) log jjajj) (2-2) r := (=p)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x
24 Alternative to Block Wiedemann Scheme to compute A b : (-) H(z) := 2m X i= ua i v :z i mod p O~(sn 2 log jjajj) (-2) compute H mod p from H(z) O~(s 2 n log jjajj) (-3) r := b for i := to k k = O~(n) (2-) x i := V H U:r mod p O~((n 2 + sn) log jjajj) (2-2) r := (=p)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x Not yet satisfying : applying matrices U and V is too costly
25 Applying block Krylov subspaces V = v Av : : : A m va 2 IF n n and v 2 IF n s can be rewrite as V = v A + A v A + : : : + A m v A Therefore, applying V to a vector corresponds to : m linear combinations of columns of v m applications of A
26 Applying block Krylov subspaces V = v Av : : : A m va 2 IF n n and v 2 IF n s can be rewrite as V = v A + A v A + : : : + A m v A Therefore, applying V to a vector corresponds to : m linear combinations of columns of v m applications of A O(m sn log jjajj) O(mn log jjajj)
27 Applying block Krylov subspaces V = v Av : : : A m va 2 IF n n and v 2 IF n s can be rewrite as V = v A + A v A + : : : + A m v A Therefore, applying V to a vector corresponds to : m linear combinations of columns of v m applications of A O(m sn log jjajj) How to improve the complexity?
28 Applying block Krylov subspaces V = v Av : : : A m va 2 IF n n and v 2 IF n s can be rewrite as V = v A + A v A + : : : + A m v A Therefore, applying V to a vector corresponds to : m linear combinations of columns of v m applications of A O(m sn log jjajj) How to improve the complexity? ) using special block projections u and v
29 Candidates as suitable block projections Considering A 2 IF n n non-singular and n = m s. Let us denote K(A; v) := v Av A m v 2 IF n n Conjecture : for any non-singular A 2 IF n n and sjn there exists a suitable block projection (R; u; v) 2 IF n n IF s n IF n s such that :. K(RA; v) and K((RA) T ; u T ) are non-singular, 2. R can be applied to a vector with O~(n) operations, 3. u, u T, v and v T can be applied to a vector with O~(n) operations.
30 A structured block projection Let v be dened as follow v T = B v : : : v m vm+ : : : v 2m... v n m+ : : : v n C A 2 IF s n where v i 's are chosen randomly from a sucient large set.
31 A structured block projection Let v be dened as follow v T = B v : : : v m vm+ : : : v 2m... v n m+ : : : v n C A 2 IF s n where v i 's are chosen randomly from a sucient large set. open question : Let R diagonal and v as dened above, is K(RA; v) necessarily non-singular? We prooved it for case s = 2 but answer is still unknown for s > 2
32 Our new algorithm Scheme to compute A b : (-) choose structured blocks u and v (-2) choose R and A := R:A, b := R:b (-3) H(z) := X 2m i= ua i v :z i mod p (-4) compute H mod p from H(z) (-5) r := b for i := to k (2-) x i := V H U:r mod p (2-2) r := (=p)(r A:x i ) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x
33 Our new algorithm Scheme to compute A b : (-) choose structured blocks u and v (-2) choose R and A := R:A, b := R:b (-3) H(z) := X 2m i= ua i v :z i mod p O~(n 2 log jjajj) (-4) compute H mod p from H(z) O~(s 2 n log jjajj) (-5) r := b for i := to k k = O~(n) (2-) x i := V H U:r mod p O~((mn + sn) log jjajj) (2-2) r := (=p)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x
34 Our new algorithm Scheme to compute A b : (-) choose structured blocks u and v (-2) choose R and A := R:A, b := R:b (-3) H(z) := X 2m i= ua i v :z i mod p O~(n 2 log jjajj) (-4) compute H mod p from H(z) O~(s 2 n log jjajj) (-5) r := b for i := to k k = O~(n) (2-) x i := V H U:r mod p O~((mn + sn) log jjajj) (2-2) r := (=p)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i :p i (3-2) rational reconstruction on x taking the optimal m = s = p n gives a complexity of O~(n 2:5 log jjajj)
35 High level implementation LinBox project (Canada-France-USA) : Our tools : BLAS-based matrix multiplication and matrix-vector product polynomial matrix arithmetic (block Hankel inversion),! FFT, Karatsuba, middle product fast application of H is needed to get O~(n 2:5 log jjajj)
36 High level implementation LinBox project (Canada-France-USA) : Our tools : BLAS-based matrix multiplication and matrix-vector product polynomial matrix arithmetic (block Hankel inversion),! FFT, Karatsuba, middle product fast application of H is needed to get O~(n 2:5 log jjajj) I Lagrange's representation of H at the beginning (Horner's scheme) I use evaluation/interpolation on polynomial vectors,! use Vandermonde matrix to have dense matrix operations
37 Is our new algorithm ecient in practice?
38 Comparing performances use of LinBox library on Itanium II -.3Ghz, 28Gb RAM random systems with 3 bits entries and elts/row (plus identity) system order extra memory Maple 849s 98s O() CRA-Wied 68s 7s 3 857s 452s 28 s O(n) P-adic-Wied 3s 693s 2 629s 8 34s 2 s O(n) Dixon s 42s 78s 429s 257s O(n 2 ) Our algo. 5s 6s 75s 426s 937s O(n :5 ) The expected p n improvement is unfortunately amortized by a high constant in the complexity.
39 Sparse solver vs Dixon's algorithm 6 sparsity = elts/row Our algo. Dixon 6 sparsity = 3elts/row Our algo. Dixon Time ( 4 s) 3 Time ( 4 s) System order System order 44 Our algorithm performances are depending on matrix sparsity
40 Practical eect of blocking factors p n blocking factor value is theoretically optimal Is this still true in practice?
41 Practical eect of blocking factors p n blocking factor value is theoretically optimal Is this still true in practice? system order =, optimal block = block size timing 723s 5264s 459s 3833s 4332s system order = 2, optimal block 4 block size timing 4472s 35967s 3854s 2852s 3738s
42 Practical eect of blocking factors p n blocking factor value is theoretically optimal Is this still true in practice? system order =, optimal block = block size timing 723s 5264s 459s 3833s 4332s system order = 2, optimal block 4 block size timing 4472s 35967s 3854s 2852s 3738s best practical blocking factor is dependent upon the ratio of sparse matrix/dense matrix operations eciency
43 Conclusions We provide a new approach for solving sparse integer linear systems : I improve the complexty by a factor p n (heuristic). I improve eciency by minimizing sparse matrix operations and maximizing BLAS use. drawback : not taking advantage of low degree minimal polynomial We propose special block projections for sparse linear algebra,! inverse of sparse matrix in O(n 2:5 ) eld op.
44 Conclusions We provide a new approach for solving sparse integer linear systems : I improve the complexty by a factor p n (heuristic). I improve eciency by minimizing sparse matrix operations and maximizing BLAS use. drawback : not taking advantage of low degree minimal polynomial We propose special block projections for sparse linear algebra,! inverse of sparse matrix in O(n 2:5 ) eld op. Ongoing work : I provide an automatic choice of block dimension (non square?) I prove conjecture for our structured block projections I handle the case of singular matrix I introduce fast matrix multiplication in the complexity
45 Sparse solver vs Dixon's algorithm 4 sparsity=.7% sparsity=.3% sparsity=.% crossover line 2 speed up System order The sparser the matrices are, the earlier the crossover appears
Solving Sparse Integer Linear Systems. Pascal Giorgi. Dali team - University of Perpignan (France)
Solving Sarse Integer Linear Systems Pascal Giorgi Dali team - University of Perignan (France) in collaboration with A. Storjohann, M. Giesbrecht (University of Waterloo), W. Eberly (University of Calgary),
More informationOn the use of polynomial matrix approximant in the block Wiedemann algorithm
On the use of polynomial matrix approximant in the block Wiedemann algorithm ECOLE NORMALE SUPERIEURE DE LYON Laboratoire LIP, ENS Lyon, France Pascal Giorgi Symbolic Computation Group, School of Computer
More informationSolving Sparse Integer Linear Systems
Solving Sparse Integer Linear Systems Wayne Eberly, Mark Giesbrecht, Giorgi Pascal, Arne Storjohann, Gilles Villard To cite this version: Wayne Eberly, Mark Giesbrecht, Giorgi Pascal, Arne Storjohann,
More informationAlgorithms for exact (dense) linear algebra
Algorithms for exact (dense) linear algebra Gilles Villard CNRS, Laboratoire LIP ENS Lyon Montagnac-Montpezat, June 3, 2005 Introduction Problem: Study of complexity estimates for basic problems in exact
More informationWhy solve integer linear systems exactly? David Saunders - U of Delaware
NSF CDI, 2007Oct30 1 Why solve integer linear systems exactly?... an enabling tool for discovery and innovation? David Saunders - U of Delaware LinBox library at linalg.org NSF CDI, 2007Oct30 2 LinBoxers
More informationTHE VECTOR RATIONAL FUNCTION RECONSTRUCTION PROBLEM
137 TH VCTOR RATIONAL FUNCTION RCONSTRUCTION PROBLM ZACH OLSH and ARN STORJOHANN David R. Cheriton School of Computer Science University of Waterloo, Ontario, Canada N2L 3G1 -mail: astorjoh@uwaterloo.ca
More informationLattice reduction of polynomial matrices
Lattice reduction of polynomial matrices Arne Storjohann David R. Cheriton School of Computer Science University of Waterloo Presented at the SIAM conference on Applied Algebraic Geometry at the Colorado
More informationToward High Performance Matrix Multiplication for Exact Computation
Toward High Performance Matrix Multiplication for Exact Computation Pascal Giorgi Joint work with Romain Lebreton (U. Waterloo) Funded by the French ANR project HPAC Séminaire CASYS - LJK, April 2014 Motivations
More informationInverting integer and polynomial matrices. Jo60. Arne Storjohann University of Waterloo
Inverting integer and polynomial matrices Jo60 Arne Storjohann University of Waterloo Integer Matrix Inverse Input: An n n matrix A filled with entries of size d digits. Output: The matrix inverse A 1.
More informationFast computation of normal forms of polynomial matrices
1/25 Fast computation of normal forms of polynomial matrices Vincent Neiger Inria AriC, École Normale Supérieure de Lyon, France University of Waterloo, Ontario, Canada Partially supported by the mobility
More informationFast computation of normal forms of polynomial matrices
1/25 Fast computation of normal forms of polynomial matrices Vincent Neiger Inria AriC, École Normale Supérieure de Lyon, France University of Waterloo, Ontario, Canada Partially supported by the mobility
More informationFraction-free Row Reduction of Matrices of Skew Polynomials
Fraction-free Row Reduction of Matrices of Skew Polynomials Bernhard Beckermann Laboratoire d Analyse Numérique et d Optimisation Université des Sciences et Technologies de Lille France bbecker@ano.univ-lille1.fr
More informationA Fast Las Vegas Algorithm for Computing the Smith Normal Form of a Polynomial Matrix
A Fast Las Vegas Algorithm for Computing the Smith Normal Form of a Polynomial Matrix Arne Storjohann and George Labahn Department of Computer Science University of Waterloo Waterloo, Ontario, Canada N2L
More informationBlack Box Linear Algebra
Black Box Linear Algebra An Introduction to Wiedemann s Approach William J. Turner Department of Mathematics & Computer Science Wabash College Symbolic Computation Sometimes called Computer Algebra Symbols
More informationBlack Box Linear Algebra with the LinBox Library. William J. Turner North Carolina State University wjturner
Black Box Linear Algebra with the LinBox Library William J. Turner North Carolina State University www.math.ncsu.edu/ wjturner Overview 1. Black Box Linear Algebra 2. Project LinBox 3. LinBox Objects 4.
More informationOutput-sensitive algorithms for sumset and sparse polynomial multiplication
Output-sensitive algorithms for sumset and sparse polynomial multiplication Andrew Arnold Cheriton School of Computer Science University of Waterloo Waterloo, Ontario, Canada Daniel S. Roche Computer Science
More informationVector Rational Number Reconstruction
Vector Rational Number Reconstruction Curtis Bright cbright@uwaterloo.ca David R. Cheriton School of Computer Science University of Waterloo, Ontario, Canada N2L 3G1 Arne Storjohann astorjoh@uwaterloo.ca
More informationParallel Computation of the Rank of Large Sparse Matrices from Algebraic K-theory
Parallel Computation of the Rank of Large Sparse Matrices from Algebraic K-theory Jean-Guillaume Dumas Laboratoire Jean Kuntzmann, UMR CNRS 5224 Université Joseph Fourier, B.P. 53 X, 38041 Grenoble, France.
More informationBlack Box Linear Algebra with the LinBox Library. William J. Turner North Carolina State University wjturner
Black Box Linear Algebra with the LinBox Library William J. Turner North Carolina State University www.math.ncsu.edu/ wjturner Overview 1. Black Box Linear Algebra 2. Project LinBox 3. LinBox Objects 4.
More informationFast algorithms for multivariate interpolation problems
Fast algorithms for multivariate interpolation problems Vincent Neiger,, Claude-Pierre Jeannerod Éric Schost Gilles Villard AriC, LIP, École Normale Supérieure de Lyon, France ORCCA, Computer Science Department,
More informationExactly Solve Integer Linear Systems Using Numerical Methods
Exactly Solve Integer Linear Systems Using Numerical Methods Zhendong Wan Department of Computer Science, University of Delaware, Newark, DE, USA Abstract We present a new fast way to exactly solve non-singular
More informationHigh Performance and Reliable Algebraic Computing
High Performance and Reliable Algebraic Computing Soutenance d habilitation à diriger des recherches Clément Pernet Université Joseph Fourier (Grenoble 1) November 25, 2014 Rapporteurs : Daniel Augot,
More informationComputing minimal interpolation bases
Computing minimal interpolation bases Vincent Neiger,, Claude-Pierre Jeannerod Éric Schost Gilles Villard AriC, LIP, École Normale Supérieure de Lyon, France University of Waterloo, Ontario, Canada Supported
More informationVector Rational Number Reconstruction
Vector Rational Number Reconstruction By Curtis Bright A research paper presented to the University of Waterloo In partial fulfillment of the requirements for the degree of Master of Mathematics in Computational
More informationHecke Operators for Arithmetic Groups via Cell Complexes. Mark McConnell. Center for Communications Research, Princeton
Hecke Operators for Arithmetic Groups via Cell Complexes 1 Hecke Operators for Arithmetic Groups via Cell Complexes Mark McConnell Center for Communications Research, Princeton Hecke Operators for Arithmetic
More informationNumeric-Symbolic Exact Rational Linear System Solver
Numeric-Symbolic Exact Rational Linear System Solver B. David Saunders, David Harlan Wood, and Bryan S. Youse Dept. of Computer and Information Sciences, University of Delaware Newark, Delaware, USA saunders@udel.edu,
More informationThe Seven Dwarfs of Symbolic Computation and the Discovery of Reduced Symbolic Models. Erich Kaltofen North Carolina State University google->kaltofen
The Seven Dwarfs of Symbolic Computation and the Discovery of Reduced Symbolic Models Erich Kaltofen North Carolina State University google->kaltofen Outline 2 The 7 Dwarfs of Symbolic Computation Discovery
More informationFast Computation of Hermite Normal Forms of Random Integer Matrices
Fast Computation of Hermite Normal Forms of Random Integer Matrices Clément Pernet 1 William Stein 2 Abstract This paper is about how to compute the Hermite normal form of a random integer matrix in practice.
More informationBlock Wiedemann likes Schirokauer maps
Block Wiedemann likes Schirokauer maps E. Thomé INRIA/CARAMEL, Nancy. /* EPI CARAMEL */ C,A, /* Cryptologie, Arithmétique : */ R,a, /* Matériel et Logiciel */ M,E, L,i= 5,e, d[5],q[999 ]={0};main(N ){for
More informationSparse Polynomial Interpolation and Berlekamp/Massey algorithms that correct Outlier Errors in Input Values
Sparse Polynomial Interpolation and Berlekamp/Massey algorithms that correct Outlier Errors in Input Values Clément PERNET joint work with Matthew T. COMER and Erich L. KALTOFEN : LIG/INRIA-MOAIS, Grenoble
More informationAlgorithms for Solving Linear Differential Equations with Rational Function Coefficients
Algorithms for Solving Linear Differential Equations with Rational Function Coefficients Erdal Imamoglu Department of Mathematics Florida State University May 25, 2017 1 / 37 1 Introduction 2 Formal Solutions,
More informationTowards High Performance Multivariate Factorization. Michael Monagan. This is joint work with Baris Tuncer.
Towards High Performance Multivariate Factorization Michael Monagan Center for Experimental and Constructive Mathematics Simon Fraser University British Columbia This is joint work with Baris Tuncer. To
More informationREDUNDANT TRINOMIALS FOR FINITE FIELDS OF CHARACTERISTIC 2
REDUNDANT TRINOMIALS FOR FINITE FIELDS OF CHARACTERISTIC 2 CHRISTOPHE DOCHE Abstract. In this paper we introduce so-called redundant trinomials to represent elements of nite elds of characteristic 2. The
More informationKaltofen s division-free determinant algorithm differentiated for matrix adjoint computation
Kaltofen s division-free determinant algorithm differentiated for matrix adjoint computation Gilles Villard Université de Lyon, CNRS, ENS de Lyon, INRIA, UCBL Laboratoire LIP, 46, Allée d Italie, 69364
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 information4. Parallel computation. a) Processor efficient NC and RNC algorithms for matrix computations with extensions. My paper with
RESEARCH HIGHLIGHTS The results of my research have been published in 3 books and over 230 research articles in journals and refereed proceedings of conferences and have been surveyed in my 18 surveys
More informationTURNER, WILLIAM JONATHAN. Black Box Linear Algebra with the LinBox Library.
Abstract TURNER, WILLIAM JONATHAN. Black Box Linear Algebra with the LinBox Library. (Under the direction of Erich Kaltofen.) Black box algorithms for exact linear algebra view a matrix as a linear operator
More informationParallel Algorithms for Computing the Smith Normal Form of Large Matrices
Parallel Algorithms for Computing the Smith Normal Form of Large Matrices Gerold Jäger Mathematisches Seminar der Christian-Albrechts-Universität zu Kiel, Christian-Albrechts-Platz 4, D-24118 Kiel, Germany
More informationFast, deterministic computation of the Hermite normal form and determinant of a polynomial matrix
Fast, deterministic computation of the Hermite normal form and determinant of a polynomial matrix George Labahn, Vincent Neiger, Wei Zhou To cite this version: George Labahn, Vincent Neiger, Wei Zhou.
More informationTowards High Performance Multivariate Factorization. Michael Monagan. This is joint work with Baris Tuncer.
Towards High Performance Multivariate Factorization Michael Monagan Center for Experimental and Constructive Mathematics Simon Fraser University British Columbia This is joint work with Baris Tuncer. The
More informationFraction-free Computation of Matrix Padé Systems
Fraction-free Computation of Matrix Padé Systems Bernhard Becermann, Laboratoire d Analyse Numérique et d Optimisation, UFR IEEA M3, USTL Flres Artois, F 59655 Villeneuve d Ascq CEDEX, France, e-mail:
More informationABSTRACT. YUHASZ, GEORGE L. Berlekamp/Massey Algorithms for Linearly Generated Matrix Sequences. (Under the direction of Professor Erich Kaltofen.
ABSTRACT YUHASZ, GEORGE L. Berlekamp/Massey Algorithms for Linearly Generated Matrix Sequences. (Under the direction of Professor Erich Kaltofen.) The Berlekamp/Massey algorithm computes the unique minimal
More informationFast Estimates of Hankel Matrix Condition Numbers and Numeric Sparse Interpolation*
Fast Estimates of Hankel Matrix Condition Numbers and Numeric Sparse Interpolation* Erich L Kaltofen Department of Mathematics North Carolina State University Raleigh, NC 769-80, USA kaltofen@mathncsuedu
More informationPopov Form Computation for Matrices of Ore Polynomials
Popov Form Computation for Matrices of Ore Polynomials Mohamed Khochtali University of Waterloo Canada mkhochta@uwaterlooca Johan Rosenkilde, né Nielsen Technical University of Denmark Denmark jsrn@jsrndk
More informationComputing the Rank and a Small Nullspace Basis of apolynomial Matrix
Computing the Rank and a Small Nullspace Basis of apolynomial Matrix Arne Storjohann School of Computer Science, U. Waterloo Waterloo Ontario N2L 3G1 Canada http://www.scg.uwaterloo.ca/~astorjoh Gilles
More information9. Numerical linear algebra background
Convex Optimization Boyd & Vandenberghe 9. Numerical linear algebra background matrix structure and algorithm complexity solving linear equations with factored matrices LU, Cholesky, LDL T factorization
More informationMITACS Project CV and Business Plan. Mathematics of Computer Algebra and Analysis (MOCAA)
Labahn and Monagan 1 MITACS Project CV and Business Plan November 1, 2006 Title of Project: Mathematics of Computer Algebra and Analysis (MOCAA) Project website: www.cecm.sfu.ca/~pborwein/mitacs/index.htm
More informationFinding Succinct. Ordered Minimal Perfect. Hash Functions. Steven S. Seiden 3 Daniel S. Hirschberg 3. September 22, Abstract
Finding Succinct Ordered Minimal Perfect Hash Functions Steven S. Seiden 3 Daniel S. Hirschberg 3 September 22, 1994 Abstract An ordered minimal perfect hash table is one in which no collisions occur among
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 informationComputer Algebra and Computer Algebra Systems L algèbre computationnelle et systèmes d algèbre computationnelle (Org: Michael Monagan (SFU))
Computer Algebra and Computer Algebra Systems L algèbre computationnelle et systèmes d algèbre computationnelle (Org: Michael Monagan (SFU)) OLEG GOLUBITSKY, Ontario Research Centre for Computer Algebra,
More informationBOUNDS ON THE COEFFICIENTS OF THE CHARACTERISTIC AND MINIMAL POLYNOMIALS
BOUNDS ON THE COEFFICIENTS OF THE CHARACTERISTIC AND MINIMAL POLYNOMIALS Received: 5 October, 006 Accepted: 3 April, 007 Communicated by: JEAN-GUILLAUME DUMAS Laboratoire Jean Kuntzmann Université Joseph
More informationComputing the sign or the value of the determinant of an integer matrix, a complexity survey 1
Computing the sign or the value of the determinant of an integer matrix, a complexity survey 1 Erich Kaltofen a and Gilles Villard b a Department of Mathematics, North Carolina State University Raleigh,
More informationSubquadratic Space Complexity Multiplication over Binary Fields with Dickson Polynomial Representation
Subquadratic Space Complexity Multiplication over Binary Fields with Dickson Polynomial Representation M A Hasan and C Negre Abstract We study Dickson bases for binary field representation Such representation
More information9. Numerical linear algebra background
Convex Optimization Boyd & Vandenberghe 9. Numerical linear algebra background matrix structure and algorithm complexity solving linear equations with factored matrices LU, Cholesky, LDL T factorization
More informationComputing Characteristic Polynomials of Matrices of Structured Polynomials
Computing Characteristic Polynomials of Matrices of Structured Polynomials Marshall Law and Michael Monagan Department of Mathematics Simon Fraser University Burnaby, British Columbia, Canada mylaw@sfu.ca
More informationEliminations and echelon forms in exact linear algebra
Eliminations and echelon forms in exact linear algebra Clément PERNET, INRIA-MOAIS, Grenoble Université, France East Coast Computer Algebra Day, University of Waterloo, ON, Canada, April 9, 2011 Clément
More informationAdaptive decoding for dense and sparse evaluation/interpolation codes
Adaptive decoding for dense and sparse evaluation/interpolation codes Clément PERNET INRIA/LIG-MOAIS, Grenoble Université joint work with M. Comer, E. Kaltofen, J-L. Roch, T. Roche Séminaire Calcul formel
More informationEfficient Algorithms for Order Bases Computation
Efficient Algorithms for Order Bases Computation Wei Zhou and George Labahn Cheriton School of Computer Science University of Waterloo, Waterloo, Ontario, Canada Abstract In this paper we present two algorithms
More informationDense Linear Algebra over Finite Fields: the FFLAS and FFPACK packages
Dense Linear Algebra over Finite Fields: the FFLAS and FFPACK packages Jean-Guillaume Dumas, Thierry Gautier, Pascal Giorgi, Clément Pernet To cite this version: Jean-Guillaume Dumas, Thierry Gautier,
More informationPrimitive sets in a lattice
Primitive sets in a lattice Spyros. S. Magliveras Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431, U.S.A spyros@fau.unl.edu Tran van Trung Institute for Experimental
More informationWORKING WITH MULTIVARIATE POLYNOMIALS IN MAPLE
WORKING WITH MULTIVARIATE POLYNOMIALS IN MAPLE JEFFREY B. FARR AND ROMAN PEARCE Abstract. We comment on the implementation of various algorithms in multivariate polynomial theory. Specifically, we describe
More informationDiophantine equations via weighted LLL algorithm
Cryptanalysis of a public key cryptosystem based on Diophantine equations via weighted LLL algorithm Momonari Kudo Graduate School of Mathematics, Kyushu University, JAPAN Kyushu University Number Theory
More informationMODEL ANSWERS TO THE FIRST QUIZ. 1. (18pts) (i) Give the definition of a m n matrix. A m n matrix with entries in a field F is a function
MODEL ANSWERS TO THE FIRST QUIZ 1. (18pts) (i) Give the definition of a m n matrix. A m n matrix with entries in a field F is a function A: I J F, where I is the set of integers between 1 and m and J is
More informationNearly Sparse Linear Algebra and application to Discrete Logarithms Computations
Nearly Sparse Linear Algebra and application to Discrete Logarithms Computations Antoine Joux 1,2,3 and Cécile Pierrot 1,4 1 Sorbonne Universités, UPMC Univ Paris 06, LIP6, 4 place Jussieu, 75005 PARIS,
More informationComputing over Z, Q, K[X]
Computing over Z, Q, K[X] Clément PERNET M2-MIA Calcul Exact Outline Introduction Chinese Remainder Theorem Rational reconstruction Problem Statement Algorithms Applications Dense CRT codes Extension to
More informationOn the Arithmetic and Shift Complexities of Inverting of Dierence Operator Matrices
On the Arithmetic and Shift Complexities of Inverting of Dierence Operator Matrices S. Abramov Moscow, Comp. Centre of RAS, MSU S. Abramov the arithmetic and shift complexities 1/12 Let K be a dierence
More informationParallel Complexity of Polynomial GCD and Toeplitz-like... 2 the precision of the computation decreases. The algorithms involve either a single random
Parallel Complexity of Computations with General and Toeplitz-like Matrices Filled with Integers and Extensions Victor Y. Pan Department of Mathematics and Computer Science Lehman College, City University
More informationAsymptotically Fast Polynomial Matrix Algorithms for Multivariable Systems
Laboratoire de l Informatique du Parallélisme École Normale Supérieure de Lyon Unité Mixte de Recherche CNRS-INRIA-ENS LYON-UCBL n o 5668 Asymptotically Fast Polynomial Matrix Algorithms for Multivariable
More informationFast Computation of the Nth Term of an Algebraic Series over a Finite Prime Field
Fast Computation of the Nth Term of an Algebraic Series over a Finite Prime Field Alin Bostan 1 Gilles Christol 2 Philippe Dumas 1 1 Inria (France) 2 Institut de mathématiques de Jussieu (France) Specfun
More informationA modified block Lanczos algorithm with fewer vectors
A modified block Lanczos algorithm with fewer vectors Emmanuel Thomé INRIA Nancy / LORIA 615 rue du jardin botanique 54600 Villers-lès-Nancy Abstract. The block Lanczos algorithm proposed by Peter Montgomery
More informationAlgorithms for Fast Linear System Solving and Rank Profile Computation
Algorithms for Fast Linear System Solving and Rank Profile Computation by Shiyun Yang A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master
More informationInteger factorization, part 1: the Q sieve. D. J. Bernstein
Integer factorization, part 1: the Q sieve D. J. Bernstein Sieving small integers 0 using primes 3 5 7: 1 3 3 4 5 5 6 3 7 7 8 9 3 3 10 5 11 1 3 13 14 7 15 3 5 16 17 18 3 3 19 0 5 etc. Sieving and 611 +
More informationThe following can also be obtained from this WWW address: the papers [8, 9], more examples, comments on the implementation and a short description of
An algorithm for computing the Weierstrass normal form Mark van Hoeij Department of mathematics University of Nijmegen 6525 ED Nijmegen The Netherlands e-mail: hoeij@sci.kun.nl April 9, 1995 Abstract This
More informationA New Polynomial Factorization Algorithm and its Implementation. Victor Shoup. Abstract
A New Polynomial Factorization Algorithm and its Implementation Victor Shoup Bellcore, 445 South St., Morristown, NJ 07960 shoup@bellcore.com Abstract We consider the problem of factoring univariate polynomials
More informationFast algorithms for factoring polynomials A selection. Erich Kaltofen North Carolina State University google->kaltofen google->han lu
Fast algorithms for factoring polynomials A selection Erich Kaltofen North Carolina State University google->kaltofen google->han lu Overview of my work 1 Theorem. Factorization in K[x] is undecidable
More informationFactoring Polynomials with Rational Coecients. Kenneth Giuliani
Factoring Polynomials with Rational Coecients Kenneth Giuliani 17 April 1998 1 Introduction Factorization is a problem well-studied in mathematics. Of particular focus is factorization within unique factorization
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 informationRevised Simplex Method
DM545 Linear and Integer Programming Lecture 7 Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. 2 Motivation Complexity of single pivot operation
More informationGRE Subject test preparation Spring 2016 Topic: Abstract Algebra, Linear Algebra, Number Theory.
GRE Subject test preparation Spring 2016 Topic: Abstract Algebra, Linear Algebra, Number Theory. Linear Algebra Standard matrix manipulation to compute the kernel, intersection of subspaces, column spaces,
More informationLinear Algebra 1: Computing canonical forms in exact linear algebra
Linear Algebra : Computing canonical forms in exact linear algebra Clément PERNET, LIG/INRIA-MOAIS, Grenoble Université, France ECRYPT II: Summer School on Tools, Mykonos, Grece, June st, 202 Introduction
More informationAlgorithms for Structured Linear Systems Solving and their Implementation
Algorithms for Structured Linear Systems Solving and their Implementation Seung Gyu Hyun Cheriton School of Computer Science University of Waterloo sghyun@edu.uwaterloo.ca Romain Lebreton LIRMM Université
More informationThree Ways to Test Irreducibility
Three Ways to Test Irreducibility Richard P. Brent Australian National University joint work with Paul Zimmermann INRIA, Nancy France 12 Feb 2009 Outline Polynomials over finite fields Irreducibility criteria
More informationJ. Symbolic Computation (1995) 11, 1{000 An algorithm for computing an integral basis in an algebraic function eld Mark van Hoeij Department of mathem
J. Symbolic Computation (1995) 11, 1{000 An algorithm for computing an integral basis in an algebraic function eld Mark van Hoeij Department of mathematics University of Nijmegen 6525 ED Nijmegen The Netherlands
More informationThree Ways to Test Irreducibility
Outline Three Ways to Test Irreducibility Richard P. Brent Australian National University joint work with Paul Zimmermann INRIA, Nancy France 8 Dec 2008 Polynomials over finite fields Irreducibility criteria
More informationOutline Introduction: Problem Description Diculties Algebraic Structure: Algebraic Varieties Rank Decient Toeplitz Matrices Constructing Lower Rank St
Structured Lower Rank Approximation by Moody T. Chu (NCSU) joint with Robert E. Funderlic (NCSU) and Robert J. Plemmons (Wake Forest) March 5, 1998 Outline Introduction: Problem Description Diculties Algebraic
More informationBon anniversaire Marc!
Bon anniversaire Marc! The full counting function of Gessel walks is algebraic Alin Bostan (INRIA) joint work with Manuel Kauers (INRIA/RISC) Journées GECKO/TERA 2008 Context: restricted lattice walks
More informationSubquadratic computation of vector generating polynomials and improvement of the block Wiedemann algorithm
Subquadratic computation of vector generating polynomials and improvement of the block Wiedemann algorithm Emmanuel Thomé To cite this version: Emmanuel Thomé. Subquadratic computation of vector generating
More informationExploiting off-diagonal rank structures in the solution of linear matrix equations
Stefano Massei Exploiting off-diagonal rank structures in the solution of linear matrix equations Based on joint works with D. Kressner (EPFL), M. Mazza (IPP of Munich), D. Palitta (IDCTS of Magdeburg)
More informationThe Orbit Problem in Zero and One Dimensions
The Orbit Problem in Zero and One Dimensions Master's dissertation of Ventsislav K. Chonev Supervised by Joël Ouaknine and James Worrell 1 Contents Contents 2 1 Introduction 3 1.1 Verication of linear
More informationhal , version 1-29 Feb 2008
Compressed Modular Matrix Multiplication Jean-Guillaume Dumas Laurent Fousse Bruno Salvy February 29, 2008 Abstract We propose to store several integers modulo a small prime into a single machine word.
More informationComputing the modular equation
Computing the modular equation Andrew V. Sutherland (MIT) Barcelona-Boston-Tokyo Number Theory Seminar in Memory of Fumiyuki Momose Andrew V. Sutherland (MIT) Computing the modular equation 1 of 8 The
More informationComputing Minimal Nullspace Bases
Computing Minimal Nullspace ases Wei Zhou, George Labahn, and Arne Storjohann Cheriton School of Computer Science University of Waterloo, Waterloo, Ontario, Canada {w2zhou,glabahn,astorjoh}@uwaterloo.ca
More information2 EBERHARD BECKER ET AL. has a real root. Thus our problem can be reduced to the problem of deciding whether or not a polynomial in one more variable
Deciding positivity of real polynomials Eberhard Becker, Victoria Powers, and Thorsten Wormann Abstract. We describe an algorithm for deciding whether or not a real polynomial is positive semidenite. The
More informationPolynomial functions over nite commutative rings
Polynomial functions over nite commutative rings Balázs Bulyovszky a, Gábor Horváth a, a Institute of Mathematics, University of Debrecen, Pf. 400, Debrecen, 4002, Hungary Abstract We prove a necessary
More informationExact Arithmetic on a Computer
Exact Arithmetic on a Computer Symbolic Computation and Computer Algebra William J. Turner Department of Mathematics & Computer Science Wabash College Crawfordsville, IN 47933 Tuesday 21 September 2010
More informationA Fraction Free Matrix Berlekamp/Massey Algorithm*
Linear Algebra and Applications, 439(9):2515-2526, November 2013. A Fraction Free Matrix Berlekamp/Massey Algorithm* Erich Kaltofen and George Yuhasz Dept. of Mathematics, North Carolina State University,
More informationFast and Small: Multiplying Polynomials without Extra Space
Fast and Small: Multiplying Polynomials without Extra Space Daniel S. Roche Symbolic Computation Group School of Computer Science University of Waterloo CECM Day SFU, Vancouver, 24 July 2009 Preliminaries
More informationMICA Milestones in Computer Algebra. Celebrating the research of Erich Kaltofen. July 16 18, 2016 University of Waterloo Canada
MICA 2016 Milestones in Computer Algebra Celebrating the research of Erich Kaltofen July 16 18, 2016 University of Waterloo Canada All talks are in DC 1302, the Davis Centre. Rooted in the boundary area
More informationFaster Inversion and Other Black Box Matrix Computations Using Efficient Block Projections
Faster Inversion and Other Black Box Matrix Comptations Using Efficient Block Projections Wayne Eberly 1, Mark Giesbrecht, Pascal Giorgi,, Arne Storjohann, Gilles Villard (1) Department of Compter Science,
More informationM4. Lecture 3. THE LLL ALGORITHM AND COPPERSMITH S METHOD
M4. Lecture 3. THE LLL ALGORITHM AND COPPERSMITH S METHOD Ha Tran, Dung H. Duong, Khuong A. Nguyen. SEAMS summer school 2015 HCM University of Science 1 / 31 1 The LLL algorithm History Applications of
More information