| 257 | //- Inner-sqr of a symmetric tensor |
| 258 | template<class Cmpt> |
| 259 | inline SymmTensor<Cmpt> |
| 260 | innerSqr(const SymmTensor<Cmpt>& st) |
| 261 | { |
| 262 | return SymmTensor<Cmpt> |
| 263 | ( |
| 264 | st.xx()*st.xx() + st.xy()*st.xy() + st.xz()*st.xz(), |
| 265 | st.xx()*st.xy() + st.xy()*st.yy() + st.xz()*st.yz(), |
| 266 | st.xx()*st.xz() + st.xy()*st.yz() + st.xz()*st.zz(), |
| 267 | |
| 268 | st.xy()*st.xy() + st.yy()*st.yy() + st.yz()*st.yz(), |
| 269 | st.xy()*st.xz() + st.yy()*st.yz() + st.yz()*st.zz(), |
| 270 | |
| 271 | st.xz()*st.xz() + st.yz()*st.yz() + st.zz()*st.zz() |
| 272 | ); |
| 273 | } |
| 274 | |
| 275 | |
| 276 | template<class Cmpt> |