[yt-dev] Issue #1256: OffAxisProjection Fails to find proper dataset limits (yt_analysis/yt)

Cameron Hummels issues-reply at bitbucket.org
Fri Aug 5 16:44:32 PDT 2016


New issue 1256: OffAxisProjection Fails to find proper dataset limits
https://bitbucket.org/yt_analysis/yt/issues/1256/offaxisprojection-fails-to-find-proper

Cameron Hummels:

In trying to generate some OffAxisProjections from outside of the simulation domain, I've stumbled on some problems with auto-determination of field extrema.  This simple script demonstrates:

```
#!python

import yt
ds = yt.load("IsolatedGalaxy/galaxy0030/galaxy0030")
yt.OffAxisProjectionPlot(ds, [1., 1., 0.5], 'density', [0.5, 0.5, 0.5], 1.5, north_vector=[0., 0., 1.]).save()
```

yields this STDERR:


```
#!python

yt : [INFO     ] 2016-08-05 16:38:13,256 Parameters: current_time              = 0.00600002000283
yt : [INFO     ] 2016-08-05 16:38:13,256 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2016-08-05 16:38:13,256 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2016-08-05 16:38:13,257 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2016-08-05 16:38:13,257 Parameters: cosmological_simulation   = 0.0
^MParsing Hierarchy :   0%|          | 0/173 [00:00<?, ?it/s]^MParsing Hierarchy : 100%|##########| 173/173 [00:00<00:00, 20876.19it/s]
yt : [INFO     ] 2016-08-05 16:38:13,287 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2016-08-05 16:38:15,246 xlim = -0.750000 0.750000
yt : [INFO     ] 2016-08-05 16:38:15,246 ylim = -0.750000 0.750000
yt : [INFO     ] 2016-08-05 16:38:15,246 zlim = -0.500000 0.500000
yt : [INFO     ] 2016-08-05 16:38:15,250 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800
/Users/chummels/src/yt-x86_64/src/yt-hg/yt/units/yt_array.py:964: RuntimeWarning: overflow encountered in multiply
  np.multiply(self, oth, out=self)
yt : [WARNING  ] 2016-08-05 16:38:16,264 Plot image for field ('gas', 'density') has both positive and negative values. Min = -inf, Max = 54704375419815931049795757413255569402681538393870135161391039774901627123390825374482250621889063954180467028564003619350195202205284149300884842122024669609440680126995877926691701507872243776797761225104772892079277541233926298501475834441217774443211945332828680339300087327793097590112761046546513920.000000.
yt : [WARNING  ] 2016-08-05 16:38:16,264 Switching to symlog colorbar scaling unless linear scaling is specified later
```

which results in this image:

![xDwRi6w.png](https://bitbucket.org/repo/BXbAb/images/2802526612-xDwRi6w.png)

Interestingly, I don't get this problem with different datasets, or with different views of *this* dataset, so I'm not entirely sure what is going on.  Here is the output of the same script using the `enzo_cosmology_plus/RD0009/RD0009` dataset:

![VTGN4cb.png](https://bitbucket.org/repo/BXbAb/images/3178170910-VTGN4cb.png)




More information about the yt-dev mailing list