On Adaptive Choice of Shifts in Rational Krylov Subspace Reduction of Evolutionary Problems

Abstract

We compute $u(t)=\exp(-tA)\varphi$ using rational Krylov subspace reduction for $0\leq t<\infty$, where $u(t),\varphi\in\mathbf{R}^N$ and $0<A=A^*\in\mathbf{R}^{N\times N}$. A priori optimization of the rational Krylov subspaces for this problem was considered in [V. Druskin, L. Knizhnerman, and M. Zaslavsky, SIAM J. Sci. Comput., 31 (2009), pp. 3760–3780]. There was suggested an algorithm generating sequences of equidistributed shifts, which are asymptotically optimal for the cases with uniform spectral distributions. Here we develop a recursive greedy algorithm for choice of shifts taking into account nonuniformity of the spectrum. The algorithm is based on an explicit formula for the residual in the frequency domain allowing adaptive shift optimization at negligible cost. The effectiveness of the developed approach is demonstrated in an example of the three-dimensional diffusion problem for Maxwell's equation arising in geophysical exploration. We compare our approach with the one using the above-mentioned equidistributed sequences of shifts. Numerical examples show that our algorithm is able to adapt to the spectral density of operator $A$. For examples with near-uniform spectral distributions, both algorithms show the same convergence rates, but the new algorithm produces superior convergence for cases with nonuniform spectra.