Publications

2017

• Galerkin approximations for the optimal control of nonlinear delay differential equations
  
  • Chekroun Mickaël D
  • Kröner Axel
  • Liu Honghu
  , 2017. Optimal control problems of nonlinear delay differential equations (DDEs) are considered for which we propose a general Galerkin approximation scheme built from Koornwinder polynomials. Error estimates for the resulting Galerkin-Koornwinder approximations to the optimal control and the value function, are derived for a broad class of cost function-als and nonlinear DDEs. The approach is illustrated on a delayed logistic equation set not far away from its Hopf bifurcation point in the parameter space. In this case, we show that low-dimensional controls for a standard quadratic cost functional can be efficiently computed from Galerkin-Koornwinder approximations to reduce at a nearly optimal cost the oscillation amplitude displayed by the DDE's solution. Optimal controls computed from the Pontryagin's maximum principle (PMP) and the Hamilton-Jacobi-Bellman equation (HJB) associated with the corresponding ODE systems, are shown to provide numerical solutions in good agreement. It is finally argued that the value function computed from the corresponding reduced HJB equation provides a good approximation of that obtained from the full HJB equation.
• On the heat diffusion for generic Riemannian and sub-Riemannian structures
  
  • Barilari Davide
  • Boscain Ugo
  • Charlot Grégoire
  • Neel Robert W.
  International Mathematics Research Notices, Oxford University Press (OUP), 2017, 15, pp.4639-4672. In this paper we provide the small-time heat kernel asymptotics at the cut locus in three relevant cases: generic low-dimensional Riemannian manifolds, generic 3D contact sub-Riemannian manifolds (close to the starting point) and generic 4D quasi-contact sub-Riemannian manifolds (close to a generic starting point). As a byproduct, we show that, for generic low-dimensional Riemannian manifolds, the only singularities of the exponential map, as a Lagragian map, that can arise along a minimizing geodesic are $A_3$ and $A_5$ (in the classification of Arnol'd's school). We show that in the non-generic case, a cornucopia of asymptotics can occur, even for Riemannian surfaces. (10.1093/imrn/rnw141)
  DOI : 10.1093/imrn/rnw141
• Some discussions on the Read Paper "Beyond subjective and objective in statistics" by A. Gelman and C. Hennig
  
  • Robert Christian
  • Celeux Gilles
  • Jewson Jack
  • Josse Julie
  • Marin Jean-Michel
  • Robert Christian P.
  , 2017. This note is a collection of several discussions of the paper "Beyond subjective and objective in statistics", read by A. Gelman and C. Hennig to the Royal Statistical Society on April 12, 2017, and to appear in the Journal of the Royal Statistical Society, Series A.
• Transport of power in random waveguides with turning points
  
  • Garnier Josselin
  • Borcea Liliana
  • Wood Derek
  Communications in Mathematical Sciences, International Press, 2017, 15 (8), pp.2327 - 2371. (10.4310/CMS.2017.v15.n8.a9)
  DOI : 10.4310/CMS.2017.v15.n8.a9
• Multi-Wave Medical Imaging: Mathematical Modelling and Imaging Reconstruction
  
  • Ammari Habib
  • Garnier Josselin
  • Kang Hyeonbae
  • Nguyen Loc Hoang
  • Seppecher Laurent
  , 2017.
• Pliability, or the whitney extension theorem for curves in carnot groups
  
  • Juillet Nicolas
  • Sigalotti Mario
  Analysis & PDE, Mathematical Sciences Publishers, 2017, 10 (7), pp.1637 - 1661. The Whitney extension theorem is a classical result in analysis giving a necessary and sufficient condition for a function defined on a closed set to be extendable to the whole space with a given class of regularity. It has been adapted to several settings, among which the one of Carnot groups. However, the target space has generally been assumed to be equal to R^d for some d ≥ 1. We focus here on the extendability problem for general ordered pairs (G_1,G_2) (with G_2 non-Abelian). We analyze in particular the case G_1 = R and characterize the groups G_2 for which the Whitney extension property holds, in terms of a newly introduced notion that we call pliability. Pliability happens to be related to rigidity as defined by Bryant an Hsu. We exploit this relation in order to provide examples of non-pliable Carnot groups, that is, Carnot groups so that the Whitney extension property does not hold. We use geometric control theory results on the accessibility of control affine systems in order to test the pliability of a Carnot group. In particular, we recover some recent results by Le Donne, Speight and Zimmermann about Lusin approximation in Carnot groups of step 2 and Whitney extension in Heisenberg groups. We extend such results to all pliable Carnot groups, and we show that the latter may be of arbitrarily large step. (10.2140/apde.2017.10.1637)
  DOI : 10.2140/apde.2017.10.1637
