Publications

2017

• Understanding the Time-Dependent Effective Diffusion Coefficient Measured by Diffusion MRI: the Intra-Cellular Case
  
  • Haddar Houssem
  • Li Jing-Rebecca
  • Schiavi Simona
  SIAM Journal on Applied Mathematics, Society for Industrial and Applied Mathematics, 2017. Diffusion Magnetic Resonance Imaging (dMRI) can be used to measure a time-dependent effective diffusion coefficient that can in turn reveal information about the tissue geometry. Recently a mathematical model for the time-dependent effective diffusion coefficient was obtained using homogenization techniques after imposing a certain scaling relationship for the time, the biological cell membrane permeability, the diffusion-encoding magnetic field gradient strength, and a periodicity length of the cellular geometry. With this choice of the scaling of the physical parameters, the effective diffusion coefficient of the medium can be computed after solving a diffusion equation subject to a time-dependent Neumann boundary condition, independently in the biological cells and in the extra-cellular space. In this paper, we analyze this new model, which we call the H-ADC model, in the case of finite domains, which is relevant to diffusion inside biological cells. We use both the eigenfunction expansion and the single layer potential representation for the solution of the above mentioned diffusion equation to obtain analytical expressions for the effective diffusion coefficient in different diffusion time regimes. These expressions are validated using numerical simulations in two dimensions.
• Parameter Estimation in Nonlinear Mixed Effect Models Using saemix, an R Implementation of the SAEM Algorithm
  
  • Comets Emmanuelle
  • Lavenu Audrey Paris
  • Lavielle Marc
  Journal of Statistical Software, University of California, Los Angeles, 2017, 80 (3), pp.i03. The saemix package for R provides maximum likelihood estimates of parameters in nonlinear mixed effect models, using a modern and efficient estimation algorithm, the stochastic approximation expectation-maximisation (SAEM) algorithm. In the present paper we describe the main features of the package, and apply it to several examples to illustrate its use. Making use of S4 classes and methods to provide user-friendly interaction, this package provides a new estimation tool to the R community. (10.18637/jss.v080.i03)
  DOI : 10.18637/jss.v080.i03
• On perturbed proximal gradient algorithms
  
  • Atchadé Yves
  • Fort Gersende
  • Moulines Éric
  Journal of Machine Learning Research, Microtome Publishing, 2017, 18 (10), pp.1-33.
• The infinitesimal model: definition, derivation, and implications
  
  • Barton Nick
  • Etheridge Alison M
  • Véber Amandine
  Theoretical Population Biology, Elsevier, 2017. Our focus here is on the infinitesimal model. In this model, one or several quantitative traits are described as the sum of a genetic and a non-genetic component, the first being distributed within families as a normal random variable centred at the average of the parental genetic components, and with a variance independent of the parental traits. Thus, the variance that segregates within families is not perturbed by selection, and can be predicted from the variance components. This does not necessarily imply that the trait distribution across the whole population should be Gaussian, and indeed selection or population structure may have a substantial effect on the overall trait distribution. One of our main aims is to identify some general conditions on the allelic effects for the infinitesimal model to be accurate. We first review the long history of the infinitesimal model in quantitative genetics. Then we formulate the model at the phenotypic level in terms of individual trait values and relationships between individuals, but including different evolutionary processes: genetic drift, recombination, selection, mutation , population structure, ... We give a range of examples of its application to evolutionary questions related to stabilising selection, assortative mating, effective population size and response to selection, habitat preference and speciation. We provide a mathematical justification of the model as the limit as the number M of underlying loci tends to infinity of a model with Mendelian inheritance, mutation and environmental noise, when the genetic component of the trait is purely additive. We also show how the model generalises to include epistatic effects. We prove in particular that, within each family, the genetic components of the individual trait values in the current generation are indeed normally distributed with a variance independent of ancestral traits, up to an error of order 1/\sqrt{M}. Simulations suggest that in some cases the convergence may be as fast as 1/\sqrt{M} .
