
In this paper we study a class of fast geometric image inpainting methods
based on the idea of filling the inpainting domain in successive shells from
its boundary inwards. Image pixels are filled by assigning them a color equal
to a weighted average of their already filled neighbors. However, there is
flexibility in terms of the order in which pixels are filled, the weights used
for averaging, and the neighborhood that is averaged over. Varying these
degrees of freedom leads to different algorithms, and indeed the literature
contains several methods falling into this general class. All of them are very
fast, but at the same time all of them leave undesirable artifacts such as
"kinking" (bending) or blurring of extrapolated isophotes. Our objective in
this paper is to build a theoretical model, based on a continuum limit and a
connection to stopped random walks, in order to understand why these artifacts
occur and what, if anything, can be done about them. At the same time, we
consider a semiimplicit extension in which pixels in a given shell are solved
for simultaneously by solving a linear system. We prove (within the continuum
limit) that this extension is able to completely eliminate kinking artifacts,
which we also prove must always be present in the direct method. Although our
analysis makes the strong assumption of a square inpainting domain, it makes
weak smoothness assumptions and is thus applicable to the low regularity
inherent in images.

We consider the problem of convergence to a saddle point of a concaveconvex
function via gradient dynamics. Since first introduced by Arrow, Hurwicz and
Uzawa in [1] such dynamics have been extensively used in diverse areas, there
are, however, features that render their analysis non trivial. These include
the lack of convergence guarantees when the function considered is not strictly
concaveconvex and also the nonsmoothness of subgradient dynamics. Our aim in
this two part paper is to provide an explicit characterization to the
asymptotic behaviour of general gradient and subgradient dynamics applied to a
general concaveconvex function. We show that despite the nonlinearity and
nonsmoothness of these dynamics their $\omega$limit set is comprised of
trajectories that solve only explicit linear ODEs that are characterized within
the paper.
More precisely, in Part I an exact characterization is provided to the
asymptotic behaviour of unconstrained gradient dynamics. We also show that when
convergence to a saddle point is not guaranteed then the system behaviour can
be problematic, with arbitrarily small noise leading to an unbounded variance.
In Part II we consider a general class of subgradient dynamics that restrict
trajectories in an arbitrary convex domain, and show that when an equilibrium
point exists their limiting trajectories are solutions of subgradient dynamics
on only affine subspaces. The latter is a smooth class of dynamics with an
asymptotic behaviour exactly characterized in Part I, as solutions to explicit
linear ODEs. These results are used to formulate corresponding convergence
criteria and are demonstrated with several examples and applications presented
in Part II.

In part I we considered the problem of convergence to a saddle point of a
concaveconvex function via gradient dynamics and an exact characterization was
given to their asymptotic behaviour. In part II we consider a general class of
subgradient dynamics that provide a restriction in an arbitrary convex domain.
We show that despite the nonlinear and nonsmooth character of these dynamics
their $\omega$limit set is comprised of solutions to only linear ODEs. In
particular, we show that the latter are solutions to subgradient dynamics on
affine subspaces which is a smooth class of dynamics the asymptotic properties
of which have been exactly characterized in part I. Various convergence
criteria are formulated using these results and several examples and
applications are also discussed throughout the manuscript.

In this paper, we establish uniqueness of the solution of the VlasovPoisson
system with spatial density belonging to a certain class of Orlicz spaces. This
extends the uniqueness result of Loeper (which holds for uniformly bounded
density) and the uniqueness result of the second author. Uniqueness is a direct
consequence of our main result, which provides a quantitative stability
estimate for the Wasserstein distance between two weak solutions with spatial
density in such Orlicz spaces, in the spirit of Dobrushin's proof of stability
for meanfield PDEs. Our proofs are built on the secondorder structure of the
underlying characteristic system associated to the equation.

We develop a new technique for establishing quantitative propagation of chaos
for systems of interacting particles. Using this technique we prove propagation
of chaos for diffusing particles whose interaction kernel is merely H\"older
continuous, even at long ranges. Moreover, we do not require specially prepared
initial data. On the way, we establish a law of large numbers for SDEs that
holds over a class of vector fields simultaneously. The proofs bring together
ideas from empirical process theory and stochastic flows.

We develop a technique of multiple scale asymptotic expansions along mean
flows and a corresponding notion of weak multiple scale convergence. These are
applied to homogenize convection dominated parabolic equations with rapidly
oscillating, locally periodic coefficients and $\mathcal{O}(\eps^{1})$ mean
convection term. Crucial to our analysis is the introduction of a fast time
variable, $\tau=\frac{t}{\eps}$, not apparent in the heterogeneous problem. The
effective diffusion coefficient is expressed in terms of the average of
Eulerian cell solutions along the orbits of the mean flow in the fast time
variable. To make this notion rigorous, we use the theory of ergodic algebras
with mean value.

The problem of approximating the discrete spectra of families of selfadjoint
operators that are merely strongly continuous is addressed. It is wellknown
that the spectrum need not vary continuously (as a set) under strong
perturbations. However, it is shown that under an additional compactness
assumption the spectrum does vary continuously, and a family of symmetric
finitedimensional approximations is constructed. An important feature of these
approximations is that they are valid for the entire family uniformly. An
application of this result to the study of plasma instabilities is illustrated.

We study convergence to equilibrium for the kinetic FokkerPlanck equation on
the torus. Solving the stochastic differential equation, we show exponential
convergence in the MongeKantorovichWasserstein $\mathcal{W}_2$ distance.
Finally, we investigate if such a coupling can be obtained by a coadapted
coupling, and show that then the bound must depend on the square root of the
initial distance.

The relativistic VlasovMaxwell system describes the evolution of a
collisionless plasma. The problem of linear instability of this system is
considered in two physical settings: the socalled "one and onehalf"
dimensional case, and the three dimensional case with cylindrical symmetry.
Sufficient conditions for instability are obtained in terms of the spectral
properties of certain Schr\"odinger operators that act on the spatial variable
alone (and not in full phase space). An important aspect of these conditions is
that they do not require any boundedness assumptions on the domains, nor do
they require monotonicity of the equilibrium.