Preconditioners for ill conditioned (block) Toeplitz systems: facts a
|
|
- Milton Harrison
- 6 years ago
- Views:
Transcription
1 Preconditioners for ill conditioned (block) Toeplitz systems: facts and ideas Department of Informatics, Athens University of Economics and Business, Athens, Greece. Joint work with D. Noutsos
2 Problem We are interested in the fast and efficient solution of nm nm BTTB systems T n,m (f)x = b where f is nonnegative real-valued belonging to C 2π,2π defined in the fundamental domain Q = ( π,π] 2 and is a priori known. The entries of the coefficient matrix are given by t j,k = 1 4π 2 Q f(x, y)e i(jx+ky) dxdy, for j = 0, ±1,...,±(n 1) and k = 0,...,±(m 1).
3 Connection between f and T nm (f) The main connection between T nm (f) and the generating function is described by the the following result: Theorem If f C[ π,π] 2 and c < C the extreme values of f(x, y) on Q. Then every eigenvalue λ of the block Toeplitz matrix T satisfies the strict inequalities c < λ < C Moreover as m, n then λ min (T nm (f)) c and λ max (T nm (f)) C
4 Crucial relationship The basis for the construction of effective preconditioners is described by the following theorem Theorem Let f, g 0 C[ π,π] 2 (f and g not identically zero). Then for every m, n the matrix Tnm 1 (g)t nm (f) has eigenvalues in the open interval (r, R), where r = inf Q f g and R = sup Q f g
5 Negative result for τ algebra preconditioners Theorem (Noutsos,Serra,Vassalos, TCS (2004)) Let f be equivalent to p k (x, y) = (2 2 cos(x)) k + (2 2 cos(y)) k with k 2 and let β be a fixed positive number independent of n. Then for every sequence {P n } with P n τ, n = (n 1, n 2 ), and such that uniformly with respect to ˆn, we have λ max (P 1 n T n (f)) β (1) (a) the minimal eigenvalue of Pn 1 T n (f) tends to zero. (b) the number #{λ(n) σ(p 1 n T n(f)) : λ(n) N(n) 0} tends to infinity as N(ˆn) tends to infinity.
6 Negative result for τ algebra preconditioners Theorem (Noutsos,Serra,Vassalos, TCS (2004)) Let f be equivalent to p k (x, y) = (2 2 cos(x)) k + (2 2 cos(y)) k with k 2 and let α be a fixed positive number independent of n = (n 1, n 2 ). Then for every sequence {P n } with P n τ and such that λ min (P 1 n T n(f)) α (2) uniformly with respect to n, we have (a) the maximal eigenvalue of P 1 n T n (f) tends to. (b) the number #{λ(n) σ(p 1 n T n(f)) : λ(n) N(n) } tends to infinity as N(n) tends to infinity.
7 How to solve: 2D- ill conditioned problem Direct methods: Levinson type methods cost O(n 2 m 3 ) ops while superfast methods O(nm 3 log 2 n). Stability problems. Not optimal. PCG method with matrix algebra preconditioners. Cost O(k(ǫ)nm log nm), where k(ǫ) is the required number of iterations and depends from the condition number of T nm. PCG method where the preconditioner is band Toeplitz matrix. Under some assumptions, the cost is optimal (O(nm log nm)). Multigrid methods: Promising but in early stages. Cost O(nm log nm) ops. G. Fiorentino and S. Serra-Capizzano (1996), T. Huckle and J. Staudacher (2002),H.W. Sun, X.Q. Jin and Q.S. Chang (2004), A. Arico, M. Donatelli and S. Serra-Capizzano (2004).
8 More on Band Preconditioners S. Serra-Capizzano (BIT, (1994)) and M. Ng (LAA, (1997)) proposed as preconditioner the band BTTB matrix generated by the minimum trigonometric polynomial g which has the same roots with f. Let f = g h with h > 0. Then D. Noutsos, S. Serra Capizzano and P. Vassalos (Numer. Math. (2006)) proposed as preconditioners the band BTTB matrix generated by g ĥ where ĥ: 1 is the trigonometric polynomial arises from the Fourier approximation on h. 2 arises from the Lagrange interpolation of h at 2D Chebyshev points, or from the interpolation of h using the 2D Fejer kernel.
9 Preliminaries We assume that the nonnegative function f has isolated zeros (x 1, y 1 ),(x 2, y 2 ),...,(x k, y k ), on Q each one of multiplicities (2µ 1, 2ν 1 ),(2µ 2, 2ν 2 ),..., (2µ k, 2ν k ). Then, f can be written as where g = f = g w k [(2 2 cos(x x i )) µ i + (2 2 cos(y y i )) ν i ] i=1 and w, is strictly positive on Q.
10 New Proposal For the system T nm (f)x = b we define and propose as a preconditioner the product of matrices K nm (f) = A nm ( w)t nm (g)a nm ( w) = A nm (h)t nm (g)a nm (h) with h = w, A nm {τ, C, H}, where {τ, C, H} is the set of matrices belonging to Block τ, Block Circulant and Block Hartley algebra, respectively.
11 Construction of 2D algebras A v [n] Q n τ v [n] i = πi 2 n+1 n+1 ( C = 2πi n F n = 1 n H v [n] i v [n] i ( sin(jv [n] i ) ) n e ijv[n] i = 2πi n Re(F n ) + Im(F n ) ) i,j=1 The matrices C(h),τ(h), H(h) can be written as ( ) A nm (h) = Q nm Diag f(v [nm] ) Qnm H where v [nm] = v [n] v [m] and Q nm = Q n Qm
12 Obviously K nm (f) has all the properties that a preconditioner must satisfied, i.e, is symmetric, is positive definite, The cost for the solution of the arbitrary system K nm (f)x = b is of order O(nm log nm). O(nm log nm) for the inversion of A nm (h) by 2D FFT, and O(nm) for the inversion of block band Toeplitz matrix T n,m (g) which can be done by multigrid methods. So, the only condition that must be fulfill, in order to be a competitive preconditioner, is the spectrum of [ K A nm (f) ] 1 Tnm (f) being bounded from above and below.
13 Useful Definition Definition We say that a function h is a ((k 1, k 2 ),(x 0, y 0 ))-smooth function if l 1+l 2 x l 1 y l 2 h(x 0, y 0 ) = 0, l 1 < k 1, l 2 < k 2, and l 1 +l 2 < max{k 1, k 2 } and l 1+l 2 x l 1 y l 2 h(x 0, y 0 ) is bounded for l 1 = k 1, l 2 = 0 and l 2 = k 2, l 1 = 0, and l 1 + l 2 = max{k 1, k 2 }, l 1 < k 1, l 2 < k 2.
14 τ Case: Clustering Theorem Let f belongs to the Wiener class. Then for every ǫ the spectrum of [K τ nm(f)] 1 T nm (f) lies in [1 ǫ, 1 + ǫ], for n, m sufficient large, except of O(m + n) outliers. Thus, we have weak clustering around unity of the spectrum of the preconditioned matrix.
15 τ Case: Bounds Theorem Let f C 2π,2π even function on Q with roots (x 1, y 1 ),(x 2, y 2 ),...,(x k, y k ), each one of multiplicities (2µ 1, 2ν 1 ), (2µ 2, 2ν 2 ),..., (2µ k, 2ν k ), respectively. If g is the trigonometric polynomial of minimal degree that rises the roots of f and h is (µ i 1,ν i 1) smooth function at the roots (x i, y i ), i = 1(1)k, then the spectrum of the preconditioned matrix P τ = [K τ nm(f)] 1 T nm (f) is bounded from above as well as bellow: c < λ min (P τ ) < λ max (P τ ) < C with c, C constants independent of n, m
16 Circulant Case: Clustering Theorem Let f belongs to the Wiener class on Q. Then for every ǫ the spectrum of [ K C nm (f) ] 1 Tnm (f) lies in [1 ǫ, 1 + ǫ] for n, m sufficient large except O(m + n) outliers. Thus, we have a weak clustering of the spectrum of the preconditioned matrix around unity.
17 Circulant Case: Bound of Spectrum Theorem Let f C 2π,2π even function on Q with roots (x 1, y 1 ),(x 2, y 2 ),...,(x k, y k ), on Q each one of multiplicities (2µ 1, 2ν 1 ),(2µ 2, 2ν 2 ),..., (2µ k, 2ν k ) respectively. If g is the trigonometric polynomial of minimal degree that rises the roots of f and h is (µ i,ν i ) smooth function at the roots (x i, y i ), i = 1(1)k then the spectrum of the preconditioned matrix P C = [ K C nm(f) ] 1 Tnm (f) is bounded above as well as bellow: c < λ min (P C ) < λ max (P C ) < C with c, C constants independent of n, m
18 Case where h is not smooth enough Let us suppose that f has roots at (x i, y i ), i = 1(1)k. If h = f g is not as smooth as the previous Theorems demand, then the spectrum of the preconditioned matrix could be unbound. For that case we make a smoothing correction and instead of h we use ĥ defined as { h(x) (x, y) Q/ Ω ĥ = i h(x i, y i ) + α i g i (x, y) (x, y) Ω i where Ω i = {(x, y) : (x i, y i ) (x, y) < ǫ i } and α i is defined such that h(ǫ i,ǫ i ) = h(x i, y i ) + α i g(ǫ i,ǫ i )
19 Figure: Smoothing of h(x, y) = ( x + y +1)(x 4 +y 4 ), by ĥ(x, y). (2 2cos(x)) 2 +(2 2cos(y)) 2
20 Numerical Experiments We compare our proposal with the already known band preconditioners: B, which is generated by the trigonometric polynomial g that raises the roots of f K which is the product g ĥ with ĥ being the trigonometric polynomial arising from the 2D Fejer kernel on the positive part h of f P and F which arise: from the Lagrange interpolation of h at 2D Chebyshev points, and from the approximation of h using the 2D Fourier expansion, respectively. For all the tests the stopping criterion was r k 2 r 0 2 < 10 5, the starting vector the zero one and the righthand side vector of the system was (1 1 1) T
21 Table: f 1 (x, y) = (x 4 + y 2 )( x + y 3 + 1) n = m B K 4,4 F 4,4 P 4,4 τ C
22 Table: f 2 (x, y) = (x 2 + y 2 )(x 4 + y 4 +.1) n = m B K 6,6 F 6,6 P 6,6 τ C * * * 24 * * *: The number of iterations exceeds 100 -: The matrix is singular.
23 Application to Toeplitz plus band systems Let f C 2π,2π even function on Q with roots (x 1, y 1 ),(x 2, y 2 ),...,(x k, y k ), on Q each one of multiplicities (2µ 1, 2ν 1 ),(2µ 2, 2ν 2 ),..., (2µ k, 2ν k ) respectively. Consider the system (T nm (f) + B)x = b, where B is a symmetric, p.d block band matrix. We can choose as preconditioner the A nm ( w)(t nm (g) + B)A nm ( w) where {τ, C, H} is the set of matrices belonging to Block τ, Block Circulant and Block Hartley algebra, respectively, and g(x, y) the trigonometric polynomial having the same roots with the same multiplicities with f.
24 Table: f(x, y) = ( x y )( x + y + 2) n = m I B τ(f) * 98 9
25 τ preconditioners in ill condition case We define as P nm the matrix Q nm Diag ( ) f(v [nm] ) Qnm H where Q nm is the orthogonal matrix diagonalize the τ algebra and v [nm] ij = ( πi n+1, πj m+1 ). Then, the following theorem holds true: Theorem Let f C 2π,2π even function on Q with root at (0, 0) with multiplicity (q, r) with q, r 2. Then, the spectrum of PnmT 1 nm (f) is clustered around unity. Moreover, for every n, m it holds that c < σ(p 1 nmt nm (f)) < C with c, C > 0 independent of n, m
26 Table: f(x, y) = ( x + y )(x 2 + y 2 + 1) n = m #I λ min (I) λ max (P) λ min (P) #P
27 Thank you very much for your attention!
Superlinear convergence for PCG using band plus algebra preconditioners for Toeplitz systems
Computers and Mathematics with Applications 56 (008) 155 170 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Superlinear
More informationThe Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In
The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed
More informationBLOCK DIAGONAL AND SCHUR COMPLEMENT PRECONDITIONERS FOR BLOCK TOEPLITZ SYSTEMS WITH SMALL SIZE BLOCKS
BLOCK DIAGONAL AND SCHUR COMPLEMENT PRECONDITIONERS FOR BLOCK TOEPLITZ SYSTEMS WITH SMALL SIZE BLOCKS WAI-KI CHING, MICHAEL K NG, AND YOU-WEI WEN Abstract In this paper we consider the solution of Hermitian
More informationToeplitz matrices: spectral properties and preconditioning in the CG method
Toeplitz matrices: spectral properties and preconditioning in the CG method Stefano Serra-Capizzano December 5, 2014 Abstract We consider multilevel Toeplitz matrices T n (f) generated by Lebesgue integrable
More informationStructured Matrices, Multigrid Methods, and the Helmholtz Equation
Structured Matrices, Multigrid Methods, and the Helmholtz Equation Rainer Fischer 1, Thomas Huckle 2 Technical University of Munich, Germany Abstract Discretization of the Helmholtz equation with certain
More informationPreconditioners for ill conditioned Toeplitz systems constructed from positive kernels
Preconditioners for ill conditioned Toeplitz systems constructed from positive kernels Daniel Potts Medizinische Universität Lübeck Institut für Mathematik Wallstr. 40 D 23560 Lübeck potts@math.mu-luebeck.de
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationApproximate Inverse-Free Preconditioners for Toeplitz Matrices
Approximate Inverse-Free Preconditioners for Toeplitz Matrices You-Wei Wen Wai-Ki Ching Michael K. Ng Abstract In this paper, we propose approximate inverse-free preconditioners for solving Toeplitz systems.
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationA Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations
A Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations Jin Yun Yuan Plamen Y. Yalamov Abstract A method is presented to make a given matrix strictly diagonally dominant
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
4 Interpolation 4.1 Polynomial interpolation Problem: LetP n (I), n ln, I := [a,b] lr, be the linear space of polynomials of degree n on I, P n (I) := { p n : I lr p n (x) = n i=0 a i x i, a i lr, 0 i
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More informationToeplitz-circulant Preconditioners for Toeplitz Systems and Their Applications to Queueing Networks with Batch Arrivals Raymond H. Chan Wai-Ki Ching y
Toeplitz-circulant Preconditioners for Toeplitz Systems and Their Applications to Queueing Networks with Batch Arrivals Raymond H. Chan Wai-Ki Ching y November 4, 994 Abstract The preconditioned conjugate
More informationarxiv: v1 [math.na] 21 Jan 2019
Approximating the Perfect Sampling Grids for Computing the Eigenvalues of Toeplitz-like Matrices Using the Spectral Symbol Sven-Erik Ekström see@2pi.se Athens University of Economics and Business January
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 informationSpectral analysis and numerical methods for fractional diffusion equations
Spectral analysis and numerical methods for fractional diffusion equations M. Dehghan M. Donatelli M. Mazza H. Moghaderi S. Serra-Capizzano Department of Science and High Technology, University of Insubria,
More informationApproximation theory
Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2
More informationDimitrios Noutsos Professor, Department of Mathematics, University of Ioannina
Dimitrios Noutsos Professor, Department of Mathematics, University of Ioannina PERSONAL INFO Name: Dimitrios Noutsos Date of birth: February 10, 1953. Position of birth: Neochori Ioannina Nationality:
More informationClarkson Inequalities With Several Operators
isid/ms/2003/23 August 14, 2003 http://www.isid.ac.in/ statmath/eprints Clarkson Inequalities With Several Operators Rajendra Bhatia Fuad Kittaneh Indian Statistical Institute, Delhi Centre 7, SJSS Marg,
More information0 1 n. (d(v) = diag (z i )) in terms of the Fourier matrix F = [ω (i 1)(j 1) ] n ij=1 / n,
Let be a n n matrix and L some space of matrices of the same order, but of lower complexity This means, for instance in case is non singular, that solving linear systems Xz = b, X L, is much easier than
More informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More information1 Introduction The conjugate gradient (CG) method is an iterative method for solving Hermitian positive denite systems Ax = b, see for instance Golub
Circulant Preconditioned Toeplitz Least Squares Iterations Raymond H. Chan, James G. Nagy y and Robert J. Plemmons z November 1991, Revised April 1992 Abstract We consider the solution of least squares
More informationSpectral features of matrix-sequences, GLT, symbol, and application in preconditioning Krylov methods, image deblurring, and multigrid algorithms
UNIVERSITÀ DEGLI STUDI DELL INSUBRIA Dipartimento di Scienza ed Alta Tecnologia Dottorato di Ricerca in Matematica del Calcolo: Modelli, Strutture, Algoritmi ed Applicazioni Spectral features of matrix-sequences,
More informationDisplacement structure approach to discrete-trigonometric-transform based preconditioners of G.Strang type and of T.Chan type
Displacement structure approach to discrete-trigonometric-transform based preconditioners of GStrang type and of TChan type Thomas Kailath Information System Laboratory, Stanford University, Stanford CA
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More information7: FOURIER SERIES STEVEN HEILMAN
7: FOURIER SERIES STEVE HEILMA Contents 1. Review 1 2. Introduction 1 3. Periodic Functions 2 4. Inner Products on Periodic Functions 3 5. Trigonometric Polynomials 5 6. Periodic Convolutions 7 7. Fourier
More informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More information4. Multi-linear algebra (MLA) with Kronecker-product data.
ect. 3. Tensor-product interpolation. Introduction to MLA. B. Khoromskij, Leipzig 2007(L3) 1 Contents of Lecture 3 1. Best polynomial approximation. 2. Error bound for tensor-product interpolants. - Polynomial
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationKOROVKIN TESTS, APPROXIMATION, AND ERGODIC THEORY
MATHEMATICS OF COMPUTATION Volume 69, Number 232, Pages 1533 1558 S 0025-5718(00)01217-5 Article electronically published on March 6, 2000 KOROVKIN TESTS, APPROXIMATION, AND ERGODIC THEORY STEFANO SERRA
More informationNOTES ON LINEAR ODES
NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary
More informationFundamentals of Unconstrained Optimization
dalmau@cimat.mx Centro de Investigación en Matemáticas CIMAT A.C. Mexico Enero 2016 Outline Introduction 1 Introduction 2 3 4 Optimization Problem min f (x) x Ω where f (x) is a real-valued function The
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 informationInternational Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994
International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994 1 PROBLEMS AND SOLUTIONS First day July 29, 1994 Problem 1. 13 points a Let A be a n n, n 2, symmetric, invertible
More informationFourier Series. ,..., e ixn ). Conversely, each 2π-periodic function φ : R n C induces a unique φ : T n C for which φ(e ix 1
Fourier Series Let {e j : 1 j n} be the standard basis in R n. We say f : R n C is π-periodic in each variable if f(x + πe j ) = f(x) x R n, 1 j n. We can identify π-periodic functions with functions on
More informationLinear algebra issues in Interior Point methods for bound-constrained least-squares problems
Linear algebra issues in Interior Point methods for bound-constrained least-squares problems Stefania Bellavia Dipartimento di Energetica S. Stecco Università degli Studi di Firenze Joint work with Jacek
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationImage restoration: numerical optimisation
Image restoration: numerical optimisation Short and partial presentation Jean-François Giovannelli Groupe Signal Image Laboratoire de l Intégration du Matériau au Système Univ. Bordeaux CNRS BINP / 6 Context
More informationIsogEometric Tearing and Interconnecting
IsogEometric Tearing and Interconnecting Christoph Hofer and Ludwig Mitter Johannes Kepler University, Linz 26.01.2017 Doctoral Program Computational Mathematics Numerical Analysis and Symbolic Computation
More informationAn iterative multigrid regularization method for Toeplitz discrete ill-posed problems
NUMERICAL MATHEMATICS: Theory, Methods and Applications Numer. Math. Theor. Meth. Appl., Vol. xx, No. x, pp. 1-18 (200x) An iterative multigrid regularization method for Toeplitz discrete ill-posed problems
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationMatrix structures and applications Luminy, November, Dario A. Bini
Matrix structures and applications Luminy, November, 3 7 2014 Dario A Bini September 30, 2014 Introduction Matrix structures are encountered in various guises in a wide variety of mathematical problems
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2017/18 Part 3: Iterative Methods PD
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationMath Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.
Math 5620 - Introduction to Numerical Analysis - Class Notes Fernando Guevara Vasquez Version 1990. Date: January 17, 2012. 3 Contents 1. Disclaimer 4 Chapter 1. Iterative methods for solving linear systems
More information4.3 Eigen-Analysis of Spectral Derivative Matrices 203
. Eigen-Analysis of Spectral Derivative Matrices 8 N = 8 N = 6 6 6 8 8 8 8 N = 6 6 N = 6 6 6 (Round off error) 6 8 8 6 8 6 6 Fig... Legendre collocation first-derivative eigenvalues computed with 6-bit
More informationNumerical Experiments for Finding Roots of the Polynomials in Chebyshev Basis
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017), pp. 988 1001 Applications and Applied Mathematics: An International Journal (AAM) Numerical Experiments
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationMAT Linear Algebra Collection of sample exams
MAT 342 - Linear Algebra Collection of sample exams A-x. (0 pts Give the precise definition of the row echelon form. 2. ( 0 pts After performing row reductions on the augmented matrix for a certain system
More informationIn English, this means that if we travel on a straight line between any two points in C, then we never leave C.
Convex sets In this section, we will be introduced to some of the mathematical fundamentals of convex sets. In order to motivate some of the definitions, we will look at the closest point problem from
More informationApproximation Theory
Approximation Theory Function approximation is the task of constructing, for a given function, a simpler function so that the difference between the two functions is small and to then provide a quantifiable
More informationRadial Basis Functions I
Radial Basis Functions I Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo November 14, 2008 Today Reformulation of natural cubic spline interpolation Scattered
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson
Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation
More informationMath 489AB Exercises for Chapter 1 Fall Section 1.0
Math 489AB Exercises for Chapter 1 Fall 2008 Section 1.0 1.0.2 We want to maximize x T Ax subject to the condition x T x = 1. We use the method of Lagrange multipliers. Let f(x) = x T Ax and g(x) = x T
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationHW1 solutions. 1. α Ef(x) β, where Ef(x) is the expected value of f(x), i.e., Ef(x) = n. i=1 p if(a i ). (The function f : R R is given.
HW1 solutions Exercise 1 (Some sets of probability distributions.) Let x be a real-valued random variable with Prob(x = a i ) = p i, i = 1,..., n, where a 1 < a 2 < < a n. Of course p R n lies in the standard
More informationIn particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with
Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient
More informationHomework 2 Foundations of Computational Math 2 Spring 2019
Homework 2 Foundations of Computational Math 2 Spring 2019 Problem 2.1 (2.1.a) Suppose (v 1,λ 1 )and(v 2,λ 2 ) are eigenpairs for a matrix A C n n. Show that if λ 1 λ 2 then v 1 and v 2 are linearly independent.
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationBOUNDARY CONDITIONS FOR SYMMETRIC BANDED TOEPLITZ MATRICES: AN APPLICATION TO TIME SERIES ANALYSIS
BOUNDARY CONDITIONS FOR SYMMETRIC BANDED TOEPLITZ MATRICES: AN APPLICATION TO TIME SERIES ANALYSIS Alessandra Luati Dip. Scienze Statistiche University of Bologna Tommaso Proietti S.E.F. e ME. Q. University
More informationAdaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers
Adaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers Jan Mandel University of Colorado at Denver Bedřich Sousedík Czech Technical University
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
More informationReproducing Kernel Hilbert Spaces Class 03, 15 February 2006 Andrea Caponnetto
Reproducing Kernel Hilbert Spaces 9.520 Class 03, 15 February 2006 Andrea Caponnetto About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationPreconditioning of elliptic problems by approximation in the transform domain
TR-CS-97-2 Preconditioning of elliptic problems by approximation in the transform domain Michael K. Ng July 997 Joint Computer Science Technical Report Series Department of Computer Science Faculty of
More informationMath 321 Final Examination April 1995 Notation used in this exam: N. (1) S N (f,x) = f(t)e int dt e inx.
Math 321 Final Examination April 1995 Notation used in this exam: N 1 π (1) S N (f,x) = f(t)e int dt e inx. 2π n= N π (2) C(X, R) is the space of bounded real-valued functions on the metric space X, equipped
More informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
More informationChapter 8: Taylor s theorem and L Hospital s rule
Chapter 8: Taylor s theorem and L Hospital s rule Theorem: [Inverse Mapping Theorem] Suppose that a < b and f : [a, b] R. Given that f (x) > 0 for all x (a, b) then f 1 is differentiable on (f(a), f(b))
More informationNumerical Linear Algebra and. Image Restoration
Numerical Linear Algebra and Image Restoration Maui High Performance Computing Center Wednesday, October 8, 2003 James G. Nagy Emory University Atlanta, GA Thanks to: AFOSR, Dave Tyler, Stuart Jefferies,
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationSpectra of Multiplication Operators as a Numerical Tool. B. Vioreanu and V. Rokhlin Technical Report YALEU/DCS/TR-1443 March 3, 2011
We introduce a numerical procedure for the construction of interpolation and quadrature formulae on bounded convex regions in the plane. The construction is based on the behavior of spectra of certain
More informationSchur Complement Matrix And Its (Elementwise) Approximation: A Spectral Analysis Based On GLT Sequences
Schur Complement Matrix And Its (Elementwise) Approximation: A Spectral Analysis Based On GLT Sequences Ali Dorostkar, Maya Neytcheva, and Stefano Serra-Capizzano 2 Department of Information Technology,
More informationReal Analysis Qualifying Exam May 14th 2016
Real Analysis Qualifying Exam May 4th 26 Solve 8 out of 2 problems. () Prove the Banach contraction principle: Let T be a mapping from a complete metric space X into itself such that d(tx,ty) apple qd(x,
More informationWeighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 2, pp. 401 412 (2013) http://campus.mst.edu/adsa Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationLeast Squares Approximation
Chapter 6 Least Squares Approximation As we saw in Chapter 5 we can interpret radial basis function interpolation as a constrained optimization problem. We now take this point of view again, but start
More information1 Singular Value Decomposition
1 Singular Value Decomposition Factorisation of rectangular matrix (generalisation of eigenvalue concept / spectral theorem): For every matrix A C m n there exists a factorisation A = UΣV U C m m, V C
More informationMATH 4211/6211 Optimization Constrained Optimization
MATH 4211/6211 Optimization Constrained Optimization Xiaojing Ye Department of Mathematics & Statistics Georgia State University Xiaojing Ye, Math & Stat, Georgia State University 0 Constrained optimization
More informationJae Heon Yun and Yu Du Han
Bull. Korean Math. Soc. 39 (2002), No. 3, pp. 495 509 MODIFIED INCOMPLETE CHOLESKY FACTORIZATION PRECONDITIONERS FOR A SYMMETRIC POSITIVE DEFINITE MATRIX Jae Heon Yun and Yu Du Han Abstract. We propose
More information2. The Concept of Convergence: Ultrafilters and Nets
2. The Concept of Convergence: Ultrafilters and Nets NOTE: AS OF 2008, SOME OF THIS STUFF IS A BIT OUT- DATED AND HAS A FEW TYPOS. I WILL REVISE THIS MATE- RIAL SOMETIME. In this lecture we discuss two
More informationUC Berkeley Department of Electrical Engineering and Computer Science. EECS 227A Nonlinear and Convex Optimization. Solutions 5 Fall 2009
UC Berkeley Department of Electrical Engineering and Computer Science EECS 227A Nonlinear and Convex Optimization Solutions 5 Fall 2009 Reading: Boyd and Vandenberghe, Chapter 5 Solution 5.1 Note that
More informationExercise Sheet 1.
Exercise Sheet 1 You can download my lecture and exercise sheets at the address http://sami.hust.edu.vn/giang-vien/?name=huynt 1) Let A, B be sets. What does the statement "A is not a subset of B " mean?
More informationIn class midterm Exam - Answer key
Fall 2013 In class midterm Exam - Answer key ARE211 Problem 1 (20 points). Metrics: Let B be the set of all sequences x = (x 1,x 2,...). Define d(x,y) = sup{ x i y i : i = 1,2,...}. a) Prove that d is
More informationReconstruction of sparse Legendre and Gegenbauer expansions
Reconstruction of sparse Legendre and Gegenbauer expansions Daniel Potts Manfred Tasche We present a new deterministic algorithm for the reconstruction of sparse Legendre expansions from a small number
More informationToeplitz matrices. Niranjan U N. May 12, NITK, Surathkal. Definition Toeplitz theory Computational aspects References
Toeplitz matrices Niranjan U N NITK, Surathkal May 12, 2010 Niranjan U N (NITK, Surathkal) Linear Algebra May 12, 2010 1 / 15 1 Definition Toeplitz matrix Circulant matrix 2 Toeplitz theory Boundedness
More informationCirculant Matrices. Ashley Lorenz
irculant Matrices Ashley Lorenz Abstract irculant matrices are a special type of Toeplitz matrix and have unique properties. irculant matrices are applicable to many areas of math and science, such as
More informationHomework sheet 4: EIGENVALUES AND EIGENVECTORS. DIAGONALIZATION (with solutions) Year ? Why or why not? 6 9
Bachelor in Statistics and Business Universidad Carlos III de Madrid Mathematical Methods II María Barbero Liñán Homework sheet 4: EIGENVALUES AND EIGENVECTORS DIAGONALIZATION (with solutions) Year - Is
More informationFinding normalized and modularity cuts by spectral clustering. Ljubjana 2010, October
Finding normalized and modularity cuts by spectral clustering Marianna Bolla Institute of Mathematics Budapest University of Technology and Economics marib@math.bme.hu Ljubjana 2010, October Outline Find
More informationMATH 117 LECTURE NOTES
MATH 117 LECTURE NOTES XIN ZHOU Abstract. This is the set of lecture notes for Math 117 during Fall quarter of 2017 at UC Santa Barbara. The lectures follow closely the textbook [1]. Contents 1. The set
More informationCHAPTER 10 Shape Preserving Properties of B-splines
CHAPTER 10 Shape Preserving Properties of B-splines In earlier chapters we have seen a number of examples of the close relationship between a spline function and its B-spline coefficients This is especially
More informationLecture4 INF-MAT : 5. Fast Direct Solution of Large Linear Systems
Lecture4 INF-MAT 4350 2010: 5. Fast Direct Solution of Large Linear Systems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo September 16, 2010 Test Matrix
More informationNEW FUNCTIONAL INEQUALITIES
1 / 29 NEW FUNCTIONAL INEQUALITIES VIA STEIN S METHOD Giovanni Peccati (Luxembourg University) IMA, Minneapolis: April 28, 2015 2 / 29 INTRODUCTION Based on two joint works: (1) Nourdin, Peccati and Swan
More informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationOn John type ellipsoids
On John type ellipsoids B. Klartag Tel Aviv University Abstract Given an arbitrary convex symmetric body K R n, we construct a natural and non-trivial continuous map u K which associates ellipsoids to
More informationPaul Schrimpf. October 18, UBC Economics 526. Unconstrained optimization. Paul Schrimpf. Notation and definitions. First order conditions
Unconstrained UBC Economics 526 October 18, 2013 .1.2.3.4.5 Section 1 Unconstrained problem x U R n F : U R. max F (x) x U Definition F = max x U F (x) is the maximum of F on U if F (x) F for all x U and
More informationMATH 5640: Functions of Diagonalizable Matrices
MATH 5640: Functions of Diagonalizable Matrices Hung Phan, UMass Lowell November 27, 208 Spectral theorem for diagonalizable matrices Definition Let V = X Y Every v V is uniquely decomposed as u = x +
More information