fn dst3d(input: &[f64], n: usize, dst: &dyn Dst1<f64>) -> Vec<f64>
3D separable DST-I, unnormalised: X[k] = Σ x[j]·sin(π·(j+1)·(k+1)/(N+1)) along each axis.