BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Yvon Maday

Yvon Maday appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.APmath.OCNumerical AnalysisApplicationsmath-phmath.FAmath.MPMethodologyphysics.comp-ph

Catalog footprint

What is connected

15works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

15 published item(s)

preprint2022arXiv

A flexible smoother adapted to censored data with outliers and its application to SARS-CoV-2 monitoring in wastewater

A sentinel network, Obépine, has been designed to monitor SARS-CoV-2 viral load in wastewaters arriving at wastewater treatment plants (WWTPs) in France as an indirect macro-epidemiological parameter. The sources of uncertainty in such monitoring system are numerous and the concentration measurements it provides are left-censored and contain outliers, which biases the results of usual smoothing methods. Hence the need for an adapted pre-processing in order to evaluate the real daily amount of virus arriving to each WWTP. We propose a method based on an auto-regressive model adapted to censored data with outliers. Inference and prediction are produced via a discretised smoother which makes it a very flexible tool. This method is both validated on simulations and on real data from Obépine. The resulting smoothed signal shows a good correlation with other epidemiological indicators and is currently used by Obépine to provide an estimate of virus circulation over the watersheds corresponding to about 200 WWTPs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An EIM-degradation free reduced basis method via over collocation and residual hyper reduction-based error estimation

The need for multiple interactive, real-time simulations using different parameter values has driven the design of fast numerical algorithms with certifiable accuracies. The reduced basis method (RBM) presents itself as such an option. RBM features a mathematically rigorous error estimator which drives the construction of a low-dimensional subspace. A surrogate solution is then sought in this low-dimensional space approximating the parameter-induced high fidelity solution manifold. However when the system is nonlinear or its parameter dependence nonaffine, this efficiency gain degrades tremendously, an inherent drawback of the application of the empirical interpolation method (EIM). In this paper, we augment and extend the EIM approach as a direct solver, as opposed to an assistant, for solving nonlinear partial differential equations on the reduced level. The resulting method, called Reduced Over-Collocation method (ROC), is stable and capable of avoiding the efficiency degradation. Two critical ingredients of the scheme are collocation at about twice as many locations as the number of basis elements for the reduced approximation space, and an efficient error indicator for the strategic building of the reduced solution space. The latter, the main contribution of this paper, results from an adaptive hyper reduction of the residuals for the reduced solution. Together, these two ingredients render the proposed R2-ROC scheme both offline- and online-efficient. A distinctive feature is that the efficiency degradation appearing in traditional RBM approaches that utilize EIM for nonlinear and nonaffine problems is circumvented, both in the offline and online stages. Numerical tests on different families of time-dependent and steady-state nonlinear problems demonstrate the high efficiency and accuracy of our R2-ROC and its superior stability performance.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Guaranteed a posteriori bounds for eigenvalues and eigenvectors: multiplicities and clusters

This paper presents a posteriori error estimates for conforming numerical approximations of eigenvalue clusters of second-order self-adjoint elliptic linear operators with compact resolvent. Given a cluster of eigenvalues, we estimate the error in the sum of the eigenvalues, as well as the error in the eigenvectors represented through the density matrix, i.e., the orthogonal projector on the associated eigenspace. This allows us to deal with degenerate (multiple) eigenvalues within the framework. All the bounds are valid under the only assumption that the cluster is separated from the surrounding smaller and larger eigenvalues; we show how this assumption can be numerically checked. Our bounds are guaranteed and converge with the same speed as the exact errors. They can be turned into fully computable bounds as soon as an estimate on the dual norm of the residual is available, which is presented in two particular cases: the Laplace eigenvalue problem discretized with conforming finite elements, and a Schr{ö}dinger operator with periodic boundary conditions of the form $--$Δ$ + V$ discretized with planewaves. For these two cases, numerical illustrations are provided on a set of test problems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Parametric Analytical Preconditioning and its Applications to the Reduced Collocation Methods

In this paper, we extend the recently developed reduced collocation method \cite{ChenGottlieb} to the nonlinear case, and propose two analytical preconditioning strategies. One is parameter independent and easy to implement, the other one has the traditional affinity with respect to the parameters which allows for efficient implementation through an offline-online decomposition. Overall, the preconditioning improves the quality of the error estimation uniformly on the parameter domain, and speeds up the convergence of the reduced solution to the truth approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model

We present a parareal in time algorithm for the simulation of neutron diffusion transient model. The method is made efficient by means of a coarse solver defined with large time steps and steady control rods model. Using finite element for the space discretization, our implementation provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch-Maurer-Werner (LMW) benchmark [1].

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Robin Schwarz algorithm for the NICEM Method: the Pq finite element case

In Gander et al. [2004] we proposed a new non-conforming domain decomposition paradigm, the New Interface Cement Equilibrated Mortar (NICEM) method, based on Schwarz type methods that allows for the use of Robin interface conditions on non-conforming grids. The error analysis was done for P1 finite elements, in 2D and 3D. In this paper, we provide new numerical analysis results that allow to extend this error analysis in 2D for piecewise polynomials of higher order and also prove the convergence of the iterative algorithm in all these cases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparison of some Reduced Representation Approximations

In the field of numerical approximation, specialists considering highly complex problems have recently proposed various ways to simplify their underlying problems. In this field, depending on the problem they were tackling and the community that are at work, different approaches have been developed with some success and have even gained some maturity, the applications can now be applied to information analysis or for numerical simulation of PDE's. At this point, a crossed analysis and effort for understanding the similarities and the differences between these approaches that found their starting points in different backgrounds is of interest. It is the purpose of this paper to contribute to this effort by comparing some constructive reduced representations of complex functions. We present here in full details the Adaptive Cross Approximation (ACA) and the Empirical Interpolation Method (EIM) together with other approaches that enter in the same category.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Locally adaptive greedy approximations for anisotropic parameter reduced basis spaces

