pub fn symgrad_inplace_f32(
sxx: &mut [f32],
sxy: &mut [f32],
sxz: &mut [f32],
syy: &mut [f32],
syz: &mut [f32],
szz: &mut [f32],
wx: &[f32],
wy: &[f32],
wz: &[f32],
nx: usize,
ny: usize,
nz: usize,
vsx: f32,
vsy: f32,
vsz: f32,
)Expand description
Symmetric gradient operator for TGV regularization (in-place, f32)
Computes the symmetric gradient tensor from a vector field w = (wx, wy, wz). The output is a 6-component symmetric tensor: q[0] = ∂wx/∂x (Sxx) q[1] = (∂wx/∂y + ∂wy/∂x) / 2 (Sxy) q[2] = (∂wx/∂z + ∂wz/∂x) / 2 (Sxz) q[3] = ∂wy/∂y (Syy) q[4] = (∂wy/∂z + ∂wz/∂y) / 2 (Syz) q[5] = ∂wz/∂z (Szz)