bugfix for (3x3) version of Cauchy

2019-10-23 14:39:52 +02:00 · 2019-10-23 14:39:52 +02:00 · 2d08eb68b1
parent 578cbf9bbf
commit 2d08eb68b1
1 changed files with 1 additions and 1 deletions
--- 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)