English

Approximation of the Scattering Amplitude using Nonsymmetric Saddle Point Matrices

Numerical Analysis 2015-03-24 v3

Abstract

In this paper we examine iterative methods for solving the forward (Ax=bA{\bf x}={\bf b}) and adjoint (ATy=gA^{T}{\bf y}={\bf g}) systems of linear equations used to approximate the scattering amplitude, defined by gTx=yTb{\bf g}^{T}{\bf x}={\bf y}^{T}{\bf b}. Based on an idea first proposed by Gene Golub, we use a conjugate gradient-like iteration for a nonsymmetric saddle point matrix that is constructed so as to have a real positive spectrum. Numerical experiments show that this method is more consistent than known methods for computing the scattering amplitude such as GLSQR or QMR. We then demonstrate that when combined with known preconditioning techniques, the proposed method exhibits more rapid convergence than state-of-the-art iterative methods for nonsymmetric systems.

Keywords

Cite

@article{arxiv.1503.00655,
  title  = {Approximation of the Scattering Amplitude using Nonsymmetric Saddle Point Matrices},
  author = {Amber S. Robertson and James V. Lambers},
  journal= {arXiv preprint arXiv:1503.00655},
  year   = {2015}
}

Comments

25 pages, 8 figures

R2 v1 2026-06-22T08:42:14.493Z