Georgia Institute of Technology College of Sciences

The discrete prolate spheroidal sequences (DPSS's) provide an efficient representation for discrete signals that are perfectly timelimited and nearly bandlimited. Due to the high computational complexity of projecting onto the DPSS basis - also known as the Slepian basis - this representation is often overlooked in favor of the fast Fourier transform (FFT). In this talk I will describe novel fast algorithms for computing approximate projections onto the leading Slepian basis elements with a complexity comparable to the FFT. I will also highlight applications of this Fast Slepian Transform in the context of compressive sensing and processing of sampled multiband signals.

There are numerous modern applications in data science that involve inference from incomplete data. Various geometric prior models such as sparse vectors or low-rank matrices have been employed to address the ill-posed inverse problems arising in these applications. Recently, similar ideas were adopted to tackle more challenging nonlinear inverse problems such as phase retrieval and blind deconvolution. In this talk, we consider the blind deconvolution problem where the desired information as a time series is accessed as indirect observations through a time-invariant system with uncertainty. The measurements in this case is given in the form of the convolution with an unknown kernel. Particularly, we study the mathematical theory of multichannel blind deconvolution where we observe the output of multiple channels that are all excited with the same unknown input source. From these observations, we wish to estimate the source and the impulse responses of each of the channels simultaneously. We show that this problem is well-posed if the channel impulse responses follow a simple geometric model. Under these models, we show how the channel estimates can be found by solving corresponding non-convex optimization problems. We analyze methods for solving these non-convex programs, and provide performance guarantees for each.

In this talk, we will present new central and central DG schemes for solving ideal magnetohydrodynamic (MHD) equations while preserving globally divergence-free magnetic field on triangular grids. These schemes incorporate the constrained transport (CT) scheme of Evans and Hawley with central schemes and central DG methods on overlapping cells which have no need for solving Riemann problems across cell edges where there are discontinuities of the numerical solution. The schemes are formally second-order accurate with major development on the reconstruction of globally divergence-free magnetic field on polygonal dual mesh. Moreover, the computational cost is reduced by solving the complete set of governing equations on the primal grid while only solving the magnetic induction equation on the polygonal dual mesh.

Weighted direct samplers, sometimes also called importance samplers, are Monte Carlo algorithms for generating independent, weighted samples from a given target probability distribution. They are used in, e.g., data assimilation, state estimation for dynamical systems, and computational statistical mechanics. One challenge in designing weighted samplers is to ensure the variance of the weights, and that of the resulting estimator, are well-behaved. Recently, Chorin, Tu, Morzfeld, and coworkers have introduced a class of novel weighted samplers called implicit samplers, which possess a number of nice empirical properties. In this talk, I will summarize an asymptotic analysis of implicit samplers in the small-noise limit and describe a simple method to obtain a higher-order accuracy. I will also discuss extensions to stochastic differential equatons. This is joint work with Jonathan Goodman, Andrew Leach, and Matthias Morzfeld.