diff --git a/processing/post/addStrainTensors.py b/processing/post/addStrainTensors.py index 8516d9be8..f35c965ba 100755 --- a/processing/post/addStrainTensors.py +++ b/processing/post/addStrainTensors.py @@ -137,8 +137,8 @@ for name in filenames: 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_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 + stretch['U'] = np.dot(R_inv,F) # F = RU for theStretch in stretches: (D,V) = np.linalg.eigh((stretch[theStretch]+stretch[theStretch].T)*0.5) # eigen decomposition (of symmetric(ed) matrix)