Source author record

Jun Kitagawa

Jun Kitagawa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2021arXiv

Quantitative stability in the geometry of semi-discrete optimal transport

We show quantitative stability results for the geometric "cells" arising in semi-discrete optimal transport problems. Our results show two types of stability, the first is stability of the associated Laguerre cells in measure, without any connectedness or regularity assumptions on the source measure. The second is stability in Hausdorff measure, under a Poincar{è}-Wirtinger inequality and a regularity assumption equivalent to the Ma-Trudinger-Wang conditions of regularity in Monge-Amp{è}re. This last result also yields stability in the uniform norm of the dual potential functions, all three stability results come with explicit quantitative bounds. Our methods utilize a combination of graph theory, convex geometry, and Monge-Amp{è}re regularity theory.

preprint2020arXiv

$\mathcal{W}_\infty$-transport with discrete target as a combinatorial matching problem

In this short note, we show that given a cost function $c$, any coupling $π$ of two probability measures where the second is a discrete measure can be associated to a certain bipartite graph containing a perfect matching, based on the value of the infinity transport cost $\norm{c}_{L^\infty(π)}$. This correspondence between couplings and bipartite graphs is explicitly constructed. We give two applications of this result to the $\mathcal{W}_\infty$ optimal transport problem when the target measure is discrete, the first is a condition to ensure existence of an optimal plan induced by a mapping, and the second is a numerical approach to approximating optimal plans.

preprint2020arXiv

A Newton algorithm for semi-discrete optimal transport with storage fees

We introduce and prove convergence of a damped Newton algorithm to approximate solutions of the semi-discrete optimal transport problem with storage fees, corresponding to a problem with hard capacity constraints. This is a variant of the optimal transport problem arising in queue penalization problems, and has applications to data clustering. Our result is novel as it is the first numerical method with proven convergence for this variant problem; additionally the algorithm applies to the classical semi-discrete optimal transport problem but does not require any connectedness assumptions on the support of the source measure, in contrast with existing results. Furthermore we find some stability results of the associated Laguerre cells. All of our results come with quantitative rates. We also present some numerical examples.

preprint2020arXiv

Inverse Iteration for the Monge-Ampère Eigenvalue Problem

We present an iterative method based on repeatedly inverting the Monge-Ampère operator with Dirichlet boundary condition and prescribed right-hand side on a bounded, convex domain $Ω\subset \mathbb{R}^n$. We prove that the iterates $u_k$ generated by this method converge as $k \to \infty$ to a solution of the Monge-Ampère eigenvalue problem $$\begin{cases} \text{det} D^2u = λ_{MA} (-u)^n & \quad \text{in } Ω,\\ u = 0 & \quad \text{on } \partial Ω. \end{cases}$$ Since the solutions of this problem are unique up to a positive multiplicative constant, the normalized iterates $\hat{u}_k := \frac{u_k}{||u_k||_{L^{\infty}(Ω)}}$ converge to the eigenfunction of unit height. In addition, we show that $\lim\limits_{k \to \infty} R(u_k) = \lim\limits_{k \to \infty} R(\hat{u}_k) = λ_{MA}$, where the Rayleigh quotient $R(u)$ is defined as $$R(u) := \frac{\int_Ω (-u) \ \text{det} D^2u}{\int_Ω (-u)^{n+1}}.$$ Our method converges for a wide class of initial choices $u_0$ that can be constructed explicitly, and does not rely on prior knowledge of the Monge-Ampère eigenvalue $λ_{MA}$.

preprint2020arXiv

Optimal transport and the Gauss curvature equation

In this short note, we consider the problem of prescribing the Gauss curvature and image of the Gauss map for the graph of a function over a domain in Euclidean space. The prescription of the image of the Gauss map turns this into a second boundary value problem. Our main observation is that this problem can be posed as an optimal transport problem where the target is a subset of the lower hemisphere of $\mathbb{S}^n$. As a result we obtain existence and regularity of solutions under mild assumptions on the curvature, as well as a quantitative version of a gradient blowup result due to Urbas, which turns out to fall within the optimal transport framework.

preprint2017arXiv

