pub fn dev2(vol: &VolTensorField) -> VolTensorField
dev2(T) = T − (2·tr(T)/3)·I of a tensor field → tensor field.
dev2(T) = T − (2·tr(T)/3)·I