Skip to main content

Numerical Effects of the Gaussian Recursive Filters in Solving Linear Systems in the 3Dvar Case Study

  • Salvatore Cuomo (a1), Ardelio Galletti (a2), Giulio Giunta (a2) and Livia Marcellino (a2)

In many applications, the Gaussian convolution is approximately computed by means of recursive filters, with a significant improvement of computational efficiency. We are interested in theoretical and numerical issues related to such an use of recursive filters in a three-dimensional variational data assimilation (3Dvar) scheme as it appears in the software OceanVar. In that context, the main numerical problem consists in solving large linear systems with high efficiency, so that an iterative solver, namely the conjugate gradient method, is equipped with a recursive filter in order to compute matrix-vector multiplications that in fact are Gaussian convolutions. Here we present an error analysis that gives effective bounds for the perturbation on the solution of such linear systems, when is computed by means of recursive filters. We first prove that such a solution can be seen as the exact solution of a perturbed linear system. Then we study the related perturbation on the solution and we demonstrate that it can be bounded in terms of the difference between the two linear operators associated to the Gaussian convolution and the recursive filter, respectively. Moreover, we show through numerical experiments that the error on the solution, which exhibits a kind of edge effect, i.e. most of the error is localized in the first and last few entries of the computed solution, is due to the structure of the difference of the two linear operators.

Corresponding author
*Corresponding author. Email addresses: (S. Cuomo), (A. Galletti), (G. Giunta), (L. Marcellino)
Hide All
[1] Witkin A., Scale-space filtering, Int. Joint Conf. Artif. , Karlsruhe, germany, (1983), pp. 10191022.
[2] Young I. T. and Van Vliet L.J., Recursive implementation of the Gaussian filter, Signal Process., 44 (1995), pp. 139151.
[3] Haglund L., Adaptive multidimensional filtering, LinkÃűuping University, Sweden, 1992.
[4] Lorenc A. C., Iterative analysis using covariance functions and filters, Q. J. Roy. Meteor. Soc., 118(118) (1992), pp. 569591.
[5] Purser R. J., Wu W. S., Parish D. F. and Roberts N. M., Numerical aspects of the application of recursive filters to variational statistical analysis. part II: spatially inhomogeneous and anisotropic general covariances, Mon. Weather Rev., 131(8) (2003), pp. 15361548.
[6] Weaver A. T. and Courtier P., Correlation modelling on the sphere using a generalized diffusion equation, Q. J. Roy. Meteor. Soc., 127 (2001), pp. 18151846.
[7] Deriche R., Separable recursive filtering for efficient multi-scale edge detection, In Proc. Int. Workshop Industrial Application, (1987), pp. 1823.
[8] Vliet L. V., Young I. and Verbeek P., Recursive Gaussian derivative filters, Int. C. Patt. Recog., (1998), pp. 509514.
[9] Cuomo S., De Pietro G., Farina R., Galletti A. and Sannino G., A novel inline-graphic numerical scheme for ECG signal denoising, Procedia Computer Science, 51 (2015), pp. 775784.
[10] Cuomo S., De Pietro G., Farina R., Galletti A. and Sannino G., A framework for ECG denoising for mobile devices, The, ACM International Conference ACM, (2015), pp. 14.
[11] Cuomo S., De Pietro G., Farina R., Galletti A. and Sannino G., A revised scheme for real time ECG Signal denoising based on recursive filtering, Biomed. Signal Proces., (2016), pp. 134144.
[12] Dobricic S. and Pinardi N., An oceanographic three-dimensional variational data assimilation scheme, Ocean Model., 22 (2008), pp. 89105.
[13] Farina R., Dobricic S., Storto A., Masina S. and Cuomo S., A revised scheme to compute horizontal covariances in an oceanographic 3D-VAR assimilation system, J. Comput. Phys., 284 (2015), pp. 631647.
[14] Cuomo S., Farina R., Galletti A. and Marcellino L., An error estimate of Gaussian Recursive Filter in 3Dvar problem, Comput. Sci. Inf. Syst. IEEE, 2 (2014), pp. 587595.
[15] D’ Amore L., Arcucci R., Marcellino L. and Murli A., HPC computation issues of the incremental 3D variational data assimilation scheme in OceanVarsoftware, J. Numer. Anal., Industrial and Applied Mathematics, 7(3-4) (2013), pp. 91105.
[16] Lorenc A.C., Development of an operational variational assimilation scheme, J. Meteorol. Soc. Jpn., 75 (1997), pp. 339346.
[17] Haben S., Lawless A. and Nicholos N., Conditioning and preconditioning of the variational data assimilation problem, Computers and Fluids, 46 (2011), pp. 252256.
[18] Haben S., Lawless A. and Nichols N., Conditioning of the 3Dvar data assimilation problem, University of Reading, Dept. of Mathematics, Math Report Series, 3 (2009).
[19] Cuomo S., Galletti A., Giunta G., Marcellino L., Some error bounds for K-iterated Gaussian recursive filters, AIP Conference Proceedings, 1776 (2016), 040008,
[20] Cuomo S., Farina R., Galletti A. and Marcellino L., A K-iterated scheme for the first-order Gaussian Recursive Filter with boundary conditions, Federated Conference on Computer Science and Information Systems, 5 (2015), pp. 641647.
[21] Galletti A. and Giunta G., Error analysis for the first-order Gaussian recursive filter operator, Proceedings of the 2016 Federated Conference on Computer Science and Information Systems, Ganzha M., Maciaszek L., Paprzycki M. (eds). ACSIS, 8 (2016), pp. 673678,
[22] Hayden C. and Purser R. J., Recursive filter objective analysis of meteorological field: applications to NESDIS operational processing, J. Appl. Meteorol., 34 (1995), pp. 315.
[23] Higham N. J., Efficient Algorithms for Computing the Condition Number of a Tridiagonal Matrix, SIAM J. Sci. Stat. Comput., 7(1) (1986), pp. 150165.
[24] Cuomo S., Galletti A., Giunta G. and Starace A., Surface Reconstruction from Scattered Point via RBF Interpolation on GPU, Comput. Sci. Inf. Syst. IEEE, (2013), pp. 433440.
[25] Cuomo S., Galletti A., Giunta G. and Marcellino L., Reconstruction of implicit curves and surfaces via RBF interpolation, Appl. Numer. Math., (2016), pp. 6063.
[26] Higham N. J., Accuracy and Stability of Numerical Algorithms, Jounal of the American Statistical Association, 16(94) (2002), pp. 285289.
[27] Triggs B. and Sdika M., Boundary conditions for Young-van Vliet recursive filtering, IEEE T. Signal Proces., 54(6) (2006), pp. 23652367.
[28] Farina R., Dobricic S. and Cuomo S., Some numerical enhancements in a data assimilation scheme, AIP Conf. Proc., 1558 (2013), pp.23692372.
Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

Numerical Mathematics: Theory, Methods and Applications
  • ISSN: 1004-8979
  • EISSN: 2079-7338
  • URL: /core/journals/numerical-mathematics-theory-methods-and-applications
Please enter your name
Please enter a valid email address
Who would you like to send this to? *



Full text views

Total number of HTML views: 0
Total number of PDF views: 14 *
Loading metrics...

Abstract views

Total abstract views: 489 *
Loading metrics...

* Views captured on Cambridge Core between 20th June 2017 - 22nd January 2018. This data will be updated every 24 hours.