From 2d08eb68b12c9506c5a03a0d255cf9a9f07943db Mon Sep 17 00:00:00 2001
From: Martin Diehl <m.diehl@mpie.de>
Date: Wed, 23 Oct 2019 14:39:52 +0200
Subject: [PATCH] bugfix for (3x3) version of Cauchy

---
 python/damask/mechanics.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/python/damask/mechanics.py b/python/damask/mechanics.py
index 1e91ccbb7..351feec50 100644
--- a/python/damask/mechanics.py
+++ b/python/damask/mechanics.py
@@ -15,7 +15,7 @@ def Cauchy(F,P):
 
   """
   if np.shape(F) == np.shape(P) == (3,3):
-    sigma = 1.0/np.linalg.det(F) * np.dot(F,P)
+    sigma = 1.0/np.linalg.det(F) * np.dot(P,F.T)
   else:
     sigma = np.einsum('i,ijk,ilk->ijl',1.0/np.linalg.det(F),P,F)
   return symmetric(sigma)