Function tangent 

pub(crate) fn tangent(
    l: i32,
    value: Complex,
    gradient: [Complex; 3],
    unit: [f64; 3],
) -> [Complex; 3]

The tangential part grad S(u) - l S(u) u of a degree-l solid harmonic’s gradient.