Unimodular Part of Tensor
Unimodular Part of a rank 2 Tensor.
\[\begin{align} \hat{\boldsymbol{C}} &= \det(\boldsymbol{C})^{-1/3} \ \boldsymbol{C} \\ (\det(\hat{\boldsymbol{C}}) &= 1) \nonumber \end{align}\]Property | Value |
---|---|
Result | Tensor-valued function |
Data Types | Tensor2 , Tensor2s |
Example
type(Tensor2) :: C, hatC
! type(Tensor2s) :: C, hatC
hatC = unimodular(C)
print *, "proof:"
print *, (hatC - (det(C)**(-1./3.)*C))