| 296 | //- Return the spherical part of a symmetric tensor |
| 297 | template<class Cmpt> |
| 298 | inline SphericalTensor<Cmpt> sph(const SymmTensor<Cmpt>& st) |
| 299 | { |
| 300 | return (1.0/3.0)*tr(st); |
| 301 | } |
| 302 | |
| 303 | |
| 304 | //- Return the symmetric part of a symmetric tensor, i.e. itself |