Sparse and Fast Fourier Transforms in Biomedical Imaging
|
|
- Geoffrey Lane
- 5 years ago
- Views:
Transcription
1 Sparse and Fast Fourier Transforms in Biomedical Imaging Stefan Kunis 1 1 joint work with Ines Melzer, TU Chemnitz and Helmholtz Center Munich, supported by DFG and Helmholtz association
2 Outline 1 Trigonometric polynomials and FFTs 2 Specific sparse FFTs 3 Butterfly approximation scheme 4 Numerical experiments 5 Applications & Summary
3 Trigonometric polynomials and FFTs torus T d = R d /Z d = [ 1 2, 1 2) d, index set trigonometric polynomials [ I N = Z d N 2, N ) d 2 f (x) = k I N ˆf k e 2πikx discrete Fourier transform (DFT) f = Fˆf, f j = k I N ˆf k e 2πikj/N, j I N FFT (Gauß; Cooley, Tukey; Frigo, Johnson) O ( N d log N )
4 Trigonometric polynomials and FFTs motivation, randomly drawn X I N, sensing matrix F X IN = (e 2πikj/N ) j X,k IN Fourier coefficient support T I N (Candes, Tao; Rauhut; Needell, Vershynin) recovery scheme C T log 4 N X N d 1 identify support T I N - hard part (Gilbert, Tropp; Iwen) 2 solve for the nonzero Fourier coefficients F X Tˆf T = f X well conditioned F X T F X T matrix vector multiplications with F X T in less than O ( X T )?
5 Specific sparse FFTs sparse DFT, T, X I N, S = T = X N d f j = k T ˆf k e 2πikj/N, j X O ( S 2) divide and conquer, compute on nonzeros - pruning ˆf 0 ˆf 1 ˆf 2 ˆf 3 ˆf 4 ˆf 3 ˆf 2 ˆf 1 -i w 8 -i -i w 3 8 f 0 f 4 f 2 f 2 f 1 f 3 f 3 f 1 ( ) O N(1 + log S2 N ) for S N
6 Specific sparse FFTs T HCFFT X HCFFT, S = T = X = O(N log d 1 N) (Baszenski, Delvos; Hallatschek) O(S log N) for arbitrary X T d, S = T = X (Döhler, Fenn, Kämmerer, Potts, K.) O(S log d 1 N log ε d )
7 Specific sparse FFTs general sparsity pattern T I N, S = T = m d, m N x j = j/m, j I m, coarse grid O(S log S) x j = Gj, G R d d, coarse lattice O(S(log S + log ε d )) reduction to O ( S log α N ) ε only for very specific situations
8 Butterfly approximation scheme Lexing Ying. Sparse Fourier transform via butterfly algorithm, SIAM J. Sci. Comput. 31 (2009) smooth and sparse T [ N 2, N 2 )d, X [ N 2, N 2 )d T = X = O(N d 1 ) nodes in frequency/spatial domain lth dyadic subdivision of [ N 2, N 2 )d and [ 1 2, 1 2 )d into 2 dl boxes have only O(2 (d 1)l ) nonempty ones O(N d 1 log N p d+1 ) like in FMM, H-matrices, nonequispaced FFT, local expansion degree p = log ε, independent of N?
9 Butterfly approximation scheme considered model problem, nonsparse, univariate, given T = {k l [0, N] : l = 0,..., N} X = {x j [0, N] : j = 0,..., N} ˆf = (ˆf k ) k T C N evaluate almost periodic function for x X f (x) = k T ˆf k e 2πikx/N FFT for nonequispaced nodes in time and frequency domain (type-3 nufft, Greengard, Lee; nnfft, Keiner, Potts, K.)
10 Butterfly approximation scheme well known low rank property, p max(2eπ, log 2 ε ) e2πikx/n p 1 s=0 (2πi) s N s s! ks x s ε, kx N 2 admissible partitions of T X = [0, N] 2 diam(t )diam(x ) N dyadic decompositions of T and X, examples for N = 8 T T T T X X X X (Edelman; Überhuber,...)
11 dyadic decompositions of T and X, examples for N = 4 T 20 T 21 T 22 T 23 X 00 X 10 X 11 T 10 T 11 X 20 X 21 X 22 X 23 T 00 butterfly graph, nodes are admissible pairs X 00 T 20 X 00 T 21 X 00 T 22 X 00 T 23 X 10 T 10 X 10 T 11 X 11 T 10 X 11 T 11 X 20 T 00 X 21 T 00 X 22 T 00 X 23 T 00
12 Butterfly approximation scheme local in T, global in X : start with f T 30 (x) = ˆf k e 2πikx k T 30 T approximate f T 30 on X 00 f X00T30 approximate f X 00T 30 + f X 00T 31 on X go on local in X, global in T : finally f X 10T 20 f X 30T 00 is an approximation to f T 00 = f on X 30 (..., Rokhlin, O Neill,...)
13 Butterfly approximation scheme approximations local in T, global in X f X 00 T 20 f X 00 T 21 f X 00 T 22 f X 00 T 23 f X 10 T 10 f X 10 T 11 f X 11 T 10 f X 11 T 11 f X 20 T 00 f X 21 T 00 f X 22 T 00 f X 23 T 00 local in X, global in T
14 Butterfly approximation scheme frequency band T = [k min, k max ], admissible X = [x min, x max ] almost periodic function g Π T, T T, g : X C, g(x) = k T ĝ k e 2πikx/N T p = k min + t p 1 (k max k min ), t = 0,...p 1 x s = x min+x max 2 + xmax x min interpolation operator Jp XT 2 cos 2s+1 2p π, s = 0,...p 1 : Π T Π Tp J XT p g(x s ) = g(x s ), s = 0,..., p 1
15 Butterfly approximation scheme Theorem (Melzer, K. 2010) For admissible T, X, g Π T, and p 3, we have with C p = g J XT p g C(X ) C p g C(X ) 2π p 4 p 1 p! π p c 0 c p 1. Sketch of proof. shift invariance Jp XT e 2πikx = e 2πik0x J X (T k 0) p e 2πi(k k 0)x polynomial interpolation of Jp XT g and g coincide Bernstein inequality for g Π T
16 Butterfly approximation scheme Corollary (Melzer, K. 2010) For N = 2 L, T, X = [0, N] d, f Π T, ε > 0, and ( ) 1 2c0 d(l + 1) p = log 3 log c 1 ε we have f f C(X ) ε ˆf 1. Sketch of proof. J XT p = d r=1 J X (r) T (r) p J X (r) T (r) p 1 + C p
17 Numerical experiments T, X [0, N], T = X = N f (x) = ˆf k e 2πikx/N k T computation time, naive O(N 2 ), 1d-butterfly nnfft O(N log N ( log ε + log log N) 2 )
18 Numerical experiments 1d-butterfly nnfft, accuracy p log ε + log log N ε = c N c p 1 ε = c p log N
19 Numerical experiments smooth and sparse T, X [0, N] d, T = X = N d 1 f (x) = ˆf k e 2πikx/N k T computation time naive O(N 2d 2 ) butterfly O(N d 1 log N ( log ε + log log N) d+1 ) d = 2 d = 3
20 Applications & Summary 2d spherical mean values, forward simulation N 2 image data N acoustic sensors, N times/radii spherical mean value operator M : C(R d ) C(R d+1 ) (Mf ) (y, r) = 1 f (y + rξ)dσ(ξ) ω d 1 S d 1 let e k (x) = e 2πikx then (Me k )(y, r) = 2 d 2 2 Γ ( ) J d 2 (2π k r ) d 2 2 (2π k r ) d 2 2 e 2πiky
21 Applications & Summary special case d = 3 J 1 (2π k r) = 2 2 sin(2π k r) π 2π k r 3d photoacoustic imaging, forward simulation N 3 volume data N 2 acoustic sensors, N times/radii number of floating point operations naive O(N 5 ) 4d butterfly sparse FFT O(N 3 log N ( log ε + log log N) 5 )
22 Applications & Summary Thank you!
23 Stability local approximation scheme computes expansion coefficients of Jp XT in a monomial basis matrix G C p p, g s,t = e 2πiktxs, s, t = 0,..., p 1 x s = 1 2 (2s + 1)π cos, k t = 1 2p 2 + t p 1 condition number cond(g) ( p 1 p π ) p 1 d = 2, N = 2 10, limits to p 9 and ε 10 8
Generated sets as sampling schemes for hyperbolic cross trigonometric polynomials
Generated sets as sampling schemes for hyperbolic cross trigonometric polynomials Lutz Kämmerer Chemnitz University of Technology joint work with S. Kunis and D. Potts supported by DFG-SPP 1324 1 / 16
More informationIntroduction to Practical FFT and NFFT
Introduction to Practical FFT and NFFT Michael Pippig and Daniel Potts Department of Mathematics Chemnitz University of Technology September 14, 211 supported by BMBF grant 1IH81B Table of Contents 1 Serial
More informationInterpolation lattices for hyperbolic cross trigonometric polynomials
Interpolation lattices for hyperbolic cross trigonometric polynomials Lutz Kämmerer a, Stefan Kunis b, Daniel Potts a a Chemnitz University of Technology, Faculty of Mathematics, 09107 Chemnitz, Germany
More informationIntroduction to Practical FFT and NFFT
Introduction to Practical FFT and NFFT Michael Pippig and Daniel Potts Faculty of Mathematics Chemnitz University of Technology 07.09.2010 supported by BMBF grant 01IH08001B Table of Contents 1 Serial
More informationSparse Fourier Transform via Butterfly Algorithm
Sparse Fourier Transform via Butterfly Algorithm Lexing Ying Department of Mathematics, University of Texas, Austin, TX 78712 December 2007, revised June 2008 Abstract This paper introduces a fast algorithm
More informationGreedy Signal Recovery and Uniform Uncertainty Principles
Greedy Signal Recovery and Uniform Uncertainty Principles SPIE - IE 2008 Deanna Needell Joint work with Roman Vershynin UC Davis, January 2008 Greedy Signal Recovery and Uniform Uncertainty Principles
More informationSparse Legendre expansions via l 1 minimization
Sparse Legendre expansions via l 1 minimization Rachel Ward, Courant Institute, NYU Joint work with Holger Rauhut, Hausdorff Center for Mathematics, Bonn, Germany. June 8, 2010 Outline Sparse recovery
More informationFast Algorithms for the Computation of Oscillatory Integrals
Fast Algorithms for the Computation of Oscillatory Integrals Emmanuel Candès California Institute of Technology EPSRC Symposium Capstone Conference Warwick Mathematics Institute, July 2009 Collaborators
More informationA Non-sparse Tutorial on Sparse FFTs
A Non-sparse Tutorial on Sparse FFTs Mark Iwen Michigan State University April 8, 214 M.A. Iwen (MSU) Fast Sparse FFTs April 8, 214 1 / 34 Problem Setup Recover f : [, 2π] C consisting of k trigonometric
More informationProblem descriptions: Non-uniform fast Fourier transform
Problem descriptions: Non-uniform fast Fourier transform Lukas Exl 1 1 University of Vienna, Institut für Mathematik, Oskar-Morgenstern-Platz 1, A-1090 Wien, Austria December 11, 2015 Contents 0.1 Software
More informationNumerical Minimization of Potential Energies on Specific Manifolds
Numerical Minimization of Potential Energies on Specific Manifolds SIAM Conference on Applied Linear Algebra 2012 22 June 2012, Valencia Manuel Gra f 1 1 Chemnitz University of Technology, Germany, supported
More informationRapidly Computing Sparse Chebyshev and Legendre Coefficient Expansions via SFTs
Rapidly Computing Sparse Chebyshev and Legendre Coefficient Expansions via SFTs Mark Iwen Michigan State University October 18, 2014 Work with Janice (Xianfeng) Hu Graduating in May! M.A. Iwen (MSU) Fast
More informationPARAMETER ESTIMATION FOR MULTIVARIATE EXPONENTIAL SUMS
PARAMETER ESTIMATION FOR MULTIVARIATE EXPONENTIAL SUMS DANIEL POTTS AND MANFRED TASCHE Abstract. The recovery of signal parameters from noisy sampled data is an essential problem in digital signal processing.
More informationNumerical stability of nonequispaced fast Fourier transforms
Numerical stability of nonequispaced fast Fourier transforms Daniel Potts and Manfred Tasche Dedicated to Franz Locher in honor of his 65th birthday This paper presents some new results on numerical stability
More informationCoSaMP: Greedy Signal Recovery and Uniform Uncertainty Principles
CoSaMP: Greedy Signal Recovery and Uniform Uncertainty Principles SIAM Student Research Conference Deanna Needell Joint work with Roman Vershynin and Joel Tropp UC Davis, May 2008 CoSaMP: Greedy Signal
More informationSparse Fourier Transform (lecture 1)
1 / 73 Sparse Fourier Transform (lecture 1) Michael Kapralov 1 1 IBM Watson EPFL St. Petersburg CS Club November 2015 2 / 73 Given x C n, compute the Discrete Fourier Transform (DFT) of x: x i = 1 x n
More informationAn Adaptive Sublinear Time Block Sparse Fourier Transform
An Adaptive Sublinear Time Block Sparse Fourier Transform Volkan Cevher Michael Kapralov Jonathan Scarlett Amir Zandieh EPFL February 8th 217 Given x C N, compute the Discrete Fourier Transform (DFT) of
More informationCoSaMP. Iterative signal recovery from incomplete and inaccurate samples. Joel A. Tropp
CoSaMP Iterative signal recovery from incomplete and inaccurate samples Joel A. Tropp Applied & Computational Mathematics California Institute of Technology jtropp@acm.caltech.edu Joint with D. Needell
More informationUsing NFFT 3 a software library for various nonequispaced fast Fourier transforms
Using NFFT 3 a software library for various nonequispaced fast Fourier transforms Jens Keiner University of Lübeck 23560 Lübeck, Germany keiner@math.uni-luebeck.de and Stefan Kunis Chemnitz University
More informationParticle Simulation Based on Nonequispaced Fast Fourier Transforms
Particle Simulation Based on Nonequispaced Fast Fourier Transforms Michael Pippig and Daniel Potts Chemnitz University of Technology Department of Mathematics 0907 Chemnitz, Germany E-mail: {michael.pippig,
More informationSuper-resolution meets algebraic geometry
Super-resolution meets algebraic geometry Stefan Kunis joint work with H. Michael Möller, Thomas Peter, Tim Römer, and Ulrich von der Ohe Outline Univariate trigonometric moment problem Multivariate trigonometric
More informationIntro to harmonic analysis on groups Risi Kondor
Risi Kondor Any (sufficiently smooth) function f on the unit circle (equivalently, any 2π periodic f ) can be decomposed into a sum of sinusoidal waves f(x) = k= c n e ikx c n = 1 2π f(x) e ikx dx 2π 0
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationRobust Uncertainty Principles: Exact Signal Reconstruction from Highly Incomplete Frequency Information
1 Robust Uncertainty Principles: Exact Signal Reconstruction from Highly Incomplete Frequency Information Emmanuel Candès, California Institute of Technology International Conference on Computational Harmonic
More informationChapter 1 Divide and Conquer Polynomial Multiplication Algorithm Theory WS 2015/16 Fabian Kuhn
Chapter 1 Divide and Conquer Polynomial Multiplication Algorithm Theory WS 2015/16 Fabian Kuhn Formulation of the D&C principle Divide-and-conquer method for solving a problem instance of size n: 1. Divide
More informationGREEDY SIGNAL RECOVERY REVIEW
GREEDY SIGNAL RECOVERY REVIEW DEANNA NEEDELL, JOEL A. TROPP, ROMAN VERSHYNIN Abstract. The two major approaches to sparse recovery are L 1-minimization and greedy methods. Recently, Needell and Vershynin
More informationSize properties of wavelet packets generated using finite filters
Rev. Mat. Iberoamericana, 18 (2002, 249 265 Size properties of wavelet packets generated using finite filters Morten Nielsen Abstract We show that asymptotic estimates for the growth in L p (R- norm of
More informationFast Multipole Methods for The Laplace Equation. Outline
Fast Multipole Methods for The Laplace Equation Ramani Duraiswami Nail Gumerov Outline 3D Laplace equation and Coulomb potentials Multipole and local expansions Special functions Legendre polynomials Associated
More informationMultiple Rank-1 Lattices as Sampling Schemes for Multivariate Trigonometric Polynomials
Multiple Rank-1 Lattices as Sampling Schemes for Multivariate Trigonometric Polynomials Lutz Kämmerer June 23, 2015 We present a new sampling metho that allows the unique reconstruction of (sparse) multivariate
More informationA Remark on the Fast Gauss Transform
Publ. RIMS, Kyoto Univ. 39 (2003), 785 796 A Remark on the Fast Gauss Transform By Kenta Kobayashi Abstract We propose an improvement on the Fast Gauss Transform which was presented by Greengard and Sun
More informationFast algorithms based on asymptotic expansions of special functions
Fast algorithms based on asymptotic expansions of special functions Alex Townsend MIT Cornell University, 23rd January 2015, CAM Colloquium Introduction Two trends in the computing era Tabulations Software
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationThe spectral transform method
The spectral transform method by Nils Wedi European Centre for Medium-Range Weather Forecasts wedi@ecmwf.int Advanced Numerical Methods for Earth-System Modelling Slide 1 Advanced Numerical Methods for
More informationFast Multipole Methods
An Introduction to Fast Multipole Methods Ramani Duraiswami Institute for Advanced Computer Studies University of Maryland, College Park http://www.umiacs.umd.edu/~ramani Joint work with Nail A. Gumerov
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 informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
More informationEMPIRICAL EVALUATION OF A SUB-LINEAR TIME SPARSE DFT ALGORITHM
EMPIRICAL EVALUATIO OF A SUB-LIEAR TIME SPARSE DFT ALGORITHM M. A. IWE, A. GILBERT, AD M. STRAUSS Abstract. In this paper we empirically evaluate a recently proposed Fast Approximate Discrete Fourier Transform
More informationSample Optimal Fourier Sampling in Any Constant Dimension
1 / 26 Sample Optimal Fourier Sampling in Any Constant Dimension Piotr Indyk Michael Kapralov MIT IBM Watson October 21, 2014 Fourier Transform and Sparsity Discrete Fourier Transform Given x C n, compute
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 informationChapter 1 Divide and Conquer Algorithm Theory WS 2016/17 Fabian Kuhn
Chapter 1 Divide and Conquer Algorithm Theory WS 2016/17 Fabian Kuhn Formulation of the D&C principle Divide-and-conquer method for solving a problem instance of size n: 1. Divide n c: Solve the problem
More informationThe discrete and fast Fourier transforms
The discrete and fast Fourier transforms Marcel Oliver Revised April 7, 1 1 Fourier series We begin by recalling the familiar definition of the Fourier series. For a periodic function u: [, π] C, we define
More informationRecovery of Compactly Supported Functions from Spectrogram Measurements via Lifting
Recovery of Compactly Supported Functions from Spectrogram Measurements via Lifting Mark Iwen markiwen@math.msu.edu 2017 Friday, July 7 th, 2017 Joint work with... Sami Merhi (Michigan State University)
More informationA Fast Algorithm for Spherical Grid Rotations and its Application to Singular Quadrature
A Fast Algorithm for Spherical Grid Rotations and its Application to Singular Quadrature Zydrunas Gimbutas Shravan Veerapaneni September 11, 2013 Abstract We present a fast and accurate algorithm for evaluating
More informationComputational Methods CMSC/AMSC/MAPL 460. Fourier transform
Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Ramani Duraiswami, Dept. of Computer Science Several slides from Prof. Healy s course at UMD Fourier Methods Fourier analysis ( harmonic analysis
More informationThe Seven Dwarfs of Symbolic Computation and the Discovery of Reduced Symbolic Models. Erich Kaltofen North Carolina State University google->kaltofen
The Seven Dwarfs of Symbolic Computation and the Discovery of Reduced Symbolic Models Erich Kaltofen North Carolina State University google->kaltofen Outline 2 The 7 Dwarfs of Symbolic Computation Discovery
More informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Matheon Workshop Compressive Sensing and Its Applications TU Berlin, December 11,
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 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 informationFourier extension and sampling on the sphere
Fourier extension and sampling on the sphere Daniel Potts Technische Universität Chemnitz Faculty of Mathematics 97 Chemnitz, Germany Email: potts@mathematiktu-chemnitzde iel Van Buggenhout KU Leuven Faculty
More information5.6 Convolution and FFT
5.6 Convolution and FFT Fast Fourier Transform: Applications Applications. Optics, acoustics, quantum physics, telecommunications, control systems, signal processing, speech recognition, data compression,
More informationIntroduction to the FFT
Introduction to the FFT 1 Introduction Assume that we have a signal fx where x denotes time. We would like to represent fx as a linear combination of functions e πiax or, equivalently, sinπax and cosπax
More informationA Framework for Discrete Integral Transformations I the Pseudo-polar Fourier transform
A Framework for Discrete Integral Transformations I the Pseudo-polar Fourier transform A. Averbuch, R.R. Coifman, D.L. Donoho, M. Israeli, and Y. Shkolnisky June 26, 2007 This paper is dedicated to the
More informationComputational Methods CMSC/AMSC/MAPL 460
Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Balaji Vasan Srinivasan Dept of Computer Science Several slides from Prof Healy s course at UMD Last time: Fourier analysis F(t) = A 0 /2 + A
More informationSparse recovery for spherical harmonic expansions
Rachel Ward 1 1 Courant Institute, New York University Workshop Sparsity and Cosmology, Nice May 31, 2011 Cosmic Microwave Background Radiation (CMB) map Temperature is measured as T (θ, ϕ) = k k=0 l=
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationAn algebraic perspective on integer sparse recovery
An algebraic perspective on integer sparse recovery Lenny Fukshansky Claremont McKenna College (joint work with Deanna Needell and Benny Sudakov) Combinatorics Seminar USC October 31, 2018 From Wikipedia:
More informationRobust Principal Component Analysis
ELE 538B: Mathematics of High-Dimensional Data Robust Principal Component Analysis Yuxin Chen Princeton University, Fall 2018 Disentangling sparse and low-rank matrices Suppose we are given a matrix M
More informationRecovering overcomplete sparse representations from structured sensing
Recovering overcomplete sparse representations from structured sensing Deanna Needell Claremont McKenna College Feb. 2015 Support: Alfred P. Sloan Foundation and NSF CAREER #1348721. Joint work with Felix
More informationSuper-resolution via Convex Programming
Super-resolution via Convex Programming Carlos Fernandez-Granda (Joint work with Emmanuel Candès) Structure and Randomness in System Identication and Learning, IPAM 1/17/2013 1/17/2013 1 / 44 Index 1 Motivation
More informationInterpolation on the unit circle
Lecture 2: The Fast Discrete Fourier Transform Interpolation on the unit circle So far, we have considered interpolation of functions at nodes on the real line. When we model periodic phenomena, it is
More informationA Unified Framework for Oscillatory Integral Transform: When to use NUFFT or Butterfly Factorization?
A Unified Framework for Oscillatory Integral Transform: When to use NUFFT or Butterfly Factorization? Haizhao Yang Department of Mathematics, National University of Singapore March 12, 2018 Abstract This
More informationProny s method on the sphere
Prony s method on the sphere Stefan Kunis H. Michael Möller Ulrich von der Ohe Eigenvalue analysis based methods are well suited for the reconstruction of finitely supported measures from their moments
More informationDimensionality reduction: Johnson-Lindenstrauss lemma for structured random matrices
Dimensionality reduction: Johnson-Lindenstrauss lemma for structured random matrices Jan Vybíral Austrian Academy of Sciences RICAM, Linz, Austria January 2011 MPI Leipzig, Germany joint work with Aicke
More informationResearch Statement. James Bremer Department of Mathematics, University of California, Davis
Research Statement James Bremer Department of Mathematics, University of California, Davis Email: bremer@math.ucdavis.edu Webpage: https.math.ucdavis.edu/ bremer I work in the field of numerical analysis,
More informationIMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM
IMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM Anne Gelb, Rodrigo B. Platte, Rosie Renaut Development and Analysis of Non-Classical Numerical Approximation
More informationDirect Methods for Reconstruction of Functions and their Edges from Non-Uniform Fourier Data
Direct Methods for Reconstruction of Functions and their Edges from Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu ICERM Research Cluster Computational Challenges in Sparse and Redundant
More informationLimitations in Approximating RIP
Alok Puranik Mentor: Adrian Vladu Fifth Annual PRIMES Conference, 2015 Outline 1 Background The Problem Motivation Construction Certification 2 Planted model Planting eigenvalues Analysis Distinguishing
More informationA Characterization of Locally Testable Affine-Invariant Properties via Decomposition Theorems
A Characterization of Locally Testable Affine-Invariant Properties via Decomposition Theorems Yuichi Yoshida National Institute of Informatics and Preferred Infrastructure, Inc June 1, 2014 Yuichi Yoshida
More informationUniform Uncertainty Principle and signal recovery via Regularized Orthogonal Matching Pursuit
Uniform Uncertainty Principle and signal recovery via Regularized Orthogonal Matching Pursuit arxiv:0707.4203v2 [math.na] 14 Aug 2007 Deanna Needell Department of Mathematics University of California,
More informationA multivariate generalization of Prony s method
A multivariate generalization of Prony s method Ulrich von der Ohe joint with Stefan Kunis Thomas Peter Tim Römer Institute of Mathematics Osnabrück University SPARS 2015 Signal Processing with Adaptive
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationInverse problems and sparse models (1/2) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France
Inverse problems and sparse models (1/2) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France remi.gribonval@inria.fr Structure of the tutorial Session 1: Introduction to inverse problems & sparse
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 16 Oct. 26, 2017
Sketching Algorithms for Big Data Fall 2017 Prof. Piotr Indyk Lecture 16 Oct. 26, 2017 Scribe: Chi-Ning Chou 1 Overview In the last lecture we constructed sparse RIP 1 matrix via expander and showed that
More informationSparse Fourier Transform (lecture 4)
1 / 5 Sparse Fourier Transform (lecture 4) Michael Kapralov 1 1 IBM Watson EPFL St. Petersburg CS Club November 215 2 / 5 Given x C n, compute the Discrete Fourier Transform of x: x i = 1 n x j ω ij, j
More informationParameter estimation for nonincreasing exponential sums by Prony like methods
Parameter estimation for nonincreasing exponential sums by Prony like methods Daniel Potts Manfred Tasche Let z j := e f j with f j (, 0] + i [ π, π) be distinct nodes for j = 1,..., M. With complex coefficients
More informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Mathematical Analysis and Applications Workshop in honor of Rupert Lasser Helmholtz
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 7, No. 6, pp. 1903 198 c 006 Society for Industrial and Applied Mathematics FAST ALGORITHMS FOR SPHERICAL HARMONIC EXPANSIONS VLADIMIR ROKHLIN AND MARK TYGERT Abstract. An algorithm
More informationCS on CS: Computer Science insights into Compresive Sensing (and vice versa) Piotr Indyk MIT
CS on CS: Computer Science insights into Compresive Sensing (and vice versa) Piotr Indyk MIT Sparse Approximations Goal: approximate a highdimensional vector x by x that is sparse, i.e., has few nonzero
More informationNumerical Programming I (for CSE)
Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let
More informationModel-Based Compressive Sensing for Signal Ensembles. Marco F. Duarte Volkan Cevher Richard G. Baraniuk
Model-Based Compressive Sensing for Signal Ensembles Marco F. Duarte Volkan Cevher Richard G. Baraniuk Concise Signal Structure Sparse signal: only K out of N coordinates nonzero model: union of K-dimensional
More informationLecture 1: Introduction
CBMS Conference on Fast Direct Solvers Dartmouth College June June 7, 4 Lecture : Introduction Gunnar Martinsson The University of Colorado at Boulder Research support by: Many thanks to everyone who made
More informationBoolean Product Polynomials and the Resonance Arrangement
Boolean Product Polynomials and the Resonance Arrangement Sara Billey University of Washington Based on joint work with: Lou Billera and Vasu Tewari FPSAC July 17, 2018 Outline Symmetric Polynomials Schur
More informationFast Ewald Summation based on NFFT with Mixed Periodicity
Fast Ewald Summation based on NFFT with Mixed Periodicity Franziska Nestler, Michael Pippig and Daniel Potts In this paper we develop new fast Fourier-based methods for the Coulomb problem. We combine
More informationRandom Coding for Fast Forward Modeling
Random Coding for Fast Forward Modeling Justin Romberg with William Mantzel, Salman Asif, Karim Sabra, Ramesh Neelamani Georgia Tech, School of ECE Workshop on Sparsity and Computation June 11, 2010 Bonn,
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationRandomized Kaczmarz Nick Freris EPFL
Randomized Kaczmarz Nick Freris EPFL (Joint work with A. Zouzias) Outline Randomized Kaczmarz algorithm for linear systems Consistent (noiseless) Inconsistent (noisy) Optimal de-noising Convergence analysis
More informationA Non-Linear Lower Bound for Constant Depth Arithmetical Circuits via the Discrete Uncertainty Principle
A Non-Linear Lower Bound for Constant Depth Arithmetical Circuits via the Discrete Uncertainty Principle Maurice J. Jansen Kenneth W.Regan November 10, 2006 Abstract We prove a non-linear lower bound on
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 informationCSE 548: Analysis of Algorithms. Lecture 4 ( Divide-and-Conquer Algorithms: Polynomial Multiplication )
CSE 548: Analysis of Algorithms Lecture 4 ( Divide-and-Conquer Algorithms: Polynomial Multiplication ) Rezaul A. Chowdhury Department of Computer Science SUNY Stony Brook Spring 2015 Coefficient Representation
More informationFast Multipole Methods: Fundamentals & Applications. Ramani Duraiswami Nail A. Gumerov
Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov Week 1. Introduction. What are multipole methods and what is this course about. Problems from physics, mathematics,
More informationThe Nonequispaced FFT and its Applications A Mini Course at Helsinki University of Technology
The Nonequispaced FFT and its Applications A Mini Course at Helsinki University of Technology Jens Keiner 1 Stefan Kunis 2 Antje Vollrath 1 1 Institut für Mathematik, Universität zu Lübeck 2 Fakultät für
More informationSignal Recovery from Permuted Observations
EE381V Course Project Signal Recovery from Permuted Observations 1 Problem Shanshan Wu (sw33323) May 8th, 2015 We start with the following problem: let s R n be an unknown n-dimensional real-valued signal,
More informationThe uniform uncertainty principle and compressed sensing Harmonic analysis and related topics, Seville December 5, 2008
The uniform uncertainty principle and compressed sensing Harmonic analysis and related topics, Seville December 5, 2008 Emmanuel Candés (Caltech), Terence Tao (UCLA) 1 Uncertainty principles A basic principle
More informationDesign and Analysis of Algorithms
Design and Analysis of Algorithms CSE 5311 Lecture 5 Divide and Conquer: Fast Fourier Transform Junzhou Huang, Ph.D. Department of Computer Science and Engineering CSE5311 Design and Analysis of Algorithms
More informationLecture 20: Discrete Fourier Transform and FFT
EE518 Digital Signal Processing University of Washington Autumn 2001 Dept of Electrical Engineering Lecture 20: Discrete Fourier Transform and FFT Dec 10, 2001 Prof: J Bilmes TA:
More informationFast Algorithms for Spherical Harmonic Expansions
Fast Algorithms for Spherical Harmonic Expansions Vladimir Rokhlin and Mark Tygert Research Report YALEU/DCS/RR-1309 December 17, 004 Abstract An algorithm is introduced for the rapid evaluation at appropriately
More informationHilbert Space Methods for Reduced-Rank Gaussian Process Regression
Hilbert Space Methods for Reduced-Rank Gaussian Process Regression Arno Solin and Simo Särkkä Aalto University, Finland Workshop on Gaussian Process Approximation Copenhagen, Denmark, May 2015 Solin &
More informationAn Introduction to Sparse Approximation
An Introduction to Sparse Approximation Anna C. Gilbert Department of Mathematics University of Michigan Basic image/signal/data compression: transform coding Approximate signals sparsely Compress images,
More informationSparse FFT for Functions with Short Frequency Support
Sparse FFT for Functions with Short Frequency Support Sina ittens May 10, 017 Abstract In this paper we derive a deterministic fast discrete Fourier transform algorithm for a π-periodic function f whose
More information