| 320 | //- Return the deviatoric part of a symmetric tensor |
| 321 | template<class Cmpt> |
| 322 | inline SymmTensor<Cmpt> dev(const SymmTensor<Cmpt>& st) |
| 323 | { |
| 324 | return st - SphericalTensor<Cmpt>::oneThirdI*tr(st); |
| 325 | } |
| 326 | |
| 327 | |
| 328 | //- Return the deviatoric part of a symmetric tensor |