Strong Lensing Reconstruction

We present a general linear algorithm for measuring the surface mass density 1-\kappa from the observable reduced shear g=\gamma/(1-\kappa) in the strong lensing regime. We show that in general, the observed polarization field can be decomposed into ``electric'' and ``magnetic'' components, which have independent and redundant solutions, but perfectly orthogonal noise properties. By combining these solutions, one can increase the signal-to-noise ratio by \sqrt{2}. The solutions allow dynamic optimization of signal and noise, both in real and Fourier space (using arbitrary smoothing windows). Boundary conditions have no effect on the reconstructions, apart from its effect on the signal-to-noise. Many existing reconstruction techniques are recovered as special cases of this framework. The magnetic solution has the added benefit of yielding the global and local parity of the reconstruction in a single step.
Comment: final accepted version for ApJL