Reduced order models, in particular the reduced basis method, rely on empirically built and problem dependent basis functions that are constructed during an off-line stage. In the on-line stage, the precomputed problem-dependent solution space, that is spanned by the basis functions, can then be used in order to reduce the size of the computational problem. For complex problems, the number of basis functions required to guarantee a certain error tolerance can become too large in order to benefit computationally from the model reduction. To overcome this, the present work introduces a framework where local approximation spaces (in parameter space) are used to define the reduced order approximation in order to have explicit control over the on-line cost. This approach also adapts the local approximation spaces to local anisotropic behavior in the parameter space. We present the algorithm and numerous numerical tests.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-consistent approximations of self-adjoint eigenproblems: Application to the supercell method

In this article, we introduce a general theoretical framework to analyze non-consistent approximations of the discrete eigenmodes of a self-adjoint operator. We focus in particular on the discrete eigenvalues laying in spectral gaps. We first provide a priori error estimates on the eigenvalues and eigenvectors in the absence of spectral pollution. We then show that the supercell method for perturbed periodic Schrödinger operators falls into the scope of our study. We prove that this method is spectral pollution free, and we derive optimal convergence rates for the planewave discretization method, taking numerical integration errors into account. Some numerical illustrations are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Parareal in time intermediate targets methods for optimal control problem

In this paper, we present a method that enables solving in parallel the Euler-Lagrange system associated with the optimal control of a parabolic equation. Our approach is based on an iterative update of a sequence of intermediate targets that gives rise to independent sub-problems that can be solved in parallel. This method can be coupled with the parareal in time algorithm. Numerical experiments show the efficiency of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stable parareal in time method for first and second order hyperbolic system

The parareal in time algorithm allows to perform parallel simulations of time dependent problems. This algorithm has been implemented on many types of time dependent problems with some success. Recent contributions have allowed to extend the domain of application of the parareal in time algorithm so as to handle long time simulations of Hamiltonian systems. This improvement has managed to avoid the fatal large lack of accuracy of the plain parareal in time algorithm consequence of the fact that the plain parareal in time algorithm does not conserve invariants. A somehow similar difficulty occurs for problems where the solution lacks regularity, either initially or in the evolution, like for the solution to hyperbolic system of conservation laws. In this paper we identify the problem of lack of stability of the parareal in time algorithm and propose a simple way to cure it. The new method is used to solve a linear wave equation and a non linear Burger's equation, the results illustrate the stability of this variant of the parareal in time algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An intermediate targets method for time parallelization in optimal control

In this paper, we present a method that enables to solve in parallel the Euler-Lagrange system associated with the optimal control of a parabolic equation. Our approach is based on an iterative update of a sequence of intermediate targets and gives rise independent sub-problems that can be solved in parallel. Numerical experiments show the efficiency of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Periodic Schrödinger operators with local defects and spectral pollution

This article deals with the numerical calculation of eigenvalues of perturbed periodic Schrödinger operators located in spectral gaps. Such operators are encountered in the modeling of the electronic structure of crystals with local defects, and of photonic crystals. The usual finite element Galerkin approximation is known to give rise to spectral pollution. In this article, we give a precise description of the corresponding spurious states. We then prove that the supercell model does not produce spectral pollution. Lastly, we extend results by Lewin and Séré on some no-pollution criteria. In particular, we prove that using approximate spectral projectors enables one to eliminate spectral pollution in a given spectral gap of the reference periodic Schödinger operator.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical analysis of the planewave discretization of some orbital-free and Kohn-Sham models

We provide a priori error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizsäcker (TFW) model and for the spectral discretization of the Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the ground state energy and density of molecular systems in the condensed phase. The TFW model is stricly convex with respect to the electronic density, and allows for a comprehensive analysis. This is not the case for the Kohn-Sham LDA model, for which the uniqueness of the ground state electronic density is not guaranteed. Under a coercivity assumption on the second order optimality condition, we prove that for large enough energy cut-offs, the discretized Kohn-Sham LDA problem has a minimizer in the vicinity of any Kohn-Sham ground state, and that this minimizer is unique up to unitary transform. We then derive optimal a priori error estimates for the spectral discretization method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reduced basis techniques for stochastic problems

We report here on the recent application of a now classical general reduction technique, the Reduced-Basis approach initiated in [C. Prud'homme, D. Rovas, K. Veroy, Y. Maday, A. T. Patera, and G. Turinici. Reliable real-time solution of parametrized partial differential equations: Reduced-basis output bounds methods. Journal of Fluids Engineering, 124(1):7080, 2002.], to the specific context of differential equations with random coefficients. After an elementary presentation of the approach, we review two contributions of the authors: [S. Boyaval, C. Le Bris, Y. Maday, N.C. Nguyen, and A.T. Patera. A reduced basis approach for variational problems with stochastic parameters: Application to heat conduction with variable Robin co-efficient. Computer Methods in Applied Mechanics and Engineering, 198(4144):3187-3206, 2009.], which presents the application of the RB approach for the discretization of a simple second order elliptic equation supplied with a random boundary condition, and [S. Boyaval and T. Lelièvre, A variance reduction method for parametrized stochastic differential equations using the reduced basis paradigm with T. Lelièvre, Commun. Math. Sci. 8, special Issue "Mathematical Issue on Complex Fluids" P. Zhang ed., to appear, 2010, ARXIV preprint arXiv:0906.3600], which uses a RB type approach to reduce the variance in the Monte-Carlo simulation of a stochastic differential equation. We conclude the review with some general comments and also discuss possible tracks for further research in the direction.

0 citations0 reviewsNo code linkedOpen access