Rapidly convergent series expansions for a class of resolvents

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteur principal: Milton, Graeme W.
Format: Preprint
Publié: 2024
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866913454703509504
author Milton, Graeme W.
author_facet Milton, Graeme W.
contents Following advances in the abstract theory of composites, we develop rapidly converging series expansions about $z=\infty$ for the resolvent ${\bf R}(z)=[z{\bf I}-{\bf P}^\dagger{\bf Q}{\bf P}]^{-1}$ where ${\bf Q}$ is an orthogonal projection and ${\bf P}$ is such that ${\bf P}{\bf P}^\dagger$ is an orthogonal projection. It is assumed that the spectrum of ${\bf P}^\dagger{\bf Q}{\bf P}$ lies within the interval $[z^-,z^+]$ for some known $z^+\leq 1$ and $z^-\geq 0$ and that the actions of the projections ${\bf Q}$ and ${\bf P}{\bf P}^\dagger$ are easy to compute. The series converges in the entire $z$-plane excluding the cut $[z^-,z^+]$. It is obtained using subspace substitution, where the desired resolvent is tied to a resolvent in a larger space and ${\bf Q}$ gets replaced by a projection $\underline{\bf Q}$ that is no longer orthogonal. When $z$ is real the rate of convergence of the new method matches that of the conjugate gradient method.
format Preprint
id arxiv_https___arxiv_org_abs_2407_15751
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Rapidly convergent series expansions for a class of resolvents
Milton, Graeme W.
Numerical Analysis
Mathematical Physics
Analysis of PDEs
Classical Analysis and ODEs
Following advances in the abstract theory of composites, we develop rapidly converging series expansions about $z=\infty$ for the resolvent ${\bf R}(z)=[z{\bf I}-{\bf P}^\dagger{\bf Q}{\bf P}]^{-1}$ where ${\bf Q}$ is an orthogonal projection and ${\bf P}$ is such that ${\bf P}{\bf P}^\dagger$ is an orthogonal projection. It is assumed that the spectrum of ${\bf P}^\dagger{\bf Q}{\bf P}$ lies within the interval $[z^-,z^+]$ for some known $z^+\leq 1$ and $z^-\geq 0$ and that the actions of the projections ${\bf Q}$ and ${\bf P}{\bf P}^\dagger$ are easy to compute. The series converges in the entire $z$-plane excluding the cut $[z^-,z^+]$. It is obtained using subspace substitution, where the desired resolvent is tied to a resolvent in a larger space and ${\bf Q}$ gets replaced by a projection $\underline{\bf Q}$ that is no longer orthogonal. When $z$ is real the rate of convergence of the new method matches that of the conjugate gradient method.
title Rapidly convergent series expansions for a class of resolvents
topic Numerical Analysis
Mathematical Physics
Analysis of PDEs
Classical Analysis and ODEs
url https://arxiv.org/abs/2407.15751