Georgia Institute of Technology College of Sciences

We consider an inverse problem for an inhomogeneous wave equation with discrete-in-time sources, modeling a seismic rupture. We assume that the sources occur along an unknown path with subsonic velocity, and that data is collected over time on some detection surface. We explore the question of uniqueness for these problems, and show how to recover the times and locations of sources microlocally first, and then the smooth part of the source assuming that it is the same at each source location. In case the sources (now all different) are (roughly speaking) non-negative and of limited oscillation in space, and sufficiently separated in space-time, which is a model for microseismicity, we present an explicit reconstruction, requiring sufficient local energy decay. (Joint research with L. Oksanen and J. Tittelfitz)

A wide variety of questions which range from social and economic sciences to physical and biological sciences lead to functions with values that are sets in finite or infinite dimensional spaces, or that are fuzzy sets. Set-valued and fuzzy-valued functions attract attention of a lot of researchers and allow them to look at numerous problems from a new point of view and provide them with new tools, ideas and results. In this talk we consider a generalized concept of such functions, that of functions with values in so-called L-space, that encompasses set-valued and fuzzy functions as special cases and allow to investigate them from the common point of view. We will discus several problems of Approximation Theory and Numerical Analysis for functions with values in L-spaces. In particular numerical methods of solution of Fredholm and Volterra integral equations for such functions will be presented.

I will discuss inverse problems involving elliptic partial differential equations with highly oscillating coefficients. The multiscale nature of such problems poses a challenge in both the mathematical formulation and the numerical modeling, which is hard even for forward computations. I will discuss uniqueness of the inverse in certain problem classes and give numerical methods for inversion that can be applied to problems in medical imaging and exploration seismology.

Given a finite set of matrices F, the Markovian Joint Spectral Radius represents the maximal rate of growth of products of matrices in F when the matrices are multiplied each other following some Markovian law. This quantity is important, for instance, in the study of the so called zero stability of variable stepsize BDF methods for the numerical integration of ordinary differential equations. Recently Kozyakin, based on a work by Dai, showed that, given a set F of N matrices of dimension d and a graph G, which represents the admissible products, it is possibile to compute the Markovian Joint Spectral Radius of the couple (F,G) as the classical Joint Spectral Radius of a new set of N matrices of dimension N*d, which are produced as a particular lifting of the matrices in F. Clearly by this approach the exact evaluation or the simple approximation of the Markovian Joint Spectral Radius becomes a challenge even for reasonably small values of N and d. In this talk we briefly review the theory of the Joint Spectral Radius, and we introduce the Markovian Joint Spectral Radius. Furthermore we address the question whether it is possible to reduce the exact calculation computational complexity of the Markovian Joint Spectral Radius. We show that the problem can be recast as the computation of N polytope norms in dimension d. We conclude the presentation with some numerical examples. This talk is based on a joint work with Nicola Guglielmi from the University of L'Aquila, Italy, and Vladimir Yu. Protasov from the Moscow State University, Russia.