Adapting the Lanczos algorithm to matrices with almost continuous spectra
Zimmerling, Druskin · math.NA,cs.NA · 2026-09-05 · 原文
We consider approximating B^T(A+sI)^{-1}B, where A\in\mathbb{R}^{n\times n} is large, symmetric positive definite, and B\in\mathbb{R}^{n\times p} with p\ll n. We focus on the case where A has an almost continuous (dense) spectrum: its eigenvalues fill one or more intervals so densely that Krylov methods cannot resolve them individually. Our target is computing multiple-input multiple-output transfer functions arising from large-scale discretizations of problems with continuous spectral measures, such as linear time-invariant PDEs on unbounded domains. Traditional Krylov methods, such as Lanczos or conjugate gradients, resolve individual eigenvalues of a dense discretization while ignoring the underlying continuous spectral measure these points approximate. We argue it is more efficient to model the operator's inherent branch cut than to exhaustively resolve the artificial point spectrum induced by discretization. We adapt the framework of Kreĭn-Nudelman semi-infinite strings to the block Lanczos algorithm, with parameters chosen adaptively by maximizing the energy absorbed at the string termination relative to the energy stored in the string. This yields a low-rank modifi
讲义
讲义·推断 依据「原文」自动生成的结构化摘要（推断），非原文表述；以原文为准。
1. 人话版
We consider approximating B^T(A+sI)^{-1}B, where A\in\mathbb{R}^{n\times n} is large, symmetric positive definite, and B\in\mathbb{R}^{n\times p} with p\ll n.
We focus on the case where A has an almost continuous (dense) spectrum: its eigenvalues fill one or more intervals so densely that Krylov methods cannot resolve them individually.
2. 领域脉络
本文类目：math.NA、cs.NA，属于其所在研究脉络的最新进展。
3. 机制拆解
Our target is computing multiple-input multiple-output transfer functions arising from large-scale discretizations of problems with continuous spectral measures, such as linear time-invariant PDEs on unbounded domains.
Traditional Krylov methods, such as Lanczos or conjugate gradients, resolve individual eigenvalues of a dense discretization while ignoring the underlying continuous spectral measure these points approximate.
We argue it is more efficient to model the operator's inherent branch cut than to exhaustively resolve the artificial point spectrum induced by discretization.
4. 证据与数字
摘要未给出量化结果——留意原文的实验与数据。
5. 反例与边界
摘要未声明局限与反例——这是需要警惕的信号，精读时先问边界。
6. 跨领域连接与意外收获
横跨 2 个类目（math.NA、cs.NA），关注其在你兴趣板块间的迁移面。
7. 可复用方法
把本文机制与你手头项目对照，找一个两周内能验证的最小实验。
8. 术语表
精读时把不熟的术语记入此处，作为下次回忆的锚点。