pub fn dev_symm(t: Tensor) -> SymmTensor
Deviatoric of the symmetric part symm(T) - (tr/3)*I.
symm(T) - (tr/3)*I