fn prox_wavelet_l1(x: &mut [f64], plan: &WaveletPlan, th: f64)
Wavelet-L1 proximal operator for an orthonormal transform: soft-threshold all coefficients (approximation + details) by th, in place.
th