223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

checked for typo

This commit is contained in:
Samad Vakili 2020-06-09 14:26:22 +02:00
parent 627e62439b
commit 824682d236
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
python/tests/test_mechanics.py
View File

@ -42,7 +42,7 @@ def Mises_stress(sigma):
def PK2(P,F): def PK2(P,F):
S = np.dot(_np.linalg.inv(F),P) S = np.dot(np.linalg.inv(F),P)
return symmetric(S) return symmetric(S)
@ -69,10 +69,10 @@ def strain_tensor(F,t,m):
w,n = np.linalg.eigh(C) w,n = np.linalg.eigh(C)
if m > 0.0: if m > 0.0:
eps = 1.0/(2.0*abs(m)) * (+ np.matmul(n,_np.einsum('ij,ikj->ijk',w**m,n)) eps = 1.0/(2.0*abs(m)) * (+ np.matmul(n,np.einsum('ij,ikj->ijk',w**m,n))
- np.einsum('ijk->ijk',np.eye(3))) - np.einsum('ijk->ijk',np.eye(3)))
elif m < 0.0: elif m < 0.0:
eps = 1.0/(2.0*abs(m)) * (- np.matmul(n,_np.einsum('ij,ikj->ijk',w**m,n)) eps = 1.0/(2.0*abs(m)) * (- np.matmul(n,np.einsum('ij,ikj->ijk',w**m,n))
+ np.einsum('ijk->ijk',np.eye(3))) + np.einsum('ijk->ijk',np.eye(3)))
else: else:
eps = np.matmul(n,np.einsum('ij,ikj->ijk',0.5*np.log(w),n)) eps = np.matmul(n,np.einsum('ij,ikj->ijk',0.5*np.log(w),n))
@ -157,7 +157,7 @@ class TestMechanics:
@pytest.mark.parametrize('vectorized,single',[(mechanics.strain_tensor,strain_tensor)]) @pytest.mark.parametrize('vectorized,single',[(mechanics.strain_tensor,strain_tensor)])
def test_vectorize_strain_tensor(self,vectroized,single): def test_vectorize_strain_tensor(self,vectorized,single):
F = np.random.rand(self.n,3,3) F = np.random.rand(self.n,3,3)
F_vec = np.random.rand(self.n//10,10,3,3) F_vec = np.random.rand(self.n//10,10,3,3)
t = ['V','U'][np.random.randint(0,2)] t = ['V','U'][np.random.randint(0,2)]