Operator Shifting for General Noisy Matrix Systems

TL;DR

This paper extends the operator shifting technique from symmetric to general nonsymmetric matrices to reduce error in noisy linear systems, with theoretical guarantees and practical demonstrations.

Contribution

It generalizes the operator shifting framework to nonsymmetric matrices and analyzes conditions under which shifting reduces error.

Findings

Shifting towards the origin often reduces error in nonsymmetric matrices.

Symmetry and residual norm assumptions ensure optimal error reduction.

Numerical experiments show significant error improvements with operator shifting.

Abstract

In the computational sciences, one must often estimate model parameters from data subject to noise and uncertainty, leading to inaccurate results. In order to improve the accuracy of models with noisy parameters, we consider the problem of reducing error in a linear system with the operator corrupted by noise. Our contribution in this paper is to extend the elliptic operator shifting framework from Etter, Ying '20 to the general nonsymmetric matrix case. Roughly, the operator shifting technique is a matrix analogue of the James-Stein estimator. The key insight is that a shift of the matrix inverse estimate in an appropriately chosen direction will reduce average error. In our extension, we interrogate a number of questions -- namely, whether or not shifting towards the origin for general matrix inverses always reduces error as it does in the elliptic case. We show that this is usually the case, but that there are three key features of the general nonsingular matrices that allow for adversarial examples not possible in the symmetric case. We prove that when these adversarial possibilities are eliminated by the assumption of noise symmetry and the use of the residual norm as the error metric, the optimal shift is always towards the origin, mirroring results from Etter, Ying '20. We also investigate behavior in the small noise regime and other scenarios. We conclude by presenting numerical experiments (with accompanying source code) inspired by Reinforcement Learning to demonstrate that operator shifting can yield substantial reductions in error.