fn dst3d(input: &[f64], n: usize, sin_table: &[Vec<f64>]) -> Vec<f64>
3D separable DST-I (Type I Discrete Sine Transform).