Skip to main content Accessibility help
×
Home
Hostname: page-component-5bf98f6d76-64tqw Total loading time: 0.266 Render date: 2021-04-21T21:56:39.552Z Has data issue: true Feature Flags: { "shouldUseShareProductTool": true, "shouldUseHypothesis": true, "isUnsiloEnabled": true, "metricsAbstractViews": false, "figures": false, "newCiteModal": false, "newCitedByModal": true }

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

Published online by Cambridge University Press:  15 August 2012

WILLIAM MCLEAN
Affiliation:
School of Mathematics and Statistics, The University of New South Wales, Sydney 2052, Australia (email: w.mclean@unsw.edu.au)
VIDAR THOMÉE
Affiliation:
Mathematical Sciences, Chalmers University of Technology and University of Gothenburg, S-412 96 Gothenburg, Sweden (email: thomee@chalmers.se)
Corresponding
Rights & Permissions[Opens in a new window]

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.

Type
Research Article
Copyright
Copyright © Australian Mathematical Society 2012

References

[1]Bell, W. N., Olson, L. N. and Schroder, J. B., “PyAMG: algebraic multigrid solvers in Python v2.0”, http://www.pyamg.org.Google Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.Google Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.Google Scholar
[14]Meerbergen, K., “The solution of parametrized symmetric linear systems”, SIAM J. Matrix Anal. Appl. 24 (2003) 10381059; doi:10.1137/S0895479800380386.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[18]Simoncini, V., “Restarted full orthogonalization method for shifted linear systems”, BIT Numer. Math. 43 (2003) 459466; doi:10.1023/A:1026000105893.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar

Full text views

Full text views reflects PDF downloads, PDFs sent to Google Drive, Dropbox and Kindle and HTML full text views.

Total number of HTML views: 0
Total number of PDF views: 101 *
View data table for this chart

* Views captured on Cambridge Core between September 2016 - 21st April 2021. This data will be updated every 24 hours.

You have Access

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

Reply to: Submit a response


Your details


Conflicting interests

Do you have any conflicting interests? *