From f5177b6d2dc3150a2eddbfa733959e2fdd221c76 Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Tue, 20 Feb 2018 10:20:35 -0500 Subject: [PATCH] [skip ci] added option for non-periodic geometries --- processing/pre/geom_vicinityOffset.py | 9 +++++++-- 1 file changed, 7 insertions(+), 2 deletions(-) diff --git a/processing/pre/geom_vicinityOffset.py b/processing/pre/geom_vicinityOffset.py index 263b4e8d9..d4eb8e097 100755 --- a/processing/pre/geom_vicinityOffset.py +++ b/processing/pre/geom_vicinityOffset.py @@ -29,9 +29,14 @@ parser.add_option('-m', '--microstructureoffset', type = 'int', metavar = 'int', help = 'offset (positive or negative) for tagged microstructure indices. '+ '"0" selects maximum microstructure index [%default]') +parser.add_option('-n', '--nonperiodic', + dest = 'mode', + action = 'store_const', const = 'nearest', + help = 'assume geometry to be non-periodic') parser.set_defaults(vicinity = 1, offset = 0, + mode = 'wrap', ) (options, filenames) = parser.parse_args() @@ -79,8 +84,8 @@ for name in filenames: if options.offset == 0: options.offset = microstructure.max() - microstructure = np.where(ndimage.filters.maximum_filter(microstructure,size=1+2*options.vicinity,mode='wrap') == - ndimage.filters.minimum_filter(microstructure,size=1+2*options.vicinity,mode='wrap'), + microstructure = np.where(ndimage.filters.maximum_filter(microstructure,size=1+2*options.vicinity,mode=options.mode) == + ndimage.filters.minimum_filter(microstructure,size=1+2*options.vicinity,mode=options.mode), microstructure, microstructure + options.offset) newInfo['microstructures'] = microstructure.max()