Source author record

Adrian Nachman

Adrian Nachman 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

9works
7topics
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

9 published item(s)

preprint2022arXiv

Wasserstein GANs with Gradient Penalty Compute Congested Transport

Wasserstein GANs with Gradient Penalty (WGAN-GP) are a very popular method for training generative models to produce high quality synthetic data. While WGAN-GP were initially developed to calculate the Wasserstein 1 distance between generated and real data, recent works (e.g. [23]) have provided empirical evidence that this does not occur, and have argued that WGAN-GP perform well not in spite of this issue, but because of it. In this paper we show for the first time that WGAN-GP compute the minimum of a different optimal transport problem, the so-called congested transport [7]. Congested transport determines the cost of moving one distribution to another under a transport model that penalizes congestion. For WGAN-GP, we find that the congestion penalty has a spatially varying component determined by the sampling strategy used in [12] which acts like a local speed limit, making congestion cost less in some regions than others. This aspect of the congested transport problem is new, in that the congestion penalty turns out to be unbounded and depends on the distributions to be transported, and so we provide the necessary mathematical proofs for this setting. One facet of our discovery is a formula connecting the gradient of solutions to the optimization problem in WGAN-GP to the time averaged momentum of the optimal mass flow. This is in contrast to the gradient of Kantorovich potentials for the Wasserstein 1 distance, which is just the normalized direction of flow. Based on this and other considerations, we speculate on how our results explain the observed performance of WGAN-GP. Beyond applications to GANs, our theorems also point to the possibility of approximately solving large scale congested transport problems using neural network techniques.

preprint2020arXiv

A Multiscale Theory for Image Registration and Nonlinear Inverse Problems

In an influential paper, Tadmor, Nezzar and Vese (Multiscale Model. Simul. (2004)) introduced a hierarchical decomposition of an image as a sum of constituents of different scales. Here we construct analogous hierarchical expansions for diffeomorphisms, in the context of image registration, with the sum replaced by composition of maps. We treat this as a special case of a general framework for multiscale decompositions, applicable to a wide range of imaging and nonlinear inverse problems. As a paradigmatic example of the latter, we consider the Calderón inverse conductivity problem. We prove that we can simultaneously perform a numerical reconstruction and a multiscale decomposition of the unknown conductivity, driven by the inverse problem itself. We provide novel convergence proofs which work in the general abstract settings, yet are sharp enough to settle an open problem on the hierarchical decompostion of Tadmor, Nezzar and Vese for arbitrary functions in $L^2$. We also give counterexamples that show the optimality of our general results.

preprint2016arXiv

A weighted minimum gradient problem with complete electrode model boundary conditions for conductivity imaging

We consider the inverse problem of recovering an isotropic electrical conductivity from interior knowledge of the magnitude of one current density field generated by applying current on a set of electrodes. The required interior data can be obtained by means of MRI measurements. On the boundary we only require knowledge of the electrodes, their impedances, and the corresponding average input currents. From the mathematical point of view, this practical question leads us to consider a new weighted minimum gradient problem for functions satisfying the boundary conditions coming from the Complete Electrode Model of Somersalo, Cheney and Isaacson. This variational problem has non-unique solutions. The surprising discovery is that the physical data is still sufficient to determine the geometry of the level sets of the minimizers. In particular, we obtain an interesting phase retrieval result: knowledge of the input current at the boundary allows determination of the full current vector field from its magnitude. We characterize the non-uniqueness in the variational problem. We also show that additional measurements of the voltage potential along one curve joining the electrodes yield unique determination of the conductivity. A nonlinear algorithm is proposed and implemented to illustrate the theoretical results.

preprint2014arXiv

Uniqueness of minimizers of weighted least gradient problems arising in conductivity imaging

We prove uniqueness for minimizers of the weighted least gradient problem \[\inf \left\lbrace \int_Ω a|Du|: \ \ u\in BV(Ω), \ \ u|_{\partial Ω}=f \right\rbrace.\] The weight function $a$ is assumed to be continuous and it is allowed to vanish in certain subsets of $Ω$. Existence is assumed a priori. Our approach is motivated by the hybrid inverse problem of imaging electric conductivity from interior knowledge (obtainable by MRI) of the magnitude of one current density vector field.

