Function bdiv_inplace 

pub fn bdiv_inplace(
    div: &mut [f64],
    gx: &[f64],
    gy: &[f64],
    gz: &[f64],
    nx: usize,
    ny: usize,
    nz: usize,
    vsx: f64,
    vsy: f64,
    vsz: f64,
)

Backward divergence operator (in-place)

Computes backward divergence with periodic boundary conditions. Writes result directly into pre-allocated output buffer.

Arguments

• div - Output divergence (must be pre-allocated to nxnynz)
• gx, gy, gz - Gradient components
• nx, ny, nz - Array dimensions
• vsx, vsy, vsz - Voxel sizes