• Correction to Black--Scholes Formula Due to Fractional Stochastic Volatility
  
  • Garnier Josselin
  • Sølna Knut
  SIAM Journal on Financial Mathematics, Society for Industrial and Applied Mathematics, 2017, 8 (1), pp.560 - 588. (10.1137/15M1036749)
  DOI : 10.1137/15M1036749
• Consistent functional cross field design for mesh quadrangulation
  
  • Azencot Omri
  • Corman Etienne
  • Ben-Chen Mirela
  • Ovsjanikov Maks
  ACM Transactions on Graphics, Association for Computing Machinery, 2017, 36 (4), pp.92. We propose a novel technique for computing consistent cross fields on a pair of triangle meshes given an input correspondence, which we use as guiding fields for approximately consistent quadrangulations. Unlike the majority of existing methods our approach does not assume that the meshes share the same connectivity or even have the same number of vertices, and furthermore does not place any restrictions on the topology (genus) of the shapes. Importantly, our method is robust with respect to small perturbations of the given correspondence, as it only relies on the transportation of real-valued functions and thus avoids the costly and error-prone estimation of the map differential. Key to this robustness is a novel formulation, which relies on the previously-proposed notion of power vectors, and we show how consistency can be enforced without pre-alignment of local basis frames, in which these power vectors are computed. We demonstrate that using the same formulation we can both compute a quadrangulation that would respect a given symmetry on the same shape or a map across a pair of shapes. We provide quantitative and qualitative comparison of our method with several baselines and show that it both provides more accurate results and allows to handle more general cases than existing techniques. (10.1145/3072959.3073696)
  DOI : 10.1145/3072959.3073696
• Cell Averaging Two-Scale Convergence: Applications to Periodic Homogenization
  
  • Alouges François
  • Di Fratta Giovanni
  Multiscale Modeling and Simulation: A SIAM Interdisciplinary Journal, Society for Industrial and Applied Mathematics, 2017, 15 (4), pp.1651-1671. (10.1137/16M1085309)
  DOI : 10.1137/16M1085309
• 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.
• An approximation of the M 2 closure: application to radiotherapy dose simulation
  
  • Pichard T
  • Alldredge G W
  • Brull Stéphane
  • Dubroca B
  • Frank M
  Journal of Scientific Computing, Springer Verlag, 2017. Particle transport in radiation therapy can be modelled by a kinetic equation which must be solved numerically. Unfortunately, the numerical solution of such equations is generally too expensive for applications in medical centers. Moment methods provide a hierarchy of models used to reduce the numerical cost of these simulations while preserving basic properties of the solutions. Moment models require a closure because they have more unknowns than equations. The entropy-based closure is based on the physical description of the particle interactions and provides desirable properties. However, computing this closure is expensive. We propose an approximation of the closure for the first two models in the hierarchy, the M 1 and M 2 models valid in one, two or three dimensions of space. Compared to other approximate closures, our method works in multiple dimensions. We obtain the approximation by a careful study of the domain of realizability and by invariance properties of the entropy minimizer. The M 2 model is shown to provide significantly better accuracy than the M 1 model for the numerical simulation of a dose computation in radiotherapy. We propose a numerical solver using those approximated closures. Numerical experiments in dose computation test cases show that the new method is more efficient compared to numerical solution of the minimum entropy problem using standard software tools.
• HOMOGENIZATION OF STOKES SYSTEM USING BLOCH WAVES
  
  • Allaire Grégoire
  • Ghosh Tuhin
  • Vanninathan Muthusamy
  Networks and Heterogeneous Media, American Institute of Mathematical Sciences, 2017, 12 (4), pp.525-550. In this work, we study the Bloch wave homogenization for the Stokes system with periodic viscosity coefficient. In particular, we obtain the spectral interpretation of the homogenized tensor. The presence of the incompressibility constraint in the model raises new issues linking the homogenized tensor and the Bloch spectral data. The main difficulty is a lack of smoothness for the bottom of the Bloch spectrum, a phenomenon which is not present in the case of the elasticity system. This issue is solved in the present work, completing the homogenization process of the Stokes system via the Bloch wave method.
