Jean-Marie Mirebeau
Jean-Marie Mirebeau
Winner of the 8th Popov Prize, for Outstanding Research Achievements in Approximation Theory !
Prize received at the 15th International Conference on Approximation Theory, May 2016.
Leader of the ANR project NS-LBR, for the development of Numerical Schemes using Lattice Basis Reduction. Members, Introductory notebooks, Publications.
Research Interests
I develop numerical schemes for Partial Differential Equations (PDEs), focusing on issues related to Anisotropy. I also investigate their applications, in particular to medical image processing.
Keywords : Anisotropic Partial Differential Equations,
Anisotropic Diffusion, Anisotropic Eikonal equations, Monge-Ampere equations, Convexity Constraint.
Anisotropic Diffusion
Nonlinear Anisotropic Diffusion is involved in state of the art image processing methods, as well as in physical phenomena. We developed a discretization for this PDE, on grids, which combines low numerical cost and strong mathematical guarantees.
Scheme description and analysis:
-Fehrenbach, M., Sparse Non-Negative Stencils for Anisotropic Diffusion, JMIV, 2014
-M., Minimal Stencils for Monotony or Causality Preserving Discretizations of Anisotropic PDEs, Preprint
Open source code available as an Insight Toolkit module :
Anisotropic diffusion in ITK, 2015 (module description), Interactive Figures
Anisotropic Eikonal Equations
The eikonal equation is a Partial Differential Equation (PDE) which characterizes the geodesic distance on a (Finsler or Riemannian) manifold. It has numerous applications, including medical image processing.
I have developed solvers for this PDE which are particularly fast and accurate in the context of large anisotropies. Slides
Design of causal stencils:
-M., Efficient Fast Marching with Finsler Metrics, Numerische Math, 2014
-M., Anisotropic Fast-Marching on cartesian grids using Lattice Basis Reduction, SINUM, 2014
Applications to medical image processing:
- Sanguinetti, Bekkers, Duits, Janssen, Mashtakov, M., Sub-Riemannian Fast Marching in SE2, Preprint, 2015
- Chen, Cohen, M.,Vessel Extraction using Anisotropic Minimal Paths and Path Score, ICIP 2014
Open source code for the Insight Toolkit :
Anisotropic Fast Marching in ITK, 2015
Monge-Ampere Equations and the Constraint of Convexity
The Monge-Ampere equation is a fully non-linear PDE, which requires specific discretizations. I study two classes of numerical schemes, based on the general method of monotone wide stencil schemes, and the more specific approach of semi-discrete optimal transport based on weighted Voronoi diagrams, implemented in CGAL.
Slides on semi discrete optimal transport, and volume preserving maps (by Q. Merigot)
Slides on the constraint of convexity.
PDE discretizations based on wide stencil schemes:
-M., Adaptive, Anisotropic and Hierarchical cones of Discrete Convex Functions, Numerische Math, 2015
-Benamou, Collino, M., Monotone and Consistent Discretization of the Monge Ampere Operator, Math of Comp, 2015
-M., Discretization of the 3D Monge-Ampere Operator, between Wide Stencils and Power Diagrams, M2AN, 2015
Applications of semi-discrete optimal transport:
-M., Numerical resolution of Euler equations, through semi-discrete optimal transport, Proceedings des journées EDP de Roscoff 2015
-Merigot, M., Minimal geodesics along volume-preserving maps, through semi-discrete optimal transport, Preprint
Adaptive and Anisotropic Finite Element Approximation
Video réalisée par EADS sur mes travaux
Mesh adaption procedures for finite elements approximation allow to adapt locally the resolution, by local refinement around the places of strong variation of the function. The use of anisotropic triangles allows to improve the efficiency of the procedure by introducing long and thin triangles fitting in particular the directions of the possible curves of discontinuity.
It is therefore natural to try to characterize, an produce optimal meshes, in the sense of the compromise between approximation error and complexity of the mesh, for a given function.
In the case of smooth functions, we have derived optimal error estimates and propose a caracterisation of the corresponding meshes (Article). We have also elaborated a hierarchical algorithm that produces such anisotropic triangulations, and which optimality has been proved for certain classes of functions.
Sharp asymptotic interpolation error:
-M., Optimally Adapted Meshes for Finite Elements of Arbitrary Order and W1p Norms, Numer Math, 2011
-Y. Babenko, T. Leskevich, M., Sharp asymptotics of the L p approximation error for interpolation on block partitions, Numerische Mathematik, 2011
-M., Optimal meshes for finite elements of arbitrary order, Constructive Approximation, 2010
Greedy refinement algorithms:
-A. Cohen, N. Dyn, F. Hecht, M., Adaptive multiresolution analysis based on anisotropic triangulations, Mathematics of Computation, 2011
-A. Cohen, M., Greedy bisection generates optimally adapted triangulations, Mathematics of Computation, 2011
-A. Cohen, M., Adaptive and anisotropic piecewise polynomial approximation, Chapter of the book Multiscale, Nonlinear and Adaptive Approximation , 2009
Link with image analysis:
-A. Cohen, M., Anisotropic smoothness classes: from finite element approximation to image models, JMIV, 2010
Design of Riemannian metrics for mesh generation:
-M., The optimal aspect ratio for piecewise quadratic anisotropic finite element approximation, proceedings of the conference SampTA 2011
Other subjects in numerical analysis
-J. Bleyer, G. Carlier, V. Duval, M., G. Peyré, A Gamma-Convergence Result for the Upper Bound Limit Analysis of Plates, M2AN, 2015
-M., Non conforming vector finite elements for H (curl) intersected with H (div), Applied Mathematic Letters, 2011
Awards :
Prix de la meilleure thèse de la fondation d'entreprise EADS (2011), en Mathématiques et leurs interactions.
Prix solennel de la chancellerie des universités de Paris (2011), Perrissin-Pirasset/Schneider, en Mathématiques Fondamentales et Appliquées, pour ma thèse.
CNRS researcher in Applied Mathematics
Laboratoire de mathématiques d’Orsay, Université Paris-Sud, CNRS,
Université Paris-Saclay, 91405 Orsay
Building 425, Office 314
jean-marie.mirebeau «at» math.u-psud.fr
Address : Département de Mathématiques Bâtiment 425
Faculté des Sciences d'Orsay Université Paris-Sud
F-91405 Orsay Cedex
Determination of the optimal aspect ratio, for anisotropic mesh generation
Anisotropic refinement by bisection
(1) Geometry of numerical scheme stencils for anisotropic eikonal equations. (2) An example of front propagation.
Three kinds of consumers in a principal agent model.
(Left) Original image. (Right) Effect of anisotropic diffusion.
Receiving the Popov Prize, with Albert Cohen (left, former PhD advisor), Pencho Petrushev (right, chairman of the selection comittee). Larger picture.