English

The Index of Invariance and its Implications for a Parameterized Least Squares Problem

Numerical Analysis 2020-09-02 v3 Numerical Analysis Optimization and Control

Abstract

We study the problem xb,ω:=arg minxS(A+ωI)1/2(bAx)2x_{b,\omega} := \text{arg min}_{x \in \mathcal{S}} \|(A + \omega I)^{-1/2} (b - Ax)\|_2, with A=AA = A^*, for a subspace S\mathcal{S} of Fn\mathbb{F}^n (F=R\mathbb{F} = \mathbb{R} or C\mathbb{C}), and ω>λmin(A)\omega > -\lambda_{min}(A). We show that there exists a subspace Y\mathcal{Y} of Fn\mathbb{F}^n, independent of bb, such that {xb,ωxb,μω,μ>λmin(A)}Y\{x_{b,\omega} - x_{b,\mu} \mid \omega,\mu > -\lambda_{min}(A)\} \subseteq \mathcal{Y}, where dim(Y)dim(S+AS)dim(S)=IndA(S)\dim(\mathcal{Y}) \leq \dim(\mathcal{S} + A\mathcal{S}) - \dim(\mathcal{S}) = \mathbf{Ind}_A(\mathcal{S}), a quantity which we call the index of invariance of S\mathcal{S} with respect to AA. In particular if S\mathcal{S} is a Krylov subspace, this implies the low dimensionality result of Hallman & Gu (2018). The problem is also such that when AA is positive and S\mathcal{S} is a Krylov subspace, it reduces to CG for ω=0\omega = 0 and to MINRES for ω\omega \to \infty. We study several properties of IndA(S)\mathbf{Ind}_A(\mathcal{S}) in relation to AA and S\mathcal{S}. We show that the dimension of the affine subspace Xb\mathcal{X}_b containing the solutions xb,ωx_{b,\omega} can be smaller than IndA(S)\mathbf{Ind}_A(\mathcal{S}) for all bb. However, we also exhibit some sufficient conditions on AA and S\mathcal{S}, under which X:=Span{xb,ωxb,μbFn,ω,μ>λmin(A)}\mathcal{X} := \text{Span}{\{x_{b,\omega} - x_{b,\mu} \mid b \in \mathbb{F}^n, \omega,\mu > -\lambda_{min}(A)\}} has dimension equal to IndA(S)\mathbf{Ind}_A(\mathcal{S}). We then study the injectivity of the map ωxb,ω\omega \mapsto x_{b,\omega}, leading us to a proof of the convexity result from Hallman & Gu (2018). We finish by showing that sets such as M(S,S)={AFn×nS+AS=S}M(\mathcal{S},\mathcal{S}') = \{A \in \mathbb{F}^{n \times n} \mid \mathcal{S} + A\mathcal{S} = \mathcal{S}'\}, for nested subspaces SSFn\mathcal{S} \subseteq \mathcal{S}' \subseteq \mathbb{F}^n, form smooth real manifolds, and explore some topological relationships between them.

Keywords

Cite

@article{arxiv.2008.11154,
  title  = {The Index of Invariance and its Implications for a Parameterized Least Squares Problem},
  author = {Léopold Cambier and Rahul Sarkar},
  journal= {arXiv preprint arXiv:2008.11154},
  year   = {2020}
}