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

typo caused problems with more than one tensor

This commit is contained in:
Martin Diehl 2015-06-06 15:42:29 +00:00
parent 83d5686995
commit d99e34d870
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
processing/post/addSpectralDecomposition.py
View File

@ -60,7 +60,7 @@ for file in files:
column[label] = table.labels.index(key) # remember columns of requested data column[label] = table.labels.index(key) # remember columns of requested data
# ------------------------------------------ assemble header --------------------------------------- # ------------------------------------------ assemble header ---------------------------------------
for labels in active: for label in active:
table.labels_append(['%i_eigval(%s)'%(i+1,label) for i in xrange(3)]) # extend ASCII header with new labels table.labels_append(['%i_eigval(%s)'%(i+1,label) for i in xrange(3)]) # extend ASCII header with new labels
table.labels_append(['%i_eigvec(%s)'%(i+1,label) for i in xrange(9)]) # extend ASCII header with new labels table.labels_append(['%i_eigvec(%s)'%(i+1,label) for i in xrange(9)]) # extend ASCII header with new labels
table.head_write() table.head_write()
@ -68,7 +68,7 @@ for file in files:
# ------------------------------------------ process data ------------------------------------------ # ------------------------------------------ process data ------------------------------------------
outputAlive = True outputAlive = True
while outputAlive and table.data_read(): # read next data line of ASCII table while outputAlive and table.data_read(): # read next data line of ASCII table
for labels in active: # loop over requested data for label in active: # loop over requested data
tensor = np.array(map(float,table.data[column[label]:column[label]+datainfo['tensor']['len']])).\ tensor = np.array(map(float,table.data[column[label]:column[label]+datainfo['tensor']['len']])).\
reshape((datainfo['tensor']['len']//3,3)) reshape((datainfo['tensor']['len']//3,3))
(u,v) = np.linalg.eig(tensor) (u,v) = np.linalg.eig(tensor)