Skip to main content
×
Home
    • Aa
    • Aa

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)
Abstract
Abstract

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.

Copyright
Corresponding author
*Corresponding author. Email addresses: salvatore.cuomo@unina.it (S. Cuomo), ardelio.galletti@uniparthenope.it (A. Galletti), giulio.giunta@uniparthenope.it (G. Giunta), livia.marcellino@uniparthenope.it (L. Marcellino)
Linked references
Hide All

This list contains references from the content that can be linked to their source. For a full set of references and notes please see the PDF or HTML where available.

[2] I. T. Young and L.J. Van Vliet , Recursive implementation of the Gaussian filter, Signal Process., 44 (1995), pp. 139151.

[5] R. J. Purser , W. S. Wu , D. F. Parish and N. M. Roberts , 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] A. T. Weaver and P. Courtier , Correlation modelling on the sphere using a generalized diffusion equation, Q. J. Roy. Meteor. Soc., 127 (2001), pp. 18151846.

[9] S. Cuomo , G. De Pietro , R. Farina , A. Galletti and G. Sannino , A novel numerical scheme for ECG signal denoising, Procedia Computer Science, 51 (2015), pp. 775784.

[11] S. Cuomo , G. De Pietro , R. Farina , A. Galletti and G. Sannino , A revised scheme for real time ECG Signal denoising based on recursive filtering, Biomed. Signal Proces., (2016), pp. 134144.

[12] S. Dobricic and N. Pinardi , An oceanographic three-dimensional variational data assimilation scheme, Ocean Model., 22 (2008), pp. 89105.

[17] S. Haben , A. Lawless and N. Nicholos , Conditioning and preconditioning of the variational data assimilation problem, Computers and Fluids, 46 (2011), pp. 252256.

[20] S. Cuomo , R. Farina , A. Galletti and L. Marcellino , 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.

[22] C. Hayden and R. J. Purser , Recursive filter objective analysis of meteorological field: applications to NESDIS operational processing, J. Appl. Meteorol., 34 (1995), pp. 315.

[23] N. J. Higham , Efficient Algorithms for Computing the Condition Number of a Tridiagonal Matrix, SIAM J. Sci. Stat. Comput., 7(1) (1986), pp. 150165.

[27] B. Triggs and M. Sdika , Boundary conditions for Young-van Vliet recursive filtering, IEEE T. Signal Proces., 54(6) (2006), pp. 23652367.

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? *
×

Keywords:

Metrics

Full text views

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

Abstract views

Total abstract views: 402 *
Loading metrics...

* Views captured on Cambridge Core between 20th June 2017 - 26th July 2017. This data will be updated every 24 hours.