fn calculate_angle_gradient( geometry: &Geometry, i: usize, j: usize, k: usize, ) -> ([f64; 3], [f64; 3], [f64; 3])