Inexact shift invert Arnoldi method for evolution equations
Linear and nonlinear evolution equations with a first order time derivative, such as the heat equation, the Burgers equation, and the reaction diffusion equation have been used to solve problems in various fields of science. Differential algebraic equations of the first order are derived after space discretization. In the simplest case, the computation of one matrix exponential with a special form is required. In the most complex case, the computation of matrix functions related to its exponentials need to be implemented repeatedly. When computing large matrix functions, the Krylov subspace methods is a viable alternative. The most well known method is the Arnoldi method, but it may require a number of iterations depending on the condition of the matrix. As a solution to this issue, we propose the Inexact Shift invert Arnoldi method to do this more efficiently. As the result, the numerical solution of evolution equations can be computed efficiently with this method. Pertinent numerical experiments establish the effectiveness of this proposed algorithm. References Benzi, M. and Boito, P. Decay properties for functions of matrices over (C^ algebras. Linear Algebra and its Applications 456(1) 174 198, 2014. doi:10.1016/j.laa.2013.11.027 Butcher, J. C. Numerical methods for ordinary differential equations, second edition John Wiley and Sons, Chichester, England, 2008. Caratheodory, C. Uber die gegenseitige Beziehung der Rander bei der konformen Abbildung des Inneren einer Jordanschen Kurve auf einen Kreis. Mathematische Annalen 73(2) 305 320, 1913. http:/gdz.sub.uni goettingen.de/pdfcache/PPN235181684_0073/PPN235181684_0073___LOG_0029.pdf Carra, E. J. Turner, I. W. and Perre, P. A variable stepsize Jacobian free exponential integrator for simulating transport in heterogeneous porous media: Application to wood drying. Journal of Computational Physics 233: 66 82, 2013. doi:10.1016/j.jcp.2012.07.024 Gang, W. Feng, T. and Yimin, W. An inexact shift and invert Arnoldi algorithm for Toeplitz matrix exponential. Numerical Linear Algebra with Applications 22(4) 777 792, 2015. doi:10.1002/nla.1992 Gallopoulos, E. and Saad, Y. Efficient solution of parabolic equations by Krylov approximation methods. SIAM Journal on Scientific Statistics 13(5):1236 1264, 1992. doi:10.1137/0913071 Gockler, T. Rational Krylov subspace methods for \phi functions in exponential integrators. em Karlsruher Instituts fur Technologie, 2014, Ph.D. thesis. http:/d nb.info/1060425408/34 Guttel, S. Rational Krylov methods for operator functions. em Technischen Universitat Bergakademie Freiberg, 2010, Ph.D. thesis. http:/www.qucosa.de/fileadmin/data/qucosa/documents/2764/diss_guettel.pdf Hashimoto, Y. and Nodera, T. Inexact shift invert Arnoldi method for linear evolution equations (Japanese) IPSJ Journal 57(10) 2250 2259, 2016. https:/ipsj.ixsq.nii.ac.jp/ej/?action=pages_view_main&active_action=repository_view_main_item_detail&item_id=175055&item_no=1&page_id=13&block_id=8 Higham, N. J. The scaling and squaring method for the matrix exponential revisited. SIAM Journal on Matrix Analysis and Applications 26(4) 1179 1193, 2005. doi:10.1137/04061101X Hochbruck, M. A short course on exponential integrators. Series in Contemporary Applied Mathematics 17: 29 49, 2015. http:/www.siam.org/students/g2s3/2013/lecturers/Hochbruck/Summary_Hochbruck.pdf Hochbruck, M. and Ostermann, A. Exponential Runge Kutta methods for parabolic problems. Applied Numerical Mathematics 53(2 4) 323 339, 2005. doi:10.1016/j.apnum.2004.08.005 Hochbruck, M. and Ostermann, A. Exponential integrators. Acta Numerica 19:209 286, 2010. doi:10.1017/S0962492910000048 Hochbruck, M. and Lubich, C. On Krylov subspace approximations to the matrix exponential Operator. SIAM Journal on Numerical Analysis 34(5) 1911 1925, 1997. doi:10.1137/S0036142995280572 Hochbruck, M. Lubich, C. and Selhofer, H. Exponential integrators for large systems of differential equations. SIAM Journal on Scientific Computing 19(5):1552 1574, 1997. doi:10.1137/S1064827595295337 Hongqing, Z. Huazhong, S. and Meiyu, D. Numerical solutions of two dimensional Burgers' equations by discrete Adomian decomposition method. Computers and Mathematics with Applications 60(3) 840 848, 2010. doi:10.1016/j.camwa.2010.05.03 Kamel, A. K. Numerical study of Fisher's reaction diffusion equation by the Sinc collocation method. Journal of Computational and Applied Mathematics 137(2) 245 255, 2001. doi:10.1016/S0377 0427(01)00356 9 Lee, S. Pang, H. and Sun, H. Shift invert Arnoldi approximation to the Toeplitz matrix exponential. SIAM Journal on Scientific Computing 32(2) 774 792, 2010. doi:10.1137/090758064 Moret, I. and Novati, P. RD rational approximations of the matrix exponential. BIT Numerical Mathematics 44(3) 595 615, 2004. doi:10.1023/B:BITN.0000046805.27551.3b Moler, C. and Van Loan, C. F. Nineteen dubious ways to compute the exponential of a matrix, Twenty Five Years Later. SIAM Review 45(1) 3 49, 2003. doi:10.1137/S00361445024180 Novati, P. Using the restricted denominator rational Arnoldi method for exponential Integrators. SIAM Journal on Matrix Analysis and Applications 32(4) 1537 1558, 2011. doi:10.1137/100814202 Svoboda, Z. The convective diffusion equation and its use in building physics. International Journal on Architectural Science 1(2) 68 79, 2000. http:/www.bse.polyu.edu.hk/researchCentre/Fire_Engineering/summary_of_output/journal/IJAS/V1/p.68 79.pdf Van der Vorst, H. A. Bi CGSTAB: A fast and smoothly converging variant of Bi CG for the solution of nonsymmetric linear systems. SlAM Journal on Scientific and Statistical Computing 13(2) 631 644, 1992. doi:10.1137/0913035