# Accurate and online-efficient evaluation of the a posteriori error bound in the reduced basis method

• Volume: 48, Issue: 1, page 207-229
• ISSN: 0764-583X

top

## Abstract

top
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.

## How to cite

top

Casenave, Fabien, Ern, Alexandre, and 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 48.1 (2014): 207-229. <http://eudml.org/doc/273118>.

@article{Casenave2014,
abstract = {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.},
author = {Casenave, Fabien, Ern, Alexandre, Lelièvre, Tony},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique},
keywords = {reduced basis method; a posteriori error bound; round-off errors; boundary element method; empirical interpolation method; acoustics; Helmholtz equation; numerical examples; diffusion problem; acoustic scattering},
language = {eng},
number = {1},
pages = {207-229},
publisher = {EDP-Sciences},
title = {Accurate and online-efficient evaluation of the a posteriori error bound in the reduced basis method},
url = {http://eudml.org/doc/273118},
volume = {48},
year = {2014},
}

TY - JOUR
AU - Casenave, Fabien
AU - Ern, Alexandre
AU - Lelièvre, Tony
TI - Accurate and online-efficient evaluation of the a posteriori error bound in the reduced basis method
JO - ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique
PY - 2014
PB - EDP-Sciences
VL - 48
IS - 1
SP - 207
EP - 229
AB - 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.
LA - eng
KW - reduced basis method; a posteriori error bound; round-off errors; boundary element method; empirical interpolation method; acoustics; Helmholtz equation; numerical examples; diffusion problem; acoustic scattering
UR - http://eudml.org/doc/273118
ER -

## References

top
1. [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. Zbl1012.65136MR1936100
2. [2] M.A. Bahayou, Sur le problème de Helmholtz. Rendiconti del Seminario matematico della Università e Politecnico di Torino (2007) 427–450. Zbl1187.35028MR2402854
3. [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. Paris339 (2004) 667–672. Zbl1061.65118MR2103208
4. [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. Zbl0747.65026MR1146660
5. [5] S. Boyaval, Mathematical modelling and numerical simulation in materials science. Ph.D. thesis, Université Paris-Est (2009).
6. [6] A. Buffa and R. Hiptmair, Regularized combined field integral equations. Numer. Math.100 (2005) 1–19. Zbl1067.65137MR2129699
7. [7] R.L. Burden and J.D. Faires, Numerical Analysis. PWS Publishing Company (1993). Zbl0788.65001
8. [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. Zbl1259.65098MR2864637
9. [9] F. Casenave, Accurate a posteriori error evaluation in the reduced basis method. C. R. Math. Acad. Sci. Paris350 (2012) 539–542. Zbl1245.65105MR2929064
10. [10] F. Casenave, Ph.D. thesis, in preparation (2013).
11. [11] F. Casenave, M. Ghattassi and R. Joubaud, A multiscale problem in thermal science. ESAIM: Proceedings 38 (2012) 202–219. Zbl1334.80003
12. [12] A. Chatterjee, An introduction to the proper orthogonal decomposition. Curr. Sci.78 (2000) 808–817.
13. [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). Zbl1253.78045MR2924023
14. [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. Zbl1181.78019MR2588434
15. [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. [16] A. Delnevo, I. Terrasse, Code ACTI3S harmonique : Justifications Mathématiques : Partie I. Technical report, EADS CCR (2001).
17. [17] A. Delnevo, I. Terrasse, Code ACTI3S, Justifications Mathématiques : Partie II, présence d’un écoulement uniforme. Technical report, EADS CCR (2002).
18. [18] A. Ern and J.L. Guermond, Theory and Practice of Finite Elements, in vol. 159 of Applied Mathematical Sciences. Springer (2004). Zbl1059.65103MR2050138
19. [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. Zbl1220.78045MR2799523
20. [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. Zbl1027.65050MR1936517
21. [21] D. Goldberg, What every computer scientist should know about floating point arithmetic. ACM Computing Surveys23 (1991) 5–48.
22. [22] G.H. Golub and C.F. Van Loan, Matrix Computations. Johns Hopkins Studies in the Mathematical Sciences. Johns Hopkins University Press (1996). Zbl1268.65037MR1417720
23. [23] R.J. Guyan, Reduction of stiffness and mass matrices. AIAA J. 3 (1965) 380.
24. [24] R. Hiptmair, Coercive combined field integral equations. J. Numer. Math.11 (2003) 115–134. Zbl1115.76356MR1987591
25. [25] R. Hiptmair and P. Meury, Stable FEM-BEM Coupling for Helmholtz Transmission Problems. ETH, Seminar für Angewandte Mathematik (2005). Zbl1221.65308MR2263042
26. [26] G.C. Hsiao and W.L. Wendland, Boundary Element Methods: Foundation and Error Analysis. John Wiley & Sons, Ltd (2004).
27. [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. Paris345 (2007) 473–478. Zbl1127.65086MR2367928
28. [28] P. Langlois, S. Graillat and N. Louvet, Compensated Horner scheme. Schloss Dagstuhl – Leibniz-Zentrum für Informatik (2006).
29. [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. Paris331 (2000) 153–158. Zbl0960.65063MR1781533
30. [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. Zbl1184.65020
31. [31] W.C.H. McLean, Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press (2000). Zbl0948.35001MR1742312
32. [32] A. Nouy and O.P. Le Maître, Generalized spectral decomposition for stochastic nonlinear problems. J. Comput. Phys.228 (2009) 202–235. Zbl1157.65009MR2464076
33. [33] A.T. Patera, Private communication (2012).
34. [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). Zbl1304.65251
35. [35] M. Paz, Dynamic condensation. AIAA J.22 (1984) 724–727.
36. [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. [37] S.A. Sauter and C. Schwab, Boundary Element Methods. Springer Series in Computational Mathematics. Springer (2010). Zbl1215.65183MR2743235
38. [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. Zbl1323.68321MR2120317
39. [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 Fluids47 (2005) 773–788. Zbl1134.76326MR2123791
40. [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. Paris337 (2003) 619–624. Zbl1036.65075MR2017737

## NotesEmbed?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.