Pith. sign in

Rational Krylov for Stieltjes matrix functions: convergence and pole selection

1 Pith paper cite this work. Polarity classification is still indexing.

1 Pith paper citing it
abstract

Evaluating the action of a matrix function on a vector, that is $x=f(\mathcal M)v$, is an ubiquitous task in applications. When $\mathcal M$ is large, one usually relies on Krylov projection methods. In this paper, we provide effective choices for the poles of the rational Krylov method for approximating $x$ when $f(z)$ is either Cauchy-Stieltjes or Laplace-Stieltjes (or, which is equivalent, completely monotonic) and $\mathcal M$ is a positive definite matrix. Relying on the same tools used to analyze the generic situation, we then focus on the case $\mathcal M=I \otimes A - B^T \otimes I$, and $v$ obtained vectorizing a low-rank matrix; this finds application, for instance, in solving fractional diffusion equation on two-dimensional tensor grids. We see how to leverage tensorized Krylov subspaces to exploit the Kronecker structure and we introduce an error analysis for the numerical approximation of $x$. Pole selection strategies with explicit convergence bounds are given also in this case.

years

2026 1

verdicts

CONDITIONAL 1

representative citing papers

citing papers explorer

Showing 1 of 1 citing paper.

  • A Shared Observation Shields Collective Fluctuations while Preserving Local Independence cond-mat.stat-mech · 2026-08-08 · conditional · none · ref 118 · internal anchor

    Conditioning an exchangeable population on its own noisy collective record produces a rank-one, nonpositive Schur shield correction to collective fluctuations while keeping pair correlations at O(1/z), giving a calculable baseline for dynamical heterogeneity in glassy systems.