Publications

2012

• Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb')
  
  • Colonna Jean-François
  , 2012. Close-up on a pseudo-octonionic Mandelbrot set (a 'Mandelbulb') (Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb'))
• Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb')
  
  • Colonna Jean-François
  , 2012. Close-up on a pseudo-octonionic Mandelbrot set (a 'Mandelbulb') (Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb'))
• Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb')
  
  • Colonna Jean-François
  , 2012. Close-up on a pseudo-octonionic Mandelbrot set (a 'Mandelbulb') (Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb'))
• Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb')
  
  • Colonna Jean-François
  , 2012. Close-up on a pseudo-octonionic Mandelbrot set (a 'Mandelbulb') (Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb'))
• Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb')
  
  • Colonna Jean-François
  , 2012. Close-up on a pseudo-octonionic Mandelbrot set (a 'Mandelbulb') (Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb'))
• Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb')
  
  • Colonna Jean-François
  , 2012. Close-up on a pseudo-octonionic Mandelbrot set (a 'Mandelbulb') (Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb'))
• Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb')
  
  • Colonna Jean-François
  , 2012. Close-up on a pseudo-octonionic Mandelbrot set (a 'Mandelbulb') (Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb'))
• Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb')
  
  • Colonna Jean-François
  , 2012. Close-up on a pseudo-octonionic Mandelbrot set (a 'Mandelbulb') (Agrandissement d'un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'Mandelbulb'))
• Noise Source Localization in an Attenuating Medium
  
  • Bretin Elie
  • Ammari Habib
  • Garnier Josselin
  • Wahab Abdul
  SIAM Journal on Applied Mathematics, Society for Industrial and Applied Mathematics, 2012, 72 (1), pp.317 - 336. (10.1137/11083191X)
  DOI : 10.1137/11083191X
• A general stochastic model for sporophytic self-incompatibility
  
  • Billiard Sylvain
  • Tran Viet Chi
  Journal of Mathematical Biology, Springer, 2012, 64 (1-2), pp.163-210. Disentangling the processes leading populations to extinction is a major topic in ecology and conservation biology. The difficulty to find a mate in many species is one of these processes. Here, we investigate the impact of self-incompatibility in flowering plants, where several inter-compatible classes of individuals exist but individuals of the same class cannot mate. We model pollen limitation through different relationships between mate availability and fertilization success. After deriving a general stochastic model, we focus on the simple case of distylous plant species where only two classes of individuals exist. We first study the dynamics of such a species in a large population limit and then, we look for an approximation of the extinction probability in small populations. This leads us to consider inhomogeneous random walks on the positive quadrant. We compare the dynamics of distylous species to self-fertile species with and without inbreeding depression, to obtain the conditions under which self-incompatible species could be less sensitive to extinction while they can suffer more pollen limitation. (10.1007/s00285-011-0410-z)
  DOI : 10.1007/s00285-011-0410-z
• Computing estimates of material properties from transmission eigenvalues
  
  • Giorgi Giovanni
  • Haddar Houssem
  Inverse Problems, IOP Publishing, 2012, 28 (5), pp.055009, 23. This work is motivated by inverse scattering problems, those problems where one is interested in reconstructing the shape and the material properties of an inclusion from electromagnetic farfield measurements. More precisely, we are interested in complementing the so-called sampling methods by providing an estimate of the material properties of the sought inclusion. We use for this purpose a measure of the first transmission eigenvalue. Our method is then based on computing the desired estimate by reformulating the so-called interior transmission eigenvalue problem as an eigenvalue problem for the material coefficients. We will restrict ourselves to the two-dimensional setting of the problem and treat the cases of both transverse electric and transverse magnetic polarizations. We present a number of numerical experiments that validate our methodology for homogeneous and inhomogeneous inclusions and backgrounds. We also treat the case of a background with absorption and the case of scatterers with multiple connected components of different refractive indices. (10.1088/0266-5611/28/5/055009)
  DOI : 10.1088/0266-5611/28/5/055009
• Contacts in dimension 2, A penalization method
  
  • Pantz Olivier
  , 2012.
• On Simultaneous Identification of the Shape and Generalized Impedance Boundary Condition in Obstacle Scattering
  
  • Bourgeois Laurent
  • Chaulet Nicolas
  • Haddar Houssem
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2012, 34 (3), pp.A1824-A1848. We consider the inverse obstacle scattering problem of determining both the shape and the "equiva- lent impedance" from far field measurements at a fixed frequency. In this work, the surface impedance is represented by a second order surface differential operator (refer to as generalized impedance boundary condition) as opposed to a scalar function. The generalized impedance boundary condition can be seen as a more accurate model for effective impedances and is widely used in the scattering problem for thin coatings. Our approach is based on a least square optimization technique. A major part of our analysis is to characterize the derivative of the cost function with respect to the boundary and this complex surface impedance configuration. In particular, we provide an extension of the notion of shape derivative to the case where the involved impedance parameters do not need to be surface traces of given functions, which leads (in general) to a non-vanishing tangential boundary perturbation. The efficiency of considering this type of derivative is illustrated by several 2D numerical experiments based on a (classical) steepest descent method. The feasibility of retrieving both the shape and the impedance parameters is also discussed in our numerical experiments. (10.1137/110850347)
  DOI : 10.1137/110850347
• Asymptotic and non asymptotic approximations for option valuation
  
  • Bompis Romain
  • Gobet Emmanuel
  , 2012, pp.80. We give a broad overview of approximation methods to derive analytical formulas for accurate and quick evaluation of option prices. We compare different approaches, from the theoretical point of view regarding the tools they require, and also from the numerical point of view regarding their performances. In the case of local volatility models with general time-dependency, we derive new formulas using the local volatility function at the mid-point between strike and spot: in general, our approximations outperform previous ones by Hagan and Henry-Labordère. We also provide approximations of the option delta.
• Analysis of exposure–response of CI-945 in patients with epilepsy: application of novel mixed hidden Markov modeling methodology
  
  • Delattre Maud
  • Savic Radojka M.
  • Miller Raymond
  • Karlsson Mats O.
  • Lavielle Marc
  Journal of Pharmacokinetics and Pharmacodynamics, Springer Verlag, 2012, 39 (3), pp.263 - 271. We propose to describe exposure–response relationship of an antiepileptic agent, using mixed hidden Markov modeling methodology, to reveal additional insights in the mode of the drug action which the novel approach offers. Daily seizure frequency data from six clinical studies including patients who received gabapentin were available for the analysis. In the model, seizure frequencies are governed by underlying unobserved disease activity states. Individual neighbouring states are dependent, like in reality and they exhibit their own dynamics with patients transitioning between low and high disease states, according to a set of transition probabilities. Our methodology enables estimation of unobserved disease dynamics and daily seizure frequencies in all disease states. Additional modes of drug action are achievable: gabapentin may influence both daily seizure frequencies and disease state dynamics. Gabapentin significantly reduced seizure frequencies in both disease activity states; however it did not significatively affect disease dynamics. Mixed hidden Markov modeling is able to mimic dynamics of seizure frequencies very well. It offers novel insights into understanding disease dynamics in epilepsy and gabapentin mode of action. (10.1007/s10928-012-9248-2)
  DOI : 10.1007/s10928-012-9248-2
• Adiabatic control of the Schrödinger equation via conical intersections of the eigenvalues
  
  • Boscain Ugo
  • Chittaro Francesca
  • Mason Paolo
  • Sigalotti Mario
  IEEE Transactions on Automatic Control, Institute of Electrical and Electronics Engineers, 2012, 57 (8), pp.1970-1983. In this paper, we present a constructive method to control the bilinear Schrödinger equation via two controls. The method is based on adiabatic techniques and works if the spectrum of the Hamiltonian admits eigenvalue intersections, and if the latter are conical (as it happens generically). In this framework, we are able to spread on several levels connected by conical intersections a state initially concentrated in a single energy level. We provide sharp estimates on the dependence of the error with respect to the controllability time. Moreover, we identify some special curves in the space of controls that improve the precision of the adiabatic approximation, when passing through conical intersections, with respect to classical adiabatic theory. (10.1109/TAC.2012.2195862)
  DOI : 10.1109/TAC.2012.2195862
• On 2-step, corank 2 nilpotent sub-Riemannian metrics
  
  • Barilari Davide
  • Boscain Ugo
  • Gauthier Jean-Paul
  SIAM Journal on Control and Optimization, Society for Industrial and Applied Mathematics, 2012, 50 (1), pp.559-582. In this paper we study the nilpotent 2-step, corank 2 sub-Riemannian metrics that are nilpotent approximations of general sub-Riemannian metrics. We exhibit optimal syntheses for these problems. It turns out that in general the cut time is not equal to the first conjugate time but has a simple explicit expression. As a byproduct of this study we get some smoothness properties of the spherical Hausdorff measure in the case of a generic 6 dimensional, 2-step corank 2 sub-Riemannian metric.
• Sensitivity analysis of energy contracts management problem by stochastic programming techniques
  
  • Cen Zhihao
  • Bonnans J. Frederic
  • Christel Thibault
  , 2012, 12 (2012), pp.447-471.. We consider a model of medium-term commodity contracts management. Randomness takes place only in the prices on which the commodities are exchanged whilst state variable is multi-dimensional. In our previous article, we proposed an algorithm to deal with such problem, based on quantization of random process and a dual dynamic programming type approach. We obtained accurate estimates of the optimal value and a suboptimal strategy from this algorithm. In this paper, we analyse the sensitivity with respect to parameters driving the price model. We discuss the estimate of marginal price based on the Danskin's theorem. Finally, some numerical results applied to realistic energy market problems have been performed. Comparisons between results obtained by our algorithm and other classical methods are provided and evidence the accuracy of the estimate of marginal prices.
• A la recherche de l'inégalité perdue
  
  • Allaire Grégoire
  Matapli, Société de Mathématiques Appliquées et Industrielles (SMAI), 2012 (98), pp.52-64. A l'occasion du centenaire de la disparition de Henri Poincaré on s'intéresse à l'origine de la célèbre inégalité de Poincaré, si utile dans l'analyse de certaines équations aux dérivées partielles. Comme souvent en mathématiques les bons résultats ne sont pas attribués aux bonnes personnes... En fait, cette inégalité est probablement due à Neumann, Schwarz ou Scheeffer et on doit plutôt à Henri Poincaré une variante, appelée communément inégalité de Poincaré-Wirtinger. Sa démonstration originale est très intéressante car elle est constructive, au contraire d'une démonstration ultérieure par contradiction, utilisant le théorème de compacité de Rellich, que l'on trouve désormais dans la plupart des ouvrages sur les équations aux dérivées partielles. La motivation de cette inégalité est l'étude des valeurs propres et fonctions propres du Laplacien.
• Convergence of time-space adaptive algorithms for nonlinear conservation laws
  
  • Coquel F.
  • Postel M.
  • Tran Q.-H.
  IMA Journal of Numerical Analysis, Oxford University Press (OUP), 2012, 32 (4), pp.1440 - 1483. A family of explicit adaptive algorithms is designed to solve nonlinear scalar one-dimensional conservation laws. Based on the Godunov scheme on a uniform grid, a first strategy uses the multiresolution analysis of the solution to design an adaptive grid that evolves in time according to the time-dependent local smoothness. The method is furthermore enhanced by a local time-stepping strategy. Both numerical schemes are shown to converge towards the unique entropy solution. (10.1093/imanum/drr054)
  DOI : 10.1093/imanum/drr054
• Palaeoenvironmental shifts drove the adaptive radiation of a noctuid stemborer tribe (Lepidoptera, Noctuidae, Apameini) in the Miocene
  
  • Toussaint Emmanuel F.A.
  • Condamine Fabien L.
  • Kergoat Gael
  • Capdevielle-Dulac Claire
  • Barbut Jérôme
  • Silvain Jean-François
  • Le Ru Bruno P.
  PLoS ONE, Public Library of Science, 2012, 7 (7), pp.np. Between the late Oligocene and the early Miocene, climatic changes have shattered the faunal and floral communities and drove the apparition of new ecological niches. Grassland biomes began to supplant forestlands, thus favouring a large-scale ecosystem turnover. The independent adaptive radiations of several mammal lineages through the evolution of key innovations are classic examples of these changes. However, little is known concerning the evolutionary history of other herbivorous groups in relation with this modified environment. It is especially the case in phytophagous insect communities, which have been rarely studied in this context despite their ecological importance. Here, we investigate the phylogenetic and evolutionary patterns of grass-specialist moths from the species-rich tribe Apameini (Lepidoptera, Noctuidae). The molecular dating analyses carried out over the corresponding phylogenetic framework reveal an origin around 29 million years ago for the Apameini. Ancestral state reconstructions indicate (i) a potential Palaearctic origin of the tribe Apameini associated with a major dispersal event in Afrotropics for the subtribe Sesamiina; (ii) a recent colonization from Palaearctic of the New World and Oriental regions by several independent lineages; and (iii) an ancestral association of the tribe Apameini over grasses (Poaceae). Diversification analyses indicate that diversification rates have not remained constant during the evolution of the group, as underlined by a significant shift in diversification rates during the early Miocene. Interestingly, this age estimate is congruent with the development of grasslands at this time. Rather than clade ages, variations in diversification rates among genera better explain the current differences in species diversity. Our results underpin a potential adaptive radiation of these phytophagous moths with the family Poaceae in relation with the major environmental shifts that have occurred in the Miocene. (10.1371/journal.pone.0041377)
  DOI : 10.1371/journal.pone.0041377
• Un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'MandelBulb')
  
  • Colonna Jean-François
  , 2012. A pseudo-octonionic Mandelbrot set (a 'MandelBulb') (Un ensemble de Mandelbrot dans l'ensemble des pseudo-octonions (un 'MandelBulb'))
• Transmission eigenvalues for inhomogeneous media containing obstacles
  
  • Cakoni Fioralba
  • Cossonnière Anne
  • Haddar Houssem
  Inverse Problems and Imaging, AIMS American Institute of Mathematical Sciences, 2012, 6 (3), pp.373-398. (10.3934/ipi.2012.6.373)
  DOI : 10.3934/ipi.2012.6.373
• Generalized fractional smoothness and Lp-variation of BSDEs with non-Lipschitz terminal condition
  
  • Geiss Christel
  • Geiss Stefan
  • Gobet Emmanuel
  Stochastic Processes and their Applications, Elsevier, 2012, 122 (5), pp.2078--2116. We relate the $L_p$-variation, $2\le p < \infty$, of a solution of a backward stochastic differential equation with a path-dependent terminal condition to a generalized notion of fractional smoothness. This concept of fractional smoothness takes into account the quantitative propagation of singularities in time.
• A Patchy Dynamic Programming Scheme for a Class of Hamilton-Jacobi-Bellman Equations
  
  • Cacace Simone
  • Cristiani Emiliano
  • Falcone Maurizio
  • Picarelli Athena
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2012, 34 (5), pp.A2625-A2649. In this paper we present a new parallel algorithm for the solution of Hamilton-Jacobi- Bellman equations related to optimal control problems. The main idea is to divide the domain of computation into subdomains following the dynamics of the control problem. This results in a rather complex geometrical subdivision, but has the advantage that every subdomain is invariant with respect to the optimal controlled vector field, so that we can compute the value function in each subdomain assigning the task to a processor and avoiding the classical transmission condition on the boundaries of the subdomains. For this specific feature the subdomains are patches in the sense introduced by Ancona and Bressan in [1]. Several examples in dimension two and three illustrate the properties of the new method. (10.1137/110841576)
  DOI : 10.1137/110841576