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))