Free discontinuities in optimal transport

We prove a nonsmooth implicit function theorem applicable to the zero set of the difference of convex functions. This theorem is explicit and global: it gives a formula representing this zero set as a difference of convex functions which holds throughout the entire domain of the original functions. As applications, we prove results on the stability of singularities of envelopes of semi-convex functions, and solutions to optimal transport problems under appropriate perturbations, along with global structure theorems on certain discontinuities arising in optimal transport maps for Ma-Trudinger-Wang costs. For targets whose components satisfy additional convexity, separation, multiplicity and affine independence assumptions we show these discontinuities occur on submanifolds of the appropriate codimension which are parameterized locally as differences of convex functions (DC, hence $C^2$ rectifiable), and --- depending on the precise assumptions --- $C^{1,α}$ smooth. In this case the highest codimension submanifolds consists of isolated points, each uniquely identified by the (affinely independent) components of the target to which it is transported.

preprint2016arXiv

Pointwise estimates and regularity in geometric optics and other Generated Jacobian Equations

The study of reflector surfaces in geometric optics necessitates the analysis of certain nonlinear equations of Monge-Ampère type known as generated Jacobian equations. These equations, whose general existence theory has been recently developed by Trudinger go beyond the framework of optimal transport. We obtain pointwise estimates for weak solutions of such equations under a condition analogous to the A3w condition of Ma, Trudinger and Wang. Estimates of this type have played an important role in the regularity theory for optimal transport maps and were previously unknown in this context, including the important case of the near field reflector problem.

preprint2015arXiv

The multi-marginal optimal partial transport problem

We introduce and study a multi-marginal optimal partial transport problem. Under a natural and sharp condition on the dominating marginals, we establish uniqueness of the optimal plan. Our strategy of proof establishes and exploits a connection with another novel problem, which we call the Monge-Kantorovich partial barycenter problem (with quadratic cost). This latter problem has a natural interpretation as a variant of the mines and factories description of optimal transport. We then turn our attention to various analytic properties of these two problems. Of particular interest, we show that monotonicity of the active marginals can fail, a surprising difference from the two marginal case.

preprint2014arXiv

Coronal magnetic field and the plasma beta determined from radio and multiple satellite observations

We derived the coronal magnetic field, plasma density, and temperature from the observation of polarization and intensity of radio thermal free-free emission using the Nobeyama Radioheliograph (NoRH) and extreme ultraviolet (EUV) observations. We observed a post-flare loop on the west limb 11 April 2013. The line-of-sight magnetic field was derived from the circularly polarized free-free emission observed by NoRH. The emission measure and temperature were derived from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The derived temperature was used to estimate the emission measure from the NoRH radio free-free emission observations. The derived density from NoRH was larger than that determined using AIA, which can be explained by the fact that the low temperature plasma is not within the temperature coverage of the AIA filters used in this study. We also discuss the other observation of the post-flare loops by the EUV Imager onboard the Solar Terrestrial Relations Observatory (STEREO), which can be used in future studies to reconstruct the coronal magnetic field strength. The derived plasma parameters and magnetic field were used to derive the plasma beta, which is a ratio between the magnetic pressure and the plasma pressure. The derived plasma beta is about 5.7*10^(-4) to 7.6*10^(-4) at the loop top region.

preprint2013arXiv

The local geometry of maps with c-convex potentials

We identify a condition for regularity of optimal transport maps that requires only three derivatives of the cost function, for measures given by densities that are only bounded above and below. This new condition is equivalent to the weak Ma-Trudinger-Wang condition when the cost is $C^4$. Moreover, we only require (non-strict) c-convexity of the support of the target measure, removing the hypothesis of strong c-convexity in a previous result of Figalli, Kim, and McCann, but at the added cost of assuming compact containment of the supports of both the source and target measures.

preprint2012arXiv

An iterative scheme for solving the optimal transportation problem

We demonstrate an iterative scheme to approximate the optimal transportation problem with a discrete target measure under certain standard conditions on the cost function. Additionally, we give a finite upper bound on the number of iterations necessary for the scheme to terminate, in terms of the error tolerance and number of points in the support of the discrete target measure.