• Adaptive multipreconditioned FETI: scalability results and robustness assessment
  
  • Bovet Christophe
  • Parret-Fréaud Augustin
  • Spillane Nicole
  • Gosselet Pierre
  Computers & Structures, Elsevier, 2017, 193, pp.1-20. The purpose of this article is to assess the adaptive multipreconditioned FETI solvers (AMPFETI) on realistic industrial problems and hardware. The multi-preconditioned FETI algorithm (first introduced as Simultaneous FETI [1]) is a non-overlapping domain decomposition method which exhibits good robust-ness properties without requiring the explicit knowledge of the original partial differential equation, or any a priori analysis of the algebraic system through eigenvalues problems. Multipreconditioned FETI solves critical problems in significantly fewer iterations than classical FETI but each iteration involves a larger computational effort. An adaptive strategy (known as the adaptive mul-tipreconditioned conjugate gradient algorithm [2]) has been proposed to achieve balance between robustness and efficiency and we will observe that it provides an efficient solver for the problems considered here. (10.1016/j.compstruc.2017.07.010)
  DOI : 10.1016/j.compstruc.2017.07.010
• EFSA-GMO-DE-2017-141 1
  
  • Angevin Frédérique
  • Bagnis Claude
  • Bar-Hen Avner
  • Barny Marie-Anne
  • Boireau Pascal
  • Brévault Thierry
  • Chauvel Bruno B.
  • Collonnier Cécile
  • Couvet Denis
  • Dassa Elie
  • de Verneuil Hubert
  • Demeneix Barbara
  • Franche Claudine
  • Guerche Philippe
  • Guillemain Joël
  • Hernandez Raquet Guillermina
  • Khalife Jamal
  • Klonjkowski Bernard
  • Lavielle Marc
  • Le Corre Valérie
  • Lefèvre François
  • Lemaire Olivier
  • Lereclus Didier D.
  • Maximilien Rémy
  • Meurs Eliane
  • Naffakh Nadia
  • Négre Didier
  • Noyer Jean-Louis
  • Ochatt Sergio
  • Pages Jean-Christophe
  • Raynaud Xavier
  • Regnault-Roger Catherine
  • Renard Michel M.
  • Renault Tristan
  • Saindrenan Patrick
  • Simonet Pascal
  • Troadec Marie-Bérengère
  • Vaissière Bernard
  • Vilotte Jean-Luc
  , 2017.
• Self-administered vitamin D status predictor : older adults are able to use a self-questionnaire for evaluating their vitamin D status
  
  • Annweiler Cédric
  • Kabeshova Anastasiia
  • Callens Alix
  • Paty Marie-Liesse
  • Holick Michael F.
  • Duval Guillaume T.
  PLoS ONE, Public Library of Science, 2017, 12 (11). The 16-item Vitamin D Status Predictor (VDSP) questionnaire helps to identify, without resorting to a blood test, older adults with low vitamin D concentrations. Our objective was to determine whether a self-administered VDSP was concordant with the VDSP administered by a physician, and to examine the concordance of every single item of the VDSP. Methods A total of 349 older in- and outpatients (mean, 83.2±7.2years; 59% female) were consecutively recruited in the geriatric ward of the University Hospital of Angers, France. All participants completed a self-administered VDSP questionnaire (self-VDSP) in paper format composed of 17 items exploring age, gender, general condition, nutrition, vision, mood, cognition, gait and falls, and osteoporosis. All participants underwent subsequently a full clinical examination by a physician exploring the same areas (rater-VDSP). Results The agreement between the self-VDSP and the rater-VDSP was almost perfect for the probability of having low vitamin D concentrations, regardless of the definition used (i.e., ≤25, ≤50 or ≤75 nmol/L). The agreements between physicians’ and patients’ responses were significant for every single VDSP item. The agreement was fair to perfect for all items, except for cognitive disorders, undernutrition and polymorbidity (poor agreement). Conclusions Older adults are able to evaluate their own probabilities of severe vitamin D deficiency, deficiency and insufficiency. A self-questionnaire may promote the use of the VDSP tool in this population, and help clinicians in decisions to supplement their patients in a reasoned way. (10.1371/journal.pone.0186578)
  DOI : 10.1371/journal.pone.0186578
