Solving Sparse Integer Linear Systems. Pascal Giorgi. Dali team - University of Perpignan (France)
|
|
- Alexandrina Tamsyn Payne
- 5 years ago
- Views:
Transcription
1 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), G. Villard (ENS Lyon) LMC - MOSAIC Seminar, June 2, 26
2 Motivations Large linear systems are involved in many mathematical alications over a eld : I integers factorization [Odlyzko 999], I discrete logarithm [Odlyzko 999 ; Thome 23], over the integers : I number theory [Cohen 993], I grou theory [Newman 972], I integer rogramming [Aardal, Hurkens, Lenstra 999]
3 Problem Let A a non-singular matrix and b a vector dened over Z. Problem : Comute x = A b over the rational numbers A = x = A b = C A, b = C A C A. Main diculty : exression swell
4 Problem Let A a non-singular matrix and b a vector dened over Z. Problem : Comute x = A b over the rational numbers A = x = A b = C A C A, b = C A. Main diculty : exression swell and take advantage of sarsity
5 Interest in linear algebra Integer linear systems are central in recent linear algebra algorithms I Determinant [Abbott, Bronstein, Mulders 999 ; Storjohann 25] I Smith Form [Eberly, Giesbrecht, Villard 2] I Nullsace, Kernel [Chen, Storjohann 25] I Diohantine solutions [Giesbrecht 997 ; Giesbrecht, Lobo, Saunders 998 ; Mulders, Storjohann 23 ; Mulders 24]
6 Algorithms for non-singular system solving Dense matrices : I Gaussian elimination and CRA,! O~(n!+ log jjajj) bit oerations I P-adic lifting [Monck, Carter 979 ; Dixon 982],! O~(n 3 log jjajj) bit oerations I High order lifting [Storjohann 25],! O~(n! log jjajj) bit oerations Sarse matrices : I P-adic lifting or CRA [Wiedemann 986 ; Kaltofen, Saunders 99],! O(n 2 (log(n) + log jjajj)) bit oerations with non-zero elts.
7 P-adic algorithm with matrix inversion Scheme to comute A b : (-) B := A mod (-2) r := b for i := to k (2-) x i := B:r mod (2-2) r := (=)(r A:x i ) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
8 P-adic algorithm with matrix inversion Scheme to comute A b : (-) B := A mod O~(n 3 log jjajj) (-2) r := b for i := to k k = O~(n) (2-) x i := B:r mod O~(n 2 log jjajj) (2-2) r := (=)(r A:x i ) O~(n 2 log jjajj) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
9 P-adic algorithm with matrix inversion Scheme to comute A b : (-) B := A mod O~(n 3 log jjajj) (-2) r := b for i := to k k = O~(n) (2-) x i := B:r mod O~(n 2 log jjajj) (2-2) r := (=)(r A:x i ) O~(n 2 log jjajj) (3-) x := P k i= x i : i (3-2) rational reconstruction on x Main oerations : matrix inversion and matrix-vector roducts
10 Dense linear system in ractice Ecient imlementations are available : LinBox. [ IML library [ Details : level 3 BLAS-based matrix inversion over rime eld with LQUP factorization [Dumas, Giorgi, Pernet 24] with Echelon form [Chen, Storjohann 25] level 2 BLAS-based matrix-vector roduct use of CRT over the integers rational number reconstruction half GCD [Schonage 97] heuristic using integer multilication [NTL library]
11 Timing for dense linear system solving use of LinBox library on Pentium 4-3.4Ghz, 2Go RAM random dense linear system with coecients over 3 bits : n time :6s 4:3s 3:s 99:6s 236:8s 449:2s random dense linear system with coecients over 2 bits : n time :8s 2:9s 9:5s 299:7s 76:4s MT erformances imrovement by a factor comare to NTL's tuned imlementation
12 what does haen when matrices are sarse? we consider sarse matrices with O(n) non zero elements,! matrix-vector roduct needs only O(n) oerations.
13 Scheme to comute A b : (-) B := A mod certainly dense (-2) r := b for i := to k (2-) x i := B:r mod dense roduct (2-2) r := (=)(r A:x i ) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
14 Sarse linear system and P-adic lifting P-adic lifting doesn't imrove comlexity as in dense case.,! comuting the modular inverse is roscribed due to ll-in Solution [Wiedemann 986 ; Kaltofen, Saunders 99] :,! use modular minimal olynomial instead of inverse
15 Sarse linear system and P-adic lifting P-adic lifting doesn't imrove comlexity as in dense case.,! comuting the modular inverse is roscribed due to ll-in Solution [Wiedemann 986 ; Kaltofen, Saunders 99] :,! use modular minimal olynomial instead of inverse Let A 2 Z n n of full rank and b 2 Z n. Then x = A b can be exressed as a linear combination of the Krylov subsace fb; Ab; :::; A n bg Let () = c + c + ::: + d 2 Z [] be the minimal olynomial of A
16 Sarse linear system and P-adic lifting P-adic lifting doesn't imrove comlexity as in dense case.,! comuting the modular inverse is roscribed due to ll-in Solution [Wiedemann 986 ; Kaltofen, Saunders 99] :,! use modular minimal olynomial instead of inverse Let A 2 Z n n of full rank and b 2 Z n. Then x = A b can be exressed as a linear combination of the Krylov subsace fb; Ab; :::; A n bg Let () = c + c + ::: + d 2 Z [] be the minimal olynomial of A A b = c (c b + c 2 Ab + ::: + A d b) {z }
17 Sarse linear system and P-adic lifting P-adic lifting doesn't imrove comlexity as in dense case.,! comuting the modular inverse is roscribed due to ll-in Solution [Wiedemann 986 ; Kaltofen, Saunders 99] :,! use modular minimal olynomial instead of inverse Let A 2 Z n n of full rank and b 2 Z n. Then x = A b can be exressed as a linear combination of the Krylov subsace fb; Ab; :::; A n bg Let () = c + c + ::: + d 2 Z [] be the minimal olynomial of A A b = c (c b + c 2 Ab + ::: + A d b) {z } x
18 P-adic algorithm for sarse systems Scheme to comute A b : (-) := minoly (A) mod (-2) r := b for i := to k X deg (2-) x i := [] [i] A i r mod i= (2-2) r := (=)(r A:x i ) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
19 P-adic algorithm for sarse systems Scheme to comute A b : (-) := minoly (A) mod O~(n 2 log jjajj) (-2) r := b for i := to k X k = O~(n) deg (2-) x i := [] [i] A i r mod O~(n 2 log jjajj) i= (2-2) r := (=)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
20 P-adic algorithm for sarse systems Scheme to comute A b : (-) := minoly (A) mod (-2) r := b for i := to k X k = O~(n) deg (2-) x i := [] [i] A i r mod O~(n 2 log jjajj) i= (2-2) r := (=)(r A:x i ) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
21 Integer sarse linear system in ractice use of LinBox library on Itanium II -.3Ghz, 28Go RAM random non-singular sarse linear system with coecients over 3 bits and non zero elements er row. system order Male 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
22 Integer sarse linear system in ractice use of LinBox library on Itanium II -.3Ghz, 28Go RAM random non-singular sarse linear system with coecients over 3 bits and non zero elements er row. system order Male 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 (Dixon) P (2-) x i := deg [] i= [i] A i r mod (P-adic-Wied) Remark : n sarse matrix alications is far from level 2 BLAS in ractice.
23 Our objectives In ratice : Integrate level 2 and 3 BLAS in integer sarse solver In theory : Imrove bit comlexity of sarse linear system solving =) O~(n ) bits oerations with < 3?
24 Integration of BLAS in sarse solver Our goals : minimize the number of sarse matrix-vector roducts. maximize the number of level 2 and 3 BLAS oerations.,! Block Wiedemann algorithm seems to be a good candidate Let s be the blocking factor of Block Wiedemann algorithm. then I the number of sarse matrix-vector roduct is divided by roughly s. I order s matrix oerations are integrated.
25 Block Wiedemann and P-adic Relace vector rojections by block of vectors rojections s f ( u A i s A@ v A f b is st column of v Let A 2 Z n n of full rank, b 2 Z n and n = m s. One can use a column of the minimal generating matrix olynomial P 2 Z [x] s s of sequence fua i vg to exress A of block krylov subsace fv ; Av ; : : : ; A m vg b as a linear combination
26 Block Wiedemann and P-adic Relace vector rojections by block of vectors rojections s f ( u A i s A@ v A f b is st column of v Let A 2 Z n n of full rank, b 2 Z n and n = m s. One can use a column of the minimal generating matrix olynomial P 2 Z [x] s s of sequence fua i vg to exress A of block krylov subsace fv ; Av ; : : : ; A m vg b as a linear combination the cost to comute P is : I O~(s 3 m) eld o. [Beckermann, Labahn 994 ; Kaltofen 995 ; Thome 22], I O~(s! m) eld o. [Giorgi, Jeannerod, Villard 23].
27 Block Wiedemann and P-adic Scheme to comute A b : (-) r := b for i := to k (2-) v ; := r (2-2) P := block minoly fua i vg mod (2-3) x i := linear combi (A i v ; P) mod (2-4) r := (=)(r A:x i ) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
28 Block Wiedemann and P-adic Scheme to comute A b : (-) r := b for i := to k (2-) v ; := r k = O~(n) (2-2) P := block minoly fua i vg mod O~(s 2 n log jjajj) (2-3) x i := linear combi (A i v ; P) mod O~(n 2 log jjajj) (2-4) r := (=)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
29 Block Wiedemann and P-adic Scheme to comute A b : (-) r := b for i := to k (2-) v ; := r k = O~(n) (2-2) P := block minoly fua i vg mod O~(s 2 n log jjajj) (2-3) x i := linear combi (A i v ; P) mod O~(n 2 log jjajj) (2-4) r := (=)(r A:x i ) (3-) x := P k i= x i : i (3-2) rational reconstruction on x Not satisfying : comutation of block minoly. at each stes How to avoid the comutation of the block minimal olynomial?
30 Alternative to Block Wiedemann Exress the inverse of the sarse matrix through a structured form,! block Hankel/Toelitz structures Let u 2 Z s n U = and v 2 Z n s u ua. ua m C s.t. following matrices are non-singular B C A, V v Av : : : Am v A 2 Zn n
31 Alternative to Block Wiedemann Exress the inverse of the sarse matrix through a structured form,! block Hankel/Toelitz structures Let u 2 Z s n U = and v 2 Z n s u ua. ua m C s.t. following matrices are non-singular B C A, V v Av : : : Am v A 2 Zn n then we can dene the block Hankel matrix H = UAV = 2 m 2 3 m+ C. A, m m 2m C i = ua i v 2 Z s s and thus we have A = VH U
32 Alternative to Block Wiedemann Nice roerty on block Hankel matrix inverse [Gohberg, Krunik 972, Labahn, Choi, Cabay 99] H = : : :.... C A : : : C A {z } {z } : : :.... C A : : : C A {z } {z } H T H 2 T 2 where H, H 2 are block Hankel matrices and T, T 2 are block Toelitz matrices
33 Alternative to Block Wiedemann Nice roerty on block Hankel matrix inverse [Gohberg, Krunik 972, Labahn, Choi, Cabay 99] H = : : :.... C A : : : C A {z } {z } : : :.... C A : : : C A {z } {z } H T H 2 T 2 where H, H 2 are block Hankel matrices and T, T 2 are block Toelitz matrices Block coecients in H, H 2, T, T 2 come from Hermite Pade aroximants of H(z) = + 2 z + : : : + 2m z 2m 2 [Labahn, Choi, Cabay 99]. Comlexity of H Jeannerod, Villard 23]. reduces to olynomial matrix multilication [Giorgi,
34 Alternative to Block Wiedemann Scheme to comute A b : (-) H(z) := 2m X i= ua i v :z i mod (-2) comute H mod from H(z) (-3) r := b for i := to k (2-) x i := V H U:r mod (2-2) r := (=)(r A:x i ) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
35 Alternative to Block Wiedemann Scheme to comute A b : (-) H(z) := 2m X i= ua i v :z i mod O~(sn 2 log jjajj) (-2) comute H mod 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 O~((n 2 + sn) log jjajj) (2-2) r := (=)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
36 Alternative to Block Wiedemann Scheme to comute A b : (-) H(z) := 2m X i= ua i v :z i mod O~(sn 2 log jjajj) (-2) comute H mod 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 O~((n 2 + sn) log jjajj) (2-2) r := (=)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i : i (3-2) rational reconstruction on x Not yet satisfying : alying matrices U and V is too costly
37 Alying block Krylov subsaces V v Av : : : A m v A 2 Z n n and v 2 Z n s can be rewrite as V v A + v A + : : : + A v A Therefore, alying V to a vector corresonds to : m linear combinations of columns of v m alications of A
38 Alying block Krylov subsaces V v Av : : : A m v A 2 Z n n and v 2 Z n s can be rewrite as V v A + v A + : : : + A v A Therefore, alying V to a vector corresonds to : m linear combinations of columns of v m alications of A O(m sn log jjajj) O(mn log jjajj)
39 Alying block Krylov subsaces V v Av : : : A m v A 2 Z n n and v 2 Z n s can be rewrite as V v A + v A + : : : + A v A Therefore, alying V to a vector corresonds to : m linear combinations of columns of v m alications of A O(m sn log jjajj) How to imrove the comlexity?
40 Alying block Krylov subsaces V v Av : : : A m v A 2 Z n n and v 2 Z n s can be rewrite as V v A + v A + : : : + A v A Therefore, alying V to a vector corresonds to : m linear combinations of columns of v m alications of A O(m sn log jjajj) How to imrove the comlexity? ) using secial block rojections u and v
41 Candidates as suitable block rojections Considering A 2 Z n n non-singular and n = m s. Let us denote K(A; v) := v Av A m v 2 Z n n A suitable block rojection is dened through the trile such that : (R; u; v) 2 Z n n Z s n Z n s. K(RA; v) and K((RA) T ; u T ) are non-singular, 2. R can be alied to a vector with O~(n) oerations, 3. u, u T, v and v T can be alied to a vector with O~(n) oerations.
42 Candidates as suitable block rojections Considering A 2 Z n n non-singular and n = m s. Let us denote K(A; v) := v Av A m v 2 Z n n A suitable block rojection is dened through the trile such that : (R; u; v) 2 Z n n Z s n Z n s. K(RA; v) and K((RA) T ; u T ) are non-singular, 2. R can be alied to a vector with O~(n) oerations, 3. u, u T, v and v T can be alied to a vector with O~(n) oerations. Conjecture : for any non-singular A 2 Z n n block rojection (R; u; v) and sjn there exists a suitable
43 A structured block rojection Let u and v be dened as follow u = v T = u : : : u m um+ : : : u 2m... v : : : v m vm+ : : : v 2m... u n m+ : : : u n v n m+ : : : v n C A 2 Z s n C A 2 Z s n where u i and v i are chosen randomly from a sucient large set.
44 A structured block rojection Let u and v be dened as follow u = v T = u : : : u m um+ : : : u 2m... v : : : v m vm+ : : : v 2m... u n m+ : : : u n v n m+ : : : v n C A 2 Z s n C A 2 Z s n where u i and v i are chosen randomly from a sucient large set. oen question : Let R diagonal and v as dened above, is K(RA; v) necessarily non-singular? We rooved it for case s = 2 but answer is still unknown for s > 2
45 Our new algorithm Scheme to comute A b : (-) choose block rojection 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 (-4) comute H mod from H(z) (-5) r := b for i := to k (2-) x i := V H U:r mod (2-2) r := (=)(r A:x i ) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
46 Our new algorithm Scheme to comute A b : (-) choose block rojection 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 O~(n 2 log jjajj) (-4) comute H mod 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 O~((mn + sn) log jjajj) (2-2) r := (=)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i : i (3-2) rational reconstruction on x
47 Our new algorithm Scheme to comute A b : (-) choose block rojection 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 O~(n 2 log jjajj) (-4) comute H mod 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 O~((mn + sn) log jjajj) (2-2) r := (=)(r A:x i ) O~(n log jjajj) (3-) x := P k i= x i : i (3-2) rational reconstruction on x taking the otimal m = s = n gives a comlexity of O~(n 2:5 log jjajj)
48 High level imlementation LinBox roject (Canada-France-USA) : Our tools : BLAS-based matrix multilication and matrix-vector roduct olynomial matrix arithmetic (block Hankel inversion),! FFT, Karatsuba, middle roduct fast alication of H is needed to get O~(n 2:5 log jjajj)
49 High level imlementation LinBox roject (Canada-France-USA) : Our tools : BLAS-based matrix multilication and matrix-vector roduct olynomial matrix arithmetic (block Hankel inversion),! FFT, Karatsuba, middle roduct fast alication of H is needed to get O~(n 2:5 log jjajj) I Lagrange's reresentation of H at the beginning (Horner's scheme) I use evaluation/interolation on olynomial vectors,! use Vandermonde matrix to have dense matrix oerations
50 Is our new algorithm ecient in ractice?
51 Performances use of LinBox library on Itanium II -.3Ghz, 28Go RAM random non-singular sarse linear system with coecients over 3 bits and non zero elements er row. system order Male 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 Our algo. 2:4s 5s 6s 75s 426s The exected n imrovement is unfortunately amortized by a high constant in the comlexity.
52 Sarse solver vs Dixon's algorithm 6 sarsity = elts/row Our algo. Dixon 6 sarsity = 3elts/row Our algo. Dixon Time ( 4 s) 3 Time ( 4 s) System order System order 44 Our algorithm erformances are deending on matrix sarsity
53 Sarse solver vs Dixon's algorithm 4 sarsity=.7% sarsity=.3% sarsity=.% crossover line 2 seed u System order The sarser the matrices are, the earlier the crossover aears
54 Practical eect on blocking factors n blocking factor value is theoretically otimal Is this still true in ractice?
55 Practical eect on blocking factors n blocking factor value is theoretically otimal Is this still true in ractice? system order =, otimal block = block size timing 723s 5264s 459s 3833s 4332s system order = 2, otimal block 4 block size timing 4472s 35967s 3854s 2852s 3738s
56 Practical eect on blocking factors n blocking factor value is theoretically otimal Is this still true in ractice? system order =, otimal block = block size timing 723s 5264s 459s 3833s 4332s system order = 2, otimal block 4 block size timing 4472s 35967s 3854s 2852s 3738s best ractical blocking factor is certainly deending on the ratio of sarse matrix/dense matrix oerations eciency
57 Conclusions We rovide a new aroach for solving sarse integer linear systems : I imrove the comlexty by a factor n (heuristic). I allow eciency by minimizing sarse matrix oerations and maximizing BLAS use. We introduce secial block rojections for sarse linear algebra,! inverse of sarse matrix in O(n 2:5 ) eld o. drawback : not taking advantage of low degree minimal olynomial
58 Conclusions We rovide a new aroach for solving sarse integer linear systems : I imrove the comlexty by a factor n (heuristic). I allow eciency by minimizing sarse matrix oerations and maximizing BLAS use. We introduce secial block rojections for sarse linear algebra,! inverse of sarse matrix in O(n 2:5 ) eld o. drawback : not taking advantage of low degree minimal olynomial On going work : I rovide an automatic choice of block dimension (non square?) I roove conjecture for secial block rojections I how to handle the case of singular matrix? I how to introduce fast matrix multilication in the comlexity?
Solving 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 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 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 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 informationRelaxed p-adic Hensel lifting for algebraic systems
Relaxed -adic Hensel lifting for algebraic systems Jérémy Berthomieu Laboratoire de Mathématiques Université de Versailles Versailles, France jeremy.berthomieu@uvsq.fr Romain Lebreton Laboratoire d Informatique
More informationGOOD MODELS FOR CUBIC SURFACES. 1. Introduction
GOOD MODELS FOR CUBIC SURFACES ANDREAS-STEPHAN ELSENHANS Abstract. This article describes an algorithm for finding a model of a hyersurface with small coefficients. It is shown that the aroach works in
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 informationOn Nonlinear Polynomial Selection and Geometric Progression (mod N) for Number Field Sieve
On onlinear Polynomial Selection and Geometric Progression (mod ) for umber Field Sieve amhun Koo, Gooc Hwa Jo, and Soonhak Kwon Email: komaton@skku.edu, achimheasal@nate.com, shkwon@skku.edu Det. of Mathematics,
More informationMATH 361: NUMBER THEORY ELEVENTH LECTURE
MATH 361: NUMBER THEORY ELEVENTH LECTURE The subjects of this lecture are characters, Gauss sums, Jacobi sums, and counting formulas for olynomial equations over finite fields. 1. Definitions, Basic Proerties
More informationRound-off Errors and Computer Arithmetic - (1.2)
Round-off Errors and Comuter Arithmetic - (.). Round-off Errors: Round-off errors is roduced when a calculator or comuter is used to erform real number calculations. That is because the arithmetic erformed
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 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 information2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized
BIT 38 :2 (998), 283{292. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES II KARL ENTACHER y Deartment of Mathematics, University of Salzburg, Hellbrunnerstr. 34 A-52 Salzburg,
More informationThe analysis and representation of random signals
The analysis and reresentation of random signals Bruno TOÉSNI Bruno.Torresani@cmi.univ-mrs.fr B. Torrésani LTP Université de Provence.1/30 Outline 1. andom signals Introduction The Karhunen-Loève Basis
More informationElliptic Curves and Cryptography
Ellitic Curves and Crytograhy Background in Ellitic Curves We'll now turn to the fascinating theory of ellitic curves. For simlicity, we'll restrict our discussion to ellitic curves over Z, where is a
More informationFor q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i
Comuting with Haar Functions Sami Khuri Deartment of Mathematics and Comuter Science San Jose State University One Washington Square San Jose, CA 9519-0103, USA khuri@juiter.sjsu.edu Fax: (40)94-500 Keywords:
More informationApplicable Analysis and Discrete Mathematics available online at HENSEL CODES OF SQUARE ROOTS OF P-ADIC NUMBERS
Alicable Analysis and Discrete Mathematics available online at htt://efmath.etf.rs Al. Anal. Discrete Math. 4 (010), 3 44. doi:10.98/aadm1000009m HENSEL CODES OF SQUARE ROOTS OF P-ADIC NUMBERS Zerzaihi
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
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 informationCubic Sieve Congruence of the Discrete Logarithm Problem, and Fractional Part Sequences
Cubic Sieve Congruence of the Discrete Logarithm Problem, and Fractional Part Sequences Srinivas Vivek University of Luxembourg, Luxembourg C. E. Veni Madhavan Deartment of Comuter Science and Automation,
More information#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS
#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS Norbert Hegyvári ELTE TTK, Eötvös University, Institute of Mathematics, Budaest, Hungary hegyvari@elte.hu François Hennecart Université
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 informationON POLYNOMIAL SELECTION FOR THE GENERAL NUMBER FIELD SIEVE
MATHEMATICS OF COMPUTATIO Volume 75, umber 256, October 26, Pages 237 247 S 25-5718(6)187-9 Article electronically ublished on June 28, 26 O POLYOMIAL SELECTIO FOR THE GEERAL UMBER FIELD SIEVE THORSTE
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More informationarxiv: v1 [math.nt] 4 Nov 2015
Wall s Conjecture and the ABC Conjecture George Grell, Wayne Peng August 0, 018 arxiv:1511.0110v1 [math.nt] 4 Nov 015 Abstract We show that the abc conjecture of Masser-Oesterlé-Sziro for number fields
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 informationOutline. EECS150 - Digital Design Lecture 26 Error Correction Codes, Linear Feedback Shift Registers (LFSRs) Simple Error Detection Coding
Outline EECS150 - Digital Design Lecture 26 Error Correction Codes, Linear Feedback Shift Registers (LFSRs) Error detection using arity Hamming code for error detection/correction Linear Feedback Shift
More informationGalois representations on torsion points of elliptic curves NATO ASI 2014 Arithmetic of Hyperelliptic Curves and Cryptography
Galois reresentations on torsion oints of ellitic curves NATO ASI 04 Arithmetic of Hyerellitic Curves and Crytograhy Francesco Paalardi Ohrid, August 5 - Setember 5, 04 Lecture - Introduction Let /Q be
More informationA Recursive Block Incomplete Factorization. Preconditioner for Adaptive Filtering Problem
Alied Mathematical Sciences, Vol. 7, 03, no. 63, 3-3 HIKARI Ltd, www.m-hiari.com A Recursive Bloc Incomlete Factorization Preconditioner for Adative Filtering Problem Shazia Javed School of Mathematical
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationHENSEL S LEMMA KEITH CONRAD
HENSEL S LEMMA KEITH CONRAD 1. Introduction In the -adic integers, congruences are aroximations: for a and b in Z, a b mod n is the same as a b 1/ n. Turning information modulo one ower of into similar
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
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 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 informationNumbers and functions. Introduction to Vojta s analogy
Numbers and functions. Introduction to Vojta s analogy Seminar talk by A. Eremenko, November 23, 1999, Purdue University. Absolute values. Let k be a field. An absolute value v is a function k R, x x v
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 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 ALGORITHM FOR MATRIX EXTENSION AND WAVELET CONSTRUCTION W. LAWTON, S. L. LEE AND ZUOWEI SHEN Abstract. This paper gives a practical method of exten
AN ALGORITHM FOR MATRIX EXTENSION AND WAVELET ONSTRUTION W. LAWTON, S. L. LEE AND ZUOWEI SHEN Abstract. This aer gives a ractical method of extending an nr matrix P (z), r n, with Laurent olynomial entries
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 informationSection 0.10: Complex Numbers from Precalculus Prerequisites a.k.a. Chapter 0 by Carl Stitz, PhD, and Jeff Zeager, PhD, is available under a Creative
Section 0.0: Comlex Numbers from Precalculus Prerequisites a.k.a. Chater 0 by Carl Stitz, PhD, and Jeff Zeager, PhD, is available under a Creative Commons Attribution-NonCommercial-ShareAlike.0 license.
More informationCryptanalysis of Pseudorandom Generators
CSE 206A: Lattice Algorithms and Alications Fall 2017 Crytanalysis of Pseudorandom Generators Instructor: Daniele Micciancio UCSD CSE As a motivating alication for the study of lattice in crytograhy we
More informationA SIMPLE AD EFFICIET PARALLEL FFT ALGORITHM USIG THE BSP MODEL MARCIA A. IDA AD ROB H. BISSELIG Abstract. In this aer, we resent a new arallel radix-4
Universiteit-Utrecht * Deartment of Mathematics A simle and ecient arallel FFT algorithm using the BSP model by Marcia A. Inda and Rob H. Bisseling Prerint nr. 3 March 2000 A SIMPLE AD EFFICIET PARALLEL
More informationMA3H1 TOPICS IN NUMBER THEORY PART III
MA3H1 TOPICS IN NUMBER THEORY PART III SAMIR SIKSEK 1. Congruences Modulo m In quadratic recirocity we studied congruences of the form x 2 a (mod ). We now turn our attention to situations where is relaced
More informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationElliptic Curves Spring 2015 Problem Set #1 Due: 02/13/2015
18.783 Ellitic Curves Sring 2015 Problem Set #1 Due: 02/13/2015 Descrition These roblems are related to the material covered in Lectures 1-2. Some of them require the use of Sage, and you will need to
More informationThe Arm Prime Factors Decomposition
The Arm Prime Factors Decomosition Arm Boris Nima arm.boris@gmail.com Abstract We introduce the Arm rime factors decomosition which is the equivalent of the Taylor formula for decomosition of integers
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 informationIDENTIFYING CONGRUENCE SUBGROUPS OF THE MODULAR GROUP
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 24, Number 5, May 996 IDENTIFYING CONGRUENCE SUBGROUPS OF THE MODULAR GROUP TIM HSU (Communicated by Ronald M. Solomon) Abstract. We exhibit a simle
More information0.6 Factoring 73. As always, the reader is encouraged to multiply out (3
0.6 Factoring 7 5. The G.C.F. of the terms in 81 16t is just 1 so there is nothing of substance to factor out from both terms. With just a difference of two terms, we are limited to fitting this olynomial
More informationQuaternionic Projective Space (Lecture 34)
Quaternionic Projective Sace (Lecture 34) July 11, 2008 The three-shere S 3 can be identified with SU(2), and therefore has the structure of a toological grou. In this lecture, we will address the question
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 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 informationDynamic Countermeasure Against the Zero Power Analysis
Dynamic Countermeasure Against the Zero Power Analysis Jean-Luc Danger 1,2, Sylvain Guilley 1,2, Philie Hoogvorst 2, Cédric Murdica 1,2, and David Naccache 3 1 Secure-IC S.A.S., 80 avenue des Buttes de
More informationSolving Cyclotomic Polynomials by Radical Expressions Andreas Weber and Michael Keckeisen
Solving Cyclotomic Polynomials by Radical Exressions Andreas Weber and Michael Keckeisen Abstract: We describe a Male ackage that allows the solution of cyclotomic olynomials by radical exressions. We
More informationResearch Article A Note on the Modified q-bernoulli Numbers and Polynomials with Weight α
Abstract and Alied Analysis Volume 20, Article ID 54534, 8 ages doi:0.55/20/54534 Research Article A Note on the Modified -Bernoulli Numbers and Polynomials with Weight α T. Kim, D. V. Dolgy, 2 S. H. Lee,
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 informationPOINTS ON CONICS MODULO p
POINTS ON CONICS MODULO TEAM 2: JONGMIN BAEK, ANAND DEOPURKAR, AND KATHERINE REDFIELD Abstract. We comute the number of integer oints on conics modulo, where is an odd rime. We extend our results to conics
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 information18.312: Algebraic Combinatorics Lionel Levine. Lecture 12
8.3: Algebraic Combinatorics Lionel Levine Lecture date: March 7, Lecture Notes by: Lou Odette This lecture: A continuation of the last lecture: comutation of µ Πn, the Möbius function over the incidence
More informationIntroduction to Group Theory Note 1
Introduction to Grou Theory Note July 7, 009 Contents INTRODUCTION. Examles OF Symmetry Grous in Physics................................. ELEMENT OF GROUP THEORY. De nition of Grou................................................
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 informationc Copyright by Helen J. Elwood December, 2011
c Coyright by Helen J. Elwood December, 2011 CONSTRUCTING COMPLEX EQUIANGULAR PARSEVAL FRAMES A Dissertation Presented to the Faculty of the Deartment of Mathematics University of Houston In Partial Fulfillment
More informationCatalan s Equation Has No New Solution with Either Exponent Less Than 10651
Catalan s Euation Has No New Solution with Either Exonent Less Than 065 Maurice Mignotte and Yves Roy CONTENTS. Introduction and Overview. Bounding One Exonent as a Function of the Other 3. An Alication
More informationAn Overview of Witt Vectors
An Overview of Witt Vectors Daniel Finkel December 7, 2007 Abstract This aer offers a brief overview of the basics of Witt vectors. As an alication, we summarize work of Bartolo and Falcone to rove that
More informationMAT 311 Solutions to Final Exam Practice
MAT 311 Solutions to Final Exam Practice Remark. If you are comfortable with all of the following roblems, you will be very well reared for the midterm. Some of the roblems below are more difficult than
More informationGalois Fields, Linear Feedback Shift Registers and their Applications
Galois Fields, Linear Feedback Shift Registers and their Alications With 85 illustrations as well as numerous tables, diagrams and examles by Ulrich Jetzek ISBN (Book): 978-3-446-45140-7 ISBN (E-Book):
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
More informationNUMBER SYSTEMS. Number theory is the study of the integers. We denote the set of integers by Z:
NUMBER SYSTEMS Number theory is the study of the integers. We denote the set of integers by Z: Z = {..., 3, 2, 1, 0, 1, 2, 3,... }. The integers have two oerations defined on them, addition and multilication,
More informationChater Matrix Norms and Singular Value Decomosition Introduction In this lecture, we introduce the notion of a norm for matrices The singular value de
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Deartment of Electrical Engineering and Comuter Science Massachuasetts Institute of Technology c Chater Matrix Norms
More informationIntroduction to Arithmetic Geometry Fall 2013 Lecture #10 10/8/2013
18.782 Introduction to Arithmetic Geometry Fall 2013 Lecture #10 10/8/2013 In this lecture we lay the groundwork needed to rove the Hasse-Minkowski theorem for Q, which states that a quadratic form over
More informationON OPTIMIZATION OF THE MEASUREMENT MATRIX FOR COMPRESSIVE SENSING
8th Euroean Signal Processing Conference (EUSIPCO-2) Aalborg, Denmark, August 23-27, 2 ON OPTIMIZATION OF THE MEASUREMENT MATRIX FOR COMPRESSIVE SENSING Vahid Abolghasemi, Saideh Ferdowsi, Bahador Makkiabadi,2,
More informationAdditive results for the generalized Drazin inverse in a Banach algebra
Additive results for the generalized Drazin inverse in a Banach algebra Dragana S. Cvetković-Ilić Dragan S. Djordjević and Yimin Wei* Abstract In this aer we investigate additive roerties of the generalized
More informationDedekind sums and continued fractions
ACTA ARITHMETICA LXIII.1 (1993 edekind sums and continued fractions by R. R. Hall (York and M. N. Huxley (Cardiff Let ϱ(t denote the row-of-teeth function ϱ(t = [t] t + 1/2. Let a b c... r be ositive integers.
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationAN IMPROVED BABY-STEP-GIANT-STEP METHOD FOR CERTAIN ELLIPTIC CURVES. 1. Introduction
J. Al. Math. & Comuting Vol. 20(2006), No. 1-2,. 485-489 AN IMPROVED BABY-STEP-GIANT-STEP METHOD FOR CERTAIN ELLIPTIC CURVES BYEONG-KWEON OH, KIL-CHAN HA AND JANGHEON OH Abstract. In this aer, we slightly
More information(Workshop on Harmonic Analysis on symmetric spaces I.S.I. Bangalore : 9th July 2004) B.Sury
Is e π 163 odd or even? (Worksho on Harmonic Analysis on symmetric saces I.S.I. Bangalore : 9th July 004) B.Sury e π 163 = 653741640768743.999999999999.... The object of this talk is to exlain this amazing
More informationDISCRIMINANTS IN TOWERS
DISCRIMINANTS IN TOWERS JOSEPH RABINOFF Let A be a Dedekind domain with fraction field F, let K/F be a finite searable extension field, and let B be the integral closure of A in K. In this note, we will
More informationComputer arithmetic. Intensive Computation. Annalisa Massini 2017/2018
Comuter arithmetic Intensive Comutation Annalisa Massini 7/8 Intensive Comutation - 7/8 References Comuter Architecture - A Quantitative Aroach Hennessy Patterson Aendix J Intensive Comutation - 7/8 3
More information#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS
#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS Ramy F. Taki ElDin Physics and Engineering Mathematics Deartment, Faculty of Engineering, Ain Shams University, Cairo, Egyt
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 informationRAMANUJAN-NAGELL CUBICS
RAMANUJAN-NAGELL CUBICS MARK BAUER AND MICHAEL A. BENNETT ABSTRACT. A well-nown result of Beuers [3] on the generalized Ramanujan-Nagell equation has, at its heart, a lower bound on the quantity x 2 2
More informationEvaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models
Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu
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 information#A45 INTEGERS 12 (2012) SUPERCONGRUENCES FOR A TRUNCATED HYPERGEOMETRIC SERIES
#A45 INTEGERS 2 (202) SUPERCONGRUENCES FOR A TRUNCATED HYPERGEOMETRIC SERIES Roberto Tauraso Diartimento di Matematica, Università di Roma Tor Vergata, Italy tauraso@mat.uniroma2.it Received: /7/, Acceted:
More informationCSE 599d - Quantum Computing When Quantum Computers Fall Apart
CSE 599d - Quantum Comuting When Quantum Comuters Fall Aart Dave Bacon Deartment of Comuter Science & Engineering, University of Washington In this lecture we are going to begin discussing what haens to
More informationRECIPROCITY LAWS JEREMY BOOHER
RECIPROCITY LAWS JEREMY BOOHER 1 Introduction The law of uadratic recirocity gives a beautiful descrition of which rimes are suares modulo Secial cases of this law going back to Fermat, and Euler and Legendre
More informationMA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES 2018
MA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES 2018 J. E. CREMONA Contents 0. Introduction: What is Number Theory? 2 Basic Notation 3 1. Factorization 4 1.1. Divisibility in Z 4 1.2. Greatest Common
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationTensor-Based Sparsity Order Estimation for Big Data Applications
Tensor-Based Sarsity Order Estimation for Big Data Alications Kefei Liu, Florian Roemer, João Paulo C. L. da Costa, Jie Xiong, Yi-Sheng Yan, Wen-Qin Wang and Giovanni Del Galdo Indiana University School
More informationPythagorean triples and sums of squares
Pythagorean triles and sums of squares Robin Chaman 16 January 2004 1 Pythagorean triles A Pythagorean trile (x, y, z) is a trile of ositive integers satisfying z 2 + y 2 = z 2. If g = gcd(x, y, z) then
More informationAKRON: An Algorithm for Approximating Sparse Kernel Reconstruction
: An Algorithm for Aroximating Sarse Kernel Reconstruction Gregory Ditzler Det. of Electrical and Comuter Engineering The University of Arizona Tucson, AZ 8572 USA ditzler@email.arizona.edu Nidhal Carla
More informationA STUDY ON THE UTILIZATION OF COMPATIBILITY METRIC IN THE AHP: APPLYING TO SOFTWARE PROCESS ASSESSMENTS
ISAHP 2005, Honolulu, Hawaii, July 8-10, 2003 A SUDY ON HE UILIZAION OF COMPAIBILIY MERIC IN HE AHP: APPLYING O SOFWARE PROCESS ASSESSMENS Min-Suk Yoon Yosu National University San 96-1 Dundeok-dong Yeosu
More informationMATH342 Practice Exam
MATH342 Practice Exam This exam is intended to be in a similar style to the examination in May/June 2012. It is not imlied that all questions on the real examination will follow the content of the ractice
More informationConstruction algorithms for good extensible lattice rules
Construction algorithms for good extensible lattice rules Harald Niederreiter and Friedrich Pillichshammer Dedicated to Ian H. Sloan on the occasion of his 70th birthday Abstract Extensible olynomial lattice
More information