From b11ad8c774e9353558e32fa15b36d9ce834f6c53 Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Thu, 23 Mar 2023 07:51:34 +0100
Subject: [PATCH 1/3] helpful information on strain calculation

---
 python/damask/_result.py | 21 +++++++++++++++++----
 1 file changed, 17 insertions(+), 4 deletions(-)

diff --git a/python/damask/_result.py b/python/damask/_result.py
index f6dddc9a0..b18d9a11c 100644
--- a/python/damask/_result.py
+++ b/python/damask/_result.py
@@ -1203,7 +1203,7 @@ class Result:
                 'label': f"epsilon_{t}^{m}({F['label']})",
                 'meta':  {
                           'unit':        F['meta']['unit'],
-                          'description': f'strain tensor of order {m} based on {side} stretch tensor '+\
+                          'description': f'Seth-Hill strain tensor of order {m} based on {side} stretch tensor '+\
                                          f"of {F['label']} ({F['meta']['description']})",
                           'creator':     'add_strain'
                           }
@@ -1212,10 +1212,11 @@ class Result:
                    F: str = 'F',
                    t: Literal['V', 'U'] = 'V',
                    m: float = 0.0):
-        """
-        Add strain tensor of a deformation gradient.
+        r"""
+        Add strain tensor (Seth-Hill family) of a deformation gradient.
 
-        For details, see damask.mechanics.strain.
+        By default, the logarithmic strain based on the
+        left stretch tensor is added.
 
         Parameters
         ----------
@@ -1249,6 +1250,18 @@ class Result:
         spatial/Eulerian strain measures (based on 'V') for elastic strains
         when calculating averages.
 
+        The strain is defined as:
+
+        .. math::
+
+            \vb*{\epsilon}_V^{(m)} = \frac{1}{2m} (\vb{V}^{2m} - \vb{I}) \\\\
+            \vb*{\epsilon}_U^{(m)} = \frac{1}{2m} (\vb{U}^{2m} - \vb{I})
+
+        References
+        ----------
+        | https://en.wikipedia.org/wiki/Finite_strain_theory
+        | https://de.wikipedia.org/wiki/Verzerrungstensor
+
         """
         self._add_generic_pointwise(self._add_strain,{'F':F},{'t':t,'m':m})
 

From 09858317d8c8077df9cbe7c719e34e554804a8ca Mon Sep 17 00:00:00 2001
From: Yi Hu <y.hu@mpie.de>
Date: Tue, 4 Jul 2023 16:44:31 +0200
Subject: [PATCH 2/3] add m=0 (default) case in docstring

---
 PRIVATE                  | 2 +-
 python/damask/_result.py | 5 +++++
 2 files changed, 6 insertions(+), 1 deletion(-)

diff --git a/PRIVATE b/PRIVATE
index 9f4ffce8b..fc04b9ef6 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 9f4ffce8b2df951191a14dc3229de1aee6e544e6
+Subproject commit fc04b9ef621161e60a2f8b72bfa8c99e77687c71
diff --git a/python/damask/_result.py b/python/damask/_result.py
index fd70ff1c0..5d04a79b7 100644
--- a/python/damask/_result.py
+++ b/python/damask/_result.py
@@ -1252,6 +1252,11 @@ class Result:
 
         .. math::
 
+            m = 0 \\\\
+            \vb*{\epsilon}_V^{(0)} = \ln (\vb{V}) \\\\
+            \vb*{\epsilon}_U^{(0)} = \ln (\vb{U}) \\\\
+
+            m \neq 0 \\\\
             \vb*{\epsilon}_V^{(m)} = \frac{1}{2m} (\vb{V}^{2m} - \vb{I}) \\\\
             \vb*{\epsilon}_U^{(m)} = \frac{1}{2m} (\vb{U}^{2m} - \vb{I})
 

From bd0a35a94648478eb2b9c7f3c100b7a07341236e Mon Sep 17 00:00:00 2001
From: Yi Hu <y.hu@mpie.de>
Date: Tue, 4 Jul 2023 18:30:06 +0200
Subject: [PATCH 3/3] upd PRIVATE

---
 PRIVATE | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/PRIVATE b/PRIVATE
index fc04b9ef6..9f4ffce8b 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit fc04b9ef621161e60a2f8b72bfa8c99e77687c71
+Subproject commit 9f4ffce8b2df951191a14dc3229de1aee6e544e6