From feaccae6520cfb97f59c8fd209f549385c3759ee Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Tue, 19 Nov 2013 19:39:14 +0000 Subject: [PATCH] set threshold for values treated as zero to 1.e-12 (was 1.e-15) to avoid flip in strain direction --- processing/post/addStrainTensors.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/processing/post/addStrainTensors.py b/processing/post/addStrainTensors.py index d6502db85..a6457fbad 100755 --- a/processing/post/addStrainTensors.py +++ b/processing/post/addStrainTensors.py @@ -147,7 +147,7 @@ for file in files: stretch['V'] = numpy.dot(F,numpy.linalg.inv(R)) for theStretch in stretches: for i in range(9): - if abs(stretch[theStretch][i%3,i//3]) < 1e-15: # kill nasty noisy data + if abs(stretch[theStretch][i%3,i//3]) < 1e-12: # kill nasty noisy data stretch[theStretch][i%3,i//3] = 0.0 (D,V) = numpy.linalg.eig(stretch[theStretch]) # eigen decomposition (of symmetric matrix) for i,eigval in enumerate(D):