fn apply_ht_h(
dker: &[f64],
mask: &[f64],
x: &[f64],
nx: usize,
ny: usize,
nz: usize,
) -> Vec<f64>Expand description
Apply H’H to a vector: ifft(DKER * fft(M * ifft(DKER * fft(x))))
H(x) = M * ifft(DKER * fft(x)) H’(y) = ifft(DKER * fft(M * y))