From 478acb96ceae26723d75c6c0be7945aa0a2ef53d Mon Sep 17 00:00:00 2001
From: Philip Eisenlohr <eisenlohr@egr.msu.edu>
Date: Sun, 24 Apr 2016 13:38:36 -0500
Subject: [PATCH] exchanged by more modern vtk_* machinery...

---
 processing/pre/geom_check.py | 108 -----------------------------------
 processing/pre/geom_check.sh |  15 +++++
 2 files changed, 15 insertions(+), 108 deletions(-)
 delete mode 100755 processing/pre/geom_check.py
 create mode 100755 processing/pre/geom_check.sh

diff --git a/processing/pre/geom_check.py b/processing/pre/geom_check.py
deleted file mode 100755
index 518e8c3b0..000000000
--- a/processing/pre/geom_check.py
+++ /dev/null
@@ -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()
diff --git a/processing/pre/geom_check.sh b/processing/pre/geom_check.sh
new file mode 100755
index 000000000..05e524049
--- /dev/null
+++ b/processing/pre/geom_check.sh
@@ -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