@conference {peters2019SEGans,
	title = {A numerical solver for least-squares sub-problems in 3D wavefield reconstruction inversion and related problem formulations},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2019},
	note = {(SEG, San Antonio)},
	month = {09},
	pages = {1536-1540},
	abstract = {Recent years saw a surge of interest in seismic waveform inversion approaches based on quadratic-penalty or augmented-Lagrangian methods, including Wavefield Reconstruction Inversion. These methods typically need to solve a least-squares sub-problem that contains a discretization of the Helmholtz equation. Memory requirements for direct solvers are often prohibitively large in three dimensions, and this limited the examples in the literature to two dimensions. We present an algorithm that uses iterative Helmholtz solvers as a black-box to solve the least-squares problem corresponding to 3D grids. This algorithm enables Wavefield Reconstruction Inversion and related formulations, in three dimensions. Our new algorithm also includes a root-finding method to convert a penalty into a constraint on the data-misfit without additional computational cost, by reusing precomputed quantities. Numerical experiments show that the cost of parallel communication and other computations are small compared to the main cost of solving one Helmholtz problem per source and one per receiver.},
	keywords = {Full-waveform inversion, least-squares, numerical linear algebra, penalty method, quadratic constraints, root finding, SEG},
	doi = {10.1190/segam2019-3216638.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/peters2019SEGans/peters2019SEGans.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/peters2019SEGans/peters2019SEGans_pres.pdf},
	author = {Bas Peters and Felix J. Herrmann}
}
@conference {peters2015PRECONasl,
	title = {An algorithm for solving least-squares problems with a Helmholtz block and multiple right-hand-sides},
	booktitle = {International Conference On Preconditioning Techniques For Scientific And Industrial Applications},
	year = {2015},
	note = {(PRECON, The Netherlands)},
	month = {06},
	keywords = {CGMN, Helmholtz, least-squares, matrix-free, randomized linear algebra},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/PRECON/2015/peters2015PRECONasl/peters2015PRECONasl_pres.pdf},
	author = {Bas Peters and Chen Greif and Felix J. Herrmann}
}
@conference {peters2015SIAMmfq,
	title = {Matrix-free quadratic-penalty methods for PDE-constrained optimization},
	booktitle = {SIAM Conference on Computational Science and Engineering},
	year = {2015},
	note = {(SIAM, Salt Lake City)},
	month = {03},
	abstract = {The large scale of seismic waveform inversion makes matrix free implementations essential. We show how to exploit the quadratic penalty structure to construct matrix free reduced-space and full-space algorithms, which have some advantages over the commonly used Lagrangian based methods for PDE-constrained optimization. This includes the construction of effective and sparse Hessian approximations and reduced sensitivity to the initial guess. A computational bottleneck is the need to solve a large least squares problem with a PDE block. When direct solvers are not available, we propose a fast matrix free iterative approach with reasonable memory requirements. It takes advantage of the structure of the least squares problem with a combination of preconditioning, low rank decomposition and deflation.},
	keywords = {least-squares, low-rank, PDE-constrained optimization, quadratic-penalty, waveform inversion},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2015/peters2015SIAMmfq/peters2015SIAMmfq.pdf},
	author = {Bas Peters and Felix J. Herrmann and Chen Greif}
}
