pub fn voigt_average(c_a: &Tensor4, c_b: &Tensor4) -> Tensor4
Arithmetic mean of two stiffness tensors (Voigt bound): C_Voigt = (C_a + C_b) / 2.