Added unimodular to calculate unimodular part of a tensor

adtzlr · May 10, 2024 · 6f5d542 · 6f5d542
1 parent 73d42b8
commit 6f5d542
Showing 1 changed file with 11 additions and 0 deletions.
diff --git a/src/matadi/math.py b/src/matadi/math.py
@@ -257,3 +257,14 @@ def astensor(A, scale=1):
 
     else:
         raise ValueError("Unknown shape of input.")
+
+def unimodular(T):
+   """
+   Compute the unimodular part of a tensor.
+
+   The unimodular part of a tensor is a modified version of the tensor where
+   the determinant is raised to the power of (-1/3) and multiplied to the tensor.
+   This operation preserves the isochoric (volume-preserving) part of the tensor
+   while removing the volumetric part.
+   """
+   return (det(T) ** (-1 / 3)) * T