pub fn dev_two_symm(t: Tensor) -> SymmTensor
Deviatoric of twice the symmetric part twoSymm(T) - (2*tr/3)*I.
twoSymm(T) - (2*tr/3)*I