Inverting integer and polynomial matrices. Jo60. Arne Storjohann University of Waterloo
|
|
- Madeleine Charla Burns
- 6 years ago
- Views:
Transcription
1 Inverting integer and polynomial matrices Jo60 Arne Storjohann University of Waterloo
2 Integer Matrix Inverse Input: An n n matrix A filled with entries of size d digits. Output: The matrix inverse A 1. Example: 6 6 with d = = input A has total size O(n 2 d) words output A 1 has total size O (n 3 d) words 1
3 Polynomial Matrix Inverse Input: An n n matrix A K[x] of degree bounded by d. Output: The matrix inverse A 1. Example: 3 3 over Z 7 [x] with d = 3 6x 2 + 4x+2 2x 2 + 5x+3 2x 2 + x 5x 2 + 3x+2 6x 2 + 6x+4 4x x 2 + 3x 5x 2 + 4x+3 5x = 2x 4 +3x 2 +6x+2 x 6 +4x 5 +5x 4 +6x 3 +6x 2 +4x+2 6x 3 +4x 2 +3 x 5 +5x 4 +3x 3 +2x 2 +x+5 6x 4 +3x 3 +4x 2 +x+4 x 6 +4x 5 +5x 4 +6x 3 +6x 2 +4x+2 6x 3 +3x 2 +5x+6 x 6 +4x 5 +5x 4 +6x 3 +6x 2 +4x+2 4x 3 +3x 2 +5x+4 x 5 +5x 4 +3x 3 +2x 2 +x+5 3x 4 +3x 3 +2x+3 x 6 +4x 5 +5x 4 +6x 3 +6x 2 +4x+2 2x 4 +6x 3 +6x 2 +4x+4 x 6 +4x 5 +5x 4 +6x 3 +6x 2 +4x+2 6x 2 +2x+5 x 5 +5x 4 +3x 3 +2x 2 +x+5 x 4 +3x 3 +x 2 +6x+6 x 6 +4x 5 +5x 4 +6x 3 +6x 2 +4x+2 input A has total size O(n 2 d) field elements from K output A 1 has total size O(n 3 d) 2
4 Notation nonsingular input A has dimension n K[x]: size of entries is d dega count field operations fromk Z: size of entries is d log A count word operations ω is a feasible exponent for matrix multiplication: 2 ω 3 Classical inversion algorithms Hensel-Newton iteration: p-adic lifting Homomorphic imaging and Chinese remaindering Classical Goal O (n ω+1 d) O (n 3 d) 3
5 Outline: Three completely different approaches for inversion A1 Sparse inverse expansion useful representation of A 1 mod p n applicable for integer and polynomial inputs A2 Balanced kernel basis nearly optimal inversion algorithm for generic polynomials inputs [Jeannerod & Villard, 2005] A3 Outer product adjoint formula nearly optimal inversion algorithm for all polynomial inputs well conditioned integer inputs 4
6 Alternative representations for the inverse of A = [ 3 2 x 5 ] A 1 = = [ 2 6+2x x 6+2x 2 6+2x 4 6+2x ] [ 1+x+x 2 + x 3 + 1xx 2 x 3 + xx 2 x x+x 2 + x 3 + ] = [ ] + [ ] x+ [ ] x 2 +. A 1 over K[[x]] modulo x 2nd sufficient to reconstruct A 1 over K(x) compute expansion over K[[x]] using quadratic Newton iteration 5
7 The sparse inverse formula for A= [ 1cx ] Explicit inverse modulo x 32 (1cx) 1 mod x 32 1+cx+c 2 x 2 + c 3 x 3 + c 4 x 4 + c 5 x 5 + +c 31 x 31 6
8 The sparse inverse formula for A= [ 1cx ] Explicit inverse modulo x 32 (1cx) 1 mod x 32 1+cx+c 2 x 2 + c 3 x 3 + c 4 x 4 + c 5 x 5 + +c 31 x 31 (1+cx)(1+c 2 x 2 )(1+c 4 x 4 )(1+c 8 x 8 )(1+c 16 x 16 ) 7
9 The sparse inverse formula for A= [ 1cx ] Explicit inverse modulo x 32 (1cx) 1 mod x 32 = 32 coefficients {}}{ 1+cx+c 2 x 2 + c 3 x 3 + c 4 x 4 + c 5 x 5 + +c 31 x 31 = (1+cx)(1+c 2 x 2 )(1+c 4 x 4 )(1+c 8 x 8 )(1+c 16 x 16 ) }{{} 5 coefficients R 4 {}}{ c 4 x 4 ) R {}}{ 1 {}}{ R 2 } (1+ c x)(1+ {{ c 2 x 2 )(1+ } (1+ (1cx) 1 mod x 8 R 8 {}}{ c 8 x 8 )(1+ R 16 {}}{ c 16 x 16 ) 8
10 Computing the residues for A=1cx R 2 x 2 = I A(1+cx) = c 2 x 2 { A 1 mod }} x 4 { R 4 x 4 = I A( 1+cx+c 2 x 2 + c 3 x 3 ) = c 4 x 4 { A 1 mod }} x 8 { R 8 x 8 = I A( 1+cx+c 2 x 2 + c 3 x 3 + c 4 x 4 + c 5 x 5 + c 6 x 6 + c 7 x 7 ) = c 8 x 8 9
11 Computing the residues for A=1cx R 2 x 2 = I A(1+cx) = c 2 x 2 { A 1 mod }} x 4 { R 4 x 4 = I A( 1+cx+c 2 x 2 + c 3 x 3 ) = c 4 x 4 { A 1 mod }} x 8 { R 8 x 8 = I A( 1+cx+c 2 x 2 + c 3 x 3 + c 4 x 4 + c 5 x 5 + c 6 x 6 + c 7 x 7 ) = c 8 x 8 If dega=d then R 2,R 4,R 8, have degree d 1. Can use reverse short product need only top d coefficients of A 1 mod x 10
12 Standard quadratic lifting High-order component lifting
13 Sparse inverse formula for integer and polynomial matrices Example for 7 1 mod x 32, x= mod x 32 = = (((52(13x 2 )+2x 4 )(15x 4 )+4x 8 )(13x 8 )+2x 16 )(15x 16 )+4x 32 General case formula can be computed in time O (n ω d) sparse inverse formula for A 1 mod x n has size O(n 2 (logn)d) can be used to compute A 1 b in time O (n ω d) fast algorithms for linear algebra over Z and K[x] randomized Las Vegas algorithm for explicit inverse? 12
14 Outline: Three completely different approaches for inversion A1 Sparse inverse expansion useful representation of A 1 mod p n applicable for integer and polynomial inputs A2 Balanced kernel basis nearly optimal inversion algorithm for generic polynomials inputs [Jeannerod & Villard, 2005] A3 Outer product adjoint formula nearly optimal inversion algorithm for all polynomial inputs well conditioned integer inputs 13
15 Minimal kernel basis of generic polynomial matrices Example 1: A K[x] 2 1,K=Z 97 [ N [ ] A ] f g f = 0 g 1 1 Example 2: A K[x] 8 4,K=Z 97 N x x 45x+51 51x+44 33x x x+92 67x+90 83x+29 47x x x+64 42x+38 12x+67 67x x+78 31x+50 4x+52 93x+28 17x+52 A 60x+16 52x20 4x89 77x+69 64x+89 16x+59 69x46 33x x+36 91x22 51x27 31x27 65x+88 10x6 80x84 26x51 88x+97 67x+58 29x+37 9x91 81x+65 12x+78 5x63 42x+9 21x27 79x22 51x+16 95x+86 97x14 83x96 8x54 =0 4 4 generic A K[x] 2n n of degree d has kernel of degree d effective algorithms, e.g., [Giorgi, Jeannerod & Villard, ISSAC 03] 14
16 From generic kernel to generic inversion split input A K[x] n n of degree d into two halves: A= [ A L A R ] compute minimal kernel U for A L and U for A R [ U U ] [ AL A R ] = [ UAL UA R ] recurse on degree 2d matrices UA L and UA R Example: n=8 and d = 1 N 3 N 2 N 1 A [ AL A R ] =(deta)a 1 15
17 From generic kernel to generic inversion N 3 N 2 N 1 A cost to compute N 1,N 2,...,N log2 n is O (n ω d)) cost to compute (deta)a 1 = N logn N 2 N 1 is O (n 3 d) [ AL A R ] =(deta)a 1 Notes requires generic input seems to require n to be a power of 2 16
18 Outline: Three completely different approaches for inversion A1 Sparse inverse expansion useful representation of A 1 mod p n applicable for integer and polynomial inputs A2 Balanced kernel basis nearly optimal inversion algorithm for generic polynomials inputs [Jeannerod & Villard, 2005] A3 Outer product adjoint formula nearly optimal inversion algorithm for all polynomial inputs well conditioned integer inputs 17
19 Special case of the outer product adjoint formula Example input matrix A= , deta= The adjoint A adj = [ ] mod
20 General case of the output product adjoint formla suppose Smith form of A is Diag(s 1,s 2,...,s n ) A adj deta v n u n + deta v n1 u n1 + + deta v 1 u 1 (mod deta) s n s n1 s 1 Integer matrices only obtain inverse modulo deta inverse of A Z n n in time O (n 3 d+ n 3 logκ(a)) nearly optimal for well conditioned input randomized Las Vegas Polynomial matrices use randomization and reversion to extend to arbitrary input inverse of any A K[x] n n in time O (n 3 d) operations fromk randomized Las Vegas 19
21 Applications of the nearly optimal inversion formulas 20
22 Example: Linear solving over K: K=Z 7, n=5 A= b= deta=3 A 1 b= Note: no expression swell over Z 7 21
23 Example input over K[x]: n=5, d = 2 = expression swell A= 6x 2 + x+5 6x+1 2x 2 x+4 3x 2 + 6x+5 3x+5 2x 2 + x+2 6x 2 + 2x+3 3x 2 + 2x+1 x 2 + 2x 5x 2 + x+4 6x 2 + 3x+6 5x 2 + 2x 5x 2 + 6x+4 3x 2 + x 6x 2 + 6x+5 2x 2 + 3x+2 3x 2 + 6x 3x 2 + 2x+2 5x 2 + 5x+4 3x 2 + 2x+5 5x 2 + 3x+1 4x+1 5x 2 + 4x 2x 2 + 3x+2 b= x 2 + x+1 2x 2 + 2x+6 x 2 + 3x+2 2x 2 + x+4 x 2 + 4x+3 deta=6x x 9 + x 8 + 3x 6 + x 5 + x 4 + 4x 2 + 2x+3 Degree is n d where d is degrees in input matrix. x x 9 + 2x 8 + x 7 + 3x 6 + x 4 + 6x 3 + 3x 5x 9 + x 8 + 6x 7 + x 5 + 4x 4 + 5x 3 + x 2 + 3x+1 A 1 b= 6x 10 + x 9 + 3x 8 + 6x 7 + 3x 6 + 2x 5 + 2x 3 + 5x x x 9 + 4x 8 + 3x 7 + 2x 6 + 3x 5 + 5x 4 + 6x 3x x 9 + 4x 7 + 5x 6 + 2x 5 + x+6 (1/detA) 22
24 Nearly optimal linear system solving with precomputation Recall situation over K Input: Nonsingular A K n n Precompute: Inverse A 1 Compute: A 1 v for any given v K n 1 A b = A 1 b size of A 1 b is Ω(n) time to compute A 1 b is O(n 2 ) 23
25 Nearly optimal linear system solving with precomputation over K[x] Input: nonsingular A K[x] n n of degree d Precompute: adjoint decomposition N logn N 2 N 1 =(deta)a 1. Compute: A 1 v for any given v K n 1 of degree d. N 3 N 2 N 1 b [1] [1] [1] [1] [1] [1] [1] [1] =(deta)a 1 b size of A 1 b is Ω(n 2 d) cost to compute A 1 b is O (n 2 d) result for nongeneric input via outer product adjoint formula 24
26 Krylov basis and matrix powers Input: matrix A K n n and vector v K n 1 Output: v,av,a 2,A 3 v,...,a n1 v Standard algorithm for matrix vector interates compute Av,A(Av),A(A 2 v),a(a 3 v),...,a(a n2 v) requires n 1 matrix-vector products total cost is O(n 3 ) operations from K Fast algorithm [Keller-Gehrig, 1985] precompute A 2,(A 2 ) 2,(A 4 ) 2,..., compute iterates via matrix multiplication total cost is O(n w logn) operations fromk 25
27 Nearly optimal computation of matrix powers Input: matrix A K n n Output: I,A,A 2,A 3,A 4,...,A n1 consider matrix I xa K[x] compute(i xa) 1 mod x n = I+ xa+x 2 A 2 + x 3 A 3 + x n1 A n1 total size of output is Θ(n 3 ) elements of K total cost is O (n 3 ) using fast inversion algorithm 26
Lattice 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 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 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 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 informationComputing with polynomials: Hensel constructions
Course Polynomials: Their Power and How to Use Them, JASS 07 Computing with polynomials: Hensel constructions Lukas Bulwahn March 28, 2007 Abstract To solve GCD calculations and factorization of polynomials
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 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 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 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 informationThe shifted number system for fast linear algebra on integer matrices
The shifted number system for fast linear algebra on integer matrices Arne Storjohann 1 School of Computer Science, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1 Abstract The shifted number
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 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 informationHASSE-MINKOWSKI THEOREM
HASSE-MINKOWSKI THEOREM KIM, SUNGJIN 1. Introduction In rough terms, a local-global principle is a statement that asserts that a certain property is true globally if and only if it is true everywhere locally.
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 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 informationAn introduction to the algorithmic of p-adic numbers
An introduction to the algorithmic of p-adic numbers David Lubicz 1 1 Universté de Rennes 1, Campus de Beaulieu, 35042 Rennes Cedex, France Outline Introduction 1 Introduction 2 3 4 5 6 7 8 When do we
More informationDeterminants Chapter 3 of Lay
Determinants Chapter of Lay Dr. Doreen De Leon Math 152, Fall 201 1 Introduction to Determinants Section.1 of Lay Given a square matrix A = [a ij, the determinant of A is denoted by det A or a 11 a 1j
More informationENGR-1100 Introduction to Engineering Analysis. Lecture 21. Lecture outline
ENGR-1100 Introduction to Engineering Analysis Lecture 21 Lecture outline Procedure (algorithm) for finding the inverse of invertible matrix. Investigate the system of linear equation and invertibility
More informationENGR-1100 Introduction to Engineering Analysis. Lecture 21
ENGR-1100 Introduction to Engineering Analysis Lecture 21 Lecture outline Procedure (algorithm) for finding the inverse of invertible matrix. Investigate the system of linear equation and invertibility
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 informationZEROES OF INTEGER LINEAR RECURRENCES. 1. Introduction. 4 ( )( 2 1) n
ZEROES OF INTEGER LINEAR RECURRENCES DANIEL LITT Consider the integer linear recurrence 1. Introduction x n = x n 1 + 2x n 2 + 3x n 3 with x 0 = x 1 = x 2 = 1. For which n is x n = 0? Answer: x n is never
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 informationTC10 / 3. Finite fields S. Xambó
TC10 / 3. Finite fields S. Xambó The ring Construction of finite fields The Frobenius automorphism Splitting field of a polynomial Structure of the multiplicative group of a finite field Structure of the
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 informationAlgebra Homework, Edition 2 9 September 2010
Algebra Homework, Edition 2 9 September 2010 Problem 6. (1) Let I and J be ideals of a commutative ring R with I + J = R. Prove that IJ = I J. (2) Let I, J, and K be ideals of a principal ideal domain.
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 informationFast Polynomial Multiplication
Fast Polynomial Multiplication Marc Moreno Maza CS 9652, October 4, 2017 Plan Primitive roots of unity The discrete Fourier transform Convolution of polynomials The fast Fourier transform Fast convolution
More informationComputational Methods CMSC/AMSC/MAPL 460. Eigenvalues and Eigenvectors. Ramani Duraiswami, Dept. of Computer Science
Computational Methods CMSC/AMSC/MAPL 460 Eigenvalues and Eigenvectors Ramani Duraiswami, Dept. of Computer Science Eigen Values of a Matrix Recap: A N N matrix A has an eigenvector x (non-zero) with corresponding
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 informationCS 4424 GCD, XGCD
CS 4424 GCD, XGCD eschost@uwo.ca GCD of polynomials First definition Let A and B be in k[x]. k[x] is the ring of polynomials with coefficients in k A Greatest Common Divisor of A and B is a polynomial
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 informationALGEBRA PH.D. QUALIFYING EXAM September 27, 2008
ALGEBRA PH.D. QUALIFYING EXAM September 27, 2008 A passing paper consists of four problems solved completely plus significant progress on two other problems; moreover, the set of problems solved completely
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 informationMath 312/ AMS 351 (Fall 17) Sample Questions for Final
Math 312/ AMS 351 (Fall 17) Sample Questions for Final 1. Solve the system of equations 2x 1 mod 3 x 2 mod 7 x 7 mod 8 First note that the inverse of 2 is 2 mod 3. Thus, the first equation becomes (multiply
More informationMathematics and Algorithms for Computer Algebra. 7: Introduction to modular and p-adic methods
Mathematics and Algorithms for Computer Algebra Part 1 c 1992 Dr Francis J. Wright CBPF, Rio de Janeiro July 9, 2003 7: Introduction to modular and p-adic methods Modular and p-adic methods are often presented
More informationHW2 Solutions
8.024 HW2 Solutions February 4, 200 Apostol.3 4,7,8 4. Show that for all x and y in real Euclidean space: (x, y) 0 x + y 2 x 2 + y 2 Proof. By definition of the norm and the linearity of the inner product,
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 informationMA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES
MA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES 2018 57 5. p-adic Numbers 5.1. Motivating examples. We all know that 2 is irrational, so that 2 is not a square in the rational field Q, but that we can
More informationHermite Forms of Polynomial Matrices
Hermite Forms of Polynomial Matrices by Somit Gupta A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master of Mathematics in Computer Science
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationarithmetic properties of weighted catalan numbers
arithmetic properties of weighted catalan numbers Jason Chen Mentor: Dmitry Kubrak May 20, 2017 MIT PRIMES Conference background: catalan numbers Definition The Catalan numbers are the sequence of integers
More informationSolving Third Order Linear Differential Equations in Terms of Second Order Equations
Solving Third Order Linear Differential Equations in Terms of Second Order Equations Mark van Hoeij (Florida State University) ISSAC 2007 Talk presented by: George Labahn (University of Waterloo) Notations.
More informationQuestion: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI?
Section 5. The Characteristic Polynomial Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Property The eigenvalues
More informationEquality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.
Introduction Matrix Operations Matrix: An m n matrix A is an m-by-n array of scalars from a field (for example real numbers) of the form a a a n a a a n A a m a m a mn The order (or size) of A is m n (read
More informationCounting points on curves: the general case
Counting points on curves: the general case Jan Tuitman, KU Leuven October 14, 2015 Jan Tuitman, KU Leuven Counting points on curves: the general case October 14, 2015 1 / 26 Introduction Algebraic curves
More informationp-adic fields Chapter 7
Chapter 7 p-adic fields In this chapter, we study completions of number fields, and their ramification (in particular in the Galois case). We then look at extensions of the p-adic numbers Q p and classify
More informationGeneralization of Hensel lemma: nding of roots of p-adic Lipschitz functions
Generalization of Hensel lemma: nding of roots of p-adic Lipschitz functions (joint talk with Andrei Khrennikov) Dr. Ekaterina Yurova Axelsson Linnaeus University, Sweden September 8, 2015 Outline Denitions
More informationRings in Coding Theory
Rings in Coding Theory Steven T. Dougherty July 3, 2013 Cyclic Codes Cyclic Codes were first studied by Prange in 1957. Prange, E. Cyclic error-correcting codes in two symbols. Technical Note TN-57-103,
More information235 Final exam review questions
5 Final exam review questions Paul Hacking December 4, 0 () Let A be an n n matrix and T : R n R n, T (x) = Ax the linear transformation with matrix A. What does it mean to say that a vector v R n is an
More informationFFT: Fast Polynomial Multiplications
FFT: Fast Polynomial Multiplications Jie Wang University of Massachusetts Lowell Department of Computer Science J. Wang (UMass Lowell) FFT: Fast Polynomial Multiplications 1 / 20 Overview So far we have
More informationA Course in Computational Algebraic Number Theory
Henri Cohen 2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. A Course in Computational Algebraic Number Theory Springer
More informationPart II. Number Theory. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 Paper 3, Section I 1G 70 Explain what is meant by an Euler pseudoprime and a strong pseudoprime. Show that 65 is an Euler
More informationCPSC 518 Introduction to Computer Algebra Schönhage and Strassen s Algorithm for Integer Multiplication
CPSC 518 Introduction to Computer Algebra Schönhage and Strassen s Algorithm for Integer Multiplication March, 2006 1 Introduction We have now seen that the Fast Fourier Transform can be applied to perform
More informationAlgebraic number theory Solutions to exercise sheet for chapter 4
Algebraic number theory Solutions to exercise sheet for chapter 4 Nicolas Mascot n.a.v.mascot@warwick.ac.uk), Aurel Page a.r.page@warwick.ac.uk) TAs: Chris Birkbeck c.d.birkbeck@warwick.ac.uk), George
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 informationSolving 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 information6.S897 Algebra and Computation February 27, Lecture 6
6.S897 Algebra and Computation February 7, 01 Lecture 6 Lecturer: Madhu Sudan Scribe: Mohmammad Bavarian 1 Overview Last lecture we saw how to use FFT to multiply f, g R[x] in nearly linear time. We also
More information= 1 2x. x 2 a ) 0 (mod p n ), (x 2 + 2a + a2. x a ) 2
8. p-adic numbers 8.1. Motivation: Solving x 2 a (mod p n ). Take an odd prime p, and ( an) integer a coprime to p. Then, as we know, x 2 a (mod p) has a solution x Z iff = 1. In this case we can suppose
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
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 informationMATH 211 review notes
MATH 211 review notes Notes written by Mark Przepiora 1 Determinants 1. Let A = 3 1 1 2 1 2 5 2 1 (a) Compute det(a). Solution: To find the determinant, we expand along the first row. det(a) = 2 1 1 2
More informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 22, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M A and MSc Scholarship Test September 22, 2018 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions that follow INSTRUCTIONS TO CANDIDATES
More informationFactoring Algorithms Pollard s p 1 Method. This method discovers a prime factor p of an integer n whenever p 1 has only small prime factors.
Factoring Algorithms Pollard s p 1 Method This method discovers a prime factor p of an integer n whenever p 1 has only small prime factors. Input: n (to factor) and a limit B Output: a proper factor of
More informationProblem Set 9 Due: In class Tuesday, Nov. 27 Late papers will be accepted until 12:00 on Thursday (at the beginning of class).
Math 3, Fall Jerry L. Kazdan Problem Set 9 Due In class Tuesday, Nov. 7 Late papers will be accepted until on Thursday (at the beginning of class).. Suppose that is an eigenvalue of an n n matrix A and
More information1 Determinants. 1.1 Determinant
1 Determinants [SB], Chapter 9, p.188-196. [SB], Chapter 26, p.719-739. Bellow w ll study the central question: which additional conditions must satisfy a quadratic matrix A to be invertible, that is to
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
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 informationComputer Algebra: General Principles
Computer Algebra: General Principles For article on related subject see SYMBOL MANIPULATION. Computer algebra is a branch of scientific computation. There are several characteristic features that distinguish
More informationFast Convolution; Strassen s Method
Fast Convolution; Strassen s Method 1 Fast Convolution reduction to subquadratic time polynomial evaluation at complex roots of unity interpolation via evaluation at complex roots of unity 2 The Master
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
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 informationZEROS OF SPARSE POLYNOMIALS OVER LOCAL FIELDS OF CHARACTERISTIC p
ZEROS OF SPARSE POLYNOMIALS OVER LOCAL FIELDS OF CHARACTERISTIC p BJORN POONEN 1. Statement of results Let K be a field of characteristic p > 0 equipped with a valuation v : K G taking values in an ordered
More informationMathematics for Cryptography
Mathematics for Cryptography Douglas R. Stinson David R. Cheriton School of Computer Science University of Waterloo Waterloo, Ontario, N2L 3G1, Canada March 15, 2016 1 Groups and Modular Arithmetic 1.1
More informationHomework 11/Solutions. (Section 6.8 Exercise 3). Which pairs of the following vector spaces are isomorphic?
MTH 9-4 Linear Algebra I F Section Exercises 6.8,4,5 7.,b 7.,, Homework /Solutions (Section 6.8 Exercise ). Which pairs of the following vector spaces are isomorphic? R 7, R, M(, ), M(, 4), M(4, ), P 6,
More informationMATH 310: Homework 7
1 MATH 310: Homework 7 Due Thursday, 12/1 in class Reading: Davenport III.1, III.2, III.3, III.4, III.5 1. Show that x is a root of unity modulo m if and only if (x, m 1. (Hint: Use Euler s theorem and
More informationMath 504 (Fall 2011) 1. (*) Consider the matrices
Math 504 (Fall 2011) Instructor: Emre Mengi Study Guide for Weeks 11-14 This homework concerns the following topics. Basic definitions and facts about eigenvalues and eigenvectors (Trefethen&Bau, Lecture
More informationCPSC 518 Introduction to Computer Algebra Asymptotically Fast Integer Multiplication
CPSC 518 Introduction to Computer Algebra Asymptotically Fast Integer Multiplication 1 Introduction We have now seen that the Fast Fourier Transform can be applied to perform polynomial multiplication
More informationFast Direct Methods for Gaussian Processes
Fast Direct Methods for Gaussian Processes Mike O Neil Departments of Mathematics New York University oneil@cims.nyu.edu December 12, 2015 1 Collaborators This is joint work with: Siva Ambikasaran Dan
More informationFast reversion of power series
Fast reversion of power series Fredrik Johansson November 2011 Overview Fast power series arithmetic Fast composition and reversion (Brent and Kung, 1978) A new algorithm for reversion Implementation results
More informationIntroduction to Arithmetic Geometry Fall 2013 Lecture #2 09/10/2013
18.78 Introduction to Arithmetic Geometry Fall 013 Lecture # 09/10/013.1 Plane conics A conic is a plane projective curve of degree. Such a curve has the form C/k : ax + by + cz + dxy + exz + fyz with
More informationSection 13 Homomorphisms
Section 13 Homomorphisms Instructor: Yifan Yang Fall 2006 Homomorphisms Definition A map φ of a group G into a group G is a homomorphism if for all a, b G. φ(ab) = φ(a)φ(b) Examples 1. Let φ : G G be defined
More informationQuasi-reducible Polynomials
Quasi-reducible Polynomials Jacques Willekens 06-Dec-2008 Abstract In this article, we investigate polynomials that are irreducible over Q, but are reducible modulo any prime number. 1 Introduction Let
More informationLecture 6: Deterministic Primality Testing
Lecture 6: Deterministic Primality Testing Topics in Pseudorandomness and Complexity (Spring 018) Rutgers University Swastik Kopparty Scribe: Justin Semonsen, Nikolas Melissaris 1 Introduction The AKS
More informationMA4H9 Modular Forms: Problem Sheet 2 Solutions
MA4H9 Modular Forms: Problem Sheet Solutions David Loeffler December 3, 010 This is the second of 3 problem sheets, each of which amounts to 5% of your final mark for the course This problem sheet will
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 25, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M. A. and M.Sc. Scholarship Test September 25, 2010 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions on the following page 1 INSTRUCTIONS
More informationEigenvalues and Eigenvectors
Contents Eigenvalues and Eigenvectors. Basic Concepts. Applications of Eigenvalues and Eigenvectors 8.3 Repeated Eigenvalues and Symmetric Matrices 3.4 Numerical Determination of Eigenvalues and Eigenvectors
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 informationDynamics of finite linear cellular automata over Z N
Dynamics of finite linear cellular automata over Z N F. Mendivil, D. Patterson September 9, 2009 Abstract We investigate the behaviour of linear cellular automata with state space Z N and only finitely
More information12x + 18y = 30? ax + by = m
Math 2201, Further Linear Algebra: a practical summary. February, 2009 There are just a few themes that were covered in the course. I. Algebra of integers and polynomials. II. Structure theory of one endomorphism.
More informationParallel Numerical Algorithms
Parallel Numerical Algorithms Chapter 13 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign CS 554 / CSE 512 Michael T. Heath Parallel Numerical Algorithms
More informationMATHEMATICS 217: HONORS LINEAR ALGEBRA Spring Term, First Week, February 4 8
MATHEMATICS 217: HONORS LINEAR ALGEBRA Spring Term, 2002 The textbook for the course is Linear Algebra by Hoffman and Kunze (second edition, Prentice-Hall, 1971). The course probably will cover most of
More information3.4. ZEROS OF POLYNOMIAL FUNCTIONS
3.4. ZEROS OF POLYNOMIAL FUNCTIONS What You Should Learn Use the Fundamental Theorem of Algebra to determine the number of zeros of polynomial functions. Find rational zeros of polynomial functions. Find
More informationax b mod m. has a solution if and only if d b. In this case, there is one solution, call it x 0, to the equation and there are d solutions x m d
10. Linear congruences In general we are going to be interested in the problem of solving polynomial equations modulo an integer m. Following Gauss, we can work in the ring Z m and find all solutions to
More informationJORDAN NORMAL FORM NOTES
18.700 JORDAN NORMAL FORM NOTES These are some supplementary notes on how to find the Jordan normal form of a small matrix. First we recall some of the facts from lecture, next we give the general algorithm
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 JORDAN-CHEVALLEY DECOMPOSITION
THE JORDAN-CHEVALLEY DECOMPOSITION JOO HEON YOO Abstract. This paper illustrates the Jordan-Chevalley decomposition through two related problems. Contents 1. Introduction 1 2. Linear Algebra Review 1 3.
More informationMATH 5524 MATRIX THEORY Problem Set 4
MATH 5524 MATRIX THEORY Problem Set 4 Posted Tuesday 28 March 217. Due Tuesday 4 April 217. [Corrected 3 April 217.] [Late work is due on Wednesday 5 April.] Complete any four problems, 25 points each.
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 information