Function diff_4d_temporal 

pub fn diff_4d_temporal(arr: &Array4D<f64>) -> Array4D<f64>

Compute temporal finite differences: output[t] = input[t+1] - input[t]. Output shape is [T-1, D, H, W].