[yt-svn] commit/yt: 5 new changesets
Bitbucket
commits-noreply at bitbucket.org
Thu Feb 9 08:27:42 PST 2012
5 new commits in yt:
https://bitbucket.org/yt_analysis/yt/changeset/365052299349/
changeset: 365052299349
branch: yt
user: jwise77
date: 2012-01-28 21:49:41
summary: Replacing 'mpc' with 'Mpc' in unit labels.
affected #: 1 file
diff -r bf342c6ac7d13fd3eb5f129e88c16414ee3bc568 -r 36505229934926db5228c466e9f0afb1330d37be yt/visualization/eps_writer.py
--- a/yt/visualization/eps_writer.py
+++ b/yt/visualization/eps_writer.py
@@ -293,6 +293,7 @@
_xlabel = ""
_ylabel = ""
else:
+ units = units.replace('mpc', 'Mpc')
_xlabel = '%s (%s)' % (x_names[plot.data.axis], units)
_ylabel = '%s (%s)' % (y_names[plot.data.axis], units)
_tickcolor = pyx.color.cmyk.white
https://bitbucket.org/yt_analysis/yt/changeset/71c2bfaa3b5f/
changeset: 71c2bfaa3b5f
branch: yt
user: jwise77
date: 2012-01-28 21:54:15
summary: Bug fixes for improved ghost zone calculation: Fixing...
* some incomplete overlap issues with covering grids and their
original data;
* round-off errors when computing the left and right edges of the box
when cascading down the refinement tree
affected #: 1 file
diff -r 36505229934926db5228c466e9f0afb1330d37be -r 71c2bfaa3b5f911ac5a88242cc36011679e4aa58 yt/data_objects/data_containers.py
--- a/yt/data_objects/data_containers.py
+++ b/yt/data_objects/data_containers.py
@@ -3487,6 +3487,7 @@
self._base_dx = (
(self.pf.domain_right_edge - self.pf.domain_left_edge) /
self.pf.domain_dimensions.astype("float64"))
+ self.global_endindex = None
AMRCoveringGridBase.__init__(self, *args, **kwargs)
self._final_start_index = self.global_startindex
@@ -3498,6 +3499,7 @@
# difference between neighboring areas.
nz = 0
buf = 0.0
+ self.min_level = 0
dl = ((self.global_startindex.astype("float64") + 1)
/ (self.pf.refine_by**self.level))
dr = ((self.global_startindex.astype("float64")
@@ -3507,29 +3509,38 @@
nz = 2 * self.pf.refine_by**self.level
buf = self._base_dx
if nz <= self.pf.refine_by**3: # delta level of 3
- cg = self.pf.h.covering_grid(self.level,
- self.left_edge - buf, self.ActiveDimensions + nz)
- cg._use_pbar = False
- count = cg.ActiveDimensions.prod()
- for g in cg._grids:
- count -= cg._get_data_from_grid(g, [])
- if count <= 0:
- min_level = g.Level
- break
+ last_buf = [None,None,None]
+ count = 0
+ # Repeat until no more grids are covered (up to a delta level of 3)
+ while na.any(buf != last_buf) or count == 3:
+ cg = self.pf.h.covering_grid(self.level,
+ self.left_edge - buf, self.ActiveDimensions + nz)
+ cg._use_pbar = False
+ count = cg.ActiveDimensions.prod()
+ for g in cg._grids:
+ count -= cg._get_data_from_grid(g, [])
+ if count <= 0:
+ self.min_level = g.Level
+ break
+ last_buf = buf
+ # Increase box by 2 cell widths at the min covering level
+ buf = 2*self._base_dx / self.pf.refine_by**self.min_level
+ nz += 4 * self.pf.refine_by**(self.level-self.min_level)
+ count += 1
else:
nz = buf = 0
- min_level = 0
+ self.min_level = 0
# This should not cost substantial additional time.
BLE = self.left_edge - buf
BRE = self.right_edge + buf
if na.any(BLE < self.pf.domain_left_edge) or \
na.any(BRE > self.pf.domain_right_edge):
grids,ind = self.pf.hierarchy.get_periodic_box_grids_below_level(
- BLE, BRE, self.level, min_level)
+ BLE, BRE, self.level, self.min_level)
else:
grids,ind = self.pf.hierarchy.get_box_grids_below_level(
BLE, BRE, self.level,
- min(self.level, min_level))
+ min(self.level, self.min_level))
sort_ind = na.argsort(self.pf.h.grid_levels.ravel()[ind])
self._grids = self.pf.hierarchy.grids[ind][(sort_ind,)]
@@ -3547,11 +3558,11 @@
# We jump-start our task here
mylog.debug("Getting fields %s from %s possible grids",
fields_to_get, len(self._grids))
- self._update_level_state(0, fields_to_get)
+ self._update_level_state(self.min_level, fields_to_get, initialize=True)
if self._use_pbar: pbar = \
get_pbar('Searching grids for values ', len(self._grids))
# The grids are assumed to be pre-sorted
- last_level = 0
+ last_level = self.min_level
for gi, grid in enumerate(self._grids):
if self._use_pbar: pbar.update(gi)
if grid.Level > last_level and grid.Level <= self.level:
@@ -3570,27 +3581,31 @@
raise KeyError(n_bad)
if self._use_pbar: pbar.finish()
- def _update_level_state(self, level, fields = None):
+ def _update_level_state(self, level, fields = None, initialize=False):
dx = self._base_dx / self.pf.refine_by**level
self.field_data['cdx'] = dx[0]
self.field_data['cdy'] = dx[1]
self.field_data['cdz'] = dx[2]
LL = self.left_edge - self.pf.domain_left_edge
+ RL = self.right_edge - self.pf.domain_left_edge
self._old_global_startindex = self.global_startindex
- self.global_startindex = na.rint(LL / dx).astype('int64') - 1
+ self._old_global_endindex = self.global_endindex
+ # We use one grid cell at LEAST, plus one buffer on all sides
+ self.global_startindex = na.floor(LL / dx).astype('int64') - 1
+ self.global_endindex = na.ceil(RL / dx).astype('int64') + 1
self.domain_width = na.rint((self.pf.domain_right_edge -
self.pf.domain_left_edge)/dx).astype('int64')
- if level == 0 and self.level > 0:
- # We use one grid cell at LEAST, plus one buffer on all sides
- idims = na.rint((self.right_edge-self.left_edge)/dx).astype('int64') + 2
+ if (level == 0 or initialize) and self.level > 0:
+ idims = self.global_endindex - self.global_startindex
fields = ensure_list(fields)
for field in fields:
self.field_data[field] = na.zeros(idims,dtype='float64')-999
self._cur_dims = idims.astype("int32")
- elif level == 0 and self.level == 0:
+ elif (level == 0 or initialize) and self.level == 0:
DLE = self.pf.domain_left_edge
self.global_startindex = na.array(na.floor(LL/ dx), dtype='int64')
idims = na.rint((self.right_edge-self.left_edge)/dx).astype('int64')
+ #idims = self.global_endindex - self.global_startindex
fields = ensure_list(fields)
for field in fields:
self.field_data[field] = na.zeros(idims,dtype='float64')-999
@@ -3599,15 +3614,16 @@
def _refine(self, dlevel, fields):
rf = float(self.pf.refine_by**dlevel)
- input_left = (self._old_global_startindex + 0.5) * rf
- dx = na.fromiter((self['cd%s' % ax] for ax in 'xyz'), count=3, dtype='float64')
- output_dims = na.rint((self.right_edge-self.left_edge)/dx+0.5).astype('int32') + 2
+ input_left = (self._old_global_startindex + 0.5) * rf
+ input_right = (self._old_global_endindex - 0.5) * rf
+ output_left = self.global_startindex + 0.5
+ output_right = self.global_endindex - 0.5
+ output_dims = (output_right - output_left + 1).astype('int32')
self._cur_dims = output_dims
for field in fields:
output_field = na.zeros(output_dims, dtype="float64")
- output_left = self.global_startindex + 0.5
ghost_zone_interpolate(rf, self[field], input_left,
output_field, output_left)
self.field_data[field] = output_field
https://bitbucket.org/yt_analysis/yt/changeset/c75bc6741e6a/
changeset: c75bc6741e6a
branch: yt
user: jwise77
date: 2012-01-29 01:16:55
summary: Merging in latest changes from the past week.
affected #: 2 files
diff -r 37fb03bc3fea8a6024b09bc29da1c6f0f48708cd -r c75bc6741e6a56d84825ddaf652024be07044f07 yt/data_objects/data_containers.py
--- a/yt/data_objects/data_containers.py
+++ b/yt/data_objects/data_containers.py
@@ -3483,16 +3483,62 @@
self._base_dx = (
(self.pf.domain_right_edge - self.pf.domain_left_edge) /
self.pf.domain_dimensions.astype("float64"))
+ self.global_endindex = None
AMRCoveringGridBase.__init__(self, *args, **kwargs)
self._final_start_index = self.global_startindex
def _get_list_of_grids(self):
if self._grids is not None: return
- buffer = ((self.pf.domain_right_edge - self.pf.domain_left_edge)
- / self.pf.domain_dimensions).max()
- AMRCoveringGridBase._get_list_of_grids(self, buffer)
- # We reverse the order to ensure that coarse grids are first
- self._grids = self._grids[::-1]
+ # Check for ill-behaved AMR schemes (Enzo) where we may have
+ # root-tile-boundary issues. This is specific to the root tiles not
+ # allowing grids to cross them and also allowing > 1 level of
+ # difference between neighboring areas.
+ nz = 0
+ buf = 0.0
+ self.min_level = 0
+ dl = ((self.global_startindex.astype("float64") + 1)
+ / (self.pf.refine_by**self.level))
+ dr = ((self.global_startindex.astype("float64")
+ + self.ActiveDimensions - 1)
+ / (self.pf.refine_by**self.level))
+ if na.any(dl == na.rint(dl)) or na.any(dr == na.rint(dr)):
+ nz = 2 * self.pf.refine_by**self.level
+ buf = self._base_dx
+ if nz <= self.pf.refine_by**3: # delta level of 3
+ last_buf = [None,None,None]
+ count = 0
+ # Repeat until no more grids are covered (up to a delta level of 3)
+ while na.any(buf != last_buf) or count == 3:
+ cg = self.pf.h.covering_grid(self.level,
+ self.left_edge - buf, self.ActiveDimensions + nz)
+ cg._use_pbar = False
+ count = cg.ActiveDimensions.prod()
+ for g in cg._grids:
+ count -= cg._get_data_from_grid(g, [])
+ if count <= 0:
+ self.min_level = g.Level
+ break
+ last_buf = buf
+ # Increase box by 2 cell widths at the min covering level
+ buf = 2*self._base_dx / self.pf.refine_by**self.min_level
+ nz += 4 * self.pf.refine_by**(self.level-self.min_level)
+ count += 1
+ else:
+ nz = buf = 0
+ self.min_level = 0
+ # This should not cost substantial additional time.
+ BLE = self.left_edge - buf
+ BRE = self.right_edge + buf
+ if na.any(BLE < self.pf.domain_left_edge) or \
+ na.any(BRE > self.pf.domain_right_edge):
+ grids,ind = self.pf.hierarchy.get_periodic_box_grids_below_level(
+ BLE, BRE, self.level, self.min_level)
+ else:
+ grids,ind = self.pf.hierarchy.get_box_grids_below_level(
+ BLE, BRE, self.level,
+ min(self.level, self.min_level))
+ sort_ind = na.argsort(self.pf.h.grid_levels.ravel()[ind])
+ self._grids = self.pf.hierarchy.grids[ind][(sort_ind,)]
def get_data(self, field=None):
self._get_list_of_grids()
@@ -3508,11 +3554,11 @@
# We jump-start our task here
mylog.debug("Getting fields %s from %s possible grids",
fields_to_get, len(self._grids))
- self._update_level_state(0, fields_to_get)
+ self._update_level_state(self.min_level, fields_to_get, initialize=True)
if self._use_pbar: pbar = \
get_pbar('Searching grids for values ', len(self._grids))
# The grids are assumed to be pre-sorted
- last_level = 0
+ last_level = self.min_level
for gi, grid in enumerate(self._grids):
if self._use_pbar: pbar.update(gi)
if grid.Level > last_level and grid.Level <= self.level:
@@ -3530,27 +3576,31 @@
raise KeyError(n_bad)
if self._use_pbar: pbar.finish()
- def _update_level_state(self, level, fields = None):
+ def _update_level_state(self, level, fields = None, initialize=False):
dx = self._base_dx / self.pf.refine_by**level
self.field_data['cdx'] = dx[0]
self.field_data['cdy'] = dx[1]
self.field_data['cdz'] = dx[2]
LL = self.left_edge - self.pf.domain_left_edge
+ RL = self.right_edge - self.pf.domain_left_edge
self._old_global_startindex = self.global_startindex
- self.global_startindex = na.rint(LL / dx).astype('int64') - 1
+ self._old_global_endindex = self.global_endindex
+ # We use one grid cell at LEAST, plus one buffer on all sides
+ self.global_startindex = na.floor(LL / dx).astype('int64') - 1
+ self.global_endindex = na.ceil(RL / dx).astype('int64') + 1
self.domain_width = na.rint((self.pf.domain_right_edge -
self.pf.domain_left_edge)/dx).astype('int64')
- if level == 0 and self.level > 0:
- # We use one grid cell at LEAST, plus one buffer on all sides
- idims = na.rint((self.right_edge-self.left_edge)/dx).astype('int64') + 2
+ if (level == 0 or initialize) and self.level > 0:
+ idims = self.global_endindex - self.global_startindex
fields = ensure_list(fields)
for field in fields:
self.field_data[field] = na.zeros(idims,dtype='float64')-999
self._cur_dims = idims.astype("int32")
- elif level == 0 and self.level == 0:
+ elif (level == 0 or initialize) and self.level == 0:
DLE = self.pf.domain_left_edge
self.global_startindex = na.array(na.floor(LL/ dx), dtype='int64')
idims = na.rint((self.right_edge-self.left_edge)/dx).astype('int64')
+ #idims = self.global_endindex - self.global_startindex
fields = ensure_list(fields)
for field in fields:
self.field_data[field] = na.zeros(idims,dtype='float64')-999
@@ -3559,15 +3609,16 @@
def _refine(self, dlevel, fields):
rf = float(self.pf.refine_by**dlevel)
- input_left = (self._old_global_startindex + 0.5) * rf
- dx = na.fromiter((self['cd%s' % ax] for ax in 'xyz'), count=3, dtype='float64')
- output_dims = na.rint((self.right_edge-self.left_edge)/dx).astype('int32') + 2
+ input_left = (self._old_global_startindex + 0.5) * rf
+ input_right = (self._old_global_endindex - 0.5) * rf
+ output_left = self.global_startindex + 0.5
+ output_right = self.global_endindex - 0.5
+ output_dims = (output_right - output_left + 1).astype('int32')
self._cur_dims = output_dims
for field in fields:
output_field = na.zeros(output_dims, dtype="float64")
- output_left = self.global_startindex + 0.5
ghost_zone_interpolate(rf, self[field], input_left,
output_field, output_left)
self.field_data[field] = output_field
diff -r 37fb03bc3fea8a6024b09bc29da1c6f0f48708cd -r c75bc6741e6a56d84825ddaf652024be07044f07 yt/visualization/eps_writer.py
--- a/yt/visualization/eps_writer.py
+++ b/yt/visualization/eps_writer.py
@@ -293,6 +293,7 @@
_xlabel = ""
_ylabel = ""
else:
+ units = units.replace('mpc', 'Mpc')
_xlabel = '%s (%s)' % (x_names[plot.data.axis], units)
_ylabel = '%s (%s)' % (y_names[plot.data.axis], units)
_tickcolor = pyx.color.cmyk.white
https://bitbucket.org/yt_analysis/yt/changeset/47db55b3a9c1/
changeset: 47db55b3a9c1
branch: yt
user: jwise77
date: 2012-02-03 16:58:34
summary: max_child was not stored, and it is needed when printing out the ASCII
merger tree.
affected #: 1 file
diff -r c75bc6741e6a56d84825ddaf652024be07044f07 -r 47db55b3a9c1b5986ac23aaad0fea4ecee100b12 yt/analysis_modules/halo_merger_tree/enzofof_merger_tree.py
--- a/yt/analysis_modules/halo_merger_tree/enzofof_merger_tree.py
+++ b/yt/analysis_modules/halo_merger_tree/enzofof_merger_tree.py
@@ -295,6 +295,7 @@
Maximum number of child halos each leaf can have.
"""
self.halonum = halonum
+ self.max_children = max_children
self.output_numbers = sorted(self.relationships, reverse=True)
self.levels = {}
trunk = self.output_numbers[0]
@@ -376,7 +377,7 @@
print "--> Most massive progenitor == Halo %d" % \
(br.progenitor)
for i,c in enumerate(br.children):
- if i > max_child: break
+ if i > self.max_children: break
print "--> Halo %8.8d :: fraction = %g" % (c[0], c[1])
def write_dot(self, filename=None):
https://bitbucket.org/yt_analysis/yt/changeset/d023f7617ecd/
changeset: d023f7617ecd
branch: yt
user: jwise77
date: 2012-02-06 18:47:24
summary: Merging latest changes from yt_analysis into my fork.
affected #: 3 files
diff -r 3987bb10d37a09c90b63f1935e555c56981bb59d -r d023f7617ecde1f6bcc60d1c5b1e2cf6d04cc6ef yt/analysis_modules/halo_merger_tree/enzofof_merger_tree.py
--- a/yt/analysis_modules/halo_merger_tree/enzofof_merger_tree.py
+++ b/yt/analysis_modules/halo_merger_tree/enzofof_merger_tree.py
@@ -295,6 +295,7 @@
Maximum number of child halos each leaf can have.
"""
self.halonum = halonum
+ self.max_children = max_children
self.output_numbers = sorted(self.relationships, reverse=True)
self.levels = {}
trunk = self.output_numbers[0]
@@ -376,7 +377,7 @@
print "--> Most massive progenitor == Halo %d" % \
(br.progenitor)
for i,c in enumerate(br.children):
- if i > max_child: break
+ if i > self.max_children: break
print "--> Halo %8.8d :: fraction = %g" % (c[0], c[1])
def write_dot(self, filename=None):
diff -r 3987bb10d37a09c90b63f1935e555c56981bb59d -r d023f7617ecde1f6bcc60d1c5b1e2cf6d04cc6ef yt/data_objects/data_containers.py
--- a/yt/data_objects/data_containers.py
+++ b/yt/data_objects/data_containers.py
@@ -3483,16 +3483,62 @@
self._base_dx = (
(self.pf.domain_right_edge - self.pf.domain_left_edge) /
self.pf.domain_dimensions.astype("float64"))
+ self.global_endindex = None
AMRCoveringGridBase.__init__(self, *args, **kwargs)
self._final_start_index = self.global_startindex
def _get_list_of_grids(self):
if self._grids is not None: return
- buffer = ((self.pf.domain_right_edge - self.pf.domain_left_edge)
- / self.pf.domain_dimensions).max()
- AMRCoveringGridBase._get_list_of_grids(self, buffer)
- # We reverse the order to ensure that coarse grids are first
- self._grids = self._grids[::-1]
+ # Check for ill-behaved AMR schemes (Enzo) where we may have
+ # root-tile-boundary issues. This is specific to the root tiles not
+ # allowing grids to cross them and also allowing > 1 level of
+ # difference between neighboring areas.
+ nz = 0
+ buf = 0.0
+ self.min_level = 0
+ dl = ((self.global_startindex.astype("float64") + 1)
+ / (self.pf.refine_by**self.level))
+ dr = ((self.global_startindex.astype("float64")
+ + self.ActiveDimensions - 1)
+ / (self.pf.refine_by**self.level))
+ if na.any(dl == na.rint(dl)) or na.any(dr == na.rint(dr)):
+ nz = 2 * self.pf.refine_by**self.level
+ buf = self._base_dx
+ if nz <= self.pf.refine_by**3: # delta level of 3
+ last_buf = [None,None,None]
+ count = 0
+ # Repeat until no more grids are covered (up to a delta level of 3)
+ while na.any(buf != last_buf) or count == 3:
+ cg = self.pf.h.covering_grid(self.level,
+ self.left_edge - buf, self.ActiveDimensions + nz)
+ cg._use_pbar = False
+ count = cg.ActiveDimensions.prod()
+ for g in cg._grids:
+ count -= cg._get_data_from_grid(g, [])
+ if count <= 0:
+ self.min_level = g.Level
+ break
+ last_buf = buf
+ # Increase box by 2 cell widths at the min covering level
+ buf = 2*self._base_dx / self.pf.refine_by**self.min_level
+ nz += 4 * self.pf.refine_by**(self.level-self.min_level)
+ count += 1
+ else:
+ nz = buf = 0
+ self.min_level = 0
+ # This should not cost substantial additional time.
+ BLE = self.left_edge - buf
+ BRE = self.right_edge + buf
+ if na.any(BLE < self.pf.domain_left_edge) or \
+ na.any(BRE > self.pf.domain_right_edge):
+ grids,ind = self.pf.hierarchy.get_periodic_box_grids_below_level(
+ BLE, BRE, self.level, self.min_level)
+ else:
+ grids,ind = self.pf.hierarchy.get_box_grids_below_level(
+ BLE, BRE, self.level,
+ min(self.level, self.min_level))
+ sort_ind = na.argsort(self.pf.h.grid_levels.ravel()[ind])
+ self._grids = self.pf.hierarchy.grids[ind][(sort_ind,)]
def get_data(self, field=None):
self._get_list_of_grids()
@@ -3508,11 +3554,11 @@
# We jump-start our task here
mylog.debug("Getting fields %s from %s possible grids",
fields_to_get, len(self._grids))
- self._update_level_state(0, fields_to_get)
+ self._update_level_state(self.min_level, fields_to_get, initialize=True)
if self._use_pbar: pbar = \
get_pbar('Searching grids for values ', len(self._grids))
# The grids are assumed to be pre-sorted
- last_level = 0
+ last_level = self.min_level
for gi, grid in enumerate(self._grids):
if self._use_pbar: pbar.update(gi)
if grid.Level > last_level and grid.Level <= self.level:
@@ -3530,27 +3576,31 @@
raise KeyError(n_bad)
if self._use_pbar: pbar.finish()
- def _update_level_state(self, level, fields = None):
+ def _update_level_state(self, level, fields = None, initialize=False):
dx = self._base_dx / self.pf.refine_by**level
self.field_data['cdx'] = dx[0]
self.field_data['cdy'] = dx[1]
self.field_data['cdz'] = dx[2]
LL = self.left_edge - self.pf.domain_left_edge
+ RL = self.right_edge - self.pf.domain_left_edge
self._old_global_startindex = self.global_startindex
- self.global_startindex = na.rint(LL / dx).astype('int64') - 1
+ self._old_global_endindex = self.global_endindex
+ # We use one grid cell at LEAST, plus one buffer on all sides
+ self.global_startindex = na.floor(LL / dx).astype('int64') - 1
+ self.global_endindex = na.ceil(RL / dx).astype('int64') + 1
self.domain_width = na.rint((self.pf.domain_right_edge -
self.pf.domain_left_edge)/dx).astype('int64')
- if level == 0 and self.level > 0:
- # We use one grid cell at LEAST, plus one buffer on all sides
- idims = na.rint((self.right_edge-self.left_edge)/dx).astype('int64') + 2
+ if (level == 0 or initialize) and self.level > 0:
+ idims = self.global_endindex - self.global_startindex
fields = ensure_list(fields)
for field in fields:
self.field_data[field] = na.zeros(idims,dtype='float64')-999
self._cur_dims = idims.astype("int32")
- elif level == 0 and self.level == 0:
+ elif (level == 0 or initialize) and self.level == 0:
DLE = self.pf.domain_left_edge
self.global_startindex = na.array(na.floor(LL/ dx), dtype='int64')
idims = na.rint((self.right_edge-self.left_edge)/dx).astype('int64')
+ #idims = self.global_endindex - self.global_startindex
fields = ensure_list(fields)
for field in fields:
self.field_data[field] = na.zeros(idims,dtype='float64')-999
@@ -3559,15 +3609,16 @@
def _refine(self, dlevel, fields):
rf = float(self.pf.refine_by**dlevel)
- input_left = (self._old_global_startindex + 0.5) * rf
- dx = na.fromiter((self['cd%s' % ax] for ax in 'xyz'), count=3, dtype='float64')
- output_dims = na.rint((self.right_edge-self.left_edge)/dx).astype('int32') + 2
+ input_left = (self._old_global_startindex + 0.5) * rf
+ input_right = (self._old_global_endindex - 0.5) * rf
+ output_left = self.global_startindex + 0.5
+ output_right = self.global_endindex - 0.5
+ output_dims = (output_right - output_left + 1).astype('int32')
self._cur_dims = output_dims
for field in fields:
output_field = na.zeros(output_dims, dtype="float64")
- output_left = self.global_startindex + 0.5
ghost_zone_interpolate(rf, self[field], input_left,
output_field, output_left)
self.field_data[field] = output_field
diff -r 3987bb10d37a09c90b63f1935e555c56981bb59d -r d023f7617ecde1f6bcc60d1c5b1e2cf6d04cc6ef yt/visualization/eps_writer.py
--- a/yt/visualization/eps_writer.py
+++ b/yt/visualization/eps_writer.py
@@ -293,6 +293,7 @@
_xlabel = ""
_ylabel = ""
else:
+ units = units.replace('mpc', 'Mpc')
_xlabel = '%s (%s)' % (x_names[plot.data.axis], units)
_ylabel = '%s (%s)' % (y_names[plot.data.axis], units)
_tickcolor = pyx.color.cmyk.white
Repository URL: https://bitbucket.org/yt_analysis/yt/
--
This is a commit notification from bitbucket.org. You are receiving
this because you have the service enabled, addressing the recipient of
this email.
More information about the yt-svn
mailing list