From 79d2432c6cfd3f8df7f6eecdd0e039c08d192d17 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 13 Sep 2019 06:33:00 -0700 Subject: [PATCH] R not needed --- processing/post/addStrainTensors.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/processing/post/addStrainTensors.py b/processing/post/addStrainTensors.py index 97c83a9c4..b7b21c92e 100755 --- a/processing/post/addStrainTensors.py +++ b/processing/post/addStrainTensors.py @@ -136,9 +136,9 @@ for name in filenames: for column in items['tensor']['column']: # loop over all requested defgrads F = np.array(list(map(float,table.data[column:column+items['tensor']['dim']])),'d').reshape(items['tensor']['shape']) (U,S,Vh) = np.linalg.svd(F) # singular value decomposition - R = np.dot(U,Vh) # rotation of polar decomposition - stretch['U'] = np.dot(R.T,F) # F = RU - stretch['V'] = np.dot(F,R.T) # F = VR + R_inv = np.dot(U,Vh).T # rotation of polar decomposition + stretch['U'] = np.dot(R_inv,F) # F = RU + stretch['V'] = np.dot(F,R_inv) # F = VR for theStretch in stretches: stretch[theStretch] = np.where(abs(stretch[theStretch]) < 1e-12, 0, stretch[theStretch]) # kill nasty noisy data