Advances in Random Matrix Theory: Let there be tools
|
|
- Raymond Price
- 5 years ago
- Views:
Transcription
1 Advances in Random Matrix Theory: Let there be tools Alan Edelman Brian Sutton, Plamen Koev, Ioana Dumitriu, Raj Rao and others MIT: Dept of Mathematics, Computer Science AI Laboratories World Congress, Bernoulli Society Barcelona, Spain Wednesday July 28, 2004
2 Message Ingredient: Take Any important mathematics Then Randomize! This will have many applications! We can t keep this in the hands of specialists anymore: Need tools! 2
3 Tools So many applications Random matrix theory: catalyst for 2 st century special functions, analytical techniques, statistical techniques In addition to mathematics and papers Need tools for the novice! Need tools for the engineers! Need tools for the specialists! 3
4 Themes of this talk Tools for general beta What is beta? Think of it as a measure of (inverse) volatility in classical random matrices. Tools for complicated derived random matrices Tools for numerical computation and simulation 4
5 Wigner s Semi-Circle The classical & most famous rand eig theorem Let S = random symmetric Gaussian MATLAB: A=randn(n); S=(A+A )/2; S known as the Gaussian Orthogonal Ensemble Normalized eigenvalue histogram is a semi-circle n=20; s=30000; d=.05; %matrix size, samples, sample dist e=[]; %gather up eigenvalues im=; %imaginary() or real(0) for i=:s, a=randn(n)+im*sqrt(-)*randn(n);a=(a+a')/(2*sqrt(2*n*(im+))); v=eig(a)'; e=[e v]; end hold off; [m x]=hist(e,-.5:d:.5); bar(x,m*pi/(2*d*n*s)); axis('square'); axis([ ]); hold on; t=-:.0:; plot(t,sqrt(-t.^2),'r'); Precise statements require n etc. 5
6 sym matrix to tridiagonal form χ 4 χ 3 Same eigenvalue distribution χ G as 3 χgoe: O(n) storage!! O(n 2 ) compute 2 G χ 6 χ 6 G χ 5 χ 5 G χ 4 G χ 2 G χ χ G 6
7 General beta G χ 6β χ 6β G beta: : reals 2: complexes 4: quaternions χ 5β χ 5β G χ Bidiagonal 4β Version corresponds χ To Wishart G 4β χmatrices of Statistics 3β χ G 3β χ 2β χ 2β G χ β χ β G 7
8 Tools Motivation: A condition number problem Jack & Hypergeometric of Matrix Argument MOPS: Ioana Dumitriu s talk The Polynomial Method The tridiagonal numerical 09 trick 8
9 Tools Motivation: A condition number problem Jack & Hypergeometric of Matrix Argument MOPS: Ioana Dumitriu s talk The Polynomial Method The tridiagonal numerical 09 trick 9
10 Numerical Analysis: Condition Numbers κ(a) = condition number of A If A=UΣV is the svd, then κ(a) = σ max /σ min. Alternatively, κ(a) = λ max (A A)/ λ min (A A) One number that measures digits lost in finite precision and general matrix badness Small=good Large=bad The condition of a random matrix??? 0
11 Von Neumann & co. Solve Ax=b via x= (A A) - A b M A - Matrix Residual: AM-I 2 AM-I 2 < 200κ 2 n ε How should we estimate κ? Assume, as a model, that the elements of A are independent standard normals!
12 Von Neumann & co. estimates (947-95) For a random matrix of order n the expectation value has been shown to be about n X Goldstine, von Neumann we choose two different values of κ, namely n and 0n P(κ<n) 0.02 P(κ< 0n) 0.44 Bargmann, Montgomery, vn With a probability ~ κ < 0n P(κ<0n) 0.80 Goldstine, von Neumann 2
13 Random cond numbers, n y Distribution of κ/n = x e 2 / x 2 / x 3 x 2 Experiment with n=200 3
14 Finite n n=0 n=25 n=50 n=00 4
15 Condition Number Distributions P(κ/n > x) 2/x P(κ/n 2 > x) 4/x 2 Real n x n, næ Generalizations: β: =real, 2=complex finite matrices rectangular: m x n Complex n x n, næ 5
16 Condition Number Distributions P(κ/n > x) 2/x P(κ/n 2 > x) 4/x 2 Square, næ : P(κ/n β > x) (2β β- /Γ(β))/x β (All Betas!!) General Real Formula: n x n, næ P(κ>x) ~Cµ/x β(n-m+), where µ= β(n-m+)/2th moment of the largest eigenvalue of W m-,n+ (β) and C is a known geometrical constant. Complex n x n, næ Density for the largest eig of W is known in terms of F ((β/2)(n+), ((β/2)(n+m-); -(x/2)i m- ) from which µ is available Tracy-Widom law applies probably all beta for large m,n. Johnstone shows at least beta=,2. 6
17 Tools Motivation: A condition number problem Jack & Hypergeometric of Matrix Argument MOPS: Ioana Dumitriu s talk The Polynomial Method The tridiagonal numerical 09 trick 7
18 Multivariate Orthogonal Polynomials & Hypergeometrics of Matrix Argument Ioana Dumitriu s talk The important special functions of the 2 st century Begin with w(x) on I p κ (x)p λ (x) (x) β i w(x i )dx i = δ κλ Jack Polynomials orthogonal for w= on the unit circle. Analogs of x m 8
19 Multivariate Hypergeometric Functions 9
20 Multivariate Hypergeometric Functions 20
21 Wishart (Laguerre) Matrices! 2
22 Plamen s clever idea 22
23 Tools Motivation: A condition number problem Jack & Hypergeometric of Matrix Argument MOPS: Ioana Dumitriu s talk The Polynomial Method The tridiagonal numerical 09 trick 23
24 Mops (Dumitriu etc.) Symbolic 24
25 Symbolic MOPS applications A=randn(n); S=(A+A )/2; trace(s^4) det(s^3) 25
26 Symbolic MOPS applications β=3; hist(eig(s)) 26
27 Smallest eigenvalue statistics A=randn(m,n); hist(min(svd(a).^2)) 27
28 Tools Motivation: A condition number problem Jack & Hypergeometric of Matrix Argument MOPS: Ioana Dumitriu s talk The Polynomial Method -- Raj! The tridiagonal numerical 09 trick 28
29 RM Tool Raj! Courtesy of the Polynomial Method 29
30 30
31 = = Γ = 0 ) ( ) ζ( k s u s k du e u s s 6 π ) ζ 2 ( = = 3 The Riemann Zeta Function On the real line with x>, for example May be analytically extended to the complex plane, with singularity only at x=.
32 ½ 2 3 = = Γ = 0 ) ( ) ζ( k x u x k du e u x x All nontrivial roots of ζ(x) satisfy Re(x)=/2. (Trivial roots at negative even integers.) The Riemann Hypothesis
33 ζ( x) ζ(x) along Re(x)=/2 Zeros =.5+i* The Riemann Hypothesis u x = du = Γ( x) e u k x k= ½ All nontrivial roots of ζ(x) satisfy Re(x)=/2. (Trivial roots at negative even integers.) 33
34 Computation of Zeros Odlyzko s fantastic computation of 0^k+ through 0^k+0,000 for k=2,2,22. See Spacings behave like the eigenvalues of A=randn(n)+i*randn(n); S=(A+A )/2; 34
35 Nearest Neighbor Spacings & Pairwise Correlation Functions 35
36 Painlevé Equations 36
37 Spacings Take a large collection of consecutive zeros/eigenvalues. Normalize so that average spacing =. Spacing Function = Histogram of consecutive differences (the (k+)st the kth) Pairwise Correlation Function = Histogram of all possible differences (the kth the jth) Conjecture: These functions are the same for random matrices and Riemann zeta 37
38 Tools Motivation: A condition number problem Jack & Hypergeometric of Matrix Argument MOPS: Ioana Dumitriu s talk The Polynomial Method The tridiagonal numerical 09 trick 38
39 Everyone s Favorite Tridiagonal n d 2 dx 2 39
40 Everyone s Favorite Tridiagonal n (βn) /2 G G G d 2 dx 2 + dw β /2 40
41 4 Stochastic Operator Limit, N(0,2) χ χ N(0,2) χ χ N(0,2) χ χ N(0,2) nβ 2 ~ H β β 2β 2)β (n )β (n )β (n β n β 2 x dx d 2 2 +, G β 2 H H n n β n + dw,
42 Largest Eigenvalue Plots 42
43 MATLAB beta=; n=e9; opts.disp=0;opts.issym=; alpha=0; k=round(alpha*n^(/3)); % cutoff parameters d=sqrt(chi2rnd( beta*(n:-:(n-k-))))'; H=spdiags( d,,k,k)+spdiags(randn(k,),0,k,k); H=(H+H')/sqrt(4*n*beta); eigs(h,,,opts) 43
44 Tricks to get O(n 9 ) speedup Sparse matrix storage (Only O(n) storage is used) Tridiagonal Ensemble Formulas (Any beta is available due to the tridiagonal ensemble) The Lanczos Algorithm for Eigenvalue Computation ( This allows the computation of the extreme eigenvalue faster than typical general purpose eigensolvers.) The shift-and-invert accelerator to Lanczos and Arnoldi (Since we know the eigenvalues are near, we can accelerate the convergence of the largest eigenvalue) The ARPACK software package as made available seamlessly in MATLAB (The Arnoldi package contains state of the art data structures and numerical choices.) The observation that if k = 0n /3, then the largest eigenvalue is determined numerically by the top k k segment of n. (This is an interesting mathematical statement related to the decay of the Airy function.) 44
45 Level Densities 45
46 Open Problems The distribution for general beta Seems to be governed by a convection-diffusion equation 46
47 Random matrix tools! 47
Alan Edelman. Why are random matrix eigenvalues cool? MAA Mathfest 2002 Thursday, August 1. MIT: Dept of Mathematics, Lab for Computer Science
Why are random matrix eigenvalues cool? Alan Edelman MIT: Dept of Mathematics, Lab for Computer Science MAA Mathfest 2002 Thursday, August 6/2/2004 Message Ingredient: Take Any important mathematics Then
More informationProgress in the method of Ghosts and Shadows for Beta Ensembles
Progress in the method of Ghosts and Shadows for Beta Ensembles Alan Edelman (MIT) Alex Dubbs (MIT) and Plamen Koev (SJS) Oct 8, 2012 1/47 Wishart Matrices (arbitrary covariance) G=mxn matrix of Gaussians
More informationNumerical Methods for Random Matrices
Numerical Methods for Random Matrices MIT 18.95 IAP Lecture Series Per-Olof Persson (persson@mit.edu) January 23, 26 1.9.8.7 Random Matrix Eigenvalue Distribution.7.6.5 β=1 β=2 β=4 Probability.6.5.4.4
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. MATRIX ANAL. APPL. Vol. XX, No. X, pp. XX XX c 005 Society for Industrial and Applied Mathematics DISTRIBUTIONS OF THE EXTREME EIGENVALUES OF THE COMPLEX JACOBI RANDOM MATRIX ENSEMBLE PLAMEN KOEV
More informationThe Random Matrix Technique of Ghosts and Shadows
The Random Matrix Technique of Ghosts and Shadows Alan Edelman November 22, 2009 Abstract We propose to abandon the notion that a random matrix has to be sampled for it to exist. Much of today's applied
More informationExponential tail inequalities for eigenvalues of random matrices
Exponential tail inequalities for eigenvalues of random matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify
More informationUniversality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium
Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium SEA 06@MIT, Workshop on Stochastic Eigen-Analysis and its Applications, MIT, Cambridge,
More informationRandom matrices: Distribution of the least singular value (via Property Testing)
Random matrices: Distribution of the least singular value (via Property Testing) Van H. Vu Department of Mathematics Rutgers vanvu@math.rutgers.edu (joint work with T. Tao, UCLA) 1 Let ξ be a real or complex-valued
More informationRandom Matrices Numerical Methods for Random Matrices. Per-Olof Persson
18.325 Random Matrices Numerical Methods for Random Matrices Per-Olof Persson (persson@mit.edu) December 19, 22 1 Largest Eigenvalue Distributions In this section, the distributions of the largest eigenvalue
More informationConcentration Inequalities for Random Matrices
Concentration Inequalities for Random Matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify the asymptotic
More informationUpdate on the beta ensembles
Update on the beta ensembles Brian Rider Temple University with M. Krishnapur IISC, J. Ramírez Universidad Costa Rica, B. Virág University of Toronto The Tracy-Widom laws Consider a random Hermitian n
More informationarxiv:math-ph/ v1 27 Jan 2005
arxiv:math-ph/5168v1 27 Jan 25 Numerical Methods for Eigenvalue Distributions of Random Matrices Alan Edelman and Per-Olof Persson September 7, 218 Abstract We present efficient numerical techniques for
More informationEigenvalue Statistics for Beta-Ensembles. Ioana Dumitriu
Eigenvalue Statistics for Beta-Ensembles by Ioana Dumitriu B.A., New York University, 1999 Submitted to the Department of Mathematics in partial fulfillment of the requirements for the degree of DOCTOR
More informationMIT Term Project: Numerical Experiments on Circular Ensembles and Jack Polynomials with Julia
MIT 18.338 Term Project: Numerical Experiments on Circular Ensembles and Jack Polynomials with Julia N. Kemal Ure May 15, 2013 Abstract This project studies the numerical experiments related to averages
More informationRandom Matrices and Multivariate Statistical Analysis
Random Matrices and Multivariate Statistical Analysis Iain Johnstone, Statistics, Stanford imj@stanford.edu SEA 06@MIT p.1 Agenda Classical multivariate techniques Principal Component Analysis Canonical
More informationLow-temperature random matrix theory at the soft edge. Abstract
Low-temperature random matrix theory at the soft edge Alan Edelman Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 239, USA Per-Olof Persson Department of Mathematics, University
More informationOrthogonal Polynomial Ensembles
Chater 11 Orthogonal Polynomial Ensembles 11.1 Orthogonal Polynomials of Scalar rgument Let wx) be a weight function on a real interval, or the unit circle, or generally on some curve in the comlex lane.
More informationPainlevé Representations for Distribution Functions for Next-Largest, Next-Next-Largest, etc., Eigenvalues of GOE, GUE and GSE
Painlevé Representations for Distribution Functions for Next-Largest, Next-Next-Largest, etc., Eigenvalues of GOE, GUE and GSE Craig A. Tracy UC Davis RHPIA 2005 SISSA, Trieste 1 Figure 1: Paul Painlevé,
More informationWhat is the Riemann Hypothesis for Zeta Functions of Irregular Graphs?
What is the Riemann Hypothesis for Zeta Functions of Irregular Graphs? Audrey Terras Banff February, 2008 Joint work with H. M. Stark, M. D. Horton, etc. Introduction The Riemann zeta function for Re(s)>1
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part II Manuela Girotti based on M. Girotti s PhD thesis, A. Kuijlaars notes from Les Houches Winter School 202 and B. Eynard s notes
More information1 Tridiagonal matrices
Lecture Notes: β-ensembles Bálint Virág Notes with Diane Holcomb 1 Tridiagonal matrices Definition 1. Suppose you have a symmetric matrix A, we can define its spectral measure (at the first coordinate
More informationRandom Matrices: Invertibility, Structure, and Applications
Random Matrices: Invertibility, Structure, and Applications Roman Vershynin University of Michigan Colloquium, October 11, 2011 Roman Vershynin (University of Michigan) Random Matrices Colloquium 1 / 37
More informationEigenvalues and Singular Values of Random Matrices: A Tutorial Introduction
Random Matrix Theory and its applications to Statistics and Wireless Communications Eigenvalues and Singular Values of Random Matrices: A Tutorial Introduction Sergio Verdú Princeton University National
More informationFree Probability Theory and Random Matrices. Roland Speicher Queen s University Kingston, Canada
Free Probability Theory and Random Matrices Roland Speicher Queen s University Kingston, Canada We are interested in the limiting eigenvalue distribution of N N random matrices for N. Usually, large N
More informationBulk scaling limits, open questions
Bulk scaling limits, open questions Based on: Continuum limits of random matrices and the Brownian carousel B. Valkó, B. Virág. Inventiones (2009). Eigenvalue statistics for CMV matrices: from Poisson
More informationRANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS
RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS David García-García May 13, 2016 Faculdade de Ciências da Universidade de Lisboa OVERVIEW Random Matrix Theory Introduction Matrix ensembles A sample computation:
More information1 Intro to RMT (Gene)
M705 Spring 2013 Summary for Week 2 1 Intro to RMT (Gene) (Also see the Anderson - Guionnet - Zeitouni book, pp.6-11(?) ) We start with two independent families of R.V.s, {Z i,j } 1 i
More informationRandom Correlation Matrices, Top Eigenvalue with Heavy Tails and Financial Applications
Random Correlation Matrices, Top Eigenvalue with Heavy Tails and Financial Applications J.P Bouchaud with: M. Potters, G. Biroli, L. Laloux, M. A. Miceli http://www.cfm.fr Portfolio theory: Basics Portfolio
More informationEigenvalues and Eigenvectors
Sec. 6.1 Eigenvalues and Eigenvectors Linear transformations L : V V that go from a vector space to itself are often called linear operators. Many linear operators can be understood geometrically by identifying
More informationUniversality of local spectral statistics of random matrices
Universality of local spectral statistics of random matrices László Erdős Ludwig-Maximilians-Universität, Munich, Germany CRM, Montreal, Mar 19, 2012 Joint with P. Bourgade, B. Schlein, H.T. Yau, and J.
More informationFluctuations from the Semicircle Law Lecture 4
Fluctuations from the Semicircle Law Lecture 4 Ioana Dumitriu University of Washington Women and Math, IAS 2014 May 23, 2014 Ioana Dumitriu (UW) Fluctuations from the Semicircle Law Lecture 4 May 23, 2014
More informationHow long does it take to compute the eigenvalues of a random symmetric matrix?
How long does it take to compute the eigenvalues of a random symmetric matrix? Govind Menon (Brown University) Joint work with: Christian Pfrang (Ph.D, Brown 2011) Percy Deift, Tom Trogdon (Courant Institute)
More informationRandom matrices: A Survey. Van H. Vu. Department of Mathematics Rutgers University
Random matrices: A Survey Van H. Vu Department of Mathematics Rutgers University Basic models of random matrices Let ξ be a real or complex-valued random variable with mean 0 and variance 1. Examples.
More informationx #{ p=prime p x }, as x logx
1 The Riemann zeta function for Re(s) > 1 ζ s -s ( ) -1 2 duality between primes & complex zeros of zeta using Hadamard product over zeros prime number theorem x #{ p=prime p x }, as x logx statistics
More informationJames Demmel UC Berkeley Math and EECS Depts. Joint work with Ioana Dumitriu, Olga Holtz, Plamen Koev
. Accurate and efficient expression evaluation and linear algebra, or Why it can be easier to compute accurate eigenvalues of a Vandermonde matrix than the accurate sum of 3 numbers James Demmel UC Berkeley
More informationInvertibility of random matrices
University of Michigan February 2011, Princeton University Origins of Random Matrix Theory Statistics (Wishart matrices) PCA of a multivariate Gaussian distribution. [Gaël Varoquaux s blog gael-varoquaux.info]
More informationTHE EFFICIENT EVALUATION OF THE HYPERGEOMETRIC FUNCTION OF A MATRIX ARGUMENT
MATHEMATICS OF COMPUTATION Volume 75, Number 54, Pages 833 846 S 5-5718(6)184- Article electronically published on January 19, 6 THE EFFICIENT EVALUATION OF THE HYPERGEOMETRIC FUNCTION OF A MATRIX ARGUMENT
More informationMaximal height of non-intersecting Brownian motions
Maximal height of non-intersecting Brownian motions G. Schehr Laboratoire de Physique Théorique et Modèles Statistiques CNRS-Université Paris Sud-XI, Orsay Collaborators: A. Comtet (LPTMS, Orsay) P. J.
More informationMOPS: Multivariate orthogonal polynomials (symbolically)
Journal of Symbolic Computation 42 2007) 587 620 www.elsevier.com/locate/jsc MOPS: Multivariate orthogonal polynomials symbolically) Ioana Dumitriu a,, Alan Edelman b, Gene Shuman c a Department of Mathematics,
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT -09 Computational and Sensitivity Aspects of Eigenvalue-Based Methods for the Large-Scale Trust-Region Subproblem Marielba Rojas, Bjørn H. Fotland, and Trond Steihaug
More informationMeixner matrix ensembles
Meixner matrix ensembles W lodek Bryc 1 Cincinnati April 12, 2011 1 Based on joint work with Gerard Letac W lodek Bryc (Cincinnati) Meixner matrix ensembles April 12, 2011 1 / 29 Outline of talk Random
More informationTOP EIGENVALUE OF CAUCHY RANDOM MATRICES
TOP EIGENVALUE OF CAUCHY RANDOM MATRICES with Satya N. Majumdar, Gregory Schehr and Dario Villamaina Pierpaolo Vivo (LPTMS - CNRS - Paris XI) Gaussian Ensembles N = 5 Semicircle Law LARGEST EIGENVALUE
More informationQuantum chaos on graphs
Baylor University Graduate seminar 6th November 07 Outline 1 What is quantum? 2 Everything you always wanted to know about quantum but were afraid to ask. 3 The trace formula. 4 The of Bohigas, Giannoni
More informationMarkov operators, classical orthogonal polynomial ensembles, and random matrices
Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random
More informationSEA s workshop- MIT - July 10-14
Matrix-valued Stochastic Processes- Eigenvalues Processes and Free Probability SEA s workshop- MIT - July 10-14 July 13, 2006 Outline Matrix-valued stochastic processes. 1- definition and examples. 2-
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationRandom Matrix: From Wigner to Quantum Chaos
Random Matrix: From Wigner to Quantum Chaos Horng-Tzer Yau Harvard University Joint work with P. Bourgade, L. Erdős, B. Schlein and J. Yin 1 Perhaps I am now too courageous when I try to guess the distribution
More informationRandom Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws. Symeon Chatzinotas February 11, 2013 Luxembourg
Random Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws Symeon Chatzinotas February 11, 2013 Luxembourg Outline 1. Random Matrix Theory 1. Definition 2. Applications 3. Asymptotics 2. Ensembles
More informationA determinantal formula for the GOE Tracy-Widom distribution
A determinantal formula for the GOE Tracy-Widom distribution Patrik L. Ferrari and Herbert Spohn Technische Universität München Zentrum Mathematik and Physik Department e-mails: ferrari@ma.tum.de, spohn@ma.tum.de
More informationStatistical Inference and Random Matrices
Statistical Inference and Random Matrices N.S. Witte Institute of Fundamental Sciences Massey University New Zealand 5-12-2017 Joint work with Peter Forrester 6 th Wellington Workshop in Probability and
More informationLimits of Interlacing Eigenvalues in the Tridiagonal β-hermite Matrix Model
Limits of Interlacing Eigenvalues in the Tridiagonal β-hermite Matrix Model Gopal K. Goel and Mentor Andrew Ahn PRIMES Conference 2017 Gopal K. Goel and Mentor Andrew Ahn About Beamer May 20, 2017 1 /
More informationThe Prime Unsolved Problem in Mathematics
The Prime Unsolved Problem in Mathematics March, 0 Augustana Math Club 859: Bernard Riemann proposes the Riemann Hypothesis. 859: Bernard Riemann proposes the Riemann Hypothesis. 866: Riemann dies at age
More informationLecture 1: Riemann, Dedekind, Selberg, and Ihara Zetas
Lecture 1: Riemann, Dedekind, Selberg, and Ihara Zetas Audrey Terras U.C.S.D. 2008 more details can be found in my webpage: www.math.ucsd.edu /~aterras/ newbook.pdf First the Riemann Zeta 1 The Riemann
More informationMOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS
MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties
More informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
More informationThe Tracy-Widom distribution is not infinitely divisible.
The Tracy-Widom distribution is not infinitely divisible. arxiv:1601.02898v1 [math.pr] 12 Jan 2016 J. Armando Domínguez-Molina Facultad de Ciencias Físico-Matemáticas Universidad Autónoma de Sinaloa, México
More informationSmall Ball Probability, Arithmetic Structure and Random Matrices
Small Ball Probability, Arithmetic Structure and Random Matrices Roman Vershynin University of California, Davis April 23, 2008 Distance Problems How far is a random vector X from a given subspace H in
More informationLecture 11: CMSC 878R/AMSC698R. Iterative Methods An introduction. Outline. Inverse, LU decomposition, Cholesky, SVD, etc.
Lecture 11: CMSC 878R/AMSC698R Iterative Methods An introduction Outline Direct Solution of Linear Systems Inverse, LU decomposition, Cholesky, SVD, etc. Iterative methods for linear systems Why? Matrix
More informationFredholm determinant with the confluent hypergeometric kernel
Fredholm determinant with the confluent hypergeometric kernel J. Vasylevska joint work with I. Krasovsky Brunel University Dec 19, 2008 Problem Statement Let K s be the integral operator acting on L 2
More informationThe circular law. Lewis Memorial Lecture / DIMACS minicourse March 19, Terence Tao (UCLA)
The circular law Lewis Memorial Lecture / DIMACS minicourse March 19, 2008 Terence Tao (UCLA) 1 Eigenvalue distributions Let M = (a ij ) 1 i n;1 j n be a square matrix. Then one has n (generalised) eigenvalues
More informationarxiv: v1 [math.pr] 22 Dec 2018
arxiv:1812.09618v1 [math.pr] 22 Dec 2018 Operator norm upper bound for sub-gaussian tailed random matrices Eric Benhamou Jamal Atif Rida Laraki December 27, 2018 Abstract This paper investigates an upper
More informationLocal semicircle law, Wegner estimate and level repulsion for Wigner random matrices
Local semicircle law, Wegner estimate and level repulsion for Wigner random matrices László Erdős University of Munich Oberwolfach, 2008 Dec Joint work with H.T. Yau (Harvard), B. Schlein (Cambrigde) Goal:
More informationMultiple orthogonal polynomials. Bessel weights
for modified Bessel weights KU Leuven, Belgium Madison WI, December 7, 2013 Classical orthogonal polynomials The (very) classical orthogonal polynomials are those of Jacobi, Laguerre and Hermite. Classical
More informationExperimental mathematics and integration
Experimental mathematics and integration David H. Bailey http://www.davidhbailey.com Lawrence Berkeley National Laboratory (retired) Computer Science Department, University of California, Davis October
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationNonlinear Integral Equation Formulation of Orthogonal Polynomials
Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,
More informationResearch Article On the Marginal Distribution of the Diagonal Blocks in a Blocked Wishart Random Matrix
Hindawi Publishing Corporation nternational Journal of Analysis Volume 6, Article D 59678, 5 pages http://dx.doi.org/.55/6/59678 Research Article On the Marginal Distribution of the Diagonal Blocks in
More informationEigenvalue PDFs. Peter Forrester, M&S, University of Melbourne
Outline Eigenvalue PDFs Peter Forrester, M&S, University of Melbourne Hermitian matrices with real, complex or real quaternion elements Circular ensembles and classical groups Products of random matrices
More informationAn Invitation to Modern Number Theory. Steven J. Miller and Ramin Takloo-Bighash PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD
An Invitation to Modern Number Theory Steven J. Miller and Ramin Takloo-Bighash PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Foreword Preface Notation xi xiii xix PART 1. BASIC NUMBER THEORY
More informationLecture 5 Singular value decomposition
Lecture 5 Singular value decomposition Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationMahler measure as special values of L-functions
Mahler measure as special values of L-functions Matilde N. Laĺın Université de Montréal mlalin@dms.umontreal.ca http://www.dms.umontreal.ca/~mlalin CRM-ISM Colloquium February 4, 2011 Diophantine equations
More informationA refined Lanczos method for computing eigenvalues and eigenvectors of unsymmetric matrices
A refined Lanczos method for computing eigenvalues and eigenvectors of unsymmetric matrices Jean Christophe Tremblay and Tucker Carrington Chemistry Department Queen s University 23 août 2007 We want to
More informationNumerical Evaluation of Standard Distributions in Random Matrix Theory
Numerical Evaluation of Standard Distributions in Random Matrix Theory A Review of Folkmar Bornemann s MATLAB Package and Paper Matt Redmond Department of Mathematics Massachusetts Institute of Technology
More informationNumerical analysis and random matrix theory. Tom Trogdon UC Irvine
Numerical analysis and random matrix theory Tom Trogdon ttrogdon@math.uci.edu UC Irvine Acknowledgements This is joint work with: Percy Deift Govind Menon Sheehan Olver Raj Rao Numerical analysis and random
More information8. Diagonalization.
8. Diagonalization 8.1. Matrix Representations of Linear Transformations Matrix of A Linear Operator with Respect to A Basis We know that every linear transformation T: R n R m has an associated standard
More informationConsider the following example of a linear system:
LINEAR SYSTEMS Consider the following example of a linear system: Its unique solution is x + 2x 2 + 3x 3 = 5 x + x 3 = 3 3x + x 2 + 3x 3 = 3 x =, x 2 = 0, x 3 = 2 In general we want to solve n equations
More informationLimits of spiked random matrices I
Limits of spiked random matrices I Alex Bloemendal Bálint Virág arxiv:1011.1877v2 [math.pr] 16 Sep 2011 September 16, 2011 Abstract Given a large, high-dimensional sample from a spiked population, the
More informationA Remark on Hypercontractivity and Tail Inequalities for the Largest Eigenvalues of Random Matrices
A Remark on Hypercontractivity and Tail Inequalities for the Largest Eigenvalues of Random Matrices Michel Ledoux Institut de Mathématiques, Université Paul Sabatier, 31062 Toulouse, France E-mail: ledoux@math.ups-tlse.fr
More informationA Review of Matrix Analysis
Matrix Notation Part Matrix Operations Matrices are simply rectangular arrays of quantities Each quantity in the array is called an element of the matrix and an element can be either a numerical value
More informationMIMO Capacities : Eigenvalue Computation through Representation Theory
MIMO Capacities : Eigenvalue Computation through Representation Theory Jayanta Kumar Pal, Donald Richards SAMSI Multivariate distributions working group Outline 1 Introduction 2 MIMO working model 3 Eigenvalue
More informationRandom transition-rate matrices for the master equation
PHYSICAL REVIEW E 8, 4 9 Random transition-rate matrices for the master equation Carsten Timm* Institute for Theoretical Physics, Technische Universität Dresden, 6 Dresden, Germany Received May 9; published
More informationSolution of eigenvalue problems. Subspace iteration, The symmetric Lanczos algorithm. Harmonic Ritz values, Jacobi-Davidson s method
Solution of eigenvalue problems Introduction motivation Projection methods for eigenvalue problems Subspace iteration, The symmetric Lanczos algorithm Nonsymmetric Lanczos procedure; Implicit restarts
More informationNumerical Methods I Eigenvalue Problems
Numerical Methods I Eigenvalue Problems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 2nd, 2014 A. Donev (Courant Institute) Lecture
More information1 Extrapolation: A Hint of Things to Come
Notes for 2017-03-24 1 Extrapolation: A Hint of Things to Come Stationary iterations are simple. Methods like Jacobi or Gauss-Seidel are easy to program, and it s (relatively) easy to analyze their convergence.
More informationNotes on Eigenvalues, Singular Values and QR
Notes on Eigenvalues, Singular Values and QR Michael Overton, Numerical Computing, Spring 2017 March 30, 2017 1 Eigenvalues Everyone who has studied linear algebra knows the definition: given a square
More informationarxiv:chao-dyn/ v1 3 Jul 1995
Chaotic Spectra of Classically Integrable Systems arxiv:chao-dyn/9506014v1 3 Jul 1995 P. Crehan Dept. of Mathematical Physics, University College Dublin, Belfield, Dublin 2, Ireland PCREH89@OLLAMH.UCD.IE
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationMeasures and Jacobians of Singular Random Matrices. José A. Díaz-Garcia. Comunicación de CIMAT No. I-07-12/ (PE/CIMAT)
Measures and Jacobians of Singular Random Matrices José A. Díaz-Garcia Comunicación de CIMAT No. I-07-12/21.08.2007 (PE/CIMAT) Measures and Jacobians of singular random matrices José A. Díaz-García Universidad
More informationEcon Slides from Lecture 8
Econ 205 Sobel Econ 205 - Slides from Lecture 8 Joel Sobel September 1, 2010 Computational Facts 1. det AB = det BA = det A det B 2. If D is a diagonal matrix, then det D is equal to the product of its
More informationEcon Slides from Lecture 7
Econ 205 Sobel Econ 205 - Slides from Lecture 7 Joel Sobel August 31, 2010 Linear Algebra: Main Theory A linear combination of a collection of vectors {x 1,..., x k } is a vector of the form k λ ix i for
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More informationAiry and Pearcey Processes
Airy and Pearcey Processes Craig A. Tracy UC Davis Probability, Geometry and Integrable Systems MSRI December 2005 1 Probability Space: (Ω, Pr, F): Random Matrix Models Gaussian Orthogonal Ensemble (GOE,
More informationARPACK. Dick Kachuma & Alex Prideaux. November 3, Oxford University Computing Laboratory
ARPACK Dick Kachuma & Alex Prideaux Oxford University Computing Laboratory November 3, 2006 What is ARPACK? ARnoldi PACKage Collection of routines to solve large scale eigenvalue problems Developed at
More informationCalculating determinants for larger matrices
Day 26 Calculating determinants for larger matrices We now proceed to define det A for n n matrices A As before, we are looking for a function of A that satisfies the product formula det(ab) = det A det
More informationLARGE SPARSE EIGENVALUE PROBLEMS. General Tools for Solving Large Eigen-Problems
LARGE SPARSE EIGENVALUE PROBLEMS Projection methods The subspace iteration Krylov subspace methods: Arnoldi and Lanczos Golub-Kahan-Lanczos bidiagonalization General Tools for Solving Large Eigen-Problems
More informationarxiv: v1 [math.pr] 22 May 2008
THE LEAST SINGULAR VALUE OF A RANDOM SQUARE MATRIX IS O(n 1/2 ) arxiv:0805.3407v1 [math.pr] 22 May 2008 MARK RUDELSON AND ROMAN VERSHYNIN Abstract. Let A be a matrix whose entries are real i.i.d. centered
More informationName: INSERT YOUR NAME HERE. Due to dropbox by 6pm PDT, Wednesday, December 14, 2011
AMath 584 Name: INSERT YOUR NAME HERE Take-home Final UWNetID: INSERT YOUR NETID Due to dropbox by 6pm PDT, Wednesday, December 14, 2011 The main part of the assignment (Problems 1 3) is worth 80 points.
More informationTuring and the Riemann zeta function
Turing and the Riemann zeta function Andrew Odlyzko School of Mathematics University of Minnesota odlyzko@umn.edu http://www.dtc.umn.edu/ odlyzko May 11, 2012 Andrew Odlyzko ( School of Mathematics University
More informationInterpolating between Wishart and inverse-wishart distributions
Interpolating between Wishart and inverse-wishart distributions Topological phase transitions in 1D multichannel disordered wires with a chiral symmetry Christophe Texier December 11, 2015 with Aurélien
More information