exchanged by more modern vtk_* machinery...
This commit is contained in:
Philip Eisenlohr
parent
9abc0dbeeb
commit
478acb96ce
|
|
@ -1,108 +0,0 @@
|
||||||
#!/usr/bin/env python
|
|
||||||
# -*- coding: UTF-8 no BOM -*-
|
|
||||||
|
|
||||||
import os,sys,vtk
|
|
||||||
import numpy as np
|
|
||||||
from optparse import OptionParser
|
|
||||||
import damask
|
|
||||||
|
|
||||||
scriptName = os.path.splitext(os.path.basename(__file__))[0]
|
|
||||||
scriptID = ' '.join([scriptName,damask.version])
|
|
||||||
|
|
||||||
#--------------------------------------------------------------------------------------------------
|
|
||||||
# MAIN
|
|
||||||
#--------------------------------------------------------------------------------------------------
|
|
||||||
|
|
||||||
parser = OptionParser(option_class=damask.extendableOption, usage='%prog [file[s]]', description = """
|
|
||||||
Produce VTK rectilinear mesh of structure data from geom description
|
|
||||||
|
|
||||||
""", version = scriptID)
|
|
||||||
|
|
||||||
parser.add_option('-m','--nodata',
|
|
||||||
dest = 'data',
|
|
||||||
action = 'store_false',
|
|
||||||
help = 'generate mesh without microstructure index data')
|
|
||||||
|
|
||||||
parser.set_defaults(data = True,
|
|
||||||
)
|
|
||||||
|
|
||||||
(options, filenames) = parser.parse_args()
|
|
||||||
|
|
||||||
# --- loop over input files -------------------------------------------------------------------------
|
|
||||||
|
|
||||||
if filenames == []: filenames = [None]
|
|
||||||
|
|
||||||
for name in filenames:
|
|
||||||
try:
|
|
||||||
table = damask.ASCIItable(name = name,
|
|
||||||
buffered = False, labeled = False, readonly = True)
|
|
||||||
except: continue
|
|
||||||
damask.util.report(scriptName,name)
|
|
||||||
|
|
||||||
# --- interpret header ----------------------------------------------------------------------------
|
|
||||||
|
|
||||||
table.head_read()
|
|
||||||
info,extra_header = table.head_getGeom()
|
|
||||||
|
|
||||||
damask.util.croak(['grid a b c: %s'%(' x '.join(map(str,info['grid']))),
|
|
||||||
'size x y z: %s'%(' x '.join(map(str,info['size']))),
|
|
||||||
'origin x y z: %s'%(' : '.join(map(str,info['origin']))),
|
|
||||||
'homogenization: %i'%info['homogenization'],
|
|
||||||
'microstructures: %i'%info['microstructures'],
|
|
||||||
])
|
|
||||||
|
|
||||||
errors = []
|
|
||||||
if np.any(info['grid'] < 1): errors.append('invalid grid a b c.')
|
|
||||||
if np.any(info['size'] <= 0.0): errors.append('invalid size x y z.')
|
|
||||||
|
|
||||||
#--- read microstructure information --------------------------------------------------------------
|
|
||||||
|
|
||||||
if options.data:
|
|
||||||
microstructure,ok = table.microstructure_read(info['grid'],strict = True) # read microstructure
|
|
||||||
|
|
||||||
if ok:
|
|
||||||
structure = vtk.vtkIntArray()
|
|
||||||
structure.SetName('Microstructures')
|
|
||||||
for idx in microstructure: structure.InsertNextValue(idx)
|
|
||||||
|
|
||||||
else: errors.append('mismatch between data and grid dimension.')
|
|
||||||
|
|
||||||
if errors != []:
|
|
||||||
damask.util.croak(errors)
|
|
||||||
table.close(dismiss = True)
|
|
||||||
continue
|
|
||||||
|
|
||||||
# --- generate VTK rectilinear grid ---------------------------------------------------------------
|
|
||||||
|
|
||||||
grid = vtk.vtkRectilinearGrid()
|
|
||||||
grid.SetDimensions([x+1 for x in info['grid']])
|
|
||||||
for i in xrange(3):
|
|
||||||
temp = vtk.vtkDoubleArray()
|
|
||||||
temp.SetNumberOfTuples(info['grid'][i]+1)
|
|
||||||
for j in xrange(info['grid'][i]+1):
|
|
||||||
temp.InsertTuple1(j,j*info['size'][i]/info['grid'][i]+info['origin'][i])
|
|
||||||
if i == 0: grid.SetXCoordinates(temp)
|
|
||||||
elif i == 1: grid.SetYCoordinates(temp)
|
|
||||||
elif i == 2: grid.SetZCoordinates(temp)
|
|
||||||
|
|
||||||
|
|
||||||
if options.data: grid.GetCellData().AddArray(structure)
|
|
||||||
|
|
||||||
# --- write data -----------------------------------------------------------------------------------
|
|
||||||
if name:
|
|
||||||
writer = vtk.vtkXMLRectilinearGridWriter()
|
|
||||||
(directory,filename) = os.path.split(name)
|
|
||||||
writer.SetDataModeToBinary()
|
|
||||||
writer.SetCompressorTypeToZLib()
|
|
||||||
writer.SetFileName(os.path.join(directory,os.path.splitext(filename)[0]+'.'+writer.GetDefaultFileExtension()))
|
|
||||||
else:
|
|
||||||
writer = vtk.vtkDataSetWriter()
|
|
||||||
writer.WriteToOutputStringOn()
|
|
||||||
writer.SetHeader('# powered by '+scriptID)
|
|
||||||
|
|
||||||
if vtk.VTK_MAJOR_VERSION <= 5: writer.SetInput(grid)
|
|
||||||
else: writer.SetInputData(grid)
|
|
||||||
writer.Write()
|
|
||||||
if name is None: sys.stdout.write(writer.GetOutputString()[0:writer.GetOutputStringLength()])
|
|
||||||
|
|
||||||
table.close()
|
|
||||||
|
|
@ -0,0 +1,15 @@
|
||||||
|
#!/bin/bash
|
||||||
|
|
||||||
|
for geom in "$@"
|
||||||
|
do
|
||||||
|
vtk_rectilinearGrid \
|
||||||
|
--geom $geom
|
||||||
|
|
||||||
|
geom_toTable \
|
||||||
|
< $geom \
|
||||||
|
| \
|
||||||
|
vtk_addRectilinearGridData \
|
||||||
|
--scalar microstructure \
|
||||||
|
--inplace \
|
||||||
|
--vtk ${geom%.*}.vtr
|
||||||
|
done
|
||||||
Loading…
Reference in New Issue