Function fgrad_periodic_inplace_f32 

pub(crate) fn fgrad_periodic_inplace_f32(
    gx: &mut [f32],
    gy: &mut [f32],
    gz: &mut [f32],
    x: &[f32],
    nx: usize,
    ny: usize,
    nz: usize,
    vsx: f32,
    vsy: f32,
    vsz: f32,
)

Forward difference gradient with periodic boundary conditions (f32, in-place) Matches MATLAB’s gradfp_mex used inside MEDI iterations. At boundaries, wraps around: dx(end) = (x(1) - x(end)) / h