Solution (source code)

= Solution

Let $(\sigma_j,u_j,v_j)$ be a singular system for $A$. The <Picard criterion> for $f\in\operatorname{dom}(A^\dagger)$ is
$$
\sum_j\frac{|\langle f,v_j\rangle|^2}{\sigma_j^2}<\infty
$$
after discarding the component in $(\operatorname{ran}A)^\perp$. For $0<\tau<\|A\|^{-2}$, the partial series acts diagonally:
$$
Q_Nf
=\sum_j
\frac{1-(1-\tau\sigma_j^2)^{N+1}}{\sigma_j}
\langle f,v_j\rangle u_j.
$$
For each $j$ the multiplier in the numerator tends to one and lies in $[0,1]$. The Picard summability condition therefore supplies an $\ell^2$ dominating sequence, so the <dominated convergence theorem> gives
$$
\boxed{Q_Nf\longrightarrow
\sum_j\frac{\langle f,v_j\rangle}{\sigma_j}u_j
=A^\dagger f}.
$$
This is the series form of <Landweber iteration>.