[cig-commits] r19068 - short/3D/PyLith/trunk/playpen/postproc

[willic3 at geodynamics.org]

Author: willic3
Date: 2011-10-13 20:23:52 -0700 (Thu, 13 Oct 2011)
New Revision: 19068

Modified:
   short/3D/PyLith/trunk/playpen/postproc/princaxes.py
Log:
Discovered problem when extracting eigenvectors.
I was getting rows rather than columns.


Modified: short/3D/PyLith/trunk/playpen/postproc/princaxes.py
===================================================================
--- short/3D/PyLith/trunk/playpen/postproc/princaxes.py	2011-10-14 01:05:54 UTC (rev 19067)
+++ short/3D/PyLith/trunk/playpen/postproc/princaxes.py	2011-10-14 03:23:52 UTC (rev 19068)
@@ -217,9 +217,9 @@
     for point in xrange(self.numTensorPoints):
       tensor = self.tensorSorted[point, :]
       tensorOrdered, eigenValuesOrdered = self._compPrincAxes(tensor)
-      self.minPrincAxis[point,:] = tensorOrdered[0]
-      self.intPrincAxis[point,:] = tensorOrdered[1]
-      self.maxPrincAxis[point,:] = tensorOrdered[2]
+      self.minPrincAxis[point,:] = tensorOrdered[:,0]
+      self.intPrincAxis[point,:] = tensorOrdered[:,1]
+      self.maxPrincAxis[point,:] = tensorOrdered[:,2]
       self.minEigenValue[point] = eigenValuesOrdered[0]
       self.intEigenValue[point] = eigenValuesOrdered[1]
       self.maxEigenValue[point] = eigenValuesOrdered[2]