Function shallow_water::stafft::forrdx4

pub fn forrdx4(
    a: ArrayView3<'_, f64>, 
    b: ArrayViewMut3<'_, f64>, 
    nv: usize, 
    lv: usize, 
    cosine: ArrayView2<'_, f64>, 
    sine: ArrayView2<'_, f64>
)

Radix four physical to Hermitian FFT with 'decimation in time'.