Skip to main content
×
×
Home

ITERATIVE SOLUTION OF SHIFTED POSITIVE-DEFINITE LINEAR SYSTEMS ARISING IN A NUMERICAL METHOD FOR THE HEAT EQUATION BASED ON LAPLACE TRANSFORMATION AND QUADRATURE

  • WILLIAM MCLEAN (a1) and VIDAR THOMÉE (a2)
Abstract

In earlier work we have studied a method for discretization in time of a parabolic problem, which consists of representing the exact solution as an integral in the complex plane and then applying a quadrature formula to this integral. In application to a spatially semidiscrete finite-element version of the parabolic problem, at each quadrature point one then needs to solve a linear algebraic system having a positive-definite matrix with a complex shift. We study iterative methods for such systems, considering the basic and preconditioned versions of first the Richardson algorithm and then a conjugate gradient method.

    • Send article to Kindle

      To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

      Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

      Find out more about the Kindle Personal Document Service.

      ITERATIVE SOLUTION OF SHIFTED POSITIVE-DEFINITE LINEAR SYSTEMS ARISING IN A NUMERICAL METHOD FOR THE HEAT EQUATION BASED ON LAPLACE TRANSFORMATION AND QUADRATURE
      Available formats
      ×
      Send article to Dropbox

      To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

      ITERATIVE SOLUTION OF SHIFTED POSITIVE-DEFINITE LINEAR SYSTEMS ARISING IN A NUMERICAL METHOD FOR THE HEAT EQUATION BASED ON LAPLACE TRANSFORMATION AND QUADRATURE
      Available formats
      ×
      Send article to Google Drive

      To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

      ITERATIVE SOLUTION OF SHIFTED POSITIVE-DEFINITE LINEAR SYSTEMS ARISING IN A NUMERICAL METHOD FOR THE HEAT EQUATION BASED ON LAPLACE TRANSFORMATION AND QUADRATURE
      Available formats
      ×
Copyright
Corresponding author
For correspondence; e-mail: w.mclean@unsw.edu.au
References
Hide All
[1]Bell, W. N., Olson, L. N. and Schroder, J. B., “PyAMG: algebraic multigrid solvers in Python v2.0”, http://www.pyamg.org.
[2]Benzi, M. and Bertaccini, D., “Block preconditioning of real-valued iterative algorithms for complex linear systems”, IMA J. Numer. Anal. 28 (2008) 598618; doi:10.1093/imanum/drm039.
[3]Davis, T. A., “Algorithm 832: UMFPACK V4.3—an unsymmetric-pattern multifrontal method”, ACM Trans. Math. Software 30 (2004) 196199; doi:10.1145/992200.992206.
[4]Freund, R., “On conjugate gradient type methods and polynomial preconditioners for a class of complex non-Hermitian matrices”, Numer. Math. 57 (1990) 285312; doi:10.1007/BF01386412.
[5]Gavrilyuk, I. P. and Makarov, V. L., “Exponentially convergent algorithms for the operator exponential with applications to inhomogeneous problems in Banach spaces”, SIAM J. Numer. Anal. 43 (2005) 21442171; doi:10.1137/040611045.
[6]Geuzaine, C. and Remacle, J.-F., “Gmsh: a three-dimensional finite element mesh generator with built-in pre- and post-processing facilities”, http://www.geuz.org/gmsh.
[7]in ’t Hout, K. J. and Weideman, J. A. C., “A contour integral method for the Black–Scholes and Heston equations”, SIAM J. Sci. Comput. 33 (2011) 763785; doi:10.1137/090776081.
[8]Jones, M. T. and Plassmann, P. E., “Algorithm 740: Fortran subroutines to compute improved incomplete Cholesky factorizations”, ACM Trans. Math. Software 21 (1995) 1819; doi:10.1145/200979.200986.
[9]McLean, W., Sloan, I. H. and Thomée, V., “Time discretization via Laplace tranformation of an integro-differential equation of parabolic type”, Numer. Math. 102 (2006) 497522; doi:10.1007/s00211-005-0657-7.
[10]McLean, W. and Thomée, V., “Time discretization of an evolution equation via Laplace transforms”, IMA J. Numer. Anal. 24 (2004) 439463; doi:10.1093/imanum/24.3.439.
[11]McLean, W. and Thomée, V., “Maximum-norm error analysis of a numerical solution via Laplace transformation and quadrature of a fractional order evolution equation”, IMA J. Numer. Anal. 30 (2010) 208230; doi:10.1093/imanum/drp004.
[12]McLean, W. and Thomée, V., “Numerical solution via Laplace transforms of a fractional-order evolution equation”, J. Integral Equations Appl. 22 (2010) 5794; doi:10.1216/JIE-2010-22-1-57.
[13]McLean, W. and Thomée, V., “Iterative methods for shifted positive definite linear systems and time discretization of the heat equation”, Preprint, 2011, http://arxiv.org/abs/1111.5105.
[14]Meerbergen, K., “The solution of parametrized symmetric linear systems”, SIAM J. Matrix Anal. Appl. 24 (2003) 10381059; doi:10.1137/S0895479800380386.
[15]Paige, C. C. and Saunders, M. A., “Solution of sparse indefinite systems of linear equations”, SIAM J. Numer. Anal. 12 (1975) 617629; doi:10.1137/0712047.
[16]Sheen, D., Sloan, I. H. and Thomée, V., “A parallel method for time-discretization of parabolic problems based on contour integral representation and quadrature”, Math. Comp. 69 (2000) 177195; doi:10.1090/S0025-5718-99-01098-4.
[17]Sheen, D., Sloan, I. H. and Thomée, V., “A parallel method for time-discretization of parabolic equations based on Laplace transformation and quadrature”, IMA J. Numer. Anal. 23 (2003) 269299; doi:10.1093/imanum/23.2.269.
[18]Simoncini, V., “Restarted full orthogonalization method for shifted linear systems”, BIT Numer. Math. 43 (2003) 459466; doi:10.1023/A:1026000105893.
[19]Simoncini, V. and Szyld, D. B., “Recent computational developments in Krylov subspace methods for linear systems”, Numer. Linear Algebra Appl. 14 (2007) 159; doi:10.1002/nla.499.
Recommend this journal

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

The ANZIAM Journal
  • ISSN: 1446-1811
  • EISSN: 1446-8735
  • URL: /core/journals/anziam-journal
Please enter your name
Please enter a valid email address
Who would you like to send this to? *
×
MathJax

Keywords:

Metrics

Full text views

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

Abstract views

Total abstract views: 171 *
Loading metrics...

* Views captured on Cambridge Core between September 2016 - 13th June 2018. This data will be updated every 24 hours.