pub fn hill_average(c_a: &Tensor4, c_b: &Tensor4) -> Option<Tensor4>
Hill average: arithmetic mean of Voigt and Reuss bounds.