• Moral hazard in dynamic risk management
  
  • Cvitanić Jakša
  • Possamaï Dylan
  • Touzi Nizar
  Management Science, INFORMS, 2017, 63 (10), pp.3328-3346. We consider a contracting problem in which a principal hires an agent to manage a risky project. When the agent chooses volatility components of the output process and the principal observes the output continuously, the principal can compute the quadratic variation of the output, but not the individual components. This leads to moral hazard with respect to the risk choices of the agent. We identify a family of admissible contracts for which the optimal agent's action is explicitly characterized, and, using the recent theory of singular changes of measures for Itô processes, we study how restrictive this family is. In particular, in the special case of the standard Homlstrom-Milgrom model with fixed volatility, the family includes all possible contracts. We solve the principal-agent problem in the case of CARA preferences, and show that the optimal contract is linear in these factors: the contractible sources of risk, including the output, the quadratic variation of the output and the cross-variations between the output and the contractible risk sources. Thus, like sample Sharpe ratios used in practice, path-dependent contracts naturally arise when there is moral hazard with respect to risk management. In a numerical example, we show that the loss of efficiency can be significant if the principal does not use the quadratic variation component of the optimal contract. (10.1287/mnsc.2016.2493)
  DOI : 10.1287/mnsc.2016.2493
• Suboptimal feedback control of PDEs by solving HJB equations on adaptive sparse grids
  
  • Garcke Jochen
  • Kröner Axel
  Journal of Scientific Computing, Springer Verlag, 2017, 70 (1), pp.1-28. An approach to solve finite time horizon suboptimal feedback control problems for partial differential equations is proposed by solving dynamic programming equations on adaptive sparse grids. The approach is illustrated for the wave equation and an extension to equations of Schrödinger type is indicated. A semi-discrete optimal control problem is introduced and the feedback control is derived from the corresponding value function. The value function can be characterized as the solution of an evolutionary Hamilton-Jacobi Bellman (HJB) equation which is defined over a state space whose dimension is equal to the dimension of the underlying semi-discrete system. Besides a low dimensional semi-discretization it is important to solve the HJB equation efficiently to address the curse of dimensionality. We propose to apply a semi-Lagrangian scheme using spatially adaptive sparse grids. Sparse grids allow the discretization of the value functions in (higher) space dimensions since the curse of dimensionality of full grid methods arises to a much smaller extent. For additional efficiency an adaptive grid refinement procedure is explored. We present several numerical examples studying the effect the parameters characterizing the sparse grid have on the accuracy of the value function and the optimal trajectory.
• Instability of dielectrics and conductors in electrostatic fields
  
  • Allaire Grégoire
  • Rauch Jeffrey
  Archive for Rational Mechanics and Analysis, Springer Verlag, 2017, 124, pp.233-268. This article proves most of the assertion in §116 of Maxwell's treatise on electromagnetism. The results go under the name Earnshaw's Theorem and assert the absence of stable equilibrium configurations of conductors and dielectrics in an external electrostatic field.
• Online Learning and Blackwell Approachability with Partial Monitoring: Optimal Convergence Rates
  
  • Kwon Joon
  • Perchet Vianney
  JMLR Papers, 2017, 54, pp.604-613. Blackwell approachability is an online learning setup generalizing the classical problem of regret minimization by allowing for instance multi-criteria optimization, global (online) optimization of a convex loss, or online linear optimization under some cumulative constraint. We consider partial monitoring where the decision maker does not necessarily observe the outcomes of his decision (unlike the traditional regret/bandit literature). Instead, he receives a random signal correlated to the decision-outcome pair, or only to the outcome. We construct, for the first time, approachability algorithms with convergence rate of order O(T −1/2) when the signal is independent of the decision and of order O(T −1/3) in the case of general signals. Those rates are optimal in the sense that they cannot be improved without further assumption on the structure of the objectives and/or the signals.