• A numerical approach to determine mutant invasion fitness and evolutionary singular strategies
  
  • Fritsch Coralie
  • Campillo Fabien
  • Ovaskainen Otso
  Theoretical Population Biology, Elsevier, 2017, 115, pp.89-99. We propose a numerical approach to study the invasion fitness of a mutant and to determine evolutionary singular strategies in evolutionary structured models in which the competitive exclusion principle holds. Our approach is based on a dual representation, which consists of the modelling of the small size mutant population by a stochastic model and the computation of its corresponding deterministic model. The use of the deterministic model greatly facilitates the numerical determination of the feasibility of invasion as well as the convergence-stability of the evolutionary singular strategy. Our approach combines standard adaptive dynamics with the link between the mutant survival criterion in the stochastic model and the sign of the eigenvalue in the corresponding deterministic model. We present our method in the context of a mass-structured individual-based chemostat model. We exploit a previously derived mathematical relationship between stochastic and deterministic representations of the mutant population in the chemostat model to derive a general numerical method for analyzing the invasion fitness in the stochastic models. Our method can be applied to the broad class of evolutionary models for which a link between the stochastic and deterministic invasion fitnesses can be established. (10.1016/j.tpb.2017.05.001)
  DOI : 10.1016/j.tpb.2017.05.001
• Sampling method for sign changing contrast
  
  • Audibert Lorenzo
  Inverse Problems and Imaging, AIMS American Institute of Mathematical Sciences, 2017. We extend the applicability of the Generalized Linear Sampling Method (GLSM) [2] and the Factorization Method (FM)[14] to the case of inhomogeneities where the contrast change sign strictly inside the obstacle. Both methods give an exact characterization of the target shapes in term of the fareld operator (at a xed frequency). One of the key ingredient to prove this exact characterization is based on a factorization of the fareld operator. This factorization involves three operators which should exhibit specic properties. This paper is concerned with the extension of the coercivity property required on one of them to the case of sign changing contrast both for isotropic and anisotropic scatters with possibly dierent supports for the isotropic and anisotropic parts. We fnally validate the method through some numerical tests 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
• Optimal control of slender microswimmers
  
  • Zoppello Marta
  • Desimone Antonio
  • Alouges François
  • Giraldi Laetitia
  • Martinon Pierre
  , 2017, pp.21. We discuss a reduced model to compute the motion of slender swimmers which propel themselves by changing the curvature of their body. Our approach is based on the use of Resistive Force Theory for the evaluation of the viscous forces and torques exerted by the surrounding fluid, and on discretizing the kinematics of the swimmer by representing its body through an articulated chain of N rigid links capable of planar deformations. The resulting system of ODEs governing the motion of the swimmer is easy to assemble and to solve, making our reduced model a valuable tool in the design and optimization of bio-inspired artificial microdevices. We prove that the swimmer is controllable in the whole plane for N is greater of equal to 3 and for almost every set of stick lengths. As a direct result, there exists an optimal swimming strategy to reach a desired configuration in minimum time. Numerical experiments for N = 3 (Purcell swimmer) suggest that the optimal strategy is periodic, namely a sequence of identical strokes. Our results indicate that this candidate for an optimal stroke indeed gives a better displacement speed than the classical Purcell stroke. (10.1007/978-3-319-73371-5_8)
  DOI : 10.1007/978-3-319-73371-5_8
• 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
• 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
• 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
• Multipoint scatterers with zero-energy bound states
  
  • Grinevich Piotr
  • Novikov Roman
  Theoretical and Mathematical Physics, Consultants bureau, 2017, 193 (2), pp.1675-1679. We study multipoint scatterers with zero-energy bound states in three dimensions. We present examples of such scatterers with multiple zero eigenvalue or with strong multipole localization of zero-energy bound states.
• Regularity for the optimal compliance problem with length penalization
  
  • Chambolle Antonin
  • Lamboley Jimmy
  • Lemenant Antoine
  • Stepanov Eugene
  SIAM Journal on Mathematical Analysis, Society for Industrial and Applied Mathematics, 2017. We prove some regularity results for a connected set S in the planar domain O, which minimizes the compliance of its complement O\S, plus its length. This problem, interpreted as to find the best location for attaching a membrane subject to a given external force f so as to minimize the compliance, can be seen as an elliptic PDE version of the average distance problem/irrigation problem (in a penalized version rather than a constrained one), which has been extensively studied in the literature. We prove that minimizers consist of a finite number of smooth curves meeting only by three at 120 degree angles, containing no loop, and possibly touching the boundary of the domain only tangentially. Several new technical tools together with the classical ones are developed for this purpose.
• 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