New issue 1319: fields as contiguous arrays for idealized galaxies
This is related to a previous bugfix: https://bitbucket.org/yt_analysis/yt/pull-requests/2498/ensure-fields-are...
If I try a related test with an idealized gadget galaxy (a test snapshot can be found here:
with the following code:
it returns an error similar to PR 2498:
jam:scratch desika$ python vel_proj_idealized.py
yt : [INFO ] 2017-02-07 11:14:47,621 Omega Lambda is 0.0, so we are turning off Cosmology.
yt : [INFO ] 2017-02-07 11:14:47,630 Parameters: current_time = 0.25
yt : [INFO ] 2017-02-07 11:14:47,630 Parameters: domain_dimensions = [2 2 2]
yt : [INFO ] 2017-02-07 11:14:47,631 Parameters: domain_left_edge = [-1000. -1000. -1000.]
yt : [INFO ] 2017-02-07 11:14:47,631 Parameters: domain_right_edge = [ 1000. 1000. 1000.]
yt : [INFO ] 2017-02-07 11:14:47,632 Parameters: cosmological_simulation = 0
yt : [INFO ] 2017-02-07 11:14:47,634 Allocating for 2.515e+07 particles (index particle type 'all')
yt : [INFO ] 2017-02-07 11:14:51,665 Identified 1.540e+06 octs
[ 0. 0. 0. ..., 0. 0. 0.] code_velocity
Traceback (most recent call last):
File "vel_proj_idealized.py", line 23, in <module>
vy = arb[ ('deposit', 'PartType0_sum_velocity_y')]
File "/Users/desika/yt/yt/data_objects/data_containers.py", line 281, in __getitem__
File "/Users/desika/yt/yt/data_objects/construction_data_containers.py", line 635, in get_data
if len(gen) > 0: self._generate_fields(gen)
File "/Users/desika/yt/yt/data_objects/data_containers.py", line 1353, in _generate_fields
fd = self._generate_field(field)
File "/Users/desika/yt/yt/data_objects/data_containers.py", line 318, in _generate_field
tr = self._generate_fluid_field(field)
File "/Users/desika/yt/yt/data_objects/data_containers.py", line 338, in _generate_fluid_field
rv = finfo(gen_obj)
File "/Users/desika/yt/yt/fields/derived_field.py", line 220, in __call__
dd = self._function(self, data)
File "/Users/desika/yt/yt/data_objects/static_output.py", line 1260, in _deposit_field
File "/Users/desika/yt/yt/data_objects/construction_data_containers.py", line 728, in deposit
op.process_grid(self, positions, fields)
File "yt/geometry/particle_deposit.pyx", line 120, in yt.geometry.particle_deposit.ParticleDepositOperation.process_grid (yt/geometry/particle_deposit.c:5181)
ValueError: Buffer and memoryview are not contiguous in the same dimension.
New issue 1318: Arepo Temperature bug?
It seems like something recently changed in the frontend that underestimates gas temperatures in Arepo file structures. I recently went back to a simulation to re-evaluate a temperature-based calculation, and found that all the ('gas','temperature') values seem too low by several orders of magnitude.
I'm not sure what's going on, but for example when I compare ('PartType0','Temperature') to ('gas','temperature') in a small volume of the simulation (cosmo. simulation centered on the galaxy), you can see that the distribution of temperatures look very different/offset from each other. (Histograms attached).
Any help is appreciated.
I just issued a pull request that backports several bugfixes to the stable
This is in preparation to do a release of yt 3.3.4. Even though there are
only a few backports this time I'd like to get it out the door since there
are some fixes that avoid crashes under the latest version of numpy as well
as annoying deprecation warnings under matplotlib 2.0.
Please let me know if I've missed any bugfixes that you'd like to see
New issue 1316: Unit conversion with equivalency
The following script fails.
ds = yt.load("IsolatedGalaxy/galaxy0030/galaxy0030")
slc = yt.SlicePlot(ds, 'z', 'temperature', width=(10,'kpc'))
#this line causes a problem
slc.set_unit('temperature', 'keV', equivalency='thermal')
yt.utilities.exceptions.YTUnitConversionError: Unit dimensionalities do not match. Tried to convert between K**2/erg (dim (temperature)**2*(time)**2/((length)**2*(mass))) and K (dim (temperature)).
This issue is related to pull request #2424. Note that calling set_buff_size() AFTER set_unit() does not trigger this issue.