One wants to build a high resolution map of the sky of properties (temperature, polarization) of the first lights in the universe (some 13.5 billions years ago).
Maths :
Build a fully algebraic 2-lvl preconditioner for a specific form of linear systems to bound the spectrum and reduce iterations.
Ones want to solve \( A x = b \) with :
\begin{equation}
A := P^T N P = \left( \sum_{i=1}^B P_{i}^T N_i P_i \right), b := P^T N d
\end{equation}
Build \(Z_i = Span\left((v_{i,l})_{l=1...k}|\lambda_{i,k} \leq \frac{1}{\tau}\right)\) where :