• Équations de Navier-Stokes incompressibles et multirésolution spatiale adaptative: sur la question des modes parasites en maillage collocalisé.
  
  • N'Guessan Marc-Arthur
  • Massot Marc
  • Tenaud Christian
  • Series Laurent
  , 2017. La simulation numérique directe (DNS) de la combustion avec chimie détaillée et transport multi-espèces représente l'un des défis les plus importants en matière de calcul scientifique dans de nombreuses applications industrielles. Un des enjeux est de coupler un solveur hydrodynamique pour la résolution des équations de Navier-Stokes, pour un mélange réactif dans la limite des faibles nombres de Mach, à une stratégie de résolution de systèmes de convection-réaction-diffusion, tout en maintenant l'efficacité algorithmique, l'adaptation temps-espace et le contrôle d'erreur. La présente communication vise à proposer une stratégie optimale pour l'élimination des modes parasites dans un contexte de maillages colocalisés en volume fini, dans un cadre de multirésolution spatiale.
• Generalized and hybrid Metropolis-Hastings overdamped Langevin algorithms
  
  • Poncet Romain
  , 2017. It has been shown that the nonreversible overdamped Langevin dynamics enjoy better convergence properties in terms of spectral gap and asymptotic variance than the reversible one. In this article we propose a variance reduction method for the Metropolis-Hastings Adjusted Langevin Algorithm (MALA) that makes use of the good behaviour of the these nonreversible dynamics. It consists in constructing a nonreversible Markov chain (with respect to the target invariant measure) by using a Generalized Metropolis-Hastings adjustment on a lifted state space. We present two variations of this method and we discuss the importance of a well-chosen proposal distribution in terms of average rejection probability. We conclude with numerical experimentations to compare our algorithms with the MALA, and show variance reduction of several orders of magnitude in some favourable toy cases.
• Exploring the complexity of the integer image problem in the max-algebra
  
  • Maccaig Marie
  Discrete Applied Mathematics, Elsevier, 2017, 217 (2), pp.261--275. We investigate the complexity of the problem of finding an integer vector in the max-algebraic column span of a matrix, which we call the integer image problem. We show some cases where we can determine in strongly polynomial time whether such an integer vector exists, and find such an integer vector if it does exist. On the other hand we also describe a group of related problems each of which we prove to be NP-hard. Our main results demonstrate that the integer image problem is equivalent to finding a special type of integer image of a matrix satisfying a property we call column typical. For a subclass of matrices this problem is polynomially solvable but if we remove the column typical assumption then it becomes NP-hard. (10.1016/j.dam.2016.09.016)
  DOI : 10.1016/j.dam.2016.09.016
• Quantitative DLA-based compressed sensing for T1-weighted acquisitions.
  
  • Svehla Pavel
  • Nguyen Khieu-Van
  • Li Jing-Rebecca
  • Ciobanu Luisa
  Journal of Magnetic Resonance, Elsevier, 2017. (10.1016/j.jmr.2017.05.002)
  DOI : 10.1016/j.jmr.2017.05.002
• Information-Geometric Optimization Algorithms: A Unifying Picture via Invariance Principles
  
  • Ollivier Yann
  • Arnold Ludovic
  • Auger Anne
  • Hansen Nikolaus
  Journal of Machine Learning Research, Microtome Publishing, 2017, 18 (18), pp.1-65. We present a canonical way to turn any smooth parametric family of probability distributions on an arbitrary search space X into a continuous-time black-box optimization method on X, the information-geometric optimization (IGO) method. Invariance as a major design principle keeps the number of arbitrary choices to a minimum. The resulting IGO flow is the flow of an ordinary differential equation conducting the natural gradient ascent of an adaptive, time-dependent transformation of the objective function. It makes no particular assumptions on the objective function to be optimized. The IGO method produces explicit IGO algorithms through time discretization. It naturally recovers versions of known algorithms and offers a systematic way to derive new ones. In continuous search spaces, IGO algorithms take a form related to natural evolution strategies (NES). The cross-entropy method is recovered in a particular case with a large time step, and can be extended into a smoothed, parametrization-independent maximum likelihood update (IGO-ML). When applied to the family of Gaussian distributions on R^d, the IGO framework recovers a version of the well-known CMA-ES algorithm and of xNES. For the family of Bernoulli distributions on {0, 1}^d, we recover the seminal PBIL algorithm and cGA. For the distributions of restricted Boltzmann machines, we naturally obtain a novel algorithm for discrete optimization on {0, 1}^d. All these algorithms are natural instances of, and unified under, the single information-geometric optimization framework. The IGO method achieves, thanks to its intrinsic formulation, maximal invariance properties: invariance under reparametrization of the search space X, under a change of parameters of the probability distribution, and under increasing transformation of the function to be optimized. The latter is achieved through an adaptive, quantile-based formulation of the objective. Theoretical considerations strongly suggest that IGO algorithms are essentially characterized by a minimal change of the distribution over time. Therefore they have minimal loss in diversity through the course of optimization, provided the initial diversity is high. First experiments using restricted Boltzmann machines confirm this insight. As a simple consequence, IGO seems to provide, from information theory, an elegant way to simultaneously explore several valleys of a fitness landscape in a single run.