preprint2013arXiv

Current Density Impedance Imaging of an Anisotropic Conductivity in a Known Conformal Class

We present a procedure for recovering the conformal factor of an anisotropic conductivity matrix in a known conformal class in a domain in Euclidean space of dimension greater than or equal to 2. The method requires one internal measurement, together with a priori knowledge of the conformal class (local orientation) of the conductivity matrix. This problem arises in the coupled-physics medical imaging modality of Current Density Impedance Imaging (CDII) and the assumptions on the data are suitable for measurements determinable from cross-property based couplings of the two imaging modalities CDII and Diffusion Tensor Imaging (DTI). We show that the corresponding electric potential is the unique solution of a constrained minimization problem with respect to a weighted total variation functional defined in terms of the physical data. Further, we show that the associated equipotential surfaces are area minimizing with respect to a Riemannian metric obtained from the data. The results are also extended to allow the presence of perfectly conducting and/or insulating inclusions.

preprint2011arXiv

A convergent algorithm for the hybrid problem of reconstructing conductivity from minimal interior data

We consider the hybrid problem of reconstructing the isotropic electric conductivity of a body $Ω$ from interior Current Density Imaging data obtainable using MRI measurements. We only require knowledge of the magnitude $|J|$ of one current generated by a given voltage $f$ on the boundary $\partialΩ$. As previously shown, the corresponding voltage potential u in $Ω$ is a minimizer of the weighted least gradient problem \[u=\hbox{argmin} \{\int_Ωa(x)|\nabla u|: u \in H^{1}(Ω), \ \ u|_{\partial Ω}=f\},\] with $a(x)= |J(x)|$. In this paper we present an alternating split Bregman algorithm for treating such least gradient problems, for $a\in L^2(Ω)$ non-negative and $f\in H^{1/2}(\partial Ω)$. We give a detailed convergence proof by focusing to a large extent on the dual problem. This leads naturally to the alternating split Bregman algorithm. The dual problem also turns out to yield a novel method to recover the full vector field $J$ from knowledge of its magnitude, and of the voltage $f$ on the boundary. We then present several numerical experiments that illustrate the convergence behavior of the proposed algorithm.

preprint2011arXiv

Conductivity imaging from one interior measurement in the presence of perfectly conducting and insulating inclusions

We consider the problem of recovering an isotropic conductivity outside some perfectly conducting or insulating inclusions from the interior measurement of the magnitude of one current density field $|J|$. We prove that the conductivity outside the inclusions, and the shape and position of the perfectly conducting and insulating inclusions are uniquely determined (except in an exceptional case) by the magnitude of the current generated by imposing a given boundary voltage. We have found an extension of the notion of admissibility to the case of possible presence of perfectly conducting and insulating inclusions. This also makes it possible to extend the results on uniqueness of the minimizers of the least gradient problem $F(u)=\int_Ωa|\nabla u|$ with $u|_{\partial Ω}=f$ to cases where $u$ has flat regions (is constant on open sets).

preprint2011arXiv

Convergence of the alternating split Bregman algorithm in infinite-dimensional Hilbert spaces

We prove results on weak convergence for the alternating split Bregman algorithm in infinite dimensional Hilbert spaces. We also show convergence of an approximate split Bregman algorithm, where errors are allowed at each step of the computation. To be able to treat the infinite dimensional case, our proofs focus mostly on the dual problem. We rely on Svaiter's theorem on weak convergence of the Douglas-Rachford splitting algorithm and on the relation between the alternating split Bregman and Douglas-Rachford splitting algorithms discovered by Setzer. Our motivation for this study is to provide a convergent algorithm for weighted least gradient problems arising in the hybrid method of imaging electric conductivity from interior knowledge (obtainable by MRI) of the magnitude of one current.

preprint2011arXiv

On the steady states of weakly reversible chemical reaction networks

A natural condition on the structure of the underlying chemical reaction network, namely weak reversibility, is shown to guarantee the existence of an equilibrium (steady state) in each positive stoichiometric compatibility class for the associated mass-action system. Furthermore, an index formula is given for the set of equilibria in a given stoichiometric compatibility class.