pub fn oamp_state_evolution(
config: &MessagePassing,
singular_values: &[f64],
n: usize,
noise_var: f64,
linear: OampLinear,
) -> Vec<StateEvolutionStep>Expand description
The empirical state evolution of OAMP with the optimal denoiser (Ma & Ping, eqs. 32
and 39): for the singular_values of $H$ ($N$ columns), from $v^{2}{0} = \mathrm{Var}[S]$,
$\tau^{2}{t} = \Phi_{t}(v^{2}{t})$ (eq. 32), the MSE of $F$ is
$\mathrm{mmse}{B}(\tau^{2}{t})$ and $v^{2}{t+1} = (1/\mathrm{mmse}{B} - 1/\tau^{2}{t})^{-1}$,
for the prior, domain and iterations of config and the linear estimator linear.