• Resolution Analysis of Passive Synthetic Aperture Imaging of Fast Moving Objects
  
  • Garnier Josselin
  • Borcea L.
  • Papanicolaou G.
  • Solna K.
  • Tsogka C.
  SIAM Journal on Imaging Sciences, Society for Industrial and Applied Mathematics, 2017, 10 (2), pp.665 - 710. (10.1137/16M109716X)
  DOI : 10.1137/16M109716X
• Functional Characterization of Intrinsic and Extrinsic Geometry
  
  • Corman Etienne
  • Solomon Justin
  • Ben-Chen Mirela
  • Guibas Leonidas
  • Ovsjanikov Maks
  ACM Transactions on Graphics, Association for Computing Machinery, 2017, 17. We propose a novel way to capture and characterize distortion between pairs of shapes by extending the recently proposed framework of shape differences built on functional maps. We modify the original definition of shape differences slightly and prove that, after this change, the discrete metric is fully encoded in two shape difference operators and can be recovered by solving two linear systems of equations. Then, we introduce an extension of the shape difference operators using offset surfaces to capture extrinsic or embedding-dependent distortion, complementing the purely intrinsic nature of the original shape differences. Finally, we demonstrate that a set of four operators is complete, capturing intrinsic and extrinsic structure and fully encoding a shape up to rigid motion in both discrete and continuous settings. We highlight the usefulness of our constructions by showing the complementary nature of our extrinsic shape differences in capturing distortion ignored by previous approaches. We additionally provide examples where we recover local shape structure from the shape difference operators, suggesting shape editing and analysis tools based on manipulating shape differences.
• An equilibrated fluxes approach to the certified descent algorithm for shape optimization using conforming finite element and discontinuous Galerkin discretizations
  
  • Giacomini Matteo
  Journal of Scientific Computing, Springer Verlag, 2017. The certified descent algorithm (CDA) is a gradient-based method for shape optimization which certifies that the direction computed using the shape gradient is a genuine descent direction for the objective functional under analysis. It relies on the computation of an upper bound of the error introduced by the finite element approximation of the shape gradient. In this paper, we present a goal-oriented error estimator which depends solely on local quantities and is fully-computable. By means of the equilibrated fluxes approach, we construct a unified strategy valid for both conforming finite element approximations and discontinuous Galerkin discretizations. The new variant of the CDA is tested on the inverse identification problem of electrical impedance tomography: both its ability to identify a genuine descent direction at each iteration and its reliable stopping criterion are confirmed. (10.1007/s10915-017-0545-1)
  DOI : 10.1007/s10915-017-0545-1
• A characterization of switched linear control systems with finite L 2 -gain
  
  • Chitour Yacine
  • Mason Paolo
  • Sigalotti Mario
  IEEE Transactions on Automatic Control, Institute of Electrical and Electronics Engineers, 2017, 62, pp.1825-1837. Motivated by an open problem posed by J.P. Hespanha, we extend the notion of Barabanov norm and extremal trajectory to classes of switching signals that are not closed under concatenation. We use these tools to prove that the finiteness of the L2-gain is equivalent, for a large set of switched linear control systems, to the condition that the generalized spectral radius associated with any minimal realization of the original switched system is smaller than one. (10.1109/tac.2016.2593678)
  DOI : 10.1109/tac.2016.2593678
