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

now aware of "origin" keyword in seeds-file

This commit is contained in:
Philip Eisenlohr 2014-06-07 18:15:08 +00:00
parent e3c4e0ecfc
commit b0ab048d7b
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
processing/pre/seeds_check.py
View File

@ -92,6 +92,7 @@ for file in files:
info = { info = {
'grid': numpy.zeros(3,'i'), 'grid': numpy.zeros(3,'i'),
'size': numpy.zeros(3,'d'), 'size': numpy.zeros(3,'d'),
'origin': numpy.zeros(3,'d'),
} }
for header in theTable.info: for header in theTable.info:
@ -108,8 +109,8 @@ for file in files:
info[headitems[0]] = mappings[headitems[0]](headitems[1]) info[headitems[0]] = mappings[headitems[0]](headitems[1])
file['croak'].write('grid a b c: %s\n'%(' x '.join(map(str,info['grid']))) + \ file['croak'].write('grid a b c: %s\n'%(' x '.join(map(str,info['grid']))) + \
'size x y z: %s\n'%(' x '.join(map(str,info['size'])))) 'size x y z: %s\n'%(' x '.join(map(str,info['size']))) + \
'origin x y z: %s\n'%(' : '.join(map(str,info['origin']))))
if 0 not in options.grid: # user-specified grid if 0 not in options.grid: # user-specified grid
info['grid'] = numpy.array(options.grid) info['grid'] = numpy.array(options.grid)
if numpy.any(info['grid'] < 1): if numpy.any(info['grid'] < 1):
@ -125,6 +126,7 @@ for file in files:
#--- read data -------------------------------------------------------------------------------- #--- read data --------------------------------------------------------------------------------
theTable.data_readArray(['x','y','z','microstructure']) theTable.data_readArray(['x','y','z','microstructure'])
theTable.data[:,0:3] *= info['size'] theTable.data[:,0:3] *= info['size']
theTable.data[:,0:3] += info['origin']
#--- generate grid -------------------------------------------------------------------------------- #--- generate grid --------------------------------------------------------------------------------
grid = vtk.vtkUnstructuredGrid() grid = vtk.vtkUnstructuredGrid()