Accurate and online-efficient evaluation of the a posteriori error bound in the reduced basis method
ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Volume 48 (2014) no. 1, p. 207-229

The reduced basis method is a model reduction technique yielding substantial savings of computational time when a solution to a parametrized equation has to be computed for many values of the parameter. Certification of the approximation is possible by means of an a posteriori error bound. Under appropriate assumptions, this error bound is computed with an algorithm of complexity independent of the size of the full problem. In practice, the evaluation of the error bound can become very sensitive to round-off errors. We propose herein an explanation of this fact. A first remedy has been proposed in [F. Casenave, Accurate a posteriori error evaluation in the reduced basis method. C. R. Math. Acad. Sci. Paris 350 (2012) 539-542.]. Herein, we improve this remedy by proposing a new approximation of the error bound using the empirical interpolation method (EIM). This method achieves higher levels of accuracy and requires potentially less precomputations than the usual formula. A version of the EIM stabilized with respect to round-off errors is also derived. The method is illustrated on a simple one-dimensional diffusion problem and a three-dimensional acoustic scattering problem solved by a boundary element method.

Classification:  65N15,  65D05,  68W25,  76Q05
Keywords: reduced basis method, a posteriori error bound, round-off errors, boundary element method, empirical interpolation method, acoustics
     author = {Casenave, Fabien and Ern, Alexandre and Leli\`evre, Tony},
     title = {Accurate and online-efficient evaluation of the a posteriori error bound in the reduced basis method},
     journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
     publisher = {EDP-Sciences},
     volume = {48},
     number = {1},
     year = {2014},
     pages = {207-229},
     doi = {10.1051/m2an/2013097},
     zbl = {1288.65157},
     language = {en},
     url = {}
Casenave, Fabien; Ern, Alexandre; Lelièvre, Tony. Accurate and online-efficient evaluation of the a posteriori error bound in the reduced basis method. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Volume 48 (2014) no. 1, pp. 207-229. doi : 10.1051/m2an/2013097.

[1] Z. Bai and D. Skoogh, Krylov subspace techniques for reduced-order modeling of large-scale dynamical systems. Appl. Numer. Math. 43 (2002) 9-44. | MR 1936100 | Zbl 1012.65136

[2] M.A. Bahayou, Sur le problème de Helmholtz. Rendiconti del Seminario matematico della Università e Politecnico di Torino (2007) 427-450. | MR 2402854 | Zbl 1187.35028

[3] M. Barrault, Y. Maday, N.C. Nguyen and A.T. Patera, An ‘empirical interpolation' method: application to efficient reduced-basis discretization of partial differential equations. C. R. Math. Acad. Sci. Paris 339 (2004) 667-672. | MR 2103208 | Zbl 1061.65118

[4] A. Björck and C.C. Paige, Loss and recapture of orthogonality in the modified Gram-Schmidt algorithm. SIAM J. Matrix Anal. Appl. 13 (1992) 176-190. | MR 1146660 | Zbl 0747.65026

[5] S. Boyaval, Mathematical modelling and numerical simulation in materials science. Ph.D. thesis, Université Paris-Est (2009).

[6] A. Buffa and R. Hiptmair, Regularized combined field integral equations. Numer. Math. 100 (2005) 1-19. | MR 2129699 | Zbl 1067.65137

[7] R.L. Burden and J.D. Faires, Numerical Analysis. PWS Publishing Company (1993). | Zbl 0788.65001

[8] E. Cancès, V. Ehrlacher and T. Lelièvre, Convergence of a greedy algorithm for high-dimensional convex nonlinear problems. Math. Models Methods Appl. Sci. 21 (2011) 2433-2467. | MR 2864637 | Zbl 1259.65098

[9] F. Casenave, Accurate a posteriori error evaluation in the reduced basis method. C. R. Math. Acad. Sci. Paris 350 (2012) 539-542. | MR 2929064 | Zbl 1245.65105

[10] F. Casenave, Ph.D. thesis, in preparation (2013).

[11] F. Casenave, M. Ghattassi and R. Joubaud, A multiscale problem in thermal science. ESAIM: Proceedings 38 (2012) 202-219.

[12] A. Chatterjee, An introduction to the proper orthogonal decomposition. Curr. Sci. 78 (2000) 808-817.

[13] Y. Chen, J.S. Hesthaven, Y. Maday, J. Rodriguez and X. Zhu, Certified reduced basis method for electromagnetic scattering and radar cross section estimation. Technical Report 2011-28, Scientific Computing Group, Brown University, Providence, RI, USA (2011). | MR 2924023 | Zbl 1253.78045

[14] Y. Chen, J.S. Hesthaven, Y. Maday and J. Rodríguez, Improved successive constraint method based a posteriori error estimate for reduced basis approximation of 2D Maxwell's problem. ESAIM: M2AN 43 (2009) 1099-1116. | Numdam | MR 2588434 | Zbl 1181.78019

[15] F. Chinesta, P. Ladeveze and C. Elías, A short review on model order reduction based on proper generalized decomposition. Arch. Comput. Methods Eng. 18 (2011) 395-404.

[16] A. Delnevo, I. Terrasse, Code Acti3S harmonique : Justifications Mathématiques : Partie I. Technical report, EADS CCR (2001).

[17] A. Delnevo, I. Terrasse, Code Acti3S, Justifications Mathématiques : Partie II, présence d'un écoulement uniforme. Technical report, EADS CCR (2002).

[18] A. Ern and J.L. Guermond, Theory and Practice of Finite Elements, in vol. 159 of Applied Mathematical Sciences. Springer (2004). | MR 2050138 | Zbl 1059.65103

[19] M. Fares, J.S. Hesthaven, Y. Maday and B. Stamm, The reduced basis method for the electric field integral equation. J. Comput. Phys. 230 (2011) 5532-5555. | MR 2799523 | Zbl 1220.78045

[20] L. Giraud and J. Langou, When modified Gram-Schmidt generates a well-conditioned set of vectors. IMA J. Numer. Anal. 22 (2002) 521-528. | MR 1936517 | Zbl 1027.65050

[21] D. Goldberg, What every computer scientist should know about floating point arithmetic. ACM Computing Surveys 23 (1991) 5-48.

[22] G.H. Golub and C.F. Van Loan, Matrix Computations. Johns Hopkins Studies in the Mathematical Sciences. Johns Hopkins University Press (1996). | MR 1417720 | Zbl 1268.65037

[23] R.J. Guyan, Reduction of stiffness and mass matrices. AIAA J. 3 (1965) 380.

[24] R. Hiptmair, Coercive combined field integral equations. J. Numer. Math. 11 (2003) 115-134. | MR 1987591 | Zbl 1115.76356

[25] R. Hiptmair and P. Meury, Stable FEM-BEM Coupling for Helmholtz Transmission Problems. ETH, Seminar für Angewandte Mathematik (2005). | MR 2263042 | Zbl 1221.65308

[26] G.C. Hsiao and W.L. Wendland, Boundary Element Methods: Foundation and Error Analysis. John Wiley & Sons, Ltd (2004).

[27] D.B.P. Huynh, G. Rozza, S. Sen and A.T. Patera, A successive constraint linear optimization method for lower bounds of parametric coercivity and inf-sup stability constants. C. R. Math. Acad. Sci. Paris 345 (2007) 473-478. | MR 2367928 | Zbl 1127.65086

[28] P. Langlois, S. Graillat and N. Louvet, Compensated Horner scheme. Schloss Dagstuhl - Leibniz-Zentrum für Informatik (2006).

[29] L. Machiels, Y. Maday, I.B. Oliveira, A.T. Patera and D.V. Rovas, Output bounds for reduced-basis approximations of symmetric positive definite eigenvalue problems. C. R. Math. Acad. Sci. Paris 331 (2000) 153-158. | MR 1781533 | Zbl 0960.65063

[30] Y. Maday, N.C. Nguyen, A.T. Patera and S. Pau, A general multipurpose interpolation procedure: the magic points. Commun. Pure Appl. Anal. 8 (2008) 383-404. | Zbl 1184.65020

[31] W.C.H. Mclean, Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press (2000). | MR 1742312 | Zbl 0948.35001

[32] A. Nouy and O.P. Le Maître, Generalized spectral decomposition for stochastic nonlinear problems. J. Comput. Phys. 228 (2009) 202-235. | MR 2464076 | Zbl 1157.65009

[33] A.T. Patera, Private communication (2012).

[34] A.T. Patera and G. Rozza, Reduced Basis Approximation and A Posteriori Error Estimation for Parametrized Partial Differential Equations. MIT Pappalardo Graduate Monographs in Mechanical Engineering (2007). | Zbl pre05344486

[35] M. Paz, Dynamic condensation. AIAA J. 22 (1984) 724-727.

[36] C. Prud'Homme, D.V. Rovas, K. Veroy, L. Machiels, Y. Maday, A.T. Patera and G. Turinici, Reliable real-time solution of parametrized partial differential equations: Reduced-basis output bound methods. J. Fluids Eng. 124 (2002) 70-80.

[37] S.A. Sauter and C. Schwab, Boundary Element Methods. Springer Series in Computational Mathematics. Springer (2010). | MR 2743235 | Zbl 1215.65183

[38] I.E. Shparlinski, Sparse polynomial approximation in finite fields. In Proceedings of the thirty-third annual ACM symposium on Theory of computing, STOC '01. ACM, New York, USA (2001) 209-215. | MR 2120317

[39] K. Veroy and A.T. Patera, Certified real-time solution of the parametrized steady incompressible Navier-Stokes equations: rigorous reduced-basis a posteriori error bounds. Int. J. Numer. Methods Fluids 47 (2005) 773-788. | MR 2123791 | Zbl 1134.76326

[40] K. Veroy, C. Prud'Homme and A.T. Patera, Reduced-basis approximation of the viscous Burgers equation: rigorous a posteriori error bounds. C. R. Math. Acad. Sci. Paris 337 (2003) 619-624. | MR 2017737 | Zbl 1036.65075