New issue 1272: Error when using load_octree https://bitbucket.org/yt_analysis/yt/issues/1272/error-when-using-load_octre...

Anonymous:

The following code:

#!python

import numpy as np

from yt.mods import load_octree, ProjectionPlot

octree_mask = np.array([8,0,0,0,0,8,0,0,0,0,0,0,0,0,0,0,8,0,0,0,0,0,0,0,0],
dtype=np.uint8)

quantities = {}
quantities[('gas', 'density')] = np.ones((22,1), dtype=float)

bbox = np.array([[-10., 10.], [-10., 10.], [-10., 10.]])

spf = load_octree(octree_mask=octree_mask,
                  data=quantities,
                  bbox=bbox,
                  over_refine_factor=0,
                  partial_coverage=0)

p = ProjectionPlot(spf, 'x', ['density'], center='c', origin='native')

produces the following error:

yt : [INFO     ] 2016-09-04 22:19:26,815 Parameters: current_time             
= 0.0
yt : [INFO     ] 2016-09-04 22:19:26,815 Parameters: domain_dimensions         = [1 1 1]
yt : [INFO     ] 2016-09-04 22:19:26,816 Parameters: domain_left_edge          = [-10.
-10. -10.]
yt : [INFO     ] 2016-09-04 22:19:26,817 Parameters: domain_right_edge         = [ 10. 
10.  10.]
yt : [INFO     ] 2016-09-04 22:19:26,817 Parameters: cosmological_simulation   = 0.0
Traceback (most recent call last):
  File "debug.py", line 18, in <module>
    p = ProjectionPlot(spf, 'x', ['density'], center='c', origin='native')
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/visualization/plot_window.py",
line 1433, in __init__
    max_level=max_level)
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/data_objects/construction_data_containers.py",
line 262, in __init__
    self.get_data(field)
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/data_objects/construction_data_containers.py",
line 336, in get_data
    self._handle_chunk(chunk, fields, tree)
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/data_objects/construction_data_containers.py",
line 436, in _handle_chunk
    d = chunk[field] * dl
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/data_objects/data_containers.py",
line 272, in __getitem__
    self.get_data(f)
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/data_objects/data_containers.py",
line 1178, in get_data
    fluids, self, self._current_chunk)
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/geometry/geometry_handler.py",
line 245, in _read_fluid_fields
    chunk_size)
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/frontends/stream/io.py",
line 252, in _read_fluid_selection
    subset.fill(field_vals, rv, selector, ind)
  File
"/Users/tom/miniconda3/envs/dev35/lib/python3.5/site-packages/yt/frontends/stream/data_structures.py",
line 1402, in fill
    selector, self.domain_id, cell_count)
  File "yt/geometry/oct_container.pyx", line 775, in
yt.geometry.oct_container.OctreeContainer.file_index_octs
(yt/geometry/oct_container.c:13230)
ValueError: Buffer dtype mismatch, expected 'uint8_t' but got 'long'

I'm using yt 3.3.1 from conda-forge and Python 3.5. This works correctly with yt 3.2.3 from conda-forge though, so I think it must be a bug in the 3.3 series?