| 541 | //- Return the skew-symmetric part of a tensor |
| 542 | template<class Cmpt> |
| 543 | inline Tensor<Cmpt> skew(const Tensor<Cmpt>& t) |
| 544 | { |
| 545 | return Tensor<Cmpt> |
| 546 | ( |
| 547 | 0.0, 0.5*(t.xy() - t.yx()), 0.5*(t.xz() - t.zx()), |
| 548 | 0.5*(t.yx() - t.xy()), 0.0, 0.5*(t.yz() - t.zy()), |
| 549 | 0.5*(t.zx() - t.xz()), 0.5*(t.zy() - t.yz()), 0.0 |
| 550 | ); |
| 551 | } |
| 552 | |
| 553 | |
| 554 | //- Return the skew-symmetric part of a symmetric tensor |