The Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation
|
|
- August Hopkins
- 6 years ago
- Views:
Transcription
1 The Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation Y. Maday (LJLL, I.U.F., Brown Univ.) Olga Mula (CEA and LJLL) G. Turinici (CEREMADE) SampTa 2013 Bremen, 01-05/07/2013 Olga MULA GEIM 1/ 26
2 Motivations of the present work Let Ω R d be a bounded spatial domain. Classical Lagrangian interpolation Classical interpolation problem: Approximate f C 0 (Ω) in the finite n-dimensional space X n = span {q i C 0 (Ω)} n i=1 where: the basis functions q i are of polynomial nature (algebraic polynomials, rational functions...). the interpolation points: d = 1: an almost optimal location is provided by the Chebychev-Gauss nodes. d > 1: more complex conditions in order for a polynomial interpolation to be well defined. Answers exist for domains of simple shapes (e.g. those obtained through tensor-product operations). Olga MULA GEIM 2/ 26
3 Motivations of the present work Extension of classical Lagrangian interpolation Let F be a compact subset of X = C 0 (Ω). Hypothesis: F is supposed to be such that any f F is approximable by linear combinations of small size we suppose that the Kolmogorov n-width d n (F,X) of F in X is small. d n (F,X) := inf sup inf x y X X n X x F y X n dim(xn)=n Extended interpolation problem: Approximate any f F by finite dimensional subspaces X n = span {q i F,i = 1,...,n} contained in the span of F. Here: the basis functions q i are very general (the only requirement is that q i F). given X n, the general conditions for which the interpolation points give an unique interpolant is an open problem. Olga MULA GEIM 3/ 26
4 Motivations of the present work Generalized Empirical Interpolation Method (GEIM) Empirical Interpolation Method (EIM) [Maday, Nguyen, Patera, Pau, 2009], [Barrault, Maday, Nguyen, Patera, 2004], [Grepl, Maday, Nguyen, Patera, 2007]: proposes a constructive approach to find the basis functions {q i F,i = 1,...,n} and the interpolation points {x i Ω, i = 1,...,n} through a Greedy algorithm. Generalized EIM: Replaces the set of interpolation points {x i Ω, i = 1,...,n} by continuous linear functions {σ i Σ, i = 1,...,n}, where the dictionary Σ L(F). A Greedy algorithm provides the basis functions and the interpolating linear functions. Unisolvence of Σ: if ϕ F is such that σ(ϕ) = 0 for any σ Σ, then ϕ = 0. Olga MULA GEIM 4/ 26
5 Outline 1 Theoretical properties of GEIM 2 Olga MULA GEIM 5/ 26
6 The generalized interpolation Assume you are given a set {q 1,...,q M } of linearly independent functions coming from a space of small Kolmogorov n-width AND a dictionary of continuous linear forms {σ}. Given a function f that you want to approximate, the problem is: find a family of scalars {α n M } 1 n M such that σ i (f) = M α M n σ i (q n ), for i = {1 n M}. n=1 where the linear forms σ i M are suitably chosen and X M = span {q n, 1 n M}. If σ i = δ( r i ) EIM. Olga MULA GEIM 6/ 26
7 The generalized interpolation Assume you are given a set {q 1,...,q M } of linearly independent functions coming from a space of small Kolmogorov n-width AND a dictionary of continuous linear forms {σ}. Given a function f that you want to approximate, the problem is: find a family of scalars {α n M } 1 n M such that σ i (f) = M α M n σ i (q n ), for i = {1 n M}. n=1 where the linear forms σ i M are suitably chosen and X M = span {q n, 1 n M}. The classical continuity assumption can be relaxed: L 2 (Ω) instead of C 0. Olga MULA GEIM 6/ 26
8 Among the questions raised by GEIM: is there an optimal selection for the linear forms σ i? is there a constructive optimal selection for the functions q i? given a set of linearly independent functions {q i } i [1,M] and a set of continuous linear forms {σ i } i [1,M], does the interpolant exist? is the interpolant unique? how does the interpolation process compares with other approximations (in particular orthogonal projections)? Under what hypothesis can we expect the GEIM approximation to converge rapidly to f? Olga MULA GEIM 7/ 26
9 The Greedy Algorithm for GEIM We are looking for a constructive way of approximating in X. From now on, we assume that X L 2. We propose a Greedy approach for selecting the interpolation linear forms and constructing the discrete spaces X M : X 1 X 2 X M X Olga MULA GEIM 8/ 26
10 The Greedy Algorithm for GEIM The Greedy algorithm for GEIM Constructive way of approximating in X: 1 Find: ϕ 0 = argsup ϕ F ϕ L 2 (Ω) and σ 0 = argsup σ Σ σ(ϕ 0 ) We define: q 0 = ϕ 0/σ 0(ϕ 0) X 1 span {q 0} J 0[ϕ] := σ 0(ϕ)q 0, ϕ F 2 By recursion: given, X N span {q j, j [0,N 1]} and {σ 0,σ 1,...,σ N 1 }, find: ϕ N = argsup ϕ F ϕ J N [ϕ] L 2 (Ω) σ N = argsup σ Σ σ(ϕ N J N [ϕ N ]) We define: q N = ϕ N J N [ϕ N ] σ N (ϕ N J N [ϕ N ]) X N+1 span {q j, j [0,N]} J N+1 : F X N+1 Olga MULA GEIM 9/ 26
11 Well-posedness It has been proven in [Maday, Mula, 2013] : Lemma For any N, the set {q j, j [0,N 1]} is linearly independent and X N is of dimension N. The generalized empirical interpolation procedure is well-posed in L 2 (Ω) and ϕ F, the interpolation error satisfies: ϕ J N [ϕ] L 2 (Ω) (1+Λ N ) inf ψ N X N ϕ ψ N L 2 (Ω) where Λ N is the Lebesgue constant in the L 2 norm: Λ N := sup ϕ F J N [ϕ] L 2 (Ω). ϕ L 2 (Ω) Olga MULA GEIM 10/ 26
12 A priory convergence As an extension of the results presented in [DeVore, Petrova, Wojtaszczyk, 2012], an optimal convergence result (modulo the Lebesgue constant) has been proven in [Maday, Mula, Turinici, 2013] Theorem 1 Assume that d n (F,L 2 (Ω)) C 0 n α for any n 1, then the interpolation error of the GEIM Greedy selection process satisfies for any ϕ F the inequality ϕ J n [ϕ] L 2 (Ω) C 0 (1+Λ n )β n n α, where the parameters γ n and β n are defined as: γ n = 1 1+Λ n β 1 = 2 n 2: β n = β 4l+k := 1 2β l1 l 2 i=1 l 1 = 2l+ 2k 3, l 2 = 2 ( l+ k 4 ) γ 1 l 2 l 1 k 4 +i (2 2) α and Olga MULA GEIM 11/ 26
13 A priory convergence As an extension of the results presented in [DeVore, Petrova, Wojtaszczyk, 2012], an optimal convergence result (modulo the Lebesgue constant) has been proven in [Maday, Mula, Turinici, 2013] Theorem 2 Assume that d n (F,L 2 (Ω)) C 0 e c 1n α for any n 1, then the interpolation error of the GEIM Greedy selection process satisfies for any ϕ F the inequality ϕ J n [ϕ] L 2 (Ω) C 0 (1+Λ n )β n e c 2n α, where β n and c 2 are defined as: c 2 := 2 1 3α c 1 β 1 = 2 for n 2: β n := n 2 i=1 γ 1 1 n 2 n 2 +i 2β n 2. Olga MULA GEIM 11/ 26
14 A priory convergence Corollary If (Λ n ) n 1 is a monotonically increasing sequence, then the sequence (γ n ) n 1 in the GEIM procedure is monotonically decreasing. The following decay rates in the generalized interpolation error can be derived: For any ϕ F, if d n (F,L 2 (Ω)) C 0 n α for any n 1, then the interpolation error of the GEIM Greedy selection process can be bounded as ϕ J n [ϕ] L 2 (Ω) C 0 2 3α+1 (1+Λ n ) 3 n α. For any ϕ F, if d n (F,L 2 (Ω)) C 0 e c 1n α for any n 1, then the interpolation error of the GEIM Greedy selection process can be bounded as ϕ J n [ϕ] L 2 (Ω) C 0 2(1+Λ n ) 3 e c 2n α. Olga MULA GEIM 12/ 26
15 A priory convergence Remark on the Lebesgue constant Our first experiments in the GEIM provide cases where it is uniformly bounded when evaluated in the L(L 2 ) norm. We do not pretend that this is universal, but only shows that the theoretical increasing upper bound is far from being optimal in a class of sets F that have a small Kolmogorov n-width Computed Lebesgue constant L 2 Lebesgue constant H 1 Lebesgue constant Dimension Figure : Estimation of the Lebesgue constant on a simple case (Laplace problem - see numerical example below) Olga MULA GEIM 13/ 26
16 Outline 1 Theoretical properties of GEIM 2 Olga MULA GEIM 14/ 26
17 Reduced basis Let E R p be a compact set and let D µ u = g µ E (1) be a family of parameter dependent PDE s representing a physical experiment. Hypothesis: {u µ, µ E} F where F X is compact and of small Kolmogorov n-width in X. Aim: For a given µ E, we wish to compute u µ in real time (i.e. at the same time as the experiment is taking place). Idea of reduced basis: Offline: We solve (1) for µ 1,...,µ n,... to build an M-dimensional reduced basis X M = span{q j F,j [1,M]}. Online: we approximate the solution by: u µ M a j q j j=1 Olga MULA GEIM 15/ 26
18 Reduced basis Suppose further that the experiment involves captor measurements we have a dictionary Σ L(F). Then, in this reduced basis approach, data from sensors can be taken into account if we apply our generalized interpolation procedure: Offline: Construction of X M = span{q j F,j [1,M]} and selection of {σ i,1 i M} by the Greedy algorithm. Online: u µ I M (u µ ). Olga MULA GEIM 16/ 26
19 Numerical application The Laplace problem Spatial domain: Ω R 2 and a non overlapping partition of it is: Ω = Ω 1 Ω2. Model: the Laplace problem. Hypothesis: only Ω 1 is driving the phenomenon. The rest (Ω 2 ) is the environment. Equation: u µ = f µ in Ω = Ω 1 Ω 2,µ f µ (x,y) = 1+(αsin(x)+βcos(γπy))χ 1 (x,y) Parameters: µ = (α,β,γ,geometry of Ω 2 ). Olga MULA GEIM 17/ 26
20 Aim We will suppose that the set of solutions {u µ Ω2, µ} F and that F is of small Kolmogorov n-width in L 2 (Ω 2 ). Given µ, we wish to: 1 Show that u µ Ω2 can be reconstructed with GEIM. 2 Run a real-time simulation: Ω 2 : Data acquisition + reconstruction with GEIM Ω 1 : Implicit calculation with the boundary conditions provided by the interpolant of Ω 2. Olga MULA GEIM 18/ 26
21 1) Reconstruction of Ω 2 with GEIM Interpolation error We have reconstructed several solutions and computed for 1 M 13: max u I M[u] ; = L 2 (Ω), H 1 (Ω) u reconstructed Interpolation error in the interpolated functions by GEIM L 2 interpolation error H 1 interpolation error Dimension Olga MULA GEIM 19/ 26
22 1) Reconstruction of Ω 2 with GEIM Comparison with POD Our GEIM approach seems to have explored the set of solutions good enough to perform as well as POD Interpolation error in the interpolated functions by GEIM L 2 interpolation error H 1 interpolation error H 1 : POD error L 2 : POD error Dimension Olga MULA GEIM 20/ 26
23 1) Reconstruction of Ω 2 with GEIM The FEM noise Going much beyond M = 13 is of no use since we reach the FEM precision ( 10 4 ) POD errors POD errors H 1 POD errors L Dimension Olga MULA GEIM 21/ 26
24 1) Reconstruction of Ω 2 with GEIM The Lebesgue constant A computational estimation of Λ M has been carried out: I M [u i ] L Λ M = max 2 (Ω) i [1,256] u i L 2 (Ω) Computed Lebesgue constant L 2 Lebesgue constant H 1 Lebesgue constant Dimension Olga MULA GEIM 22/ 26
25 2) Data acquisition coupled with simulation Application to the Laplace problem The interpolant in the environment Ω 2 can give the appropriate boudary conditions to the subdomain that has a big Kolmogorov n-width (Ω 1 ). We compute Ω 1 by FEM with the boundary conditions provided by the GEIM procedure Reconstructed analysis H 1 error in Ω H 1 error in Ω Olga MULA GEIM 23/ 26
26 Current work Treatment of noisy measurements: some developments in [Maday, Mula, 2013] Future work More complex applications (e.g. Stokes, Navier-Stokes, transport equations) Extension to time-dependent problems Olga MULA GEIM 24/ 26
27 For further reading [Maday, Nguyen, Patera, Pau, 2009] Y. Maday, N. Nguyen, A. Patera, and G. Pau, A general multipurpose interpolation procedure: the magic points Commun. Pure Appl. Anal., pp , [Barrault, Maday, Nguyen, Patera, 2004] M. Barrault, Y. Maday, N. Y. Nguyen, and A. Patera, An empirical interpolation method: Application to efficient reduced-basis discretization of partial differential equations C. R. Acad. Sci. Paris, Serie I., vol. 339, pp , [Grepl, Maday, Nguyen, Patera, 2007] M. Grepl, Y. Maday, N. Nguyen, and A. Patera, Efficient reduced basis treatment of nonaffine and nonlinear partial differential equations Math. Model. Numer. Anal., vol. 41(3), pp , Olga MULA GEIM 25/ 26
28 For further reading [DeVore, Petrova, Wojtaszczyk, 2012] R. DeVore, G. Petrova, and P. Wojtaszczyk Greedy algorithms for reduced bases in Banach spaces Constructive Approximation, pp. 1-12, 2012 [Maday, Mula, 2013] Y. Maday and O. Mula A generalized empirical interpolation method: application of reduced basis techniques to data assimilation Analysis and Numerics of Partial Differential Equations, vol. XIII, pp , 2013 [Maday, Mula, Turinici, 2013] Y. Maday, O. Mula and G. Turinici A priori convergence of the Generalized Empirical Interpolation Method accepted in the proceedings of the 10th International Conference on Sampling Theory and Applications, Olga MULA GEIM 26/ 26
Empirical Interpolation Methods
Empirical Interpolation Methods Yvon Maday Laboratoire Jacques-Louis Lions - UPMC, Paris, France IUF and Division of Applied Maths Brown University, Providence USA Doctoral Workshop on Model Reduction
More informationREDUCED BASIS METHOD
REDUCED BASIS METHOD approximation of PDE s, interpolation, a posteriori estimates Yvon Maday Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie, Paris, France, Institut Universitaire de
More informationThe Generalized Empirical Interpolation Method: Stability theory on Hilbert spaces with an application to the Stokes equation
The Generalized Empirical Interpolation Method: Stability theory on Hilbert spaces with an application to the Stokes equation The MIT Faculty has made this article openly available. Please share how this
More informationrecent developments of approximation theory and greedy algorithms
recent developments of approximation theory and greedy algorithms Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Reduced Order Modeling in
More informationGreedy control. Martin Lazar University of Dubrovnik. Opatija, th Najman Conference. Joint work with E: Zuazua, UAM, Madrid
Greedy control Martin Lazar University of Dubrovnik Opatija, 2015 4th Najman Conference Joint work with E: Zuazua, UAM, Madrid Outline Parametric dependent systems Reduced basis methods Greedy control
More informationEmpirical Interpolation of Nonlinear Parametrized Evolution Operators
MÜNSTER Empirical Interpolation of Nonlinear Parametrized Evolution Operators Martin Drohmann, Bernard Haasdonk, Mario Ohlberger 03/12/2010 MÜNSTER 2/20 > Motivation: Reduced Basis Method RB Scenario:
More informationReduced Modeling in Data Assimilation
Reduced Modeling in Data Assimilation Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Challenges in high dimensional analysis and computation
More informationMATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 05.2013 February 2013 (NEW 10.2013 A weighted empirical
More informationGreedy Control. Enrique Zuazua 1 2
Greedy Control Enrique Zuazua 1 2 DeustoTech - Bilbao, Basque Country, Spain Universidad Autónoma de Madrid, Spain Visiting Fellow of LJLL-UPMC, Paris enrique.zuazua@deusto.es http://enzuazua.net X ENAMA,
More informationModel reduction for multiscale problems
Model reduction for multiscale problems Mario Ohlberger Dec. 12-16, 2011 RICAM, Linz , > Outline Motivation: Multi-Scale and Multi-Physics Problems Model Reduction: The Reduced Basis Approach A new Reduced
More informationPOD-DEIM APPROACH ON DIMENSION REDUCTION OF A MULTI-SPECIES HOST-PARASITOID SYSTEM
POD-DEIM APPROACH ON DIMENSION REDUCTION OF A MULTI-SPECIES HOST-PARASITOID SYSTEM G. Dimitriu Razvan Stefanescu Ionel M. Navon Abstract In this study, we implement the DEIM algorithm (Discrete Empirical
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationContinuous Functions on Metric Spaces
Continuous Functions on Metric Spaces Math 201A, Fall 2016 1 Continuous functions Definition 1. Let (X, d X ) and (Y, d Y ) be metric spaces. A function f : X Y is continuous at a X if for every ɛ > 0
More informationUNITEXT La Matematica per il 3+2
UNITEXT La Matematica per il 3+2 Volume 92 Editor-in-chief A. Quarteroni Series editors L. Ambrosio P. Biscari C. Ciliberto M. Ledoux W.J. Runggaldier More information about this series at http://www.springer.com/series/5418
More informationJOINT SOURCE AND SENSOR PLACEMENT FOR SOUND FIELD CONTROL BASED ON EMPIRICAL INTERPOLATION METHOD Hongo, Bunkyo-ku, Tokyo , Japan
JOINT SOURCE AND SENSOR PLACEMENT FOR SOUND FIELD CONTROL BASED ON EMPIRICAL INTERPOLATION METHOD Shoichi Koyama 1,2, Gilles Chardon 3, and Laurent Daudet 2 1 The University of Tokyo, Graduate School of
More informationReduced-Order Greedy Controllability of Finite Dimensional Linear Systems. Giulia Fabrini Laura Iapichino Stefan Volkwein
Universität Konstanz Reduced-Order Greedy Controllability of Finite Dimensional Linear Systems Giulia Fabrini Laura Iapichino Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 364, Oktober 2017 ISSN
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationNotions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy
Banach Spaces These notes provide an introduction to Banach spaces, which are complete normed vector spaces. For the purposes of these notes, all vector spaces are assumed to be over the real numbers.
More informationPeng Chen 1, Alfio Quarteroni 1,2 and Gianluigi Rozza 3
ESAI: 2AN 48 (214 943 953 DOI: 1.151/m2an/213128 ESAI: athematical odelling and Numerical Analysis www.esaim-m2an.org A WEIGHTED EPIRICAL INTERPOLATION ETHOD: APRIORI CONVERGENCE ANALYSIS AND APPLICATIONS
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationList of Symbols, Notations and Data
List of Symbols, Notations and Data, : Binomial distribution with trials and success probability ; 1,2, and 0, 1, : Uniform distribution on the interval,,, : Normal distribution with mean and variance,,,
More informationProblem Set 5: Solutions Math 201A: Fall 2016
Problem Set 5: s Math 21A: Fall 216 Problem 1. Define f : [1, ) [1, ) by f(x) = x + 1/x. Show that f(x) f(y) < x y for all x, y [1, ) with x y, but f has no fixed point. Why doesn t this example contradict
More informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationl(y j ) = 0 for all y j (1)
Problem 1. The closed linear span of a subset {y j } of a normed vector space is defined as the intersection of all closed subspaces containing all y j and thus the smallest such subspace. 1 Show that
More informationWavelets and modular inequalities in variable L p spaces
Wavelets and modular inequalities in variable L p spaces Mitsuo Izuki July 14, 2007 Abstract The aim of this paper is to characterize variable L p spaces L p( ) (R n ) using wavelets with proper smoothness
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationPROBLEMS. (b) (Polarization Identity) Show that in any inner product space
1 Professor Carl Cowen Math 54600 Fall 09 PROBLEMS 1. (Geometry in Inner Product Spaces) (a) (Parallelogram Law) Show that in any inner product space x + y 2 + x y 2 = 2( x 2 + y 2 ). (b) (Polarization
More informationREDUCED ORDER METHODS
REDUCED ORDER METHODS Elisa SCHENONE Legato Team Computational Mechanics Elisa SCHENONE RUES Seminar April, 6th 5 team at and of Numerical simulation of biological flows Respiration modelling Blood flow
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationThe Reduced Basis Method:
The Reduced Basis Method: Application to Problems in Fluid and Solid Mechanics ++ Aachen Institute for Advanced Study in Computational Engineering Science K. Veroy-Grepl 31 January 2014 AICES & MOR I The
More informationA NOTE ON LINEAR FUNCTIONAL NORMS
A NOTE ON LINEAR FUNCTIONAL NORMS YIFEI PAN AND MEI WANG Abstract. For a vector u in a normed linear space, Hahn-Banach Theorem provides the existence of a linear functional f, f(u) = u such that f = 1.
More informationMath 550 Notes. Chapter 2. Jesse Crawford. Department of Mathematics Tarleton State University. Fall 2010
Math 550 Notes Chapter 2 Jesse Crawford Department of Mathematics Tarleton State University Fall 2010 (Tarleton State University) Math 550 Chapter 2 Fall 2010 1 / 20 Linear algebra deals with finite dimensional
More informationA-Posteriori Error Estimation for Parameterized Kernel-Based Systems
A-Posteriori Error Estimation for Parameterized Kernel-Based Systems Daniel Wirtz Bernard Haasdonk Institute for Applied Analysis and Numerical Simulation University of Stuttgart 70569 Stuttgart, Pfaffenwaldring
More informationConvergence Rates for Greedy Algorithms in Reduced Basis Methods
Convergence Rates for Greedy Algorithms in Reduced Basis Methods Peter Binev, Albert Cohen, Wolfgang Dahmen, Ronald DeVore, Guergana Petrova, and Przemyslaw Wojtaszczyk October 4, 2010 Abstract The reduced
More informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationPBDW: a non-intrusive Reduced Basis Data Assimilation Method and its application to outdoor Air Quality Models
PBDW: a non-intrusive Reduced Basis Data Assimilation Method and its application to outdoor Air Quality Models Janelle K Hammond, Rachida Chakir, Frédéric Bourquin, Yvon Maday To cite this version: Janelle
More informationReal Analysis Notes. Thomas Goller
Real Analysis Notes Thomas Goller September 4, 2011 Contents 1 Abstract Measure Spaces 2 1.1 Basic Definitions........................... 2 1.2 Measurable Functions........................ 2 1.3 Integration..............................
More informationProblem Set 6: Solutions Math 201A: Fall a n x n,
Problem Set 6: Solutions Math 201A: Fall 2016 Problem 1. Is (x n ) n=0 a Schauder basis of C([0, 1])? No. If f(x) = a n x n, n=0 where the series converges uniformly on [0, 1], then f has a power series
More informationExistence Theorem for Abstract Measure. Differential Equations Involving. the Distributional Henstock-Kurzweil Integral
Journal of Applied Mathematics & Bioinformatics, vol.4, no.1, 2014, 11-20 ISSN: 1792-6602 (print), 1792-6939 (online) Scienpress Ltd, 2014 Existence Theorem for Abstract Measure Differential Equations
More informationOutline of Fourier Series: Math 201B
Outline of Fourier Series: Math 201B February 24, 2011 1 Functions and convolutions 1.1 Periodic functions Periodic functions. Let = R/(2πZ) denote the circle, or onedimensional torus. A function f : C
More informationRolle s Theorem for Polynomials of Degree Four in a Hilbert Space 1
Journal of Mathematical Analysis and Applications 265, 322 33 (2002) doi:0.006/jmaa.200.7708, available online at http://www.idealibrary.com on Rolle s Theorem for Polynomials of Degree Four in a Hilbert
More informationA general, multipurpose interpolation procedure: the magic points
A general, multipurpose interpolation procedure: the magic points Yvon aday, Ngoc Cuong Nguyen, Anthony T. Patera, George S.H. Pau To cite this version: Yvon aday, Ngoc Cuong Nguyen, Anthony T. Patera,
More informationFREMLIN TENSOR PRODUCTS OF CONCAVIFICATIONS OF BANACH LATTICES
FREMLIN TENSOR PRODUCTS OF CONCAVIFICATIONS OF BANACH LATTICES VLADIMIR G. TROITSKY AND OMID ZABETI Abstract. Suppose that E is a uniformly complete vector lattice and p 1,..., p n are positive reals.
More information1 Definition of the Riemann integral
MAT337H1, Introduction to Real Analysis: notes on Riemann integration 1 Definition of the Riemann integral Definition 1.1. Let [a, b] R be a closed interval. A partition P of [a, b] is a finite set of
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Methods for Nonlinear Systems Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu 1 / 65 Outline 1 Nested Approximations 2 Trajectory PieceWise Linear (TPWL)
More informationFUNCTIONAL ANALYSIS-NORMED SPACE
MAT641- MSC Mathematics, MNIT Jaipur FUNCTIONAL ANALYSIS-NORMED SPACE DR. RITU AGARWAL MALAVIYA NATIONAL INSTITUTE OF TECHNOLOGY JAIPUR 1. Normed space Norm generalizes the concept of length in an arbitrary
More informationAnother Riesz Representation Theorem
Another Riesz Representation Theorem In these notes we prove (one version of) a theorem known as the Riesz Representation Theorem. Some people also call it the Riesz Markov Theorem. It expresses positive
More information2.3 Variational form of boundary value problems
2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationBounding Stability Constants for Affinely Parameter-Dependent Operators
Bounding Stability Constants for Affinely Parameter-Dependent Operators Robert O Connor a a RWTH Aachen University, Aachen, Germany Received *****; accepted after revision +++++ Presented by Abstract In
More informationarxiv: v2 [math.na] 15 May 2014
A nonintrusive Reduced Basis Method applied to aeroacoustic simulations arxiv:1401.3684v2 [math.na] 15 May 2014 March 26, 2018 Fabien Casenave 1, Alexandre Ern 1, and Tony Lelièvre 1,2 1 Université Paris-Est,
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationRank Determination for Low-Rank Data Completion
Journal of Machine Learning Research 18 017) 1-9 Submitted 7/17; Revised 8/17; Published 9/17 Rank Determination for Low-Rank Data Completion Morteza Ashraphijuo Columbia University New York, NY 1007,
More informationApprentice Linear Algebra, 1st day, 6/27/05
Apprentice Linear Algebra, 1st day, 6/7/05 REU 005 Instructor: László Babai Scribe: Eric Patterson Definitions 1.1. An abelian group is a set G with the following properties: (i) ( a, b G)(!a + b G) (ii)
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationSurrogate gravitational waveform models
Department of Physics, University of Maryland Theoretical astrophysics seminar Cornell University, November 30, 2013 Previous works - Field, Galley, Herrmann, Hesthaven, Ochsner, Tiglio (Phys. Rev. Lett.
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationNonlinear elliptic systems with exponential nonlinearities
22-Fez conference on Partial Differential Equations, Electronic Journal of Differential Equations, Conference 9, 22, pp 139 147. http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationApproximation of Parametric Derivatives by the Empirical Interpolation Method
Approximation of Parametric Derivatives by the Empirical Interpolation ethod Jens L. Eftang, artin A. Grepl, Anthony T. Patera and Einar. Rønquist Bericht Nr. 330 September 2011 Key words: function approximation;
More informationTWO MAPPINGS RELATED TO SEMI-INNER PRODUCTS AND THEIR APPLICATIONS IN GEOMETRY OF NORMED LINEAR SPACES. S.S. Dragomir and J.J.
RGMIA Research Report Collection, Vol. 2, No. 1, 1999 http://sci.vu.edu.au/ rgmia TWO MAPPINGS RELATED TO SEMI-INNER PRODUCTS AND THEIR APPLICATIONS IN GEOMETRY OF NORMED LINEAR SPACES S.S. Dragomir and
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationMath 328 Course Notes
Math 328 Course Notes Ian Robertson March 3, 2006 3 Properties of C[0, 1]: Sup-norm and Completeness In this chapter we are going to examine the vector space of all continuous functions defined on the
More informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More informationPhilippe. Functional Analysis with Applications. Linear and Nonlinear. G. Ciarlet. City University of Hong Kong. siajtl
Philippe G Ciarlet City University of Hong Kong Linear and Nonlinear Functional Analysis with Applications with 401 Problems and 52 Figures siajtl Society for Industrial and Applied Mathematics Philadelphia
More informationRegularity and compactness for the DiPerna Lions flow
Regularity and compactness for the DiPerna Lions flow Gianluca Crippa 1 and Camillo De Lellis 2 1 Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy. g.crippa@sns.it 2 Institut für Mathematik,
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationDefinition A.1. We say a set of functions A C(X) separates points if for every x, y X, there is a function f A so f(x) f(y).
The Stone-Weierstrass theorem Throughout this section, X denotes a compact Hausdorff space, for example a compact metric space. In what follows, we take C(X) to denote the algebra of real-valued continuous
More informationSTRONG CONVERGENCE OF AN IMPLICIT ITERATION PROCESS FOR ASYMPTOTICALLY NONEXPANSIVE IN THE INTERMEDIATE SENSE MAPPINGS IN BANACH SPACES
Kragujevac Journal of Mathematics Volume 36 Number 2 (2012), Pages 237 249. STRONG CONVERGENCE OF AN IMPLICIT ITERATION PROCESS FOR ASYMPTOTICALLY NONEXPANSIVE IN THE INTERMEDIATE SENSE MAPPINGS IN BANACH
More informationOn the bang-bang property of time optimal controls for infinite dimensional linear systems
On the bang-bang property of time optimal controls for infinite dimensional linear systems Marius Tucsnak Université de Lorraine Paris, 6 janvier 2012 Notation and problem statement (I) Notation: X (the
More informationGeneralized Monotonicities and Its Applications to the System of General Variational Inequalities
Generalized Monotonicities and Its Applications to the System of General Variational Inequalities Khushbu 1, Zubair Khan 2 Research Scholar, Department of Mathematics, Integral University, Lucknow, Uttar
More informationGAKUTO International Series
GAKUTO International Series Mathematical Sciences and Applications, Vol.28(2008) Proceedings of Fourth JSIAM-SIMMAI Seminar on Industrial and Applied Mathematics, pp.139-148 A COMPUTATIONAL APPROACH TO
More informationSummary of Real Analysis by Royden
Summary of Real Analysis by Royden Dan Hathaway May 2010 This document is a summary of the theorems and definitions and theorems from Part 1 of the book Real Analysis by Royden. In some areas, such as
More informationBases in Banach spaces
Chapter 3 Bases in Banach spaces Like every vector space a Banach space X admits an algebraic or Hamel basis, i.e. a subset B X, so that every x X is in a unique way the (finite) linear combination of
More informationLecture 5: Peak sets for operator algebras and quantum cardinals. David Blecher. December 2016
Lecture 5: Peak sets for operator algebras and quantum cardinals David Blecher University of Houston December 2016 Abstract for Lecture 5 Abstract: We discuss noncommutative peak sets for operator algebras,
More informationNote on the Chen-Lin Result with the Li-Zhang Method
J. Math. Sci. Univ. Tokyo 18 (2011), 429 439. Note on the Chen-Lin Result with the Li-Zhang Method By Samy Skander Bahoura Abstract. We give a new proof of the Chen-Lin result with the method of moving
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationProgress in Several Complex Variables KIAS 2018
for Progress in Several Complex Variables KIAS 2018 Outline 1 for 2 3 super-potentials for 4 real for Let X be a real manifold of dimension n. Let 0 p n and k R +. D c := { C k (differential) (n p)-forms
More informationReminder Notes for the Course on Measures on Topological Spaces
Reminder Notes for the Course on Measures on Topological Spaces T. C. Dorlas Dublin Institute for Advanced Studies School of Theoretical Physics 10 Burlington Road, Dublin 4, Ireland. Email: dorlas@stp.dias.ie
More informationBEST APPROXIMATIONS AND ORTHOGONALITIES IN 2k-INNER PRODUCT SPACES. Seong Sik Kim* and Mircea Crâşmăreanu. 1. Introduction
Bull Korean Math Soc 43 (2006), No 2, pp 377 387 BEST APPROXIMATIONS AND ORTHOGONALITIES IN -INNER PRODUCT SPACES Seong Sik Kim* and Mircea Crâşmăreanu Abstract In this paper, some characterizations of
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationFinal Year M.Sc., Degree Examinations
QP CODE 569 Page No Final Year MSc, Degree Examinations September / October 5 (Directorate of Distance Education) MATHEMATICS Paper PM 5: DPB 5: COMPLEX ANALYSIS Time: 3hrs] [Max Marks: 7/8 Instructions
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationOn the simplest expression of the perturbed Moore Penrose metric generalized inverse
Annals of the University of Bucharest (mathematical series) 4 (LXII) (2013), 433 446 On the simplest expression of the perturbed Moore Penrose metric generalized inverse Jianbing Cao and Yifeng Xue Communicated
More informationA functional model for commuting pairs of contractions and the symmetrized bidisc
A functional model for commuting pairs of contractions and the symmetrized bidisc Nicholas Young Leeds and Newcastle Universities Lecture 2 The symmetrized bidisc Γ and Γ-contractions St Petersburg, June
More informationOptimal Polynomial Admissible Meshes on the Closure of C 1,1 Bounded Domains
Optimal Polynomial Admissible Meshes on the Closure of C 1,1 Bounded Domains Constructive Theory of Functions Sozopol, June 9-15, 2013 F. Piazzon, joint work with M. Vianello Department of Mathematics.
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. Research Scholarships Screening Test. Saturday, February 2, Time Allowed: Two Hours Maximum Marks: 40
NATIONAL BOARD FOR HIGHER MATHEMATICS Research Scholarships Screening Test Saturday, February 2, 2008 Time Allowed: Two Hours Maximum Marks: 40 Please read, carefully, the instructions on the following
More informationIf Y and Y 0 satisfy (1-2), then Y = Y 0 a.s.
20 6. CONDITIONAL EXPECTATION Having discussed at length the limit theory for sums of independent random variables we will now move on to deal with dependent random variables. An important tool in this
More informationApplication of POD-DEIM Approach for Dimension Reduction of a Diffusive Predator-Prey System with Allee Effect
Application of POD-DEIM Approach for Dimension Reduction of a Diffusive Predator-Prey System with Allee Effect Gabriel Dimitriu 1, Ionel M. Navon 2 and Răzvan Ştefănescu 2 1 The Grigore T. Popa University
More informationTHE STONE-WEIERSTRASS THEOREM AND ITS APPLICATIONS TO L 2 SPACES
THE STONE-WEIERSTRASS THEOREM AND ITS APPLICATIONS TO L 2 SPACES PHILIP GADDY Abstract. Throughout the course of this paper, we will first prove the Stone- Weierstrass Theroem, after providing some initial
More informationThe Way of Analysis. Robert S. Strichartz. Jones and Bartlett Publishers. Mathematics Department Cornell University Ithaca, New York
The Way of Analysis Robert S. Strichartz Mathematics Department Cornell University Ithaca, New York Jones and Bartlett Publishers Boston London Contents Preface xiii 1 Preliminaries 1 1.1 The Logic of
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationModel reduction of parametrized aerodynamic flows: discontinuous Galerkin RB empirical quadrature procedure. Masayuki Yano
Model reduction of parametrized aerodynamic flows: discontinuous Galerkin RB empirical quadrature procedure Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony
More information