• Dependence of tropical eigenspaces
  
  • Niv Adi
  • Rowen Louis
  Communications in Algebra, Taylor & Francis, 2017, 45 (3), pp.924-942. We study the pathology that causes tropical eigenspaces of distinct su-pertropical eigenvalues of a non-singular matrix A, to be dependent. We show that in lower dimensions the eigenvectors of distinct eigenvalues are independent, as desired. The index set that differentiates between subsequent essential monomials of the characteristic polynomial, yields an eigenvalue λ, and corresponds to the columns of adj(A + λI) from which the eigenvectors are taken. We ascertain the cause for failure in higher dimensions, and prove that independence of the eigenvectors is recovered in case the " difference criterion " holds, defined in terms of disjoint differences between index sets of subsequent coefficients. We conclude by considering the eigenvectors of the matrix A ∇ := 1 det(A) adj(A) and the connection of the independence question to generalized eigenvectors. (10.1080/00927872.2016.1172603)
  DOI : 10.1080/00927872.2016.1172603
• Log-majorization of the moduli of the eigenvalues of a matrix polynomial by tropical roots
  
  • Akian Marianne
  • Gaubert Stéphane
  • Sharify Meisam
  Linear Algebra and its Applications, Elsevier, 2017, 528, pp.394--435. We show that the sequence of moduli of the eigenvalues of a matrix polynomial is log-majorized, up to universal constants, by a sequence of "tropical roots" depending only on the norms of the matrix coefficients. These tropical roots are the non-differentiability points of an auxiliary tropical polynomial, or equivalently, the opposites of the slopes of its Newton polygon. This extends to the case of matrix polynomials some bounds obtained by Hadamard, Ostrowski and Pólya for the roots of scalar polynomials. We also obtain new bounds in the scalar case, which are accurate for "fewnomials" or when the tropical roots are well separated. (10.1016/j.laa.2016.11.004)
  DOI : 10.1016/j.laa.2016.11.004
• Matched-Filter and Correlation-Based Imaging for Fast Moving Objects Using a Sparse Network of Receivers
  
  • Garnier Josselin
  • Fournier J.
  • Papanicolaou G.
  • Tsogka C.
  SIAM Journal on Imaging Sciences, Society for Industrial and Applied Mathematics, 2017, 10 (4), pp.2165 - 2216. (10.1137/17M112364X)
  DOI : 10.1137/17M112364X
• Certified Descent Algorithm for shape optimization driven by fully-computable a posteriori error estimators
  
  • Giacomini Matteo
  • Pantz Olivier
  • Trabelsi Karim
  ESAIM: Control, Optimisation and Calculus of Variations, EDP Sciences, 2017, 23 (3), pp.977-1001. In this paper we introduce a novel certified shape optimization strategy-named Certified Descent Algorithm (CDA)-to account for the numerical error introduced by the Finite Element approximation of the shape gradient. We present a goal-oriented procedure to derive a certified upper bound of the error in the shape gradient and we construct a fully-computable, constant-free a posteriori error estimator inspired by the complementary energy principle. The resulting CDA is able to identify a genuine descent direction at each iteration and features a reliable stopping criterion. After validating the error estimator, some numerical simulations of the resulting certified shape optimization strategy are presented for the well-known inverse identification problem of Electrical Impedance Tomography. (10.1051/cocv/2016021)
  DOI : 10.1051/cocv/2016021
• Pulse Reflection in a Random Waveguide with a Turning Point
  
  • Garnier Josselin
  • Borcea Liliana
  Multiscale Modeling and Simulation: A SIAM Interdisciplinary Journal, Society for Industrial and Applied Mathematics, 2017, 15 (4), pp.1472 - 1501. (10.1137/16M1094154)
  DOI : 10.